HowTo Add a Coriolis Force And Make a Tornado

From OpenFOAMWiki
Revision as of 14:37, 10 June 2016 by Shor-ty (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

1 Copy Some Source Code

if you do not yet have a solver directory in $WM_PROJECT_USER_DIR , create it:

mkdir $WM_PROJECT_USER_DIR/applications/solvers

Copy the original tutorial code for buoyantBoussinesqPimpleFoam into your solver directory, renaming it to myTornado:

cd heatTransfer
cp -r buoyantBoussinesqPimpleFoam  $WM_PROJECT_USER_DIR/applications/solvers/myTornado

With this code in its new location, we need to copy a file that was shared in the the previous location:

cp buoyantBoussinesqSimpleFoam/readTransportProperties.H  $WM_PROJECT_USER_DIR/applications/solvers/myTornado/.

Now let's move on to check that this copied code will compile in its new location, with its new name:

cd $WM_PROJECT_USER_DIR/applications/solvers
ls -l
cd myTornado

if you have this file, delete it:

rm buoyantBoussinesqPimpleFoam.dep
cd Make

If you have this directory, then delete it:

rm -rf linux64GccDPOpt

and edit the file named "files" to this:

EXE = $(FOAM_APPBIN)/myTornado

In addition you have to rename the source file:

 mv buoyantBoussinesqSimpleFoam.C myTornado.C

Though not required, you may want to delete this line in the file named "options":

  -I../buoyantBoussinesqSimpleFoam \
cd ..

You may get a few warnings. But hopefully you will see a new app myTornado (though without any ability yet to make a tornado) :


At this point, the app myTornado will be identical to buoyantBoussinesqPimpleFoam. Next we add the capability to make the tornado.

2 Add the Coriolis Force

Add a single line UEqn.H, so that the top of the file is:

    // Solve the momentum equation

    fvVectorMatrix UEqn
      + fvm::div(phi, U)
      + turbulence->divDevReff(U)
      + (2*myVector ^ U) // Coriolis force

Add one line to the end of readTransportProperties.H

dimensionedVector myVector(laminarTransport.lookup("myVector"));

Then recompile


3 Prepare the Run Directory

cd $FOAM_RUN/tutorials/heatTransfer
cp -r buoyantBoussinesqPimpleFoam/hotRoom $FOAM_RUN/myTornado

We consolidate some code that were separate in the original tutorial structure:

cp -r buoyantPimpleFoam/hotRoom/setHotRoom $FOAM_RUN/myTornado/.
cd $FOAM_RUN/myTornado

You should see something like this:

0/  Allclean*  Allrun*  constant/  setHotRoom/  system/

edit system/controlDict so that the application line reads (remove the original line rather than attempting to comment it out) :

application     myTornado;

Edit Allrun so that the compileApplication line is simply:

compileApplication setHotRoom

Now we need to specify myVector, which we will use as the rotation vector in the calculation of the Coriolis force:

cd constant

Edit transportProperties near the end to include this line:

myVector myVector [0 0 -1 0 0 0 0] (0 0.001 0);


cd ..

A quick way to proceed is to do:


After about 20 seconds of CPU time, the run should finish. Now you should be apple to view the results with paraFoam.

clip though center, show pressure and isosurface

Note the log files.

If you want to experience the details that were scripted by Allrun, do this instead:

cd setHotRoom

This should reveal a new app setHotRoom:

cd ..