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
23922f31
Commit
23922f31
authored
Aug 25, 2016
by
Douglas Serson
Committed by
Chris Cantwell
Sep 08, 2016
Browse files
Fix problem with time accuracy of VCSMapping
(cherry picked from commit
cd38f0e5
)
parent
119c3c6b
Changes
3
Hide whitespace changes
Inline
Side-by-side
solvers/IncNavierStokesSolver/EquationSystems/StandardExtrapolate.cpp
View file @
23922f31
...
...
@@ -75,10 +75,9 @@ namespace Nektar
const
Array
<
OneD
,
const
Array
<
OneD
,
NekDouble
>
>
&
N
,
NekDouble
kinvis
)
{
m_pressureCalls
++
;
if
(
m_HBCdata
.
num_elements
()
>
0
)
{
m_pressureCalls
++
;
// Rotate HOPBCs storage
RollOver
(
m_pressureHBCs
);
...
...
solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.cpp
View file @
23922f31
...
...
@@ -221,7 +221,7 @@ namespace Nektar
}
// Add mapping terms
ApplyIncNSMappingForcing
(
outarray
);
ApplyIncNSMappingForcing
(
inarray
,
outarray
);
// Calculate High-Order pressure boundary conditions
m_extrapolation
->
EvaluatePressureBCs
(
inarray
,
outarray
,
m_kinvis
);
...
...
@@ -739,6 +739,7 @@ namespace Nektar
* Explicit terms of the mapping
*/
void
VCSMapping
::
ApplyIncNSMappingForcing
(
const
Array
<
OneD
,
Array
<
OneD
,
NekDouble
>
>
&
inarray
,
Array
<
OneD
,
Array
<
OneD
,
NekDouble
>
>
&
outarray
)
{
int
physTot
=
m_fields
[
0
]
->
GetTotPoints
();
...
...
@@ -758,7 +759,7 @@ namespace Nektar
{
for
(
int
i
=
0
;
i
<
m_nConvectiveFields
;
++
i
)
{
vel
[
i
]
=
m_fields
[
i
]
->
GetPhys
()
;
vel
[
i
]
=
inarray
[
i
]
;
m_fields
[
0
]
->
HomogeneousBwdTrans
(
vel
[
i
],
velPhys
[
i
]);
}
}
...
...
@@ -766,8 +767,8 @@ namespace Nektar
{
for
(
int
i
=
0
;
i
<
m_nConvectiveFields
;
++
i
)
{
vel
[
i
]
=
m_fields
[
i
]
->
GetPhys
()
;
Vmath
::
Vcopy
(
physTot
,
m_fields
[
i
]
->
GetPhys
()
,
1
,
velPhys
[
i
],
1
);
vel
[
i
]
=
inarray
[
i
]
;
Vmath
::
Vcopy
(
physTot
,
inarray
[
i
]
,
1
,
velPhys
[
i
],
1
);
}
}
...
...
solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.h
View file @
23922f31
...
...
@@ -64,6 +64,7 @@ namespace Nektar
//
void
ApplyIncNSMappingForcing
(
const
Array
<
OneD
,
Array
<
OneD
,
NekDouble
>
>
&
inarray
,
Array
<
OneD
,
Array
<
OneD
,
NekDouble
>
>
&
outarray
);
virtual
~
VCSMapping
();
...
...
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