How the boundary conditions are enforced in Openfoam
I have a question about how the boundary condition is used in the code. For example, when the momentum equation is solved, I did not find any statements to call the boundary conditions. The codes are follows:
00001 fvVectorMatrix UEqn
00003 fvm::ddt(rho, U)
00004 + fvm::div(phi, U)
00005 + turbulence->divDevRhoReff(U)
00008 + fvOptions(rho, U)
00015 if (pimple.momentumPredictor())
00024 - ghf*fvc::snGrad(rho)
00025 - fvc::snGrad(p_rgh)
00031 K = 0.5*magSqr(U);
It seems that from the source code for pressure and energy equations, the boundary treatment does not appear in the application source files. Does anybody know where the quantities at the boundary conditions are treated?
Hi, the boundary conditions are enforced within the solve method, prior to call the sparse solvers, check the fvMatrixSolve.C file.
|All times are GMT -4. The time now is 21:56.|