# Which procedure to solve Ueqn without momentum predictor.

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

 July 2, 2013, 21:04 Which procedure to solve Ueqn without momentum predictor. #1 Senior Member   Dongyue Li Join Date: Jun 2012 Location: Torino, Italy Posts: 752 Rep Power: 10 Sponsored Links Hi, This is icoFoam's code,: Code: fvVectorMatrix UEqn ( fvm::ddt(U) + fvm::div(phi, U) - fvm::laplacian(nu, U) ); solve(UEqn == -fvc::grad(p));//this term is only executed when the momentumPredictor flag is set to true // --- PISO loop for (int corr=0; corr

 July 2, 2013, 23:37 #2 Senior Member     Santiago Marquez Damian Join Date: Aug 2009 Location: Santa Fe, Santa Fe, Argentina Posts: 430 Rep Power: 17 Hi, in the PISO method, U is predicted and then corrected in the pressure-velocity coupling loop, so that, if you don't predict the velocity, the velocity used to start the coupling loop is a derivative of the velocity in the last time step as is done in interFoam and other solvers. In the code you've cited the matrix for the momentum equation is assembled so that you have the A and H methods available for the use in the PISO loop, the line: HbyA = rAU*UEqn.H(); calculates a first approximation for U using the H operator applied to the velocity of the last-timestep or which was predicted when you solved the UEqn (if you did it) and the A operator. If you need it, please check the definition of these operators in chapter 2.7 of my Ph.D. thesis The link to my PhD thesis is provided in my user page at openfoamwiki.net: http://openfoamwiki.net/index.php/User:Santiagomarquezd and Jasak's Ph.D. and other's thesis. Regards. Tushar@cfd, sharonyue and Nicole 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 Last edited by wyldckat; May 24, 2016 at 17:39. Reason: rectified link according to request

 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 Liam CFX 28 July 16, 2013 08:24 matteoL OpenFOAM Running, Solving & CFD 0 November 18, 2009 07:58 sek OpenFOAM Running, Solving & CFD 5 March 6, 2008 17:27 ABS FLUENT 0 October 4, 2005 17:58 west_wing FLUENT 0 August 25, 2003 10:00