Hello,
I try to do CFD with
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) #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 ( ); // ************************************************** *********************** // |
See http://www.cfd-online.com/
|
Hi Mattijs
Thanks a lot for y
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 |
All times are GMT -4. The time now is 21:47. |