help me for a DNS simulation
 

I have a problem with DNS, I want a guide, the DNSFoam
my problem is a water injector, with a torpedo

Can you tell us what is your problem?

I think we have to guess, otherwise would be too easy!

omar@omar-laptop:~/Bureau/essaiDNS$ boxTurb
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.0.1-51f1de99a4bc
Exec : boxTurb
Date : Jun 14 2012
Time : 22:28:25
Host : omar-laptop
PID : 2214
Case : /home/omar/Bureau/essaiDNS
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field U

Reading boxTurbDict



--> FOAM FATAL ERROR:
calculated number of cells is incorrect

From function Kmesh::Kmesh(const fvMesh& mesh)
in file Kmesh/Kmesh.C at line 84.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::Kmesh::Kmesh(Foam::fvMesh const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/librandomProcesses.so"
#3
in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/boxTurb"
#4 __libc_start_main in "/lib/libc.so.6"
#5
in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/boxTurb"
Abandon

I have no idea, for dnsFoam, when I my problem, does not calculate the solever

I have no experience with dnsFoam, but maybe there is something wrong with your mesh? Can you share the output of 'checkMesh'?

you are trying to use boxTurb with your case? As you can see from its name that is only for a box, it makes initial condition for a homogeneous turbulence box, in your case you should only think of inlet condition, turbulence inside will be made itself.
btw DNSFoam adds forced turbulence during the simulation which you don't need either!!! just take another solver and with a fine mesh you are doing dns

omar@omar-laptop:~/Bureau/essaiDNS$ blockMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.0.1-51f1de99a4bc
Exec : blockMesh
Date : Jun 16 2012
Time : 00:09:02
Host : omar-laptop
PID : 2386
Case : /home/omar/Bureau/essaiDNS
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Creating block mesh from
"/home/omar/Bureau/essaiDNS/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Creating block mesh topology

Check topology

Basic statistics
Number of internal faces : 28
Number of boundary faces : 40
Number of defined boundary faces : 40
Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 0.01

Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (0 0 -0.05) (0.4 0.05 0.05)
nPoints: 137781
nCells: 128000
nFaces: 393600
nInternalFaces: 374400
----------------
Patches
----------------
patch 0 (start: 374400 size: 1600) name: inlet
patch 1 (start: 376000 size: 1600) name: outlet
patch 2 (start: 377600 size: 6400) name: upperWal
patch 3 (start: 384000 size: 3200) name: lowerWall
patch 4 (start: 387200 size: 3200) name: axeWall
patch 5 (start: 390400 size: 3200) name: front

End

omar@omar-laptop:~/Bureau/essaiDNS$

mar@omar-laptop:~/Bureau/essaiDNS$ checkMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.0.1-51f1de99a4bc
Exec : checkMesh
Date : Jun 16 2012
Time : 00:09:52
Host : omar-laptop
PID : 2565
Case : /home/omar/Bureau/essaiDNS
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 137781
faces: 393600
internal faces: 374400
cells: 128000
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 128000
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
inlet 1600 1701 ok (non-closed singly connected)
outlet 1600 1701 ok (non-closed singly connected)
upperWal 6400 6561 ok (non-closed singly connected)
lowerWall 3200 3321 ok (non-closed singly connected)
axeWall 3200 3321 ok (non-closed singly connected)
front 3200 3402 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (0 0 -0.05) (0.4 0.05 0.05)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-2.54409e-17 -6.19965e-15 1.35568e-17) OK.
Max cell openness = 2.04348e-15 OK.
Max aspect ratio = 118.807 OK.
Minumum face area = 1.07985e-08. Maximum face area = 3.56404e-05. Face area magnitudes OK.
Min volume = 6.23376e-11. Max volume = 3.84201e-08. Total volume = 0.00148204. Cell volumes OK.
Mesh non-orthogonality Max: 22.4534 average: 4.15421
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.00178 OK.

Mesh OK.

End

omar@omar-laptop:~/Bureau/essaiDNS$ ^C
omar@omar-laptop:~/Bureau/essaiDNS$

omar@omar-laptop:~/Bureau/essaiDNS$ boxTurb
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.0.1-51f1de99a4bc
Exec : boxTurb
Date : Jun 16 2012
Time : 00:29:02
Host : omar-laptop
PID : 2949
Case : /home/omar/Bureau/essaiDNS
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field U

Reading boxTurbDict



--> FOAM FATAL ERROR:
calculated number of cells is incorrect

From function Kmesh::Kmesh(const fvMesh& mesh)
in file Kmesh/Kmesh.C at line 84.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::Kmesh::Kmesh(Foam::fvMesh const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/librandomProcesses.so"
#3
in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/boxTurb"
#4 __libc_start_main in "/lib/libc.so.6"
#5
in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/boxTurb"
Abandon
omar@omar-laptop:~/Bureau/essaiDNS$

omar@omar-laptop:~/Bureau/essaiDNS$ dnsFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.0.1-51f1de99a4bc
Exec : dnsFoam
Date : Jun 16 2012
Time : 00:30:09
Host : omar-laptop
PID : 2960
Case : /home/omar/Bureau/essaiDNS
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh, no clear-out for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

Reading turbulenceProperties



--> FOAM FATAL ERROR:
calculated number of cells is incorrect

From function Kmesh::Kmesh(const fvMesh& mesh)
in file Kmesh/Kmesh.C at line 84.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::Kmesh::Kmesh(Foam::fvMesh const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/librandomProcesses.so"
#3
in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/dnsFoam"
#4 __libc_start_main in "/lib/libc.so.6"
#5
in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/dnsFoam"
Abandon
omar@omar-laptop:~/Bureau/essaiDNS$ ^C
omar@omar-laptop:~/Bureau/essaiDNS$

As you see it says calculating number of cells is incorrect, this is because it is using fft to energize the fluid in the fourier space which does not apply to your case. your grid number should be 2^n in any direction, even in that case you won't get a reasonable velocity field from this solver