[Orgogozo]

Dear Foamers,


Please find attached a homemade OpenFOAM solver for Richards equation, an equation of flow of water within variably saturated porous media (e.g. : soils). You may find additional details related to RichardsFOAM in the following paper :
http://authors.elsevier.com/a/1PsPx2OInFfaV

A careful attention is paid to the parallel behaviour of the solver.
Key words: geosciences, hydrosciences, transfer into soils, flows in porous media,
parallel computing, strong scaling, weak scaling,
impact of I/O on computational performances.

Best regards,

Laurent Orgogozo Universtity Of Toulouse Geosciences Environment Toulouse Laboratory

[Orgogozo]

Dear Foamers,

One can find additionnal information on RichardsFOAM (e.g.: conference presentation) on this web page :
www.richardsfoam.info

All the best,

Laurent

[wyldckat]

Greetings Laurent,

I've created a basic wiki page for your solver at openfoamwiki.net, so that it's a bit easier for people to find it: http://openfoamwiki.net/index.php/Contrib/RichardsFOAM

Best regards,
Bruno

[Orgogozo]

Dear foamers,

A new version of RichardsFoam, RichardsFoam2, has been recently published in Computer Physics Communications.
http://authors.elsevier.com/a/1Rkjj2OInJTkZ

The main new features are the direct handling of the evapotranpiration and the possibility of deal easily with fully heterogeneous porous media. Thus this new version is more suitable for modelling of water transfers in the vadose zone.

You will find additionnal details (e.g.: source code, conferences and seminar slides) on
www.richardsfoam.info

Cheers,

Laurent

[Orgogozo]

Dear Foamers,

The solver RichardsFoam2 in its former version is not compatible with OF-v1606+. In order to fix that a new version has been done, and you may find this new version at the following url :
www.richardsfoam.info

All the best,

Laurent

[xirtam]

Dear Laurent, I need to simulate the infiltration of water into the soil. I was able to compile RihcardsFoam2 using the cfdsupport.com version for Windows. But it's showing 'nan' when I run. Initially I will not need the effects of evapotranspiration. Would you help me?

[Orgogozo]

Dear Xirtam,
I never used OF under a windows OS, so I fear I can't help you.
Regards,
Laurent

[rakesh927]

Dear Laurent,

I could compile RichardsFoam2 in OF-v1712, however while running the validation case, it's showing the error

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
#4  ? at ??:?
#5  ? at ??:?
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? at ??:?
Floating point exception (core dumped)

Kindly help me, where it's getting wrong. Is it related to setting up of initial conditions of the domain?

Thanks,
Rakesh

[Orgogozo]

Dear Rakesh,
Did you try the updated version of RichardsFoam2 avalaible at :
www.richardsfoam.info
?

[rakesh927]

Dear Laurent,
Yes, I have tried RichardsFoam2 for OF-v1606+ mentioned in the link.

[Orgogozo]

OK, then there are clearly somethings to be done to solve the problem. I have no time left for that now, but I will do it soon, likely by the end of april, surely before the end of may. Until that, please stay in of-v1606+ when you want to use RichardsFoam2.

All the best,

Laurent

[rakesh927]

I got the solution for the error I was facing with.

Simply replace the function

pos()

with

pos0()

in the original code of Orgogozo and recompile the solver.


Rakesh

[Orgogozo]

Dear Rakesh,

Thanks a lot for the hint ! I will deliver an updated RichardsFoam2 package for OF-v1712 in the following weeks.

All the best,

Laurent

[wavefunction]

Hello. Could anyone please let me know what is the function of the pos() and pos0() functions?

In interMixingFoam, a volScalarField called nearInterface is defined that takes the value of 1 if it is near the interface and 0 if it is not.

Foam::tmp<Foam::volScalarField>
Foam::threePhaseInterfaceProperties::nearInterface() const
{
    return max
    (
        pos0(mixture_.alpha1() - 0.01)*pos0(0.99 - mixture_.alpha1()),
        pos0(mixture_.alpha2() - 0.01)*pos0(0.99 - mixture_.alpha2())
    );
}

How is this done?

[Orgogozo]

Dear Saichara,
pos(x), where x is a volScalarField, is a volScalarField that is equal to 0 wherever x is negative and equal to 1 wherever x is positive.
Regards,
Laurent

[rakesh927]

pos(x)=
(x>0)?1:0

pos0(x)=
(x>=0)?1:0

[wavefunction]

Thanks a lot for your replies. You helped me out a lot. Sorry I couldn't find it out myself. I tried using the grep command as well but I was unable to find what this means. Anyways, now I know.

[Orgogozo]

Dear foamers,


Please find in ResearchGate the updated RichardsFoam2 package, compatible with the version v1712 of OpenFOAM.


https://www.researchgate.net/publica...OpenFOAM-v1712


All the best,


Laurent

[Orgogozo]

Dear Foamers,


I kindly drive your attention on a new born OpenFOAM solver for Geosciences: the permaFoam solver (Orgogozo et al., 2019), which deals with coupled heat and water transfers in porous media with freeze/thaw of the pore water. The main related field is permafrost modelling, but there are also other fields of applications such as geotechnics or environmental engineering. Please find below the web-sites in which you can find the solver and the related paper.


All the best,


Laurent Orgogozo


https://www.researchgate.net/publica...rSjAvFKgNhOvVQ



https://www.researchgate.net/publica...ozo_et_al_2019

https://www.researchgate.net/publica..._rooting_depth



https://onlinelibrary.wiley.com/doi/....1002/ppp.1995

[Orgogozo]

Dear foamers,

Please find through the link below the new version announcement for RichardsFoam3, an update of RichardsFoam that encompasses new boundary conditions that are handier for hydrogeology, more flexible use through on-the-fly handling of linearisation processes, and other improvements.

https://doi.org/10.1016/j.cpc.2021.108182

The RichardsFoam3 package may be found both in the Computer Physics Communications Code Library and in the develop.openfoam.com/community interface.

https://doi.org/10.17632/vkr7sd6fhb.1

https://develop.openfoam.com/Community/hydrology/

This new version has been developped in the framework of the HiPerBorea project (hiperborea.omp.eu).

Best wishes,

Laurent Orgogozo