Solver for system of hyperbolic equations

treima · October 17, 2012, 04:29

Hi,

I derived the adjoint equations from the twodimensional Eulerequation (with boundary conditions) in theory. Now I would like to solve these with OpenFOAM. It´s a system of hyperbolic equations, so the normal discretisation of OF doesn´t work (I´ve tried

).

The system is given by

$a_{21} \Lambda_{2,x} + a_{31} \Lambda_{3,x} + a_{41} \Lambda_{4,x} + b_{21} \Lambda_{2,y} + b_{31} \Lambda_{3,y} + b_{41} \Lambda_{4,y} = 0$ ,
$\Lambda_{1,x} + a_{22} \Lambda_{2,x} + a_{32} \Lambda_{3,x} + a_{42} \Lambda_{4,x} + b_{22} \Lambda_{2,y} + b_{32} \Lambda_{3,y} + b_{42} \Lambda_{4,y} = 0$ ,
$a_{23} \Lambda_{2,x} + a_{33} \Lambda_{3,x} + a_{43} \Lambda_{4,x} + \Lambda_{1,y} + b_{23} \Lambda_{2,y} + b_{33} \Lambda_{3,y} + b_{43} \Lambda_{4,y} = 0$ ,
$a_{24} \Lambda_{2,x} + a_{44} \Lambda_{4,x} + b_{34} \Lambda_{3,y} + b_{44} \Lambda_{4,y} = 0$ .

where the indices x and y mean the x and y derivation and the a´s and b´s are factors, dependent from the solution of the Eulerequation.

According to the solution of the onedimensional advection equation (you can have a look at this thread: 1D advection equation) I tried to use a special discretisation with interpolation to the cell faces. You can find the code at the bottom.

For the solution I use the same timesteps like for the solution of the Eulerequation, but the adjoint solver is unstable and shows oscillations. My problem is, that I can´t explane me why

.

Do you have any suggestions for a correct implementation for this system of hyperbolic equations or some hints, why my solver doesn´t work?

Thank you!

regards

treima

praveen · October 21, 2012, 09:34

You need to use an upwind scheme or add some dissipation.

In the scalar case, you are using a central scheme with implicit time scheme. This is L2 stable so the computations do not blow up. But you do not get monotone solution, your solution is not between [0,1] but is going out of this interval.

Thread Tools
Show Printable Version Email this Page
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
thobois class engineTopoChangerMesh error	Peter_600	OpenFOAM	3	May 4, 2013 08:19
Is it possible to develop a generic solver for hyp. equations??	shyamdsundar	Lounge	1	January 21, 2010 08:45
Getting too many iterations by velocity solving (aborting). Changing U - Solver?	suitup	OpenFOAM Running, Solving & CFD	0	January 20, 2010 07:45
smoothSolver diverges - solution in using PBiCG solver?	makaveli_lcf	OpenFOAM Running, Solving & CFD	0	December 15, 2009 14:49
free " ordinary differential equations solver"	ztdep	Main CFD Forum	2	March 21, 2006 01:58

October 21, 2012, 09:34		#2
praveen Super Moderator Praveen. C Join Date: Mar 2009 Location: Bangalore Posts: 221 Blog Entries: 6 Rep Power: 7	You need to use an upwind scheme or add some dissipation. In the scalar case, you are using a central scheme with implicit time scheme. This is L2 stable so the computations do not blow up. But you do not get monotone solution, your solution is not between [0,1] but is going out of this interval.