# Attempt to implement the Chorin Projection Method

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

 October 26, 2013, 06:46 Attempt to implement the Chorin Projection Method #1 New Member   Martin Goddard Join Date: Mar 2013 Location: Melbourne, AU Posts: 22 Rep Power: 5 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 Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Petrov Main CFD Forum 2 June 10, 2011 09:50 Scott2 Main CFD Forum 1 August 16, 2010 20:24 misabel Main CFD Forum 3 March 11, 2010 12:19 ilker basol Main CFD Forum 8 January 28, 2003 15:51 Y.M. Wu Main CFD Forum 3 April 28, 1999 17:59

All times are GMT -4. The time now is 12:07.