Minus liquid pressure occurs when adding boundary condition in a volumetric mesh in TH2M model
Description
I am running an example about the stress distribution during water injection into a reservoir using TH2M model in OGS6.
Steps to Reproduce the Problem
I tried to use the temperature boundary condition in a sphere volumetric mesh representing the injection well instead of the heat source term to represent the cold water injection, in order to test thermal response in the reservoir. Under this condition, the temperature of the whole injection sphere is a certain value, I guess. The initial temperature condition is 353.15K. For testing, the initial temperature difference is 1K between the injection volume and the reservoir. At the beginning, the pressure source term is only 1e-13kg/(m3·s), and increases with time to 5.96e-3 kg/(m3·s) at last.
Expected behavior: [What you expected to happen] Under 1K temperature difference condition, the liquid pressure is plus and quite small(about 1e-4Pa) in 100s result and will be larger with the increase of the pressure source term.
Actual behavior: [What actually happened] In the actual result after 100s’ injection the liquid pressure is minus around the injection well and quite large, about -1e6Pa.
Specifications
- Version: (ogs6.4.1)cold_water_inj_test.prj
- Platform: (Windows10)