I am working on a poroelastic model in which I have initialized pore pressure in the entire domain (solid cylinder) under “initial values” in the Darcy’s law. On all the boundaries of this cylinder, I have applied a different pressure (which is less than the pore pressure) as the boundary condition.

I ran the transient solution for this model and the results are not as expected. For t=0, the pressure in the cylinder (slice plot) shows that the pressure is equal to the pressure applied at the boundary and not the initial value of the pore pressure.

After disabling the poroelasticity, I ran the model again (solid mechanics and Dary’s law) and then the pressure in the cylinder (slice plot) at t=0 is equal to the initialized value of the pore pressure.

I am not sure why there is a discrepancy in the result when the poroelasticity is enabled in the model.

I look forward to hearing your comments.

Thanks,

Shashank