How to "turn on" gravity for magma-2d example?

PyLith
1

Hello developers and users,

I want to incorporate gravity and body forces in the magma-2d: step01_inflation example. I set the gravity_field in the problem section of my step01_inflation.cfg file like so:

image
image1323×1059 73.2 KB

I also included “spatialdata.spatialdb.GravityField” in the features section of the pylithapp.cfg file.

Is there anywhere else where I need to specify gravity or body forces? Do I need to include gravity or body forces in the materials section of the pylithapp.cfg file in the db_auxiliary_field.values/ db_auxiliary_field.data section?

2

You may want to balance the gravitational stresses to prevent the domain from collapsing when you turn on gravity. See the reverse-2d example. You can do this either with reference stresses or incompressible material.

Cheers,
Charles

3

Hello,

I am having trouble implementing a reference stress/strain.

While reading through the manual I see that in the materials section, there is a subfield for reference stress and reference strain that can be included in the auxiliary field spatial database as shown here:

image
image849×1214 116 KB

How do I include the reference stress and reference strain in the pylithapp.cfg file?

Best,
Eddie

4

There is another issue with the reference stress/strain for poro that is being updated and tested at the moment.

I’ll make a post on the forum when it is merged, and perhaps post an example problem.