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

Boundary Conditions for Euler

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 3, 2001, 01:46
Default Boundary Conditions for Euler
  #1
Kevin
Guest
 
Posts: n/a
Hello!

I'm writing a 2nd-order accurate FVM Euler code to compute a flow around an airfoil, on triangular grids. I have two questions concerning the boundary conditions.

(1)Solid body: I store solutions at a cell center. So, on the airfoil, I need to compute the flux. How can I compute this flux? Perhaps I need to set up a ghost cell inside the body, and use the interior scheme to find the flux. But what values should I put in the ghost cells? Or is there any other ways to compute the boundary flux?

(2)Far-Field: To compute a lifting flow, I need to place a point vortex at a quarter chord of the airfoil. But I don't know how to determine the strength of the vortex. This seems impossible to me because I wouldn't know the value of the lift, which is directly related to the circulation(or vortex strength), until I finish computing the solution. Does any body know how to determine the strength?

Any comments are welcome. Thank you!

Kevin
  Reply With Quote

Old   March 3, 2001, 16:15
Default Re: Boundary Conditions for Euler
  #2
O
Guest
 
Posts: n/a
Concerning (1):

If I have this correctly you want to compute a flux using collocated FV at a solid boundary?

First the boundary cells usually have zero volume so the velocity in this cell is zero.

For the u-,v-, and w-momentum the velocity normal to the wall will (in most cases) be zero, thus a zero flux. If you are ignoring diffusion then just set your coefficient (A_i) normal to the wall equal to zero.

  Reply With Quote

Old   March 4, 2001, 23:53
Default Re: Boundary Conditions for Euler
  #3
Kevin
Guest
 
Posts: n/a
Thank you for your comments! But I'm not sure waht you mean by "the boundary cells usually have zero volume so the velocity in this cell is zero". Are these cells are those having a boundary face? Then, they all have a nonzero volume. I think they usually are, for it wouldn't make sense to introduce cells with zero volume around an airfoil.

Right, the fluxes across the boundary must be zero. But pressure remains in the flux even on the boundary. Maybe I could extrapolate the pressure to the boundary from the interior. I think this is worth trying. Thank you!

Kevin
  Reply With Quote

Old   March 6, 2001, 12:27
Default Re: Boundary Conditions for Euler
  #4
Frank
Guest
 
Posts: n/a
Hi Kevin

I actually don't understand, why the volume at boundary cells should be zero ?!

If you use a cell centered finite volume scheme with unstructured meshes I would actually not use any ghost cells. For a second order scheme, you have to reconstruct the variables anyhow. (in my experience it is better to reconstruct the primitive variables, than the conservative ones. Matter of choice ?!)

That means you will determine the gradients in the cell centers of each variable and then extrapolate the solution on the faces, where you solve the riemann problems for the convective terms. I would also not use a central scheme with artificial dissipation but use riemann solvers. At boundaries, all that changes is the stencil for the reconstruction, which becomes more one sides. One can simply extrapolate the solution on the boundary and then the flux is imply given exactly by the pressure in the momentum equations. If using riemann solvers, you can also extrapolate all variables on the boundary, compute the pressure (put it as left and right state of the riemann problem) and then before solving the riemann problem set the velocity normal to the solid wall it zero. When solving the riemann problem, this wil nothing but result in a contribution of the pressure again. This formulation of boundary conditions is only a weak one, since it is just via fluxes.

For the vortex: Look at a paper of Thomas and Salas: Thomas, J.L., Salas, M.D., Far--Field Boundary Conditions for Transonic Lifting Solutions to the Euler Equations AIAA J., 24, 1986, pp. 1074--1080

There is a formulation, how to obtain the circulation as a function of the lift. So, first you compute the lift by integrating the pressure around the airfoil (considering the normal vector on the surface and the angle of attack) , and from that you get the circulation. Of course, it is an iterative process. And the lift as well as the circulation will get its true value, when the solution is converged. At the beginning of the computation, both can get any value...

Hope it helps somewhat.

Frank ;-)

  Reply With Quote

Old   March 6, 2001, 21:54
Default Re: Boundary Conditions for Euler
  #5
Kevin
Guest
 
Posts: n/a
Thank you very much! I will try solving Riemann problem at boundary. For vortex, I didn't know it was an iterative process. I'll look at the paper.

Thank you! Your comment is very helpful!!

Kevin
  Reply With Quote

Old   March 7, 2001, 08:15
Default Re: Boundary Conditions for Euler
  #6
Frank
Guest
 
Posts: n/a
Hi Kevin,

solving Riemann problems at boundaries is of course only useful, if you use riemann solvers anyway.

I personally favor the HLLC scheme by Batten and Leschziner (Journal of computational physics.. end of 97 or 98). It is nicer and easier than Roe. The AUSMDV also needs entropy and shock fix, which the HLLC does not require. But this a personal choice.

I just find it easier to program, if I do not need to replace the flux at the wall boundary by the flux given by the pressure. I mean: In this case one has to use a special treatment for these walls. What I do is just to copy the extrapolated values on the side, which is interior of the field, to the solution exterior of the field (one always needs 2 states for the riemann problem). And then set both normal velocities (of the left and right state of the riemann problem) to zero. But AFTER I have computed the pressure ! (In case you extrapolate the conservative variables. I first get the pressure and then set the velocity to zero). And then one just do one loop over all boundary faces, including far field boundaries and wall boundaries.

Since the whole procedure is an iteration, I would suppose that the vorticity should get out correctly at the end of the computation.

In my opinion, Euler wall conditions are harder to get working really nicely than Navier Stokes. And everybody does something different.

Cheers,

Frank

  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
symmetry boundary conditions in cfx lost.identity CFX 41 May 22, 2013 07:21
Impinging Jet Boundary Conditions Anindya Main CFD Forum 24 January 11, 2012 14:40
OpenFOAM Variable Velocity Boundary Conditions NickolasPl OpenFOAM Programming & Development 2 May 19, 2011 05:37
boundary conditions and mesh exporting vaina74 Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 May 27, 2010 09:38
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15


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