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

Axisymmetric bodies wedgetype BCbs

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

Like Tree3Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   February 16, 2007, 04:32
Default Have you had a look at those f
  #21
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 17
mattijs is on a distinguished road
Have you had a look at those faces with zero area? (using e.g. foamToVTK) Are they at the axis?

Have a look at the collapseEdges utility. I think it will remove your 0 area face at the axis. (if they are quads)
mattijs is offline   Reply With Quote

Old   March 6, 2007, 15:18
Default Hello Mattijs: What is the
  #22
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 9
vatant is on a distinguished road
Hello Mattijs:
What is the fundamental difference betweeen an wedge type and cyclic type? If I consider an axi-symmetric grid (one cell thick) and impose wedge type condition, does OF compute in the rotational direction? For the same case, why doesnt a cyclic condition work?

Regards,

Vatant
vatant is offline   Reply With Quote

Old   March 16, 2007, 15:54
Default Hi all, I am new to OF so m
  #23
New Member
 
gianfilippo coppola
Join Date: Mar 2009
Posts: 14
Rep Power: 9
gfcoppola is on a distinguished road
Hi all,

I am new to OF so my question may be stupid, anyway, I am trying to simulate an axisymmetric opposed jet configuration, where each jet has a circular flange flat with the exit section. I am using icoFoam from FoamX and blockMesh to generate the mesh. Below is the blockMeshDict file:

arguments "/home/gcoppola/OpenFOAM/gcoppola-1.3/run/tutorials/icoFoam/counterflow" off;

convertToMeters 0.001;

vertices
(
(-14.5852 0 0)
(-0.297656 0 0)
(-0.297656 1.19063 -0.051984)
(-14.5852 1.19063 -0.051984)
(0.297656 0 0)
(0.297656 1.19063 -0.051984)
(14.5852 0 0)
(14.5852 1.19063 -0.051984)
(0.297656 100 -4.366094)
(-0.297656 100 -4.366094)
(-0.297656 1.19063 0.051984)
(-14.5852 1.19063 0.051984)
(0.297656 1.19063 0.051984)
(14.5852 1.19063 0.051984)
(0.297656 100 4.366094)
(-0.297656 100 4.366094)
);

blocks
(
hex (0 1 2 3 0 1 10 11) (39 20 1) simpleGrading (0.04 0.3 1)
hex (1 4 5 2 1 4 12 10) (10 20 1) simpleGrading (1 0.3 1)
hex (4 6 7 5 4 6 13 12) (39 20 1) simpleGrading (25 0.3 1)
hex (2 5 8 9 10 12 14 15) (10 395 1) simpleGrading (1 25 1)
);

edges
(
);
patches
(
empty axis
(
(0 1 1 0)
(1 4 4 1)
(4 6 6 4)
)

wedge back
(
(0 1 2 3)
(0 1 10 11)
(1 4 5 2)
(1 4 12 10)
)

wedge front
(
(4 6 7 5)
(4 6 13 12)
(2 5 8 9)
(10 12 14 15)
)

patch lx_inlet
(
(0 0 11 3)
)

patch rx_inlet
(
(6 6 7 13)
)

wall top_flange
(
(2 3 11 10)
(9 2 10 15)
)

wall btm_flange
(
(7 5 12 13)
(5 8 14 12)
)
patch outlet
(
(8 9 15 14)
)

);

mergePatchPairs
(
);


The axis of symmetry is along X, the radius along Y and the one cell thickness along Z. The mesh is quite corse and I do get the 'Severe skewness' warning.
When I run the simulation, the result seem to be much like the flow between 2D flanged jets, as the velocity between the flanges seem to be constant along the radius.
I also used another file, which gives the same result, but the FRONT and BACK (wedge) palnes are parallel for the flanges and at an angle for the jets.
I am very likely doing some very stupid mistake, and would appreciate if you could help.

