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