euler2d icon indicating copy to clipboard operation
euler2d copied to clipboard

2-D inviscid flow and adjoint solver


Data structure:

elem: Three vertices forming a triangle. Given in counter-clockwise order

edge: two vertices forming an edge. For a boundary edge, when we go from vertex one to vertex two, the computational domain lies to the left of the edge

tedge: Triangles adjacent to an edge; left is first and then right.

esue: Elements surrounding an element. Required only for implicit scheme

vedge: The vertices opposite to an edge; left first then right. Used for reconstruction.

spts: Solid wall points

bdedge: For each solid wall point, the two boundary edges sharing that point. Used for calculating wall normal, which is required for inviscid bc.

esubp: Elements surrounding a boundary point. Required for vertex averaging.

ptype: Type of vertex.


Files

flo: output: FLO.DAT

adj: input : FLO.DAT output: ADJ.DAT, ADJ.DXY

adjoint: input : ADJ.DXY (gradient wrt coordinates) output: DEF.OUT (gradient wrt boundary coordinates)

shape: input : DEF.OUT output: gradient wrt shape parameters, maybe HicksHenne


Changelog

Following functions have been changed solid_flux farfield_flux clcd costfunc

Linear interpolation has been replaced with first order treatment. This would be more robust while results dont seem to be affected.

  • Added entropy fix for |u| eigenvalues in roe flux