Commit 389214e1 authored by Andrea Cassinelli's avatar Andrea Cassinelli

Various changes and tidying of user guide - mostly IncNavierStokesSolver

parent e5702e26
......@@ -170,10 +170,10 @@ Under this section it is possible to set the parameters of the simulation.
\subsection{Solver Info}
\begin{lstlisting}[style=XmlStyle]
<SOLVERINFO>
<I PROPERTY="EQType" VALUE="APE"/>
<I PROPERTY="Projection" VALUE="DisContinuous"/>
<I PROPERTY="TimeIntegrationMethod" VALUE="ClassicalRungeKutta4"/>
<I PROPERTY="UpwindType" VALUE="LaxFriedrichs"/>
<I PROPERTY="EQType" VALUE="APE" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="TimeIntegrationMethod" VALUE="ClassicalRungeKutta4" />
<I PROPERTY="UpwindType" VALUE="LaxFriedrichs" />
</SOLVERINFO>
\end{lstlisting}
\begin{itemize}
......@@ -283,10 +283,10 @@ We require a discontinuous Galerkin projection and use an explicit
fourth-order Runge-Kutta time integration scheme. We therefore set the following
solver information:
\begin{lstlisting}[style=XmlStyle]
<I PROPERTY="EQType" VALUE="APE"/>
<I PROPERTY="Projection" VALUE="DisContinuous"/>
<I PROPERTY="EQType" VALUE="APE"/>
<I PROPERTY="Projection" VALUE="DisContinuous"/>
<I PROPERTY="TimeIntegrationMethod" VALUE="ClassicalRungeKutta4"/>
<I PROPERTY="UpwindType" VALUE="LaxFriedrichs"/>
<I PROPERTY="UpwindType" VALUE="LaxFriedrichs"/>
\end{lstlisting}
To maintain numerical stability we must use a small time-step.
......
......@@ -184,11 +184,13 @@ Since we are solving the unsteady advection problem, we must specify this in the
solver information. We also choose to use a discontinuous flux-reconstruction
projection and use a Runge-Kutta order 4 time-integration scheme.
\begin{lstlisting}[style=XMLStyle]
<I PROPERTY="EQTYPE" VALUE="UnsteadyAdvection" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="AdvectionType" VALUE="FRDG" />
<I PROPERTY="UpwindType" VALUE="Upwind" />
<I PROPERTY="TimeIntegrationMethod" VALUE="ClassicalRungeKutta4"/>
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="UnsteadyAdvection" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="AdvectionType" VALUE="FRDG" />
<I PROPERTY="UpwindType" VALUE="Upwind" />
<I PROPERTY="TimeIntegrationMethod" VALUE="ClassicalRungeKutta4" />
</SOLVERINFO>
\end{lstlisting}
We choose to advect our solution for $20$ time units with a time-step of $0.01$
......@@ -307,8 +309,8 @@ information.
</PARAMETERS>
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="Helmholtz" />
<I PROPERTY="Projection" VALUE="Continuous" />
<I PROPERTY="EQTYPE" VALUE="Helmholtz" />
<I PROPERTY="Projection" VALUE="Continuous" />
</SOLVERINFO>
\end{lstlisting}
......@@ -468,12 +470,14 @@ We choose to use a continuous projection and an first-order implicit-explicit
time-integration scheme. The \inltt{DiffusionAdvancement} and
\inltt{AdvectionAdvancement} parameters specify how these terms are treated.
\begin{lstlisting}[style=XMLStyle]
<I PROPERTY="EQTYPE" VALUE="UnsteadyAdvectionDiffusion" />
<I PROPERTY="Projection" VALUE="Continuous" />
<I PROPERTY="DiffusionAdvancement" VALUE="Implicit" />
<I PROPERTY="AdvectionAdvancement" VALUE="Explicit" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder1" />
<I PROPERTY="GlobalSysSoln" VALUE="IterativeStaticCond" />
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="UnsteadyAdvectionDiffusion" />
<I PROPERTY="Projection" VALUE="Continuous" />
<I PROPERTY="DiffusionAdvancement" VALUE="Implicit" />
<I PROPERTY="AdvectionAdvancement" VALUE="Explicit" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder1" />
<I PROPERTY="GlobalSysSoln" VALUE="IterativeStaticCond" />
</SOLVERINFO>
\end{lstlisting}
We integrate for a total of $30$ time units with a time-step of $0.0005$,
......
......@@ -235,15 +235,15 @@ conditions are employed (i.e. $T_{w}$). Default value = 300.15$K$;
Under this section it is possible to set the solver information.
\begin{lstlisting}[style=XmlStyle]
<SOLVERINFO>
<I PROPERTY="EQType" VALUE="NavierStokesCFE" />
<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" />
<I PROPERTY="EquationOfState" VALUE="IdealGas" />
<I PROPERTY="EQType" VALUE="NavierStokesCFE" />
<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" />
<I PROPERTY="EquationOfState" VALUE="IdealGas" />
</SOLVERINFO>
\end{lstlisting}
\begin{itemize}
......@@ -481,7 +481,7 @@ where $s_0 = s_\kappa - 4.25\;log_{10}(p)$.
To enable the non-smooth viscosity model, the following line has to be added to the \inltt{SOLVERINFO} section:
\begin{lstlisting}[style=XmlStyle]
<SOLVERINFO>
<I PROPERTY="ShockCaptureType" VALUE="NonSmooth" />
<I PROPERTY="ShockCaptureType" VALUE="NonSmooth" />
<SOLVERINFO>
\end{lstlisting}
The diffusivity and the sensor can be controlled by the following parameters:
......@@ -514,7 +514,7 @@ where $S_\kappa$ is a normalised sensor value and serves as a forcing term for t
To enable the smooth viscosity model, the following line has to be added to the \inltt{SOLVERINFO} section:
\begin{lstlisting}[style=XmlStyle]
<SOLVERINFO>
<I PROPERTY="ShockCaptureType" VALUE="Smooth" />
<I PROPERTY="ShockCaptureType" VALUE="Smooth" />
<SOLVERINFO>
\end{lstlisting}
Furthermore, the extra viscosity variable \inltt{eps} has to be added to the variable list:
......
......@@ -326,10 +326,10 @@ and other parameters are specified.
</PARAMETERS>
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="PulseWavePropagation" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="EQTYPE" VALUE="PulseWavePropagation" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="TimeIntegrationMethod" VALUE="RungeKutta2_ImprovedEuler" />
<I PROPERTY="UpwindTypePulse" VALUE="UpwindPulse"/>
<I PROPERTY="UpwindTypePulse" VALUE="UpwindPulse" />
</SOLVERINFO>
<VARIABLES>
......@@ -423,7 +423,7 @@ Our simulation is started as described before and the results show the time
history for the conservative variables A and u, as well as for the
characteristic variables W1 and W2 at the beginning of the ascending aorta
(Artery 1). We can see that physically correct the shape of the inflow boundary
condition appears in the forward traveling characteristic W1. As we do not have
condition appears in the forward travelling characteristic W1. As we do not have
a terminal resistance at the outflow, one would normally expect W2 to be
constant. However this is not the case, as bifurcations cause reflections if the
radii of parent and daughter vessels are not well matching, leading to changes
......@@ -521,10 +521,10 @@ Finally the domain is specified by the first composite by
\paragraph{Solver Information:~}The Discontinuous Galerkin Method is used as projection scheme and the time-integration is performed by a simple Forward Euler scheme. A full list of possible time integration scheme is given in the parameter section of the \hyperref[PulseWaveSolver]{Pulse Wave Solver}
\begin{lstlisting}[style=XMLStyle]
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="PulseWavePropagation" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="TimeIntegrationMethod" VALUE="ForwardEuler" />
<I PROPERTY="UpwindTypePulse" VALUE="UpwindPulse"/>
<I PROPERTY="EQTYPE" VALUE="PulseWavePropagation" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="TimeIntegrationMethod" VALUE="ForwardEuler" />
<I PROPERTY="UpwindTypePulse" VALUE="UpwindPulse" />
</SOLVERINFO>
\end{lstlisting}
......
......@@ -97,10 +97,10 @@ use and (in the case the discontinuous Galerkin method is used)
the choice of numerical flux. A typical example would be:
\begin{lstlisting}[style=XmlStyle]
<SOLVERINFO>
<I PROPERTY="EqType" VALUE="NonlinearSWE">
<I PROPERTY="Projection" VALUE="DisContinuous">
<I PROPERTY="TimeIntegrationMethod" VALUE="ClassicalRungeKutta4">
<I PROPERTY="UpwindType" VALUE="HLLC">
<I PROPERTY="EqType" VALUE="NonlinearSWE" />
<I PROPERTY="Projection" VALUE="DisContinuous" />
<I PROPERTY="TimeIntegrationMethod" VALUE="ClassicalRungeKutta4" />
<I PROPERTY="UpwindType" VALUE="HLLC" />
</SOLVERINFO>
\end{lstlisting}
......
......@@ -93,7 +93,7 @@ where \inltt{FieldConvert} is the executable associated to the utility
FieldConvert, \inltt{test.xml} is the session file and
\inltt{test.vtu}, \inltt{test.dat}, \inltt{test.plt} are the desired
format outputs, either Paraview, VisIt or Tecplot formats. \\
When converting to \inltt{.dat} or \inltt{plt} format, it is possible to
When converting to \inltt{.dat} or \inltt{.plt} format, it is possible to
enable output with double precision, which is more accurate but requires larger
disk space. For example, double precision output in plt. format can be produced
with the command:
......
......@@ -65,11 +65,11 @@ In the following example, the driver \inltt{Standard} is used to manage the
execution of the incompressible Navier-Stokes equations:
\begin{lstlisting}[style=XMLStyle]
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" />
<I PROPERTY="SolverType" VALUE="VelocityCorrectionScheme" />
<I PROPERTY="Projection" VALUE="Galerkin" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2" />
<I PROPERTY="Driver" VALUE="Standard" />
<I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" />
<I PROPERTY="SolverType" VALUE="VelocityCorrectionScheme" />
<I PROPERTY="Projection" VALUE="Galerkin" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2" />
<I PROPERTY="Driver" VALUE="Standard" />
</SOLVERINFO>
\end{lstlisting}
......@@ -618,12 +618,12 @@ homogeneouns dimension; here an example
\begin{lstlisting}[style=XMLStyle]
<SOLVERINFO>
<I PROPERTY="SolverType" VALUE="VelocityCorrectionScheme"/>
<I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes"/>
<I PROPERTY="AdvectionForm" VALUE="Convective"/>
<I PROPERTY="Projection" VALUE="Galerkin"/>
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2"/>
<I PROPERTY="HOMOGENEOUS" VALUE="1D"/>
<I PROPERTY="SolverType" VALUE="VelocityCorrectionScheme" />
<I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" />
<I PROPERTY="AdvectionForm" VALUE="Convective" />
<I PROPERTY="Projection" VALUE="Galerkin" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2" />
<I PROPERTY="HOMOGENEOUS" VALUE="1D" />
</SOLVERINFO>
\end{lstlisting}
......@@ -652,11 +652,11 @@ for two harmonic directions (in Y and Z direction). For Example,
\begin{lstlisting}[style=XMLStyle]
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="UnsteadyAdvectionDiffusion" />
<I PROPERTY="Projection" VALUE="Continuous"/>
<I PROPERTY="HOMOGENEOUS" VALUE="2D"/>
<I PROPERTY="DiffusionAdvancement" VALUE="Implicit"/>
<I PROPERTY="AdvectionAdvancement" VALUE="Explicit"/>
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2"/>
<I PROPERTY="Projection" VALUE="Continuous" />
<I PROPERTY="HOMOGENEOUS" VALUE="2D" />
<I PROPERTY="DiffusionAdvancement" VALUE="Implicit" />
<I PROPERTY="AdvectionAdvancement" VALUE="Explicit" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2" />
</SOLVERINFO>
<PARAMETERS>
<P> TimeStep = 0.001 </P>
......@@ -683,16 +683,16 @@ below:
\begin{lstlisting}[style=XMLStyle]
<SOLVERINFO>
<I PROPERTY="EQTYPE" VALUE="UnsteadyAdvectionDiffusion" />
<I PROPERTY="Projection" VALUE="Continuous"/>
<I PROPERTY="HOMOGENEOUS" VALUE="2D"/>
<I PROPERTY="DiffusionAdvancement" VALUE="Implicit"/>
<I PROPERTY="AdvectionAdvancement" VALUE="Explicit"/>
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2"/>
<I PROPERTY="USEFFT" VALUE="FFTW"/>
<I PROPERTY="Projection" VALUE="Continuous" />
<I PROPERTY="HOMOGENEOUS" VALUE="2D" />
<I PROPERTY="DiffusionAdvancement" VALUE="Implicit" />
<I PROPERTY="AdvectionAdvancement" VALUE="Explicit" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2" />
<I PROPERTY="USEFFT" VALUE="FFTW" />
</SOLVERINFO>
\end{lstlisting}
The number of homogenenous modes has to be even. The Quasi-3D apporach can be
The number of homogenenous modes has to be even. The Quasi-3D approach can be
created starting from a 2D mesh and adding one homogenous expansion or starting
form a 1D mesh and adding two homogeneous expansions. Not other options
available. In case of a 1D homogeneous extension, the homogeneous direction will
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment