CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

help me for a DNS simulation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 13, 2012, 13:18
Exclamation help me for a DNS simulation
  #1
New Member
 
omar_beg
Join Date: Mar 2012
Posts: 17
Rep Power: 6
dima13 is on a distinguished road
I have a problem with DNS, I want a guide, the DNSFoam
my problem is a water injector, with a torpedo
Attached Images
File Type: jpg anim.0013.jpg (16.7 KB, 23 views)
dima13 is offline   Reply With Quote

Old   June 14, 2012, 06:17
Default
  #2
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 261
Rep Power: 10
anishtain4 is on a distinguished road
Can you tell us what is your problem?
anishtain4 is offline   Reply With Quote

Old   June 14, 2012, 16:36
Default
  #3
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 464
Rep Power: 10
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
I think we have to guess, otherwise would be too easy!
lovecraft22 is offline   Reply With Quote

Old   June 14, 2012, 17:29
Default
  #4
New Member
 
omar_beg
Join Date: Mar 2012
Posts: 17
Rep Power: 6
dima13 is on a distinguished road
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:rintStack(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
dima13 is offline   Reply With Quote

Old   June 14, 2012, 17:32
Default
  #5
New Member
 
omar_beg
Join Date: Mar 2012
Posts: 17
Rep Power: 6
dima13 is on a distinguished road
I have no idea, for dnsFoam, when I my problem, does not calculate the solever
Attached Images
File Type: jpg hh.0000.jpg (52.0 KB, 20 views)
dima13 is offline   Reply With Quote

Old   June 15, 2012, 03:13
Default
  #6
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 14
Bernhard is on a distinguished road
I have no experience with dnsFoam, but maybe there is something wrong with your mesh? Can you share the output of 'checkMesh'?
Bernhard is offline   Reply With Quote

Old   June 15, 2012, 04:23
Default
  #7
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 261
Rep Power: 10
anishtain4 is on a distinguished road
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
anishtain4 is offline   Reply With Quote

Old   June 15, 2012, 19:09
Default
  #8
New Member
 
omar_beg
Join Date: Mar 2012
Posts: 17
Rep Power: 6
dima13 is on a distinguished road
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$
dima13 is offline   Reply With Quote

Old   June 15, 2012, 19:10
Default
  #9
New Member
 
omar_beg
Join Date: Mar 2012
Posts: 17
Rep Power: 6
dima13 is on a distinguished road
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$
dima13 is offline   Reply With Quote

Old   June 15, 2012, 19:29
Default
  #10
New Member
 
omar_beg
Join Date: Mar 2012
Posts: 17
Rep Power: 6
dima13 is on a distinguished road
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:rintStack(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$
dima13 is offline   Reply With Quote

Old   June 15, 2012, 19:30
Default
  #11
New Member
 
omar_beg
Join Date: Mar 2012
Posts: 17
Rep Power: 6
dima13 is on a distinguished road
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:rintStack(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$
dima13 is offline   Reply With Quote

Old   June 15, 2012, 23:21
Default
  #12
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Location: NB, Canada
Posts: 261
Rep Power: 10
anishtain4 is on a distinguished road
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
anishtain4 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Simulation of a complex wing in solidworks flow simulation niels1900 FloEFD, FloWorks & FloTHERM 6 April 20, 2011 10:44
GUI crash and simulation engine still running RPJones FLOW-3D 2 November 9, 2010 09:18
finding noise of flow by dns simulation m2montazari OpenFOAM Running, Solving & CFD 0 October 22, 2010 11:54
velocity profile export from a simulation onto another sudhirlv STAR-CCM+ 1 September 12, 2010 18:57
FSI TWO-WAY SIMULATION Smagmon CFX 1 March 6, 2009 14:24


All times are GMT -4. The time now is 13:31.