Regards,
Gianfilippo
gfcoppola is offline   Reply With Quote

Old   December 28, 2007, 12:14
Default In my setup, the mesh is cylin
  #24
sek
Member
 
Sung-Eun Kim
Join Date: Mar 2009
Posts: 76
Rep Power: 9
sek is on a distinguished road
In my setup, the mesh is cylindrical, covering the 180 deg sector in the azimuthal direction. As a result, I have wedge elements touching the axis of symmetry at r = 0. The axis is not a boundary in my setup. The solutions seems to converge fine. But the solutions near the axis are not quite what I expected, showing high velocity in the body wake where I expect typical wake velocity profiles. I used starToFoam to convert the mesh generated using GridGen. Has anyone had similar experience?
sek is offline   Reply With Quote

Old   January 7, 2008, 07:42
Default Sung-Eun, Did you solve you
  #25
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 198
Blog Entries: 1
Rep Power: 10
egp is on a distinguished road
Sung-Eun,

Did you solve your problem? If your domain is 180 deg, I don't think the wedge BC is appropriate. I would use the cyclic condition instead.

Eric
egp is offline   Reply With Quote

Old   January 7, 2008, 09:13
Default I haven;t reolved the problem
  #26
sek
Member
 
Sung-Eun Kim
Join Date: Mar 2009
Posts: 76
Rep Power: 9
sek is on a distinguished road
I haven;t reolved the problem yet. The solution converges and looks fine in the boundary layer. The velocity field in the wake at and near the axis looks wrong.

I'm using symmetry boundary for 0 and 180 deg plane. My question was about the treatments at r = 0. Apparently, the axis (r=0) is not a boundary in cell-centered unstructured FV . Yet, the solution behavior at and near the axis made me suspicious of whether it has anything to do with the treatments of the wedge-shaped cells touching the axis.
sek is offline   Reply With Quote

Old   October 22, 2010, 05:55
Default 5 wedge
  #27
Senior Member
 
Join Date: Sep 2010
Location: France
Posts: 224
Rep Power: 9
T.D. is on a distinguished road
Hi guys,
Anyone has ideas about how to setup flow between two parallel disks, as an axi-symmetric flow problem.

i tried 5 wedge geometry, and applied a rotatingWallVelocity for Upper disk and zero fixedValue for lower and for the Outer side i applied a zero Gradient.
For pressure i applied all zeroGradient except for Outer i put fixedValue 0.
for the empty axis i put in U, p a zeroGradient also.

when i run with icoFoam, i got extraordinary nonlogic results. it didn't work!

Any Ideas?
help please

thanks a lot
T.D. is offline   Reply With Quote

Old   February 8, 2013, 00:35
Default Errorneous velocity while using wedge type BC for axisymmetric flow
  #28
Member
 
Anant Diwakar
Join Date: Jan 2013
Posts: 68
Rep Power: 5
diwakaranant is on a distinguished road
Hi

I am solving axisymmetric flow problem of flow between two rotating cylinders using rhoCentralFoam. I generated the mesh in GMSH and converted it into openFoam using gmshToFoam utility.
But when I am running the simulation, I am getting non-zero radial velocity
near the outer rotating cylinder, which ideally should be zero.

The Uy and the velocity vector plots near the outer rotating wall
are attached.

The rotation is about X axis. Hence Uy is the radial velocity component.

Can some one tell why is this happening ?

Thanks
Anant
Attached Images
File Type: jpg Uy_near_wall.jpg (18.2 KB, 45 views)
File Type: jpg velocity_vectors.jpg (63.9 KB, 35 views)
diwakaranant is offline   Reply With Quote

Old   February 8, 2013, 00:48
Default Errorneous velocity while using wedge type BC for axisymmetric flow
  #29
Member
 
Anant Diwakar
Join Date: Jan 2013
Posts: 68
Rep Power: 5
diwakaranant is on a distinguished road
Hi

