Commit 86c71f01 authored by Kilian Lackhove's avatar Kilian Lackhove

APE: Added sponge layer support

parent 4c83a458
......@@ -372,7 +372,7 @@ void APE::DoOdeRhs(const Array<OneD, const Array<OneD, NekDouble> >&inarray,
std::vector<SolverUtils::ForcingSharedPtr>::const_iterator x;
for (x = m_forcing.begin(); x != m_forcing.end(); ++x)
(*x)->Apply(m_fields, outarray, outarray, m_time);
(*x)->Apply(m_fields, inarray, outarray, m_time);
......@@ -8,14 +8,14 @@
<metric type="L2" id="1">
<value variable="p" tolerance="1e-12">7.93225</value>
<value variable="u" tolerance="1e-12">0.0189841</value>
<value variable="v" tolerance="1e-12">0.0109285</value>
<value variable="p" tolerance="1e-4">7.56584</value>
<value variable="u" tolerance="1e-7">0.0174419</value>
<value variable="v" tolerance="1e-7">0.00907637</value>
<metric type="Linf" id="2">
<value variable="p" tolerance="1e-12">18.7643</value>
<value variable="u" tolerance="1e-12">0.0322549</value>
<value variable="v" tolerance="1e-12">0.0460461</value>
<value variable="p" tolerance="1e-4">18.7643</value>
<value variable="u" tolerance="1e-7">0.0322549</value>
<value variable="v" tolerance="1e-7">0.0460461</value>
......@@ -337,9 +337,9 @@
<D VAR="p" VALUE="0"/>
<D VAR="u" VALUE="0"/>
<D VAR="v" VALUE="0"/>
<D VAR="p" USERDEFINEDTYPE="RiemannInvariantBC" VALUE="0" />
<D VAR="u" USERDEFINEDTYPE="RiemannInvariantBC" VALUE="0" />
<D VAR="v" USERDEFINEDTYPE="RiemannInvariantBC" VALUE="0" />
......@@ -358,10 +358,25 @@
<E VAR="u" VALUE="0"/>
<E VAR="v" VALUE="0"/>
<FUNCTION NAME="ExactSolution"> <!-- Not really the exact solution -->
<E VAR="p" VALUE="0"/>
<E VAR="u" VALUE="0"/>
<E VAR="v" VALUE="0"/>
<FUNCTION NAME="InitialConditions">
<E VAR="p" VALUE="100*exp(-32*(x^2+y^2))"/> <!-- Gaussian pulse located at the origin -->
<E VAR="u" VALUE="0"/>
<E VAR="v" VALUE="0"/>
<FUNCTION NAME="SpongeCoefficient">
<E VAR="p" VALUE="(y>0.5) * -5000 * ((1-cos(PI*((y-0.5)/(1-0.5))))/2)^3" />
<E VAR="u" VALUE="(y>0.5) * -5000 * ((1-cos(PI*((y-0.5)/(1-0.5))))/2)^3" />
<E VAR="v" VALUE="(y>0.5) * -5000 * ((1-cos(PI*((y-0.5)/(1-0.5))))/2)^3" />
<FORCE TYPE="Absorption">
<COEFF> SpongeCoefficient </COEFF>
Markdown is supported
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment