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

Axisymmetric cases Test with blockmesh and icoFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 11, 2007, 16:39
Default Hello, I try to do CFD with
  #1
frackowi
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
(
);
// ************************************************** *********************** //
  Reply With Quote

Old   July 11, 2007, 17:47
Default See http://www.cfd-online.com/
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
See http://www.cfd-online.com/cgi-bin/Op...4944#POST14944
mattijs is offline   Reply With Quote

Old   July 11, 2007, 19:02
Default Hi Mattijs Thanks a lot for y
  #3
frackowi
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
  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
IcoFoam Not Installed and blockMesh command not found cesarbz OpenFOAM Installation 6 July 4, 2008 04:44
Test Cases: 2D or 3D? jasond CFD-Wiki 6 September 13, 2007 05: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


All times are GMT -4. The time now is 19:41.