I am solving axisymmetric flow problem of flow between two rotating cylinders using rhoCentralFoam. I generated the mesh in GMSH and converted it into openFoam using gmshToFoam utility.
But when I am running the simulation, I am getting non-zero radial velocity
near the outer rotating cylinder, which ideally should be zero.

The Uy and the velocity vector plots near the outer rotating wall
are attached.

The rotation is about X axis. Hence Uy is the radial velocity component.

Can some one tell why is this happening ?

Thanks
Anant
Attached Images
File Type: jpg Uy_near_wall.jpg (18.2 KB, 5 views)
File Type: jpg velocity_vectors.jpg (63.9 KB, 4 views)
diwakaranant is offline   Reply With Quote

Old   September 13, 2013, 04:20
Default
  #30
New Member
 
anshul bansal
Join Date: Jun 2013
Posts: 22
Rep Power: 5
anshul bansal is on a distinguished road
anyone please tell me how i put my wedge patch for back and front face in 0/p,0/T,0/U files.

means i put like this
back
{
type wedge;
}

front
{
type wedge;
}


is it correct.?


actually i find this error when i run my c-d nozzle with rhoCentralFoam solver


--> FOAM FATAL IO ERROR:

patch type 'genericPatch' not constraint type 'wedge'
for patch back of field p in file "/home/anshu/OpenFOAM/anshu-2.2.1/run/tutorials/compressible/rhoCentralFoam/with_axis/0/p"

file: /home/anshu/OpenFOAM/anshu-2.2.1/run/tutorials/compressible/rhoCentralFoam/with_axis/0/p.boundaryField.back from line 41 to line 41.

From function wedgeFvPatchField<Type>::wedgeFvPatchField
(
const fvPatch& p,
const Field<Type>& field,
dictionary& dict
)

in file fields/fvPatchFields/constraint/wedge/wedgeFvPatchField.C at line 99.

FOAM exiting
anshul bansal is offline   Reply With Quote

Old   September 13, 2013, 04:28
Default
  #31
Senior Member
 
Join Date: Sep 2010
Location: France
Posts: 224
Rep Power: 9
T.D. is on a distinguished road
Hi,

This should help you:

http://openfoamwiki.net/index.php/Ma...s/AxiSymmetric

Regards,

T.D.
T.D. is offline   Reply With Quote

Old   September 16, 2013, 06:59
Default
  #32
New Member
 
anshul bansal
Join Date: Jun 2013
Posts: 22
Rep Power: 5
anshul bansal is on a distinguished road
thanx for ur reply.
first i have problem in checkMesh of my geometry.
i got this problem-

***Number of edges not aligned with or perpendicular to non-empty directions: 48277
<<Writing 20200 points on non-aligned edges to set nonAlignedEdges

what should i do?
i can send u my geometry file in your email .
*(i made my geometry in gambit and import in openfoam)
anshul bansal is offline   Reply With Quote

Old   September 23, 2013, 05:51
Default
  #33
New Member
 
anshul bansal
Join Date: Jun 2013
Posts: 22
Rep Power: 5
anshul bansal is on a distinguished road
Can anybody help me to put my boundary condition and patch field in c-d nozzle problem.

my boundary file is-
(
top
{
type wall;
nFaces 1000;
startFace 49825;
}
inlet
{
type patch;
nFaces 175;
startFace 50825;
}
outlet
{
type patch;
nFaces 175;
startFace 51000;
}
front
{
type wedge;
nFaces 2000;
startFace 51175;
}
back
{
type wedge;
nFaces 2000;
startFace 53175;
}
)


and intial (0 file) for pressure is --
internalField uniform 101325;

boundaryField
{
inlet
{
type fixedValue;
value uniform 235000;
}

outlet
{
type zeroGradient;
}

top
{
type zeroGradient;
}


back
{
type wedge;
}
front
{
type wedge;
}


}
Attached Files
File Type: gz U.tar.gz (489 Bytes, 0 views)
File Type: gz T.tar.gz (490 Bytes, 0 views)
anshul bansal is offline   Reply With Quote

