|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Join Date: Mar 2012
Location: Munich, Germany
Posts: 67
Rep Power: 14 ![]() |
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 ![]() ![]() ![]() ![]() 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: http://www.cfd-online.com/Forums/ope...-equation.html) 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 |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Super Moderator
|
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 | Thread Starter | Forum | Replies | Last Post |
thobois class engineTopoChangerMesh error | Peter_600 | OpenFOAM | 4 | August 2, 2014 09:52 |
smoothSolver diverges - solution in using PBiCG solver? | makaveli_lcf | OpenFOAM Running, Solving & CFD | 3 | September 11, 2013 12:44 |
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 |
free " ordinary differential equations solver" | ztdep | Main CFD Forum | 2 | March 21, 2006 01:58 |