|
[Sponsors] |
How the boundary conditions are enforced in Openfoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 26, 2013, 11:28 |
How the boundary conditions are enforced in Openfoam
|
#1 |
Senior Member
Join Date: Nov 2012
Posts: 171
Rep Power: 13 |
Hi All,
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 00002 ( 00003 fvm::ddt(rho, U) 00004 + fvm::div(phi, U) 00005 + turbulence->divDevRhoReff(U) 00006 == 00007 parcels.SU(U) 00008 + fvOptions(rho, U) 00009 ); 00010 00011 UEqn.relax(); 00012 00013 fvOptions.constrain(UEqn); 00014 00015 if (pimple.momentumPredictor()) 00016 { 00017 solve 00018 ( 00019 UEqn 00020 == 00021 fvc::reconstruct 00022 ( 00023 ( 00024 - ghf*fvc::snGrad(rho) 00025 - fvc::snGrad(p_rgh) 00026 )*mesh.magSf() 00027 ) 00028 ); 00029 00030 fvOptions.correct(U); 00031 K = 0.5*magSqr(U); 00032 } 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? Thanks bestH |
|
April 1, 2013, 11:58 |
|
#2 |
Senior Member
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 452
Rep Power: 23 |
Hi, the boundary conditions are enforced within the solve method, prior to call the sparse solvers, check the fvMatrixSolve.C file.
Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC) - CONICET/UNL Tel: 54-342-4511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe - Argentina. http://www.cimec.org.ar |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
non reflective boundary conditions for incompresible flow | Pascal_doran | OpenFOAM Programming & Development | 16 | August 25, 2015 06:35 |
inlet boundary conditions | newOFuser | OpenFOAM | 1 | January 10, 2013 09:08 |
Boundary conditions for 2D Navie-Stockes simulation | L1011 | OpenFOAM | 5 | December 13, 2012 09:17 |
[Commercial meshers] exporting boundary conditions with pointwise -> openFoam | ebah6 | OpenFOAM Meshing & Mesh Conversion | 2 | September 16, 2012 17:57 |
A problem about setting boundary conditions | lyang | Main CFD Forum | 0 | September 19, 1999 19:29 |