Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in / Register
Toggle navigation
Menu
Open sidebar
Nektar
Nektar
Commits
ab29d869
Commit
ab29d869
authored
Mar 02, 2016
by
Daniele de Grazia
Browse files
Added regression test for Euler 1D
parent
24f1d307
Changes
3
Hide whitespace changes
Inline
Side-by-side
solvers/CompressibleFlowSolver/CMakeLists.txt
View file @
ab29d869
...
...
@@ -34,6 +34,7 @@ IF( NEKTAR_SOLVER_COMPRESSIBLE_FLOW )
ADD_NEKTAR_TEST
(
CylinderSubsonicMix
)
ADD_NEKTAR_TEST
(
CylinderSubsonic_P3
)
ADD_NEKTAR_TEST_LENGTHY
(
CylinderSubsonic_P8
)
ADD_NEKTAR_TEST
(
Euler1D
)
ADD_NEKTAR_TEST
(
IsentropicVortex16_P3
)
ADD_NEKTAR_TEST
(
IsentropicVortex_FRDG_SEM
)
ADD_NEKTAR_TEST_LENGTHY
(
IsentropicVortex_FRSD_SEM
)
...
...
solvers/CompressibleFlowSolver/Tests/Euler1D.tst
0 → 100644
View file @
ab29d869
<?xml version="1.0" encoding="utf-8"?>
<test>
<description>
Euler 1D P=3, WeakDG, MODIFIED
</description>
<executable>
CompressibleFlowSolver
</executable>
<parameters>
Euler1D.xml
</parameters>
<files>
<file
description=
"Session File"
>
Euler1D.xml
</file>
</files>
<metrics>
<metric
type=
"L2"
id=
"1"
>
<value
variable=
"rho"
tolerance=
"1e-12"
>
1.98838e-06
</value>
<value
variable=
"rhou"
tolerance=
"1e-12"
>
0.00067684
</value>
<value
variable=
"E"
tolerance=
"1e-12"
>
0.575708
</value>
</metric>
<metric
type=
"Linf"
id=
"2"
>
<value
variable=
"rho"
tolerance=
"1e-12"
>
1.98524e-05
</value>
<value
variable=
"rhou"
tolerance=
"1e-12"
>
0.00675771
</value>
<value
variable=
"E"
tolerance=
"1e-12"
>
5.74799
</value>
</metric>
</metrics>
</test>
solvers/CompressibleFlowSolver/Tests/Euler1D.xml
0 → 100644
View file @
ab29d869
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<GEOMETRY
DIM=
"1"
SPACE=
"1"
>
<VERTEX>
<V
ID=
"0"
>
0.000e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"1"
>
0.020e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"2"
>
0.040e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"3"
>
0.060e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"4"
>
0.080e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"5"
>
0.100e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"6"
>
0.120e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"7"
>
0.140e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"8"
>
0.160e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"9"
>
0.180e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"10"
>
0.200e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"11"
>
0.220e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"12"
>
0.240e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"13"
>
0.260e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"14"
>
0.280e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"15"
>
0.300e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"16"
>
0.320e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"17"
>
0.340e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"18"
>
0.360e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"19"
>
0.380e+00 0.000e+00 0.000e+00
</V>
<V
ID=
"20"
>
0.400e+00 0.000e+00 0.000e+00
</V>
</VERTEX>
<ELEMENT>
<S
ID=
"0"
>
0 1
</S>
<S
ID=
"1"
>
1 2
</S>
<S
ID=
"2"
>
2 3
</S>
<S
ID=
"3"
>
3 4
</S>
<S
ID=
"4"
>
4 5
</S>
<S
ID=
"5"
>
5 6
</S>
<S
ID=
"6"
>
6 7
</S>
<S
ID=
"7"
>
7 8
</S>
<S
ID=
"8"
>
8 9
</S>
<S
ID=
"9"
>
9 10
</S>
<S
ID=
"10"
>
10 11
</S>
<S
ID=
"11"
>
11 12
</S>
<S
ID=
"12"
>
12 13
</S>
<S
ID=
"13"
>
13 14
</S>
<S
ID=
"14"
>
14 15
</S>
<S
ID=
"15"
>
15 16
</S>
<S
ID=
"16"
>
16 17
</S>
<S
ID=
"17"
>
17 18
</S>
<S
ID=
"18"
>
18 19
</S>
<S
ID=
"19"
>
19 20
</S>
</ELEMENT>
<COMPOSITE>
<C
ID=
"0"
>
S[0-19]
</C>
<C
ID=
"1"
>
V[0]
</C>
<C
ID=
"2"
>
V[20]
</C>
</COMPOSITE>
<DOMAIN>
<D
ID=
"0"
>
C[0]
</D>
</DOMAIN>
</GEOMETRY>
<EXPANSIONS>
<E
COMPOSITE=
"C[0]"
NUMMODES=
"3"
FIELDS=
"rho,rhou,E"
TYPE=
"MODIFIED"
/>
</EXPANSIONS>
<CONDITIONS>
<PARAMETERS>
<P>
TimeStep = 1e-6
</P>
<P>
FinTime = 40e-6
</P>
<P>
NumSteps = FinTime/TimeStep
</P>
<P>
IO_CheckSteps = 100
</P>
<P>
IO_InfoSteps = 10
</P>
<P>
Gamma = 1.4
</P>
<P>
pInf = 101325
</P>
<P>
rhoInf = 1.225
</P>
<P>
u0 = 0.1
</P>
<P>
uInf = u0
</P>
<P>
uInfL = u0 + 0.01
</P>
<P>
uInfR = u0
</P>
</PARAMETERS>
<SOLVERINFO>
<I
PROPERTY=
"EQTYPE"
VALUE=
"EulerCFE"
/>
<I
PROPERTY=
"Projection"
VALUE=
"DisContinuous"
/>
<I
PROPERTY=
"AdvectionType"
VALUE=
"WeakDG"
/>
<I
PROPERTY=
"DiffusionType"
VALUE=
"LDGNS"
/>
<I
PROPERTY=
"TimeIntegrationMethod"
VALUE=
"ClassicalRungeKutta4"
/>
<I
PROPERTY=
"UpwindType"
VALUE=
"ExactToro"
/>
<I
PROPERTY=
"ProblemType"
VALUE=
"General"
/>
<I
PROPERTY=
"ViscosityType"
VALUE=
"Constant"
/>
</SOLVERINFO>
<VARIABLES>
<V
ID=
"0"
>
rho
</V>
<V
ID=
"1"
>
rhou
</V>
<V
ID=
"2"
>
E
</V>
</VARIABLES>
<BOUNDARYREGIONS>
<B
ID=
"0"
>
C[1]
</B>
<B
ID=
"1"
>
C[2]
</B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION
REF=
"0"
>
<D
VAR=
"rho"
VALUE=
"rhoInf"
/>
<D
VAR=
"rhou"
VALUE=
"rhoInf*uInfL"
/>
<D
VAR=
"E"
VALUE=
"pInf/(Gamma-1)+0.5*rhoInf*(uInfL*uInfL)"
/>
</REGION>
<REGION
REF=
"1"
>
<D
VAR=
"rho"
VALUE=
"rhoInf"
/>
<D
VAR=
"rhou"
VALUE=
"rhoInf*uInfR"
/>
<D
VAR=
"E"
VALUE=
"pInf/(Gamma-1)+0.5*rhoInf*(uInfR*uInfR)"
/>
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION
NAME=
"InitialConditions"
>
<E
VAR=
"rho"
DOMAIN=
"0"
VALUE=
"rhoInf"
/>
<E
VAR=
"rhou"
DOMAIN=
"0"
VALUE=
"rhoInf*u0"
/>
<E
VAR=
"E"
DOMAIN=
"0"
VALUE=
"pInf/(Gamma-1)+0.5*rhoInf*(u0*u0)"
/>
</FUNCTION>
<FUNCTION
NAME=
"ExactSolution"
>
<E
VAR=
"rho"
DOMAIN=
"0"
VALUE=
"rhoInf"
/>
<E
VAR=
"rhou"
DOMAIN=
"0"
VALUE=
"rhoInf*uInf"
/>
<E
VAR=
"E"
DOMAIN=
"0"
VALUE=
"pInf/(Gamma-1)+0.5*rhoInf*(uInf*uInf)"
/>
</FUNCTION>
</CONDITIONS>
</NEKTAR>
Dave Moxey
@dmoxey
mentioned in commit
e390b800
·
Jul 12, 2016
mentioned in commit
e390b800
mentioned in commit e390b800b03e825d2c09cc51345c180bc82fa07e
Toggle commit list
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment