Commit da29541a authored by Douglas Serson's avatar Douglas Serson Committed by Chris Cantwell

Fix another issue with time integration of VCSMapping

(cherry picked from commit 32daa119)
parent 23922f31
...@@ -190,18 +190,6 @@ namespace Nektar ...@@ -190,18 +190,6 @@ namespace Nektar
Array<OneD, Array<OneD, NekDouble> > &outarray, Array<OneD, Array<OneD, NekDouble> > &outarray,
const NekDouble time) const NekDouble time)
{ {
// Update mapping and Deal with Dirichlet boundary conditions
if (m_mapping->IsTimeDependent())
{
if (m_mapping->IsFromFunction())
{
// If the transformation is explicitly defined, update it here
// Otherwise, it will be done somewhere else (ForcingMovingBody)
m_mapping->UpdateMapping(time);
}
m_mapping->UpdateBCs(time);
}
EvaluateAdvectionTerms(inarray, outarray); EvaluateAdvectionTerms(inarray, outarray);
// Smooth advection // Smooth advection
...@@ -222,9 +210,21 @@ namespace Nektar ...@@ -222,9 +210,21 @@ namespace Nektar
// Add mapping terms // Add mapping terms
ApplyIncNSMappingForcing( inarray, outarray); ApplyIncNSMappingForcing( inarray, outarray);
// Calculate High-Order pressure boundary conditions // Calculate High-Order pressure boundary conditions
m_extrapolation->EvaluatePressureBCs(inarray,outarray,m_kinvis); m_extrapolation->EvaluatePressureBCs(inarray,outarray,m_kinvis);
// Update mapping and deal with Dirichlet boundary conditions
if (m_mapping->IsTimeDependent())
{
if (m_mapping->IsFromFunction())
{
// If the transformation is explicitly defined, update it here
// Otherwise, it will be done somewhere else (ForcingMovingBody)
m_mapping->UpdateMapping(time+m_timestep);
}
m_mapping->UpdateBCs(time+m_timestep);
}
} }
......
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