Darcy-Forchheimer OpenFOAM

MiguelMValero · March 21, 2023, 07:27

Hello everybody,

I have a concern about how the Darcy law is implemented in OpenFOAM. More specifically, I don't really understand the meaning of the vectors "d" and "f". In my case, I want:

d [0 -2 0 0 0 0 0] (2.71508e6 70891 69.2805)
f [0 -1 0 0 0 0 0] (0 0 0)

where I thought that "d" and "f" were "D" and "F", respectively, in the Darcy formula:

$S_i = -\nu (D_{ij}+\frac{1}{2}\rho |u_{kk}| F_{ij})u_i$ .

However, when I tried to manually implement Darcy-Forchheimer in OpenFOAM by removing the file "porosityProperties" file and changing the incompressible solver "pimpleFoam" with the following expression in those cells affected by the porosity condition, I'm not getting the same results.

tmp<fvVectorMatrix> tUEqn
(
fvm::ddt(U) + fvm::div(phi, U)
+ MRF.DDt(U)
+ turbulence->divDegSigma(U)
+ nu*cmptMultiply(D, U)
==
fvOptions(U)
);

I appreciate any help you can give me. I think I misunderstood how the multiplication of D and U is performed. I already checked on this document, but the operation is not specified:

https://www.tfd.chalmers.se/~hani/ku...ukurReport.pdf

Many thanks for considering my request.

SurajGangani · December 28, 2023, 11:34

Hey Valero,

I know you posted this a long time ago, but I am also doing similar stuff and want to know whether you succeeded or not?

Kind regards,
Suraj

Miguel M Valero · January 3, 2024, 08:58

Hello Suraj,

Honestly, I cannot provide a good answer to you. When I wrote the post, I was using version 8 of OpenFOAM. Then, I moved to version 9, and the fvOptions library no longer exists, but it has been split into two other libraries: fvConstraints and fvModels.

In my particular test case (I only used the linear relationship with the velocity), both methods seem to agree when coefficients constituting D are up to 1e5. When I try larger orders of magnitude, OpenFOAM appears to give good results whereas implementing the forcing term directly inside the solver gradually diverges the simulation. So, I think that OpenFOAM employs some diffusive schemes inside the fvOptions library that might stabilize the results.

I know it is not a very good reply, but unfortunately I cannot tell you further about this.

March 21, 2023, 07:27	Darcy-Forchheimer OpenFOAM	#1
MiguelMValero New Member Miguel Martínez Valero Join Date: Mar 2023 Posts: 1 Rep Power: 0	Hello everybody, I have a concern about how the Darcy law is implemented in OpenFOAM. More specifically, I don't really understand the meaning of the vectors "d" and "f". In my case, I want: d [0 -2 0 0 0 0 0] (2.71508e6 70891 69.2805) f [0 -1 0 0 0 0 0] (0 0 0) where I thought that "d" and "f" were "D" and "F", respectively, in the Darcy formula: $S_i = -\nu (D_{ij}+\frac{1}{2}\rho \|u_{kk}\| F_{ij})u_i$ . However, when I tried to manually implement Darcy-Forchheimer in OpenFOAM by removing the file "porosityProperties" file and changing the incompressible solver "pimpleFoam" with the following expression in those cells affected by the porosity condition, I'm not getting the same results. tmp<fvVectorMatrix> tUEqn ( fvm::ddt(U) + fvm::div(phi, U) + MRF.DDt(U) + turbulence->divDegSigma(U) + nucmptMultiply(D, U) == fvOptions(U) ); I appreciate any help you can give me. I think I misunderstood how the multiplication of D and U is performed. I already checked on this document, but the operation is not specified: https://www.tfd.chalmers.se/~hani/ku...ukurReport.pdf Many thanks for considering my request. Last edited by MiguelMValero; March 21, 2023 at 07:29. Reason: LaTex misunderstanding*

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Frequently Asked Questions about Installing OpenFOAM	wyldckat	OpenFOAM Installation	3	November 14, 2023 11:58
Map of the OpenFOAM Forum - Understanding where to post your questions!	wyldckat	OpenFOAM	10	September 2, 2021 05:29
OpenFOAM Training Jan-Jul 2017, Virtual, London, Houston, Berlin	CFDFoundation	OpenFOAM Announcements from Other Sources	0	January 4, 2017 06:15
UNIGE February 13th-17th - 2107. OpenFOAM advaced training days	joegi.geo	OpenFOAM Announcements from Other Sources	0	October 1, 2016 19:20
OpenFOAM Foundation releases OpenFOAM 2.2.2	opencfd	OpenFOAM Announcements from ESI-OpenCFD	0	October 14, 2013 07:18

December 28, 2023, 11:34		#2
SurajGangani New Member Suraj Gangani Join Date: Feb 2021 Location: Germany Posts: 2 Rep Power: 0	Hey Valero, I know you posted this a long time ago, but I am also doing similar stuff and want to know whether you succeeded or not? Kind regards, Suraj

January 3, 2024, 08:58		#3
Miguel M Valero New Member Miguel Martínez Valero Join Date: Dec 2020 Posts: 1 Rep Power: 0	Hello Suraj, Honestly, I cannot provide a good answer to you. When I wrote the post, I was using version 8 of OpenFOAM. Then, I moved to version 9, and the fvOptions library no longer exists, but it has been split into two other libraries: fvConstraints and fvModels. In my particular test case (I only used the linear relationship with the velocity), both methods seem to agree when coefficients constituting D are up to 1e5. When I try larger orders of magnitude, OpenFOAM appears to give good results whereas implementing the forcing term directly inside the solver gradually diverges the simulation. So, I think that OpenFOAM employs some diffusive schemes inside the fvOptions library that might stabilize the results. I know it is not a very good reply, but unfortunately I cannot tell you further about this.