CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Attempt to implement the Chorin Projection Method

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

LinkBack Thread Tools Display Modes
Old   October 26, 2013, 06:46
Default Attempt to implement the Chorin Projection Method
New Member
Martin Goddard
Join Date: Mar 2013
Location: Melbourne, AU
Posts: 22
Rep Power: 4
McFly is on a distinguished road

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.



    phi = fvc::interpolate(U) & mesh.Sf();

    dU = (-fvc::div(phi)*U + nu*fvc::laplacian(U))*runTime.deltaT();

    U = U + dU;


    phi = fvc::interpolate(U) & mesh.Sf();

    solve(fvm::laplacian(p) == fvc::div(phi)/runTime.deltaT());

    U = U - runTime.deltaT()*fvc::grad(p);

McFly is offline   Reply With Quote


chorin, openfoam, projection method

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
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 12:19
projection method ilker basol Main CFD Forum 8 January 28, 2003 15:51
How to implement LU-SGS method Y.M. Wu Main CFD Forum 3 April 28, 1999 17:59

All times are GMT -4. The time now is 14:34.