|
[Sponsors] |
Axisymmetric cases Test with blockmesh and icoFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 11, 2007, 17:39 |
Hello,
I try to do CFD with
|
#1 |
Guest
Posts: n/a
|
Hello,
I try to do CFD with axisymmetric meshes. My test case is very simple, The domain is pentahedric (inlet, upperwall, outlet, the axis and the 2 wedges.) Blockmesh works and generate the mesh (See the script at the end) However, when I run icoFoam, I get the following errormessage, whereas I don't get it if the same mesh is 2D (inlet instead of wedges, and no axis) Starting time loop Time = 2e-08 Courant Number mean: 0 max: 0 #0 Foam::error::printStack(Foam:stream&) #1 Foam::sigFpe::sigFpeHandler(int) #2 Uninterpreted: [0xffffe420] #3 Foam::polyMesh::calcDirections() const #4 Foam::polyMesh::directions() const #5 Foam::fvMatrix<foam::vector<double> >::solve(Foam::Istream&) #6 Foam::lduMatrix::solverPerformance Foam::solve<foam::vector<double> >(Foam::tmp<foam::fvmatrix<foam::vector<double> > > const&) #7 main #8 __libc_start_main #9 __gxx_personality_v0 at ../sysdeps/i386/elf/start.S:122 Exception en point flottant I try to change the axis position, but it changes nothing. What can I do? Thank you in advance Bruno convertToMeters 0.001; vertices ( (0 10 -0.436) (0 0 -1) (50 0 -1) (50 10 -0.436) (0 10 0.436) (50 10 0.436) ); blocks ( hex (0 1 2 3 4 1 2 5) (10 50 1) simpleGrading (1 1 1) ); edges ( ); patches ( empty axis ( (1 2 2 1) ) wedge back ( (0 1 2 3) ) wedge front ( (4 1 2 5) ) patch inlet ( (0 1 1 4) ) wall top_flange ( (0 3 5 4) ) patch outlet ( (3 2 2 5) ) ); mergePatchPairs ( ); // ************************************************** *********************** // |
|
July 11, 2007, 18:47 |
See http://www.cfd-online.com/
|
#2 |
Senior Member
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 26 |
||
July 11, 2007, 20:02 |
Hi Mattijs
Thanks a lot for y
|
#3 |
Guest
Posts: n/a
|
Hi Mattijs
Thanks a lot for your advice, the test seems to work now and provide good results. However, is the mesh still considered as axisymmetric from the solver, without an axis definition in the boundary list? How does the solver know, where is the axis and where is the non varying direction? Thank you in advance Bruno |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
IcoFoam Not Installed and blockMesh command not found | cesarbz | OpenFOAM Installation | 6 | July 4, 2008 05:44 |
Test Cases: 2D or 3D? | jasond | CFD-Wiki | 6 | September 13, 2007 06:58 |
FSI test cases | Kirikou | Main CFD Forum | 8 | November 6, 2006 10:49 |
test cases | Maciej Matyka | Main CFD Forum | 3 | November 24, 2004 09:27 |
NACA test cases | Ben | Main CFD Forum | 1 | February 28, 2000 15:45 |