CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Attempt to implement the Chorin Projection Method

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 26, 2013, 07:46
Default Attempt to implement the Chorin Projection Method
  #1
New Member
 
Martin Goddard
Join Date: Mar 2013
Location: Melbourne, AU
Posts: 22
Rep Power: 13
McFly is on a distinguished road
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();
McFly is offline   Reply With Quote

Reply

Tags
chorin, openfoam, projection method


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 Off
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 10:50
Chorin's Projection Method for Finite Volume Scott2 Main CFD Forum 1 August 16, 2010 21: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 18:59


All times are GMT -4. The time now is 02:40.