|
[Sponsors] |
Attempt to implement the Chorin Projection Method |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Martin Goddard
Join Date: Mar 2013
Location: Melbourne, AU
Posts: 22
Rep Power: 13 ![]() |
Hello,
As above, I'm attempting to implement the Chorin Projection Method in OpenFOAM. This is my first time playing with solvers, I'm taking most of my information from the paper, "On the Implementation of Low-Dissipative, Fractional Step Methods for Transient Flows Using OpenFOAM" V Vuorinen, J.-P Keskinen, C Duwig, O Kaario, M Larmi, B J Boersma. I have code that compiles but doesn't work. It blows up. Although I can see the results are clearly wrong in paraFoam, the equations seem sound based on the above paper and other reading. Can you help me find what I'm missing? I've been searching for days and haven't found a complete implementation of the code to look at. Cheers, Marty Code:
phi = fvc::interpolate(U) & mesh.Sf(); dU = (-fvc::div(phi)*U + nu*fvc::laplacian(U))*runTime.deltaT(); U = U + dU; U.correctBoundaryConditions(); phi = fvc::interpolate(U) & mesh.Sf(); solve(fvm::laplacian(p) == fvc::div(phi)/runTime.deltaT()); U = U - runTime.deltaT()*fvc::grad(p); U.correctBoundaryConditions(); |
|
![]() |
![]() |
![]() |
Tags |
chorin, openfoam, projection method |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Projection Methods | Petrov | Main CFD Forum | 2 | June 10, 2011 09:50 |
Chorin's Projection Method for Finite Volume | Scott2 | Main CFD Forum | 1 | August 16, 2010 20:24 |
Inner non linear iterations and Newton method | misabel | Main CFD Forum | 3 | March 11, 2010 11:19 |
projection method | ilker basol | Main CFD Forum | 8 | January 28, 2003 14:51 |
How to implement LU-SGS method | Y.M. Wu | Main CFD Forum | 3 | April 28, 1999 17:59 |