|
[Sponsors] |
May 7, 2020, 11:48 |
Porous Media with rhoCentralFoam
|
#1 |
Senior Member
Join Date: Dec 2019
Posts: 215
Rep Power: 7 |
Hi,
I have noticed that rhoCentralFoam has no porosity treatment unlike sonicFoam or rhoPimpleFoam. Since I have the best results with rhoCentralFoam I would like to add the porosity treatment. I want to use the Darcy-Forchheimer model explicitly by using fvOptions, which works well for other solvers except rhoCentralFoam. I have noticed, that in rhoPimpleFoam UEqn.H fvOptions is mentioned: Code:
// Solve the Momentum equation MRF.correctBoundaryVelocity(U); tmp<fvVectorMatrix> tUEqn ( fvm::ddt(rho, U) + fvm::div(phi, U) + MRF.DDt(rho, U) + turbulence->divDevRhoReff(U) == fvOptions(rho, U) ); fvVectorMatrix& UEqn = tUEqn.ref(); UEqn.relax(); fvOptions.constrain(UEqn); if (pimple.momentumPredictor()) { solve(UEqn == -fvc::grad(p)); fvOptions.correct(U); K = 0.5*magSqr(U); } Code:
// --- Solve momentum solve(fvm::ddt(rhoU) + fvc::div(phiUp)); U.ref() = rhoU() /rho(); U.correctBoundaryConditions(); rhoU.boundaryFieldRef() == rho.boundaryField()*U.boundaryField(); if (!inviscid) { solve ( fvm::ddt(rho, U) - fvc::ddt(rho, U) - fvm::laplacian(muEff, U) - fvc::div(tauMC) ); rhoU = rho*U; } Since my flow is viscous, I should concentrate on the if statement. Replacing it by Code:
if (!inviscid) { solve ( fvm::ddt(rho, U) - fvc::ddt(rho, U) - fvm::laplacian(muEff, U) - fvc::div(tauMC) == fvOptions(rho, U) ); rhoU = rho*U; } I would be glad to hear your advice! Regards, shock77 |
|
May 8, 2020, 08:00 |
|
#2 |
Senior Member
Join Date: Dec 2019
Posts: 215
Rep Power: 7 |
So what I want to to is to add the additional sink term of the darcy model or better the darcy forchheimer model. Unfortunately I dont know how to do it.
https://openfoamwiki.net/index.php/DarcyForchheimer I hope someone has more experience than me. |
|
May 9, 2020, 18:54 |
|
#3 |
Senior Member
Join Date: Dec 2019
Posts: 215
Rep Power: 7 |
I was mistaken, it works well.
For anyone with the same issue: Just add the additional sink term with fvOptions(rho, U) like I did. Of course you have to include fvOption.H in rhoCentralFoam.C or your new solver! Compile it and use it like the other solvers. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Thermal non-equilibrium porous media model with conjugate heat transfer | Hexahedron | FLUENT | 9 | February 22, 2023 02:55 |
use of porous media with sonicFoam or rhoCentralFoam | Calmly | OpenFOAM Running, Solving & CFD | 6 | May 7, 2020 11:12 |
Porous media setup issues in Fluent | Bernard Van | FLUENT | 29 | January 26, 2017 04:09 |
How to model granular flow through porous media | Axius | FLUENT | 2 | August 7, 2014 10:34 |
porous media: Fluent or Star-CD? | Igor | Main CFD Forum | 0 | December 5, 2002 15:16 |