Old   January 13, 2015, 07:25
Default
  #34
New Member
 
Rok
Join Date: Sep 2014
Posts: 1
Rep Power: 0
Rokerda is on a distinguished road
Quote:
Originally Posted by hjasak View Post
Please make sure that the front and back of the wedge are in two separate patches.

Enjoy,

Hrv
Thank you very much for this information. I have been struggling with my code for a few days now
Rokerda is offline   Reply With Quote

Old   July 4, 2016, 04:55
Default
  #35
Senior Member
 
Join Date: Jan 2013
Posts: 311
Rep Power: 6
openfoammaofnepo is on a distinguished road
Hi All,

For the boundary condition "wedge" in OpenFOAM, how is it treated when the numerical fluxes are computed? Are they actually assumed that there are no any fluxes (i.e. convection, molecular and turbulent diffusions) across the wedge face?

I checked the source file as follows, but did not find any information in it about how the numerical flues on them.

https://github.com/OpenFOAM/OpenFOAM...FvPatchField.C

Do you kow where I can find the relevant information in the source files? Thank you.

OFFO
openfoammaofnepo is offline   Reply With Quote

Old   August 22, 2016, 07:20
Default
  #36
New Member
 
Join Date: May 2016
Posts: 14
Rep Power: 2
dzordz is on a distinguished road
Quote:
Originally Posted by anshul bansal View Post
thanx for ur reply.
first i have problem in checkMesh of my geometry.
i got this problem-

***Number of edges not aligned with or perpendicular to non-empty directions: 48277
<<Writing 20200 points on non-aligned edges to set nonAlignedEdges

what should i do?
i can send u my geometry file in your email .
*(i made my geometry in gambit and import in openfoam)
Hi Anshul,

what I found out is that when I create a geometry with snappyHexMesh (don't know if it is the same with Gambit) the front and back are not identical! So I get the same error when using extrude wedge. So what I do first is a linear extrude in one direction, so the both front and back patches are identical. Then you do wedge extrude and this should no longer be a problem.

Cheers.
dzordz is offline   Reply With Quote

Old   October 11, 2016, 19:39
Default
  #37
Senior Member
 
Join Date: Jan 2013
Posts: 311
Rep Power: 6
openfoammaofnepo is on a distinguished road
Dear Eric,

Have you got the idea about the wedge boundary conditions implemented in OpenFOAM? I mean the governing equation in Cartersian coordinate system is solver, but the transformation for the BC is applied and then this small wedge can be used to represent the axisymmetric problem. Could you please comment the transformation implemented in OF? Thank you so much.


Quote:
Originally Posted by egp View Post
OF community,

I've recently begun experimenting with OF, and am posting my first question to the message board.

When simulating flow over an axisymmetric body, the User Manual (U-145, U-146) specifies that the wedge-type cells must be 5 degrees, and must straddle one of the coordinate planes.

What is the reason for this restriction?

So that I can understand the implementation, where would I look in the source?

Thank you,

Eric Paterson
Penn State Univ
State College, PA USA
openfoammaofnepo 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
Gmsh problem with cyclic bcbs sripplinger Open Source Meshers: Gmsh, Netgen, CGNS, ... 1 September 28, 2010 09:29
Implementing new bcbs srinath OpenFOAM 3 September 3, 2008 05:11
Basics of setting up BCbs question shawn OpenFOAM Running, Solving & CFD 0 December 15, 2006 14:23
BCbs for high Renumbers gjesing OpenFOAM Running, Solving & CFD 10 May 7, 2006 09:51
CFD code about axisymmetric bodies at hypersonic speeds Afshin Azari Main CFD Forum 5 October 26, 1998 18:09


All times are GMT -4. The time now is 00:16.