# 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: Beijing, China Posts: 841 Rep Power: 18 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: 452 Rep Power: 23 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, Nicole and 2 others 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