The Porothermoelastic Stress Evolution of Inhomogeneous Solid due to Near-Surface Reservoir

Figure 2

Flowchart describing the proposed solution: (1) define inclusion shape and location, (2) specify matrix elastic properties, (3) calculate porothermal eigenstrains, (4) determine the elliptic and potential functions, (5) solve the displacement field, (6) compute total strain, (7) check the point location inside or outside the region, (8) calculate the stress at the interior or exterior field, and (9) output stress distribution.