Am I setting up my PyLith simulation correctly?

Hello,

My goal is to model the mechanical response of a wastewater injection site.

I would like to know if I am setting up my problem correctly

My starting point was Step01_inflation.cfg from the magma-2d example.

This is a schematic of what I am trying to model:

I’m trying to model a 2D slice at depth of a wastewater injection site.
The dimensions of the mesh are the following:

I’m trying to fix all the boundaries so that they don’t move, and I hope I’ve done that correctly here:


In terms of the timestep, I’m trying to model a monthly timestep because I only have monthly injection pressure averages:

One of my questions regarding the initial conditions above is the following: instead of explicitly including gravity and reference stresses in the initial conditions, my thought is it implicitly include gravity and stresses by assigning the domain pressure to be 53 MPa. I got that value by simply calculating the lithostatic weight (pressure) as rho * g * z.

Can I do this? Would my 53 MPa pressure value be considered my lithostatic weight based on my model?

Lastly, in my pylithapp.cfg file, I didn’t update the first part of the script and only updated the material properties. However, I have a question regarding the scales for nondimensionalization:

How do I go about choosing the correct values for these?

Thank you,
Eddie

The displacement boundary conditions look correct.

In the current release, turning on gravitational body forces does not appear to be working based on reports from users. However, you can use the body force that matches the overburden pressure. I think this is what you want.

Sorry for the delayed response.