Skip to content
Snippets Groups Projects
Commit f46ef543 authored by tnagel's avatar tnagel Committed by Dmitri Naumov
Browse files

Correction to bulk modulus computation

parent 26ef137b
No related branches found
No related tags found
No related merge requests found
...@@ -668,7 +668,8 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement, ...@@ -668,7 +668,8 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement,
t, x_position, dt, temperature); t, x_position, dt, temperature);
auto const K_SR = auto const K_SR =
_ip_data[ip].solid_material.getBulkModulus(t, x_position, &C_el); _ip_data[ip].solid_material.getBulkModulus(t, x_position, &C_el) /
(1 - alpha);
auto const K_LR = auto const K_LR =
liquid_phase.property(MPL::PropertyType::bulk_modulus) liquid_phase.property(MPL::PropertyType::bulk_modulus)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment