CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   movingConeTopo modifications (https://www.cfd-online.com/Forums/openfoam-solving/143993-movingconetopo-modifications.html)

MarinaSS88 November 6, 2014 15:20

movingConeTopo modifications
 
Hi!

I'm trying to simulate a box moving on a fluid addig/removing mesh layers. But I am having several problems, as I am very beginner with openFoam.

I have started with the following tutorial:
www.cfd-online.com/OpenFOAM_Discus/messages/1/movingConeTopo-5336.unk

I have noticed that the parameters in blockMeshDict are defined in mm. I would like to run a similar geometry but changing the parameters into meters and modifying the dimensions of the box that I want to move. The modifications I have performed are:
- Change convertToMeters in blockMeshDict into 1; and point coordinates (not a very big change, for example instead of -7.5, I have writen -10, but this time in meters).
- Change the mesh size so all the mesh cells have 1m x 1m x 1m dimensions.
- In dinamicMeshDict: have changed leftEdge, leftObstacleEdge and rightObstacleEdge to values corresponding to my geometry; and min and max thicknesses to 0.5 and 1.5 meters at both sides.

When I try to run this new case, openFoam keeps moving the box defined in the tutorial (where: leftObstacleEdge= -0.007 and rightObstacleEdge= -0.0035) and not the box that I have defined after modifying dinamicMeshDict. It seems that the program maintains the old values without refreshing them. Maybe I am missing to change this values in any other file when it is necessary? I thought that these parameters were only defined here? or maybe my computer has gone mad?? I always delete all the files generated in a case-running before running a new case, so I donīt know why new values are not considered by the program :confused:

Moreover, depending on the new points coordinates, sometimes when I run moveDynamicMesh I get the error:

Face extrusion zone contains no faces. Plese check your mesh definition.

Is there any rule that you have to follow to define dinamicMeshDict correctly depending on the geometry and the background mesh?

I would strongly appreciate if anybody could help me!

Many thanks in advance!!!!

Marina

alexB March 26, 2015 05:27

Hi Marina,

tested movingConeTopo by myself, but also couldn't perform any great changes in the mesh.
I also checked the file movingConeTopoFvMesh.C in my folder /opt/openfoam230/src/topoChangerFvMesh/movingConeTopoFvMesh and noticed fixed values in the code.

Code:

            fc[faceI].x() > -0.003501
        && fc[faceI].x() < -0.003499

Code:

            fc[faceI].x() > -0.00701
        && fc[faceI].x() < -0.00699

Even when I changed them to the values in my blockmeshdict and ran the Allwmake script* in /opt/openfoam230/src/ , I still get the error "Face extrusion zone contains no faces..." .

Did you made any further progress by your problem?



*after spending several hours during the first run I wrote a reduced wmake script.
Code:

#!/bin/sh
cd ${0%/*} || exit 1    # run from this directory
makeType=${1:-libso}

wmakeCheckPwd "$WM_PROJECT_DIR/src" || {
    echo "Error: Current directory is not \$WM_PROJECT_DIR/src"
    echo "    The environment variables are inconsistent with the installation."
    echo "    Check the OpenFOAM entries in your dot-files and source them."
    exit 1
}

[ -n "$FOAM_EXT_LIBBIN" ] || {
    echo "Error: FOAM_EXT_LIBBIN not set"
    echo "    Check the OpenFOAM entries in your dot-files and source them."
    exit 1
}

set -x

# update OpenFOAM version strings if required
wmakePrintBuild -check || /bin/rm -f OpenFOAM/Make/*/global.? 2>/dev/null

wmakeLnInclude OpenFOAM
wmakeLnInclude OSspecific/${WM_OSTYPE:-POSIX}
Pstream/Allwmake $*

OSspecific/${WM_OSTYPE:-POSIX}/Allwmake $*

wmake $makeType topoChangerFvMesh

# ----------------------------------------------------------------- end-of-file



regards
Alex


All times are GMT -4. The time now is 06:08.