CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Problems with convergence and initial guessed fields in compressible steady flows

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 19, 2000, 15:51
Default Problems with convergence and initial guessed fields in compressible steady flows
Pedro Gil
Posts: n/a

I try to solve a supersonic steady problem neglecting the temporal term, in the compresible N-S ecuations with all terms term (i.e time step =inf). The geomtry is the typical test problem (channel with a circular arc "bump") My boundary conditions are density,pressure,end velocity prescribed at inlet , and extrapolate at outlet.

I start with a uniform guessed velocity field with the same value that the inlet velocity. idem for temperature and idem for pressure.

I solve Momentum and pressure correction (simple) since convergence in velocity is reached. (criterium is mean((unew-uold)/uold)).

I use strong underrelaxation in momentum and pressure correction. (p=p+alpha*p');

The two reflecting shock waves begin to develop in the channel and when iterations converge , I solve energy ecuation and actualize properties. Then begin to iterate again.

When I do this two or three times the fiels diverge and the method shut down.

Why this happens?.Is the grid incapable of support great magnitude gradient values? ,Might I regrid for more acuratte?.

The problem ussually ocurrs in the cross of two shock waves.

If anyone has any Solution or explanation to the Problem please contact me!.

You Can help me a lot with my studies-end proyect!!!!.

  Reply With Quote

Old   April 19, 2000, 17:58
Default Re: Problems with convergence and initial guessed fields in compressible steady flows
John C. Chien
Posts: n/a
(1). You need to find out first whether your method works for supersonic flow with shocks. (2). Try simple test cases to address this issue. (3). What is your Mach number? (4). Is the flow transonic or supersonic? These will affect your solutions. good luck.
  Reply With Quote

Old   April 20, 2000, 00:02
Default Re: Problems with convergence and initial guessed fields in compressible steady flows
Bernard Parent
Posts: n/a
Compressible supersonic flows at high-speed are much stiffer and more 'non-linear' that incompressible flows. In order to obtain a solution, you must march in time slowly, by retaining as much of the time accuracy of the scheme as possible. A CFL constant not much greater than one is generally required to obtain a solution when strong shocks are present, unless the solution is already close to being converged. Just write your scheme so that it preserves time accuracy and your problem should converge more smoothly.

Also, when dealing with shocks, a TVD scheme based on a monotonic first order scheme is generally a must to discretize the convective terms. If the shocks are weak (transonic flow for eg) you might get away with a Jameson-type artificial dissipation scheme but such a method will fail at higher shock Mach numbers..

In short, compressible flows cannot generally be solved using an extension of methods developped for incompressible flows, the opposite being also true: you need different schemes for different problems.

  Reply With Quote

Old   April 25, 2000, 11:34
Default Re: Problems with convergence and initial guessed fields in compressible steady flows
Dietmar Giebert
Posts: n/a
Dear Petro,

From your e-mail I can see you are using the SIMPLE-algorithm to solve the problem. Basically, this approach (originally developed for incompressible flow) should be extended for its use to solve high Mach number flows with shocks. So you have to take into account that a variation in pressure (p'in SIMPLE) leads to a variation/correction in both velocity and density. The influence upon temperature is achieved via the equation of state in the outer iterations process. How this can be done as well as the derivation of a pressure correction equation in its compressible formulation is given in 'Karki and Patankar (1989)ressure based ......, AIAA Journal, Vol 27, pp. 1167-1174'

Another important fact which strongly influences convergence of the method is the treatment of density at the surfaces of the control volumes. UPWIND scheme provides sufficient numerical diffusion to get a stabil algorithm but leads to bad shock resolution. Here, a mixed approach (e.g. 80% higher order scheme plus 20% UPWIND) improves resolution of the shock significantly while retaining stability. To derive the compressible pressure correction equation the correction of density at the surfaces can be approximated by using the UPWIND scheme. This has no influence upon the accuracy of the final solution since corrections tend to zero then.

As far as I know from my own calculations one can start to solve the whole set of equations (momentum, pressure-correction, turbulence, energy etc...) right from the beginning. But the assumption of a constant total enthalpy should be justified at a first step.

Good luck,


  Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
problem with Min/max rho tH3f0rC3 OpenFOAM 8 July 31, 2019 09:48
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
convergence problem when use pisoFoam, LES for wind tunnel case Forrest_Lei OpenFOAM 3 July 19, 2011 06:00
Compressible epsilon blows up swahono OpenFOAM 10 November 26, 2010 05:38
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 14:09

All times are GMT -4. The time now is 09:05.