Hi Mohammad,
when you changed the position of the cement, did you also adjust the
positions in other places such as x_min_cement in porosityAtPos?
In the first two zones, for cementPorosity_ and fc_porosity, you use
polynomials in x and fc respectively. This is not translation invariant,
should this perhaps be using globalPos[dimWorld-2] - x_min_cement_ and
fc1 respectively?
Both of these polynomials use coefficients gives as quadratics in time.
Do these give sensible values for the entirety of the large timeframe
you are testing?
Kind regards,
Mathis
On 10/13/22 11:17, Mohammad Hodroj (Student) wrote:
Dear DuMux Community,
I hope this email finds you well and safe.
I am writing to kindly ask for your support in the following:
a) I am working on porousmedium flow model (injection problem)
b)I have a domain of (10m x 10m) and it includes an aquifer, cement
layer, and overburden (as you can see in the attached image1)
c) all layers have a constant porosity and permeability unless the
cement layer
d) the cement layer is defined as following
bool isInCement_(const GlobalPosition &globalPos) const
{
return globalPos[dimWorld-2] > 5 + eps_ &&
globalPos[dimWorld-2] < 5.2+ eps_ && globalPos[dimWorld-1] > 0 + eps_
&& globalPos[dimWorld-1] < 7+ eps_; }
e) the cement layer has a length of 20 cm and the porosity changes as
function of globalPos[2] and time
f) the cement layer permeability is dependent on porosity
g) I am using cctpfa for discretization
h) the grid as following:
Positions0 = 0 2 5 5.2 10
Positions1 = 0 10
Cells0 = 100 100 100 100
Cells1 = 200
Grading0 = 1.0 1.0 1.0 1.0
Grading1 = 1.0
I) the resolution in the cement layer (x-direction) should be small
(e.g. 0.001).
*Issue:
*
*
*
****a) at the beginning of the simulation I get correct porosity
values of the cement but after few minutes I start to get very large
porosity values (either positive or negative) and I am not able to
understand the reason behind that.
*note: *when I change the position of the cement layer from the middle
of the domain (5-5.2) to the left of the domain (0-0.2) as you can see
below, I got the right porosity values of the cement layer (image 2).
bool isInCement_(const GlobalPosition &globalPos) const
{
return globalPos[dimWorld-2] > 0 + eps_ &&
globalPos[dimWorld-2] < 0.2+ eps_ && globalPos[dimWorld-1] > 0 + eps_
&& globalPos[dimWorld-1] < 7+ eps_; }
Positions0 = 0 0.2 5 10
Positions1 = 0 10
Cells0 = 100 100 100 100
Cells1 = 200
Grading0 = 1.0 1.0 1.0 1.0
Grading1 = 1.0
I attached the files for things to be clear.
Looking forward to your help.
Best Regards,
Mohammad Hodroj
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
********************************************************
Mathis Kelm
Universitaet Stuttgart
Institut fuer Wasser- und Umweltsystemmodellierung
Lehrstuhl fuer Hydromechanik und Hydrosystemmodellierung
Pfaffenwaldring 61, 70569 Stuttgart
Tel.: 0711/685-60146
[email protected]
https://www.iws.uni-stuttgart.de/lh2/
********************************************************
_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
