Nektar++
|
#include <cstdio>
#include <cstdlib>
#include <string>
#include <iostream>
#include <iomanip>
#include <cmath>
#include <MultiRegions/ExpList.h>
#include <MultiRegions/ExpList1D.h>
#include <MultiRegions/ExpList2D.h>
#include <MultiRegions/ExpList2DHomogeneous1D.h>
#include <MultiRegions/ExpList3DHomogeneous1D.h>
#include <MultiRegions/ExpList3DHomogeneous2D.h>
#include <MultiRegions/ExpList3D.h>
#include <MultiRegions/AssemblyMap/AssemblyMapDG.h>
#include <MultiRegions/ContField2D.h>
#include <MultiRegions/ContField3D.h>
#include <MultiRegions/DisContField2D.h>
#include <MultiRegions/DisContField3D.h>
#include <LocalRegions/MatrixKey.h>
#include <LocalRegions/Expansion2D.h>
#include <LocalRegions/Expansion3D.h>
#include <LocalRegions/Expansion.h>
#include <LibUtilities/BasicUtils/FieldIO.h>
#include <LibUtilities/BasicUtils/NekFactory.hpp>
#include <LibUtilities/BasicUtils/SessionReader.h>
#include <LibUtilities/BasicUtils/SharedArray.hpp>
#include <LibUtilities/Communication/Comm.h>
#include <LibUtilities/Memory/NekMemoryManager.hpp>
#include <SpatialDomains/MeshGraph2D.h>
#include <SpatialDomains/MeshGraph3D.h>
#include <SolverUtils/SolverUtilsDeclspec.h>
Go to the source code of this file.
Functions | |
void | COMPBL (Array< OneD, NekDouble > v, Array< OneD, NekDouble > dv) |
void | RK4 (Array< OneD, NekDouble > y, Array< OneD, NekDouble > dydx, int n, NekDouble x, NekDouble h, Array< OneD, NekDouble > yout) |
void | RKDUMB (Array< OneD, NekDouble > vstart, int nvar, NekDouble x1, NekDouble x2, int m_xpoints, Array< OneD, NekDouble > xx, Array< OneD, Array< OneD, NekDouble > > y) |
void | OUTPUT (int m_xpoints, Array< OneD, NekDouble > xx, Array< OneD, Array< OneD, NekDouble > > ff, int nQuadraturePts, Array< OneD, NekDouble > x_QuadraturePts, Array< OneD, NekDouble > y_QuadraturePts, Array< OneD, NekDouble > u_QuadraturePts, Array< OneD, NekDouble > v_QuadraturePts, Array< OneD, NekDouble > rho_QuadraturePts, Array< OneD, NekDouble > T_QuadraturePts) |
int | main (int argc, char *argv[]) |
Variables | |
NekDouble | m_Re |
NekDouble | m_Mach |
NekDouble | L |
NekDouble | m_Tinf |
NekDouble | m_Suth |
NekDouble | m_Tw |
NekDouble | m_Twall |
NekDouble | m_Gamma |
NekDouble | m_Pr |
NekDouble | m_long |
NekDouble | m_uInf |
NekDouble | m_rhoInf |
NekDouble | m_R |
NekDouble | m_vInf |
NekDouble | m_mu |
NekDouble | m_To = 273.11 |
const int | m_xpoints = 1000001 |
const NekDouble | Nvisc = 1 |
const NekDouble | Omega = 1 |
const NekDouble | etamax = 10.0 |
const NekDouble | errtol = 1e-5 |
int main | ( | int | argc, |
char * | argv[] | ||
) |
Calculate the compressible boundary layer solution for a given 3D mesh and dump the solution into a .rst file.
Definition at line 363 of file CompressibleBL.cpp.
References ASSERTL0, errtol, etamax, m_Gamma, m_long, m_Mach, m_mu, m_Pr, m_R, m_Re, m_rhoInf, m_Suth, m_Tinf, m_Tw, m_Twall, m_uInf, m_vInf, m_xpoints, OUTPUT(), RKDUMB(), Vmath::Vcopy(), and Nektar::LibUtilities::Write().
void OUTPUT | ( | int | m_xpoints, |
Array< OneD, NekDouble > | xx, | ||
Array< OneD, Array< OneD, NekDouble > > | ff, | ||
int | nQuadraturePts, | ||
Array< OneD, NekDouble > | x_QuadraturePts, | ||
Array< OneD, NekDouble > | y_QuadraturePts, | ||
Array< OneD, NekDouble > | u_QuadraturePts, | ||
Array< OneD, NekDouble > | v_QuadraturePts, | ||
Array< OneD, NekDouble > | rho_QuadraturePts, | ||
Array< OneD, NekDouble > | T_QuadraturePts | ||
) |
void RK4 | ( | Array< OneD, NekDouble > | y, |
Array< OneD, NekDouble > | dydx, | ||
int | n, | ||
NekDouble | x, | ||
NekDouble | h, | ||
Array< OneD, NekDouble > | yout | ||
) |
Perform the RK4 integration
Definition at line 141 of file CompressibleBL.cpp.
References COMPBL().
Referenced by RKDUMB().
void RKDUMB | ( | Array< OneD, NekDouble > | vstart, |
int | nvar, | ||
NekDouble | x1, | ||
NekDouble | x2, | ||
int | m_xpoints, | ||
Array< OneD, NekDouble > | xx, | ||
Array< OneD, Array< OneD, NekDouble > > | y | ||
) |
Calculate initial guess for RK4
Definition at line 188 of file CompressibleBL.cpp.
References COMPBL(), m_xpoints, and RK4().
Referenced by main().
const NekDouble errtol = 1e-5 |
Definition at line 99 of file CompressibleBL.cpp.
Referenced by main().
const NekDouble etamax = 10.0 |
Definition at line 98 of file CompressibleBL.cpp.
Definition at line 79 of file CompressibleBL.cpp.
Referenced by main(), Polylib::optdiff(), and Nektar::Utilities::InputMCF::ParseFile().
NekDouble m_Gamma |
Definition at line 84 of file CompressibleBL.cpp.
NekDouble m_long |
Definition at line 86 of file CompressibleBL.cpp.
Referenced by main().
NekDouble m_Mach |
Definition at line 78 of file CompressibleBL.cpp.
NekDouble m_mu |
Definition at line 91 of file CompressibleBL.cpp.
NekDouble m_Pr |
Definition at line 85 of file CompressibleBL.cpp.
NekDouble m_R |
Definition at line 89 of file CompressibleBL.cpp.
Referenced by main().
NekDouble m_Re |
Definition at line 77 of file CompressibleBL.cpp.
NekDouble m_rhoInf |
Definition at line 88 of file CompressibleBL.cpp.
NekDouble m_Suth |
Definition at line 81 of file CompressibleBL.cpp.
NekDouble m_Tinf |
Definition at line 80 of file CompressibleBL.cpp.
Referenced by main().
NekDouble m_To = 273.11 |
Definition at line 92 of file CompressibleBL.cpp.
NekDouble m_Tw |
Definition at line 82 of file CompressibleBL.cpp.
Referenced by main().
NekDouble m_Twall |
Definition at line 83 of file CompressibleBL.cpp.
NekDouble m_uInf |
Definition at line 87 of file CompressibleBL.cpp.
NekDouble m_vInf |
Definition at line 90 of file CompressibleBL.cpp.
Referenced by main().
const int m_xpoints = 1000001 |
Definition at line 94 of file CompressibleBL.cpp.
const NekDouble Nvisc = 1 |
Definition at line 96 of file CompressibleBL.cpp.
Referenced by COMPBL().
const NekDouble Omega = 1 |
Definition at line 97 of file CompressibleBL.cpp.
Referenced by COMPBL().