
[Sponsors] 
how to deal with sharp corner points for compressible simulations? 

LinkBack  Thread Tools  Display Modes 
January 3, 2013, 23:24 
how to deal with sharp corner points for compressible simulations?

#1 
New Member
Lipo Wang
Join Date: Dec 2012
Posts: 14
Rep Power: 5 
Dear friends,
Happy new year! I am doing a viscous compressible simulation with simple (square block) geometry. The boundary temperature is fixed. The headache now is for the pressure (or density) at the corner points. I have tried several different ones, like use the continuity equation to update density, to use the mean of two adjacent surface points' pressure etc. But after about 6000 steps, the code always diverges at the corner points (with very large velocity and negative density). Many thanks for your good ideas! 

January 4, 2013, 10:38 

#2 
Super Moderator

Hi
I am assuming you have a grid point at the corner , right ? Usually people just update density, update velocity, update energy Now you can calculate the new pressure. Then set the velocity to zero. Are you using a compressible, finite volume approach ? If yes, what is your numerical flux function ? I have a scheme in which I implement boundary conditions in a weak manner even for boundary grid points. I would be interested to see how it performs on your problem. Can you describe your geometry and flow conditions ?
__________________
http://twitter.com/cfdlab 

January 4, 2013, 10:53 

#3  
New Member
Lipo Wang
Join Date: Dec 2012
Posts: 14
Rep Power: 5 
Quote:
Thanks. The question I raised is how to update the pressure (or density) at the corner points. For boundary points on flat surface, I can use, for instance the continuity equation to update density. But if goes to the corner point, we may see multiple choices (such as the average of two adjacent surface points etc). But different methods have been tried and no one can work well. The simulation I am doing has quite simple geometry: a 2D channel with a square block inside. The flow conditions are: compressible, subsonic inlet, if the outlet is supersonic or not depends on the size of the block. Temperature on the channel surface and block surface are constant. Many thanks if you can provide good ideas. 

January 4, 2013, 11:15 

#4 
Super Moderator
Oliver Gloth
Join Date: Mar 2009
Location: Todtnau, Germany
Posts: 116
Rep Power: 9 
Hello,
do you have a nodecentred or cellcentred scheme? Usually, if the Mach number is moderate, you should be able to get decent results by updating the equations with the proper fluxes across the boundary  without any explicitly enforced values (e.g. zero velocity). The problem is that the values (assuming a finite volume scheme) are in the centre of the control volume. Even if you have a nodecentred solver, the values are not really on the nodes but in the centres of the corresponding control volumes (sort of a half cell on the boundary). If you enforce noslip here, you might run into problems. You can still have a noslip condition, but you need to assume zero velocity on the actual domain border and use a corresponding shear flux for your momentum equations. Having said this, in the past I have used the approach that Praveen sketched in his post and it worked  problems only appeared when wall functions needed to be implemented. If you have a structured code, however, I would strongly recommend to use cellcentred variables. For an unstructured code, nodecentred might be better  depending on your cell types. Cheers, Oliver 

January 4, 2013, 21:16 

#5  
New Member
Lipo Wang
Join Date: Dec 2012
Posts: 14
Rep Power: 5 
Quote:
Thanks for your continuous help. For the problem I mentioned I used the nodecentered MacCormack scheme. The boundary points are located on the solid surface and thus I must specify the values, especially the density (or pressure). The case is still laminar (no wall function). I found very bad zigzag profiles. Then try to cure such problem by decreasing time step. But the problem is still there, especially at the corner points. After about 5000 step the density around the corner points (few points away from the corner) became negative. Do you have any good suggestion? Or if you can provide me some simple codes for reference (just the Cartesian grid case) for the viscous compressible flows, that will be really appreciated. Best, Lipo 

January 5, 2013, 06:53 

#6 
Super Moderator
Oliver Gloth
Join Date: Mar 2009
Location: Todtnau, Germany
Posts: 116
Rep Power: 9 
Hi Lipo,
I might have an idea what causes your problems: If you use the MacCormack scheme in the standard finite difference formulation, you might have a somewhat degenerated stencil (see picture below). The red nodes illustrate the stencil for your corner node. The nodes on the boundary (blue circles) will give wrong first derivatives due to the enforced noslip condition. As a result your viscous flux will be too small and might lead to some weird expansion around the corner... A finite volume scheme would use the green control volume and differentiate between two different boundary contributions for this node. Maybe you could only use differences with the healthy nodes  this means do not switch between forward and backward difference, but only use the differences into the field. Of course this can lead to stability issues, but it might still work if it is done for those corners only. If I had to write a code for this problem it would be a cell centred finite volume scheme with some kind of upwind flux (Roe, AUSM, van Leer, ...) and Runge Kutta for time integration. 'hope this helps a little. Cheers, Oliver 

January 6, 2013, 01:23 

#7 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 474
Rep Power: 11 
I don't think a central difference method can converge without artificial dissipation.
That being said, if your flow is supercritical you probably need to solve the Euler equations coupled rather than uncoupled. Also, if your input subsonic flow's Mach number is high enough, your 2D channel may choke, which requires the correct input and output boundary conditions. Since the flow is viscous, I'm skeptical it really matters how you compute the pressure at the corner since the corner point is imbedded in a boundary layer, (and therefor subsonic). Of course this requires that you throw enough grid points at it. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
.STL: nonclosed manifold surface  giulio.topazio  OpenFOAM Native Meshers: snappyHexMesh and Others  29  May 1, 2016 19:04 
No layers in a small gap  bobburnquist  OpenFOAM Native Meshers: snappyHexMesh and Others  6  August 26, 2015 09:38 
[DesignModeler] DM's JScript: FPoint()'s GetPoint(i) function picks points backwards?  ANT  ANSYS Meshing & Geometry  2  July 23, 2012 15:25 
BlockMeshmergePatchPairs  hjasak  OpenFOAM Native Meshers: blockMesh  11  August 15, 2008 07:36 
ICEM meshing of sharp corner  Sans  CFX  3  January 17, 2008 08:36 