|
[Sponsors] |
March 22, 2017, 19:15 |
How to update phi in Fractional Step Method
|
#1 |
Member
Join Date: Jul 2012
Posts: 67
Rep Power: 13 |
Hi Foamers,
I intend to realise a semi-implicit fractional step in OpenFOAM but got stuck at the final step where the flux phi is updated. Please see the code here: Code:
// Create Poisson Matrix fvScalarMatrix pEqn ( fvm::laplacian(p) ); pEqn.setReference(pRefCell, pRefValue); // PressureCorrection.H solve(pEqn==fvc::div(U)/runTime.deltaT()); phi -= pEqn.flux();// This is where I have problem. U -= runTime.deltaT()*fvc::grad(p); U.correctBoundaryConditions(); Code:
phi -= pEqn.flux(); Code:
phi += pEqn.flux(); Thanks for your help, Carlen |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
courant number increases to rather large values | 6863523 | OpenFOAM Running, Solving & CFD | 22 | July 5, 2023 23:48 |
Fractional Step Method and SIMPLE | Ha Lim Choi | Main CFD Forum | 14 | June 14, 2017 11:17 |
Convergence issue with Fractional Step Method using RK2 method for time integration | mihirmakwana6 | Main CFD Forum | 23 | June 6, 2015 06:58 |
Comparing between the Fractional step method and the SIMPLE method | ghlee | Main CFD Forum | 1 | April 10, 2012 16:59 |
IcoFoam parallel woes | msrinath80 | OpenFOAM Running, Solving & CFD | 9 | July 22, 2007 02:58 |