# Singularity

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 July 20, 2000, 22:10 Singularity #1 Jim Guest   Posts: n/a Hi, Why are sharp corners considered as points of singularity in CFD. Is it because conflicting boundary conditions exist? And what are the different ways to resolve the issue ? Thanks !

 July 20, 2000, 23:04 Re: Singularity #2 John C. Chien Guest   Posts: n/a (1). There are at least two types of sharp corners: one is the 90 degree corner, such as the flow over a backward-facing step, and the other is the sharp trailing edge or lip of a nozzle. (2). In both cases, the flow always separate from the sharp corner. In other words, there is always a separating streamlines attached to the corner and divide the flow into two regions. (3). Local mesh refinement around the corner is a simple way to handle the difficulty. (4). Or you can replace the sharp corner by a rounded corner. (5). Or you can assume that there is a segment of separating streamline attached to the corner, and is parallel to the main flow direction. In this way, the sharp corner property is controlled by the main flow boundary layer. (that is the upstream side of the flow, instead of the recirculating part of the flow. both flows meet at the point of separation)

 July 21, 2000, 13:31 Re: Singularity #3 Kalyan Guest   Posts: n/a Yes, it is difficult to implement neumann boundary conditions at corners due to the difficulty associated with defining a unique surface normal. In inviscid flows, an infinite pressure gradient can be produced at the corners. There isn't a mathematically sound way of dealing with corners in such cases. IN viscous flow solvers, a neumann BC is used for the pressure at the walls. The conflicting BC are handled by defining a multi-valued pressure variable at the corners. Eg : if a 90 degrees corner exist in the x-y plane, there are many possible surface normals. dp/dx is used to get pressure at the walls and this pressure is used in the x-momentum equation. dp/dy is used to get another value of pressure at the same point and is used only in the y-momentum equation. If you use staggered grids, you do need pressure at the walls and hence there are no conflicting BC at the corners. In any case, staggered or non-staggered, corners require good resolution. This is because the pressure is relatively independent of the viscous forces (at high Re) but at these corners pressure and viscous forces are strongly coupled. This coupling should reach grid independence at the corners for an accurate solution.

 July 21, 2000, 16:41 Re: Singularity #4 Adrin Gharakhani Guest   Posts: n/a >> IN viscous flow solvers, a neumann BC is used for the pressure at the walls. The conflicting BC are handled by defining a multi-valued pressure variable at the corners. Eg : if a 90 degrees corner exist in the x-y plane, there are many possible surface normals. The pressure is single-valued even if there is a sharp corner. However, the (normal) pressure gradient is multi-valued at the corners. If the pressure were multi-valued then you'd have a pressure jump (and thus flow leakage?) at the corners! Adrin Gharakhani

 July 21, 2000, 17:50 Re: Singularity #5 Kalyan Guest   Posts: n/a From my experience, if the same value of pressure at the corner (however it is obtained) is used to compute both (dp/dx) and (dp/dy), then the solution diverges. Since both (dp/dx) and (dp/dy) qualify as normal gradients, pressure is computed using both. The values of pressure thus obtained are used in x-momentum and y-momentum equations respectively. This is OK for incompressible flows since actual value of pressure is of no consequence. Also, a multi-valued pressure is defined merely to compute normal pressure gradients along different directions. The actual value of the pressure is perhaps unique at the corner. In the limit of a infinite grid resolution, pressures computed at the corner from both direction ( x and y) become equal. In compressible flows, actual value of pressure appears in the equation of state and so the corner problem is more acute. It may be better to use staggered grids for such cases.

 July 22, 2000, 04:02 Re: Singularity #7 li yuan Guest   Posts: n/a The corner singularity, as well as axis singularity in 3D may be avoided by adopting a cell-centered finite volume method

 July 22, 2000, 11:49 Re: Singularity #8 John C. Chien Guest   Posts: n/a (1). Smart idea indeed. (2). If we use one cell centered finite volume to cover each problem domain, then each will become nice and perfect square. (3). When there is no longer internal resolution, everything becomes ONE. (4). I think, the pressure is always uniquely defined, whether in the interior or at the boundary. (5). In most cases, the real problem is modified to simulate a desired result. But it can have negative empact of the solution. (6). For the backward-facing step problem, there are typically two approaches to create the mesh. One is to use cartesian H-type mesh, in which one can model the corner flow field as if the flow separate from the corner point and moves in the direction parallel to the main stream direction. This will force the corner to be dominated by the upstream flow only. (7). The other approach is to use the body aligned curvilinear mesh, in which the wall represents a continuous streamline. This allows the flow to go around the corner. And the fan shaped mesh around the corner will allow the flow to turn. The mesh density will be high automatically around the corner, thus avoid the accuracy problem. Inside the boundary layer, near the corner, the Re is low, thus the problem is not serious. (8). For supersonic, inviscid flow around the corner, there will be corner centered expansion waves generated. But then the method of characteristics can provide the exact solutions, even though the corner is still a point. There, like any solution to the Euler equation, discontinuous solution is allowed.

 July 24, 2000, 11:48 Re: Singularity #9 TOT KTO 3HAET Guest   Posts: n/a I'm really glad to see how impressive is your way to circumvent the singularity in traction! I've come up with an impression that many people post here for the sake of posting itself, often providing completely misleading concepts, solutions and advices. Please try not to do so. There are not many persons in this forum who actually gives good answers, listen to them. As for singularity, ignoring what you encounter is not always the best solution. The singularity cannot simply be disregarded. Instead, it needs a carefull attention to get what's called "mesh independent solution". Cheers to all.

 July 24, 2000, 13:14 Re: Singularity #10 Kalyan Guest   Posts: n/a I agree with the point that the gradients at the corners may just be too high to resolve using a FD or FV scheme which typically use a 2nd or 4th order schemes for the Poisson equation. That is why I pointed out that the discretization scheme becomes more consistent with BC when the grid is refined enough to resolve these gradients and then the the two different values of pressure become equal and there is a unique pressure. Most LES of the backward facing step and cavity problems are performed using solvers with 2nd order Poisson solvers and I have not come across the "leakage" issue anywhere in the LES literature. In fact, I ensure global mass conservation in my own code and did not notice any leakage. However, my primary CFD research is not code/scheme development. But I sure would appreciate any information as to where the leakage issues are addressed.

 July 25, 2000, 15:02 Re: Singularity #11 Adrin Gharakhani Guest   Posts: n/a > The corner singularity, as well as axis singularity in 3D may be avoided by adopting a cell-centered finite volume method With a cell-centered method one "avoids" the singularity issue at the corners by (practically) assuming an averaged value of the variable within the control volume. That is essentially equivalent to numerically rounding the corners. So, you "avoid" singularity by solving another problem - a rounded-cornered one. Also, most cell-centered solution methodologies I know are of the lowest order of accuracy (which is not the issue here of course) Adrin Gharakhani

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58 litonx OpenFOAM Running, Solving & CFD 1 February 21, 2007 02:32 tangd OpenFOAM Running, Solving & CFD 0 May 29, 2006 09:07 liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07 Bigarella Main CFD Forum 16 August 6, 2002 18:23

All times are GMT -4. The time now is 01:44.

 Contact Us - CFD Online - Top