Running MixerVesselAMI2D

This is basically a note to myself. At the the first try of running MixerVesselAmi2D, had some weird error messages from openfoam.

FOAM FATAL ERROR: bool IPstream::init(int& argc, 
char**& argv) : attempt to run parallel on 1 processor

When I was expecting tutorials to run smoothly as before but this was not the case. And as always, one finds the solution, applies it and then forgets about it I decided to make a small post. This time it was VERY easy, thanks to this other blogger that had the same problem:

http://diffusionht.blogspot.de/2012/08/foamparallel-ubuntu-mpirun-error.html

The changes worked like a charm and the tutorial run without problem. It was only necessary to point mpirun to use mpirun.openmpi when running in parallel.  I listed the steps here but one can also go to the original blog.

The file to modify is the RunFunctions one:

$WM_PROJECT_DIR/bin/tools/RunFunctions
runParallel()
{
    APP_RUN=$1
    APP_NAME=${1##*/}
    shift

    if [ -f log.$APP_NAME ]
    then
        echo "$APP_NAME already run on $PWD: remove log file to re-run"
    else
        nProcs=$1
        shift
        echo "Running $APP_RUN in parallel on $PWD using $nProcs processes"

        #if [ "$WM_SCHEDULER" ]
        #then
        #    echo "$PWD: $WM_SCHEDULER -np $nProcs" 1>&2
        #    $WM_SCHEDULER -np $nProcs "( mpirun -np $nProcs $APP_RUN -parallel $* < /dev/null > log.$APP_NAME 2>&1 )"
        #else
            ( mpirun -np $nProcs $APP_RUN -parallel $* < /dev/null > log.$APP_NAME 2>&1 )
        #fi
    fi
}

And change mpirun to mpirun.openmpi:

runParallel()
{
    APP_RUN=$1
    APP_NAME=${1##*/}
    shift

    if [ -f log.$APP_NAME ]
    then
        echo "$APP_NAME already run on $PWD: remove log file to re-run"
    else
        nProcs=$1
        shift
        echo "Running $APP_RUN in parallel on $PWD using $nProcs processes"

        #if [ "$WM_SCHEDULER" ]
        #then
        #    echo "$PWD: $WM_SCHEDULER -np $nProcs" 1>&2
        #    $WM_SCHEDULER -np $nProcs "( mpirun -np $nProcs $APP_RUN -parallel $* < /dev/null > log.$APP_NAME 2>&1 )"
        #else
            ( mpirun.openmpi -np $nProcs $APP_RUN -parallel $* < /dev/null > log.$APP_NAME 2>&1 )
        #fi
    fi
}

Velocity Solution for OpenFOAM MixerAMI2D in parallel

Velocity Solution for OpenFOAM MixerAMI2D in parallel

Posted in CFD, OpenFOAM | Tagged , , , , | Leave a comment

Concentric Inclined Disk 3D plots

 

Posted in 3D | Leave a comment

HE-3 30 ml scale with TIP1

These plot are from a 30 ml vessel stirred by a HE-3 like impeller at 250 RPM.

Posted in 3D | Leave a comment

HE-3 134 liter Scale 3D particle plot

These plots are from a vessel being stirred by a HE-3 impeller. Particles were seeded at position number 1, TIP1, which is in the upper section of the vessel near the surface.

Posted in 3D | Leave a comment

Interactive 3D webGL plot

I saw the RGL WebGL interactive plot in R-bloggers and I am linking to it if someone is interested in it.

http://geolabs.wordpress.com/2012/04/24/r-interactive-3d-webgl-plot-of-time-space-cube-with-rgl/

The interactive plots of my works are created using these tools. I recommend you to subscribe to R-bloggers if you dig R.

http://www.r-bloggers.com/

 

 

 

 

Posted in 3D | Leave a comment

First 3D representation with R

Well this is one of the first 3D images I was able to get from R.  It was created with random data and bin number was set to 10 on x,y and z coordinates. Resolution was too coarse to  obtain a decent “mixing” time.

First3Drepresentation

First 3D representation of particle count in the studied vessel. Color indicates particle number, legend was still not created.

 

 

This image was created with scatterplot3D library, with symbols set to cube, Color scale denotes particle count in that “cube”. Legend was still not created. I’ll post the script later.

Posted in 3D | 4 Comments

Eccentric Inclined Disk Impeller E=0.42

This gallery contains 6 photos.

Poincaré maps created by an inclined disk impeller with eccentricity value (E) 0.42 at 250 RPM from TIP2.  

More Galleries | Leave a comment

Using R to visualize models in 3D

This was not available in the print document as it is an interactive resource. Using R it is possible to show the node position in 3D from one of the models used for the CFD simulations interactively. Each point represents 1 node.  One can Zoom in / out and rotate.

Notice the more densely populated volume around corresponding to the volume of fluid in the rotating volume element in the simulation.  A more refined mesh was created there to resolve the impeller area more clearly.

Please to be able to visualize correctly use Firefox ™ or Chrome ™

 Concentric Inclined Disk Impeller

Other geometries will be posted further on.

Posted in 3D | Leave a comment

3D particle counting and tracing with R

This gallery contains 6 photos.

3D particle plotting and counting with R, scripts will follow.

More Galleries | Leave a comment

Poincaré Maps!

This gallery contains 6 photos.

More Galleries | Leave a comment