Behaviour of Numerical Schemes

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

 September 7, 2001, 10:52 Behaviour of Numerical Schemes #1 Ravi. B. R Guest   Posts: n/a Dear Sir, I have a question regarding the behaviour of numerical schemes, in particular Roe's scheme for evaluating Inviscid fluxes. I have a NS-code which works for supersonic and hypersonic flows with implicit time marching (euler). I have noticed that though for High Mach numbers the residue oscillates and I have to reduce the courant number, though implicit schemes are unconditionaly stable and that it should converge for any courant/CFD number. I would like this point to be clarified. thanks, ravi.

 September 8, 2001, 02:55 Re: Behaviour of Numerical Schemes #2 Kang, Seok Koo Guest   Posts: n/a Unconditional stability does not mean the scheme is always stable under any numerical or flow conditions, since it is based on linear stability analysis. First-order Euler method is known as fairly stable method, however, it does not allow very large CFL number also. If your problem is really stiff in time, you can use other implicit scheme. But, in my opinion, they need exessive computation time. There can be several reasons for that your solution does not converge for high CFL number : choice of grid system, approximation of inviscid/viscous flux jacobian, using factorised scheme or not, etc.

 September 17, 2001, 07:52 Re: Behaviour of Numerical Schemes #3 Rajani KUmar Guest   Posts: n/a Hi, As for the theory it is ok to say that Implicit schems are not having the stability problem and also no CFL restriction. but the question is whether you are appling boundary conditions implicitly or not. boundary conditions in the case of euler equations is very easy. but for Navier strokes it is difficult(MacCormak paper 1982). as u are saying there are oscilations in residual u have to check whether you are appling sufficient numerical diffusion. when u deal with NS equations u have to add second order diffusion implicitly and 4th order diffusion explicitly. these conditions may not be nessary for every case. these are required for flows which are having shocks etc. so there are lot of things u have to take care. Godunov schems are good choice. by Rajani Kumar

 September 22, 2001, 10:18 Re: Behaviour of Numerical Schemes #4 Ravi Guest   Posts: n/a Hi Rajani Kumar, Thanks for the wealth of info. on schemes. Indeed, I haven't applied the boundary condition IMPLICITLY and it's one of the reasons that it behaves this way. Apart from that I have blowing at the plate surface and residues oscillate more vigourously in this case. Shocks do are present in this case as I am trying to solve for supersonic and hypersonic flows with blowing/gas injection at the surface. When there is no blowing the residue fall is quite smooth. By the way let me know your email ID. Are you a student studying in any Univ. in US. Thanks for the info again. ravi.

 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 botp OpenFOAM 2 March 11, 2011 16:27 Valerio FLUENT 2 July 28, 2008 04:02 Praveen Main CFD Forum 0 September 6, 2002 07:09 Ravi B R Main CFD Forum 2 September 9, 2001 04:23 X. Ye Main CFD Forum 7 August 31, 1999 17:05

All times are GMT -4. The time now is 15:31.