How the boundary conditions are enforced in Openfoam

 Register Blogs Members List Search Today's Posts Mark Forums Read

 March 26, 2013, 11:28 How the boundary conditions are enforced in Openfoam #1 Senior Member   Join Date: Nov 2012 Posts: 168 Rep Power: 6 Sponsored Links 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, 10:58 #2 Senior Member     Santiago Marquez Damian Join Date: Aug 2009 Location: Santa Fe, Santa Fe, Argentina Posts: 430 Rep Power: 17 Hi, the boundary conditions are enforced within the solve method, prior to call the sparse solvers, check the fvMatrixSolve.C file. Regards. wyldckat and hz283 like this. __________________ 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 Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Pascal_doran OpenFOAM Programming & Development 16 August 25, 2015 05:35 newOFuser OpenFOAM 1 January 10, 2013 09:08 L1011 OpenFOAM 5 December 13, 2012 09:17 ebah6 OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 September 16, 2012 16:57 lyang Main CFD Forum 0 September 19, 1999 18:29