Nektar++
Public Member Functions | Static Public Member Functions | Static Public Attributes | Protected Member Functions | Protected Attributes | Friends | List of all members
Nektar::PulseWavePropagation Class Reference

#include <PulseWavePropagation.h>

Inheritance diagram for Nektar::PulseWavePropagation:
[legend]

Public Member Functions

Array< OneD, NekDouble > & GetA0 ()
 
Array< OneD, NekDouble > & GetBeta ()
 
Array< OneD, NekDouble > & GetAlpha ()
 
Array< OneD, NekDouble > & GetN ()
 
NekDouble GetRho ()
 
NekDouble GetDomains ()
 
- Public Member Functions inherited from Nektar::PulseWaveSystem
int GetNdomains ()
 
Array< OneD, MultiRegions::ExpListSharedPtrUpdateVessels (void)
 
- Public Member Functions inherited from Nektar::SolverUtils::UnsteadySystem
SOLVER_UTILS_EXPORT ~UnsteadySystem () override=default
 Destructor. More...
 
SOLVER_UTILS_EXPORT NekDouble GetTimeStep (const Array< OneD, const Array< OneD, NekDouble > > &inarray)
 Calculate the larger time-step mantaining the problem stable. More...
 
SOLVER_UTILS_EXPORT NekDouble GetTimeStep ()
 
SOLVER_UTILS_EXPORT void SetTimeStep (const NekDouble timestep)
 
SOLVER_UTILS_EXPORT void SteadyStateResidual (int step, Array< OneD, NekDouble > &L2)
 
SOLVER_UTILS_EXPORT LibUtilities::TimeIntegrationSchemeSharedPtrGetTimeIntegrationScheme ()
 Returns the time integration scheme. More...
 
SOLVER_UTILS_EXPORT LibUtilities::TimeIntegrationSchemeOperatorsGetTimeIntegrationSchemeOperators ()
 Returns the time integration scheme operators. More...
 
- Public Member Functions inherited from Nektar::SolverUtils::EquationSystem
virtual SOLVER_UTILS_EXPORT ~EquationSystem ()
 Destructor. More...
 
SOLVER_UTILS_EXPORT void InitObject (bool DeclareField=true)
 Initialises the members of this object. More...
 
SOLVER_UTILS_EXPORT void DoInitialise (bool dumpInitialConditions=true)
 Perform any initialisation necessary before solving the problem. More...
 
SOLVER_UTILS_EXPORT void DoSolve ()
 Solve the problem. More...
 
SOLVER_UTILS_EXPORT void TransCoeffToPhys ()
 Transform from coefficient to physical space. More...
 
SOLVER_UTILS_EXPORT void TransPhysToCoeff ()
 Transform from physical to coefficient space. More...
 
SOLVER_UTILS_EXPORT void Output ()
 Perform output operations after solve. More...
 
SOLVER_UTILS_EXPORT std::string GetSessionName ()
 Get Session name. More...
 
template<class T >
std::shared_ptr< T > as ()
 
SOLVER_UTILS_EXPORT void ResetSessionName (std::string newname)
 Reset Session name. More...
 
SOLVER_UTILS_EXPORT LibUtilities::SessionReaderSharedPtr GetSession ()
 Get Session name. More...
 
SOLVER_UTILS_EXPORT MultiRegions::ExpListSharedPtr GetPressure ()
 Get pressure field if available. More...
 
SOLVER_UTILS_EXPORT void ExtraFldOutput (std::vector< Array< OneD, NekDouble > > &fieldcoeffs, std::vector< std::string > &variables)
 
SOLVER_UTILS_EXPORT void PrintSummary (std::ostream &out)
 Print a summary of parameters and solver characteristics. More...
 
SOLVER_UTILS_EXPORT void SetLambda (NekDouble lambda)
 Set parameter m_lambda. More...
 
SOLVER_UTILS_EXPORT SessionFunctionSharedPtr GetFunction (std::string name, const MultiRegions::ExpListSharedPtr &field=MultiRegions::NullExpListSharedPtr, bool cache=false)
 Get a SessionFunction by name. More...
 
SOLVER_UTILS_EXPORT void SetInitialConditions (NekDouble initialtime=0.0, bool dumpInitialConditions=true, const int domain=0)
 Initialise the data in the dependent fields. More...
 
SOLVER_UTILS_EXPORT void EvaluateExactSolution (int field, Array< OneD, NekDouble > &outfield, const NekDouble time)
 Evaluates an exact solution. More...
 
SOLVER_UTILS_EXPORT NekDouble L2Error (unsigned int field, const Array< OneD, NekDouble > &exactsoln, bool Normalised=false)
 Compute the L2 error between fields and a given exact solution. More...
 
SOLVER_UTILS_EXPORT NekDouble L2Error (unsigned int field, bool Normalised=false)
 Compute the L2 error of the fields. More...
 
SOLVER_UTILS_EXPORT NekDouble LinfError (unsigned int field, const Array< OneD, NekDouble > &exactsoln=NullNekDouble1DArray)
 Linf error computation. More...
 
SOLVER_UTILS_EXPORT NekDouble H1Error (unsigned int field, const Array< OneD, NekDouble > &exactsoln, bool Normalised=false)
 Compute the H1 error between fields and a given exact solution. More...
 
SOLVER_UTILS_EXPORT Array< OneD, NekDoubleErrorExtraPoints (unsigned int field)
 Compute error (L2 and L_inf) over an larger set of quadrature points return [L2 Linf]. More...
 
SOLVER_UTILS_EXPORT void Checkpoint_Output (const int n)
 Write checkpoint file of m_fields. More...
 
SOLVER_UTILS_EXPORT void Checkpoint_Output (const int n, MultiRegions::ExpListSharedPtr &field, std::vector< Array< OneD, NekDouble > > &fieldcoeffs, std::vector< std::string > &variables)
 Write checkpoint file of custom data fields. More...
 
SOLVER_UTILS_EXPORT void Checkpoint_BaseFlow (const int n)
 Write base flow file of m_fields. More...
 
SOLVER_UTILS_EXPORT void WriteFld (const std::string &outname)
 Write field data to the given filename. More...
 
SOLVER_UTILS_EXPORT void WriteFld (const std::string &outname, MultiRegions::ExpListSharedPtr &field, std::vector< Array< OneD, NekDouble > > &fieldcoeffs, std::vector< std::string > &variables)
 Write input fields to the given filename. More...
 
SOLVER_UTILS_EXPORT void ImportFld (const std::string &infile, Array< OneD, MultiRegions::ExpListSharedPtr > &pFields)
 Input field data from the given file. More...
 
SOLVER_UTILS_EXPORT void ImportFldToMultiDomains (const std::string &infile, Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int ndomains)
 Input field data from the given file to multiple domains. More...
 
SOLVER_UTILS_EXPORT void ImportFld (const std::string &infile, std::vector< std::string > &fieldStr, Array< OneD, Array< OneD, NekDouble > > &coeffs)
 Output a field. Input field data into array from the given file. More...
 
SOLVER_UTILS_EXPORT void ImportFld (const std::string &infile, MultiRegions::ExpListSharedPtr &pField, std::string &pFieldName)
 Output a field. Input field data into ExpList from the given file. More...
 
SOLVER_UTILS_EXPORT void SessionSummary (SummaryList &vSummary)
 Write out a session summary. More...
 
SOLVER_UTILS_EXPORT Array< OneD, MultiRegions::ExpListSharedPtr > & UpdateFields ()
 
SOLVER_UTILS_EXPORT LibUtilities::FieldMetaDataMapUpdateFieldMetaDataMap ()
 Get hold of FieldInfoMap so it can be updated. More...
 
SOLVER_UTILS_EXPORT NekDouble GetTime ()
 Return final time. More...
 
SOLVER_UTILS_EXPORT int GetNcoeffs ()
 
SOLVER_UTILS_EXPORT int GetNcoeffs (const int eid)
 
SOLVER_UTILS_EXPORT int GetNumExpModes ()
 
SOLVER_UTILS_EXPORT const Array< OneD, int > GetNumExpModesPerExp ()
 
SOLVER_UTILS_EXPORT int GetNvariables ()
 
SOLVER_UTILS_EXPORT const std::string GetVariable (unsigned int i)
 
SOLVER_UTILS_EXPORT int GetTraceTotPoints ()
 
SOLVER_UTILS_EXPORT int GetTraceNpoints ()
 
SOLVER_UTILS_EXPORT int GetExpSize ()
 
SOLVER_UTILS_EXPORT int GetPhys_Offset (int n)
 
SOLVER_UTILS_EXPORT int GetCoeff_Offset (int n)
 
SOLVER_UTILS_EXPORT int GetTotPoints ()
 
SOLVER_UTILS_EXPORT int GetTotPoints (int n)
 
SOLVER_UTILS_EXPORT int GetNpoints ()
 
SOLVER_UTILS_EXPORT int GetSteps ()
 
SOLVER_UTILS_EXPORT void SetSteps (const int steps)
 
SOLVER_UTILS_EXPORT NekDouble GetTimeStep ()
 
SOLVER_UTILS_EXPORT void CopyFromPhysField (const int i, Array< OneD, NekDouble > &output)
 
SOLVER_UTILS_EXPORT void CopyToPhysField (const int i, const Array< OneD, const NekDouble > &input)
 
SOLVER_UTILS_EXPORT Array< OneD, NekDouble > & UpdatePhysField (const int i)
 
SOLVER_UTILS_EXPORT void ZeroPhysFields ()
 
SOLVER_UTILS_EXPORT void FwdTransFields ()
 
SOLVER_UTILS_EXPORT void SetModifiedBasis (const bool modbasis)
 
SOLVER_UTILS_EXPORT int GetCheckpointNumber ()
 
SOLVER_UTILS_EXPORT void SetCheckpointNumber (int num)
 
SOLVER_UTILS_EXPORT int GetCheckpointSteps ()
 
SOLVER_UTILS_EXPORT void SetCheckpointSteps (int num)
 
SOLVER_UTILS_EXPORT int GetInfoSteps ()
 
SOLVER_UTILS_EXPORT void SetInfoSteps (int num)
 
SOLVER_UTILS_EXPORT void SetIterationNumberPIT (int num)
 
SOLVER_UTILS_EXPORT void SetWindowNumberPIT (int num)
 
SOLVER_UTILS_EXPORT Array< OneD, const Array< OneD, NekDouble > > GetTraceNormals ()
 
SOLVER_UTILS_EXPORT void SetTime (const NekDouble time)
 
SOLVER_UTILS_EXPORT void SetTimeStep (const NekDouble timestep)
 
SOLVER_UTILS_EXPORT void SetInitialStep (const int step)
 
SOLVER_UTILS_EXPORT void SetBoundaryConditions (NekDouble time)
 Evaluates the boundary conditions at the given time. More...
 
SOLVER_UTILS_EXPORT bool NegatedOp ()
 Identify if operator is negated in DoSolve. More...
 
- Public Member Functions inherited from Nektar::SolverUtils::ALEHelper
virtual ~ALEHelper ()=default
 
virtual SOLVER_UTILS_EXPORT void v_ALEInitObject (int spaceDim, Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
 
SOLVER_UTILS_EXPORT void InitObject (int spaceDim, Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
 
virtual SOLVER_UTILS_EXPORT void v_UpdateGridVelocity (const NekDouble &time)
 
virtual SOLVER_UTILS_EXPORT void v_ALEPreMultiplyMass (Array< OneD, Array< OneD, NekDouble > > &fields)
 
SOLVER_UTILS_EXPORT void ALEDoElmtInvMass (Array< OneD, Array< OneD, NekDouble > > &traceNormals, Array< OneD, Array< OneD, NekDouble > > &fields, NekDouble time)
 Update m_fields with u^n by multiplying by inverse mass matrix. That's then used in e.g. checkpoint output and L^2 error calculation. More...
 
SOLVER_UTILS_EXPORT void ALEDoElmtInvMassBwdTrans (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray)
 
SOLVER_UTILS_EXPORT void MoveMesh (const NekDouble &time, Array< OneD, Array< OneD, NekDouble > > &traceNormals)
 
const Array< OneD, const Array< OneD, NekDouble > > & GetGridVelocity ()
 
SOLVER_UTILS_EXPORT const Array< OneD, const Array< OneD, NekDouble > > & GetGridVelocityTrace ()
 
SOLVER_UTILS_EXPORT void ExtraFldOutputGridVelocity (std::vector< Array< OneD, NekDouble > > &fieldcoeffs, std::vector< std::string > &variables)
 

Static Public Member Functions

static EquationSystemSharedPtr create (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
 Creates an instance of this class. More...
 

Static Public Attributes

static std::string className
 Name of class. More...
 
- Static Public Attributes inherited from Nektar::SolverUtils::UnsteadySystem
static std::string cmdSetStartTime
 
static std::string cmdSetStartChkNum
 

Protected Member Functions

 PulseWavePropagation (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
 
 ~PulseWavePropagation () override=default
 
void v_InitObject (bool DeclareField=false) override
 Init object for UnsteadySystem class. More...
 
void DoOdeRhs (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
 
void DoOdeProjection (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
 
void SetPulseWaveBoundaryConditions (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
 
void GetFluxVector (const Array< OneD, Array< OneD, NekDouble > > &physfield, Array< OneD, Array< OneD, Array< OneD, NekDouble > > > &flux)
 DG Pulse Wave Propagation routines: More...
 
void v_GenerateSummary (SolverUtils::SummaryList &s) override
 
- Protected Member Functions inherited from Nektar::PulseWaveSystem
 PulseWaveSystem (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
 Initialises PulseWaveSystem class members. More...
 
 ~PulseWaveSystem () override=default
 
void v_InitObject (bool DeclareField=false) override
 
void v_DoInitialise (bool dumpInitialConditions=false) override
 Sets up initial conditions. More...
 
void v_DoSolve () override
 Solves an unsteady problem. More...
 
void LinkSubdomains (Array< OneD, Array< OneD, Array< OneD, NekDouble > > > &fields)
 Links the subdomains. More...
 
void BifurcationRiemann (Array< OneD, NekDouble > &Au, Array< OneD, NekDouble > &uu, Array< OneD, NekDouble > &beta, Array< OneD, NekDouble > &A_0, Array< OneD, NekDouble > &alpha)
 Riemann Problem for Bifurcation. More...
 
void MergingRiemann (Array< OneD, NekDouble > &Au, Array< OneD, NekDouble > &uu, Array< OneD, NekDouble > &beta, Array< OneD, NekDouble > &A_0, Array< OneD, NekDouble > &alpha)
 Riemann Problem for Merging Flow. More...
 
void InterfaceRiemann (Array< OneD, NekDouble > &Au, Array< OneD, NekDouble > &uu, Array< OneD, NekDouble > &beta, Array< OneD, NekDouble > &A_0, Array< OneD, NekDouble > &alpha)
 Riemann Problem for Interface/Junction. More...
 
void v_Output (void) override
 
void CheckPoint_Output (const int n)
 
NekDouble v_L2Error (unsigned int field, const Array< OneD, NekDouble > &exactsoln=NullNekDouble1DArray, bool Normalised=false) override
 Compute the L2 error between fields and a given exact solution. More...
 
NekDouble v_LinfError (unsigned int field, const Array< OneD, NekDouble > &exactsoln=NullNekDouble1DArray) override
 Compute the L_inf error between fields and a given exact solution. More...
 
void WriteVessels (const std::string &outname)
 Write input fields to the given filename. More...
 
void EnforceInterfaceConditions (const Array< OneD, const Array< OneD, NekDouble > > &fields)
 
- Protected Member Functions inherited from Nektar::SolverUtils::UnsteadySystem
SOLVER_UTILS_EXPORT UnsteadySystem (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
 Initialises UnsteadySystem class members. More...
 
SOLVER_UTILS_EXPORT void v_InitObject (bool DeclareField=true) override
 Init object for UnsteadySystem class. More...
 
SOLVER_UTILS_EXPORT void v_DoSolve () override
 Solves an unsteady problem. More...
 
virtual SOLVER_UTILS_EXPORT void v_PrintStatusInformation (const int step, const NekDouble cpuTime)
 Print Status Information. More...
 
virtual SOLVER_UTILS_EXPORT void v_PrintSummaryStatistics (const NekDouble intTime)
 Print Summary Statistics. More...
 
SOLVER_UTILS_EXPORT void v_DoInitialise (bool dumpInitialConditions=true) override
 Sets up initial conditions. More...
 
SOLVER_UTILS_EXPORT void v_GenerateSummary (SummaryList &s) override
 Print a summary of time stepping parameters. More...
 
virtual SOLVER_UTILS_EXPORT NekDouble v_GetTimeStep (const Array< OneD, const Array< OneD, NekDouble > > &inarray)
 Return the timestep to be used for the next step in the time-marching loop. More...
 
virtual SOLVER_UTILS_EXPORT bool v_PreIntegrate (int step)
 
virtual SOLVER_UTILS_EXPORT bool v_PostIntegrate (int step)
 
virtual SOLVER_UTILS_EXPORT bool v_RequireFwdTrans ()
 
virtual SOLVER_UTILS_EXPORT void v_SteadyStateResidual (int step, Array< OneD, NekDouble > &L2)
 
virtual SOLVER_UTILS_EXPORT bool v_UpdateTimeStepCheck ()
 
SOLVER_UTILS_EXPORT NekDouble MaxTimeStepEstimator ()
 Get the maximum timestep estimator for cfl control. More...
 
SOLVER_UTILS_EXPORT void CheckForRestartTime (NekDouble &time, int &nchk)
 
SOLVER_UTILS_EXPORT void SVVVarDiffCoeff (const Array< OneD, Array< OneD, NekDouble > > vel, StdRegions::VarCoeffMap &varCoeffMap)
 Evaluate the SVV diffusion coefficient according to Moura's paper where it should proportional to h time velocity. More...
 
SOLVER_UTILS_EXPORT void DoDummyProjection (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
 Perform dummy projection. More...
 
- Protected Member Functions inherited from Nektar::SolverUtils::EquationSystem
SOLVER_UTILS_EXPORT EquationSystem (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
 Initialises EquationSystem class members. More...
 
virtual SOLVER_UTILS_EXPORT void v_InitObject (bool DeclareFeld=true)
 Initialisation object for EquationSystem. More...
 
virtual SOLVER_UTILS_EXPORT void v_DoInitialise (bool dumpInitialConditions=true)
 Virtual function for initialisation implementation. More...
 
virtual SOLVER_UTILS_EXPORT void v_DoSolve ()
 Virtual function for solve implementation. More...
 
virtual SOLVER_UTILS_EXPORT NekDouble v_LinfError (unsigned int field, const Array< OneD, NekDouble > &exactsoln=NullNekDouble1DArray)
 Virtual function for the L_inf error computation between fields and a given exact solution. More...
 
virtual SOLVER_UTILS_EXPORT NekDouble v_L2Error (unsigned int field, const Array< OneD, NekDouble > &exactsoln=NullNekDouble1DArray, bool Normalised=false)
 Virtual function for the L_2 error computation between fields and a given exact solution. More...
 
virtual SOLVER_UTILS_EXPORT NekDouble v_H1Error (unsigned int field, const Array< OneD, NekDouble > &exactsoln=NullNekDouble1DArray, bool Normalised=false)
 Virtual function for the H_1 error computation between fields and a given exact solution. More...
 
virtual SOLVER_UTILS_EXPORT void v_TransCoeffToPhys ()
 Virtual function for transformation to physical space. More...
 
virtual SOLVER_UTILS_EXPORT void v_TransPhysToCoeff ()
 Virtual function for transformation to coefficient space. More...
 
virtual SOLVER_UTILS_EXPORT void v_GenerateSummary (SummaryList &l)
 Virtual function for generating summary information. More...
 
virtual SOLVER_UTILS_EXPORT void v_SetInitialConditions (NekDouble initialtime=0.0, bool dumpInitialConditions=true, const int domain=0)
 
virtual SOLVER_UTILS_EXPORT void v_EvaluateExactSolution (unsigned int field, Array< OneD, NekDouble > &outfield, const NekDouble time)
 
virtual SOLVER_UTILS_EXPORT void v_Output (void)
 
virtual SOLVER_UTILS_EXPORT MultiRegions::ExpListSharedPtr v_GetPressure (void)
 
virtual SOLVER_UTILS_EXPORT bool v_NegatedOp (void)
 Virtual function to identify if operator is negated in DoSolve. More...
 
virtual SOLVER_UTILS_EXPORT void v_ExtraFldOutput (std::vector< Array< OneD, NekDouble > > &fieldcoeffs, std::vector< std::string > &variables)
 

Protected Attributes

SolverUtils::RiemannSolverSharedPtr m_riemannSolver
 
SolverUtils::AdvectionSharedPtr m_advObject
 
Array< OneD, PulseWaveBoundarySharedPtrm_Boundary
 
- Protected Attributes inherited from Nektar::PulseWaveSystem
Array< OneD, MultiRegions::ExpListSharedPtrm_vessels
 
size_t m_nDomains
 
size_t m_currentDomain
 
size_t m_nVariables
 
UpwindTypePulse m_upwindTypePulse
 
Array< OneD, int > m_fieldPhysOffset
 
NekDouble m_rho
 
NekDouble m_pext
 
NekDouble m_C
 
NekDouble m_RT
 
NekDouble m_pout
 
Array< OneD, Array< OneD, NekDouble > > m_A_0
 
Array< OneD, Array< OneD, NekDouble > > m_A_0_trace
 
Array< OneD, Array< OneD, NekDouble > > m_beta
 
Array< OneD, Array< OneD, NekDouble > > m_beta_trace
 
Array< OneD, Array< OneD, NekDouble > > m_gamma
 
Array< OneD, Array< OneD, NekDouble > > m_alpha
 
Array< OneD, Array< OneD, NekDouble > > m_alpha_trace
 
Array< OneD, Array< OneD, NekDouble > > m_trace_fwd_normal
 
std::map< int, SpatialDomains::CompositeMapm_domain
 
std::vector< int > m_domOrder
 
std::map< int, std::vector< int > > m_domainToFilterIDs
 
std::map< int, int > m_filterToVesselID
 
Array< OneD, Array< OneD, NekDouble > > m_pressure
 
PulseWavePressureAreaSharedPtr m_pressureArea
 
bool extraFields = false
 
Array< OneD, Array< OneD, NekDouble > > m_PWV
 
Array< OneD, Array< OneD, NekDouble > > m_W1
 
Array< OneD, Array< OneD, NekDouble > > m_W2
 
std::vector< std::vector< InterfacePointShPtr > > m_vesselIntfcs
 
std::vector< std::vector< InterfacePointShPtr > > m_bifurcations
 
std::vector< std::vector< InterfacePointShPtr > > m_mergingJcts
 
- Protected Attributes inherited from Nektar::SolverUtils::UnsteadySystem
LibUtilities::TimeIntegrationSchemeSharedPtr m_intScheme
 Wrapper to the time integration scheme. More...
 
LibUtilities::TimeIntegrationSchemeOperators m_ode
 The time integration scheme operators to use. More...
 
Array< OneD, Array< OneD, NekDouble > > m_previousSolution
 Storage for previous solution for steady-state check. More...
 
std::vector< int > m_intVariables
 
NekDouble m_cflSafetyFactor
 CFL safety factor (comprise between 0 to 1). More...
 
NekDouble m_CFLGrowth
 CFL growth rate. More...
 
NekDouble m_CFLEnd
 Maximun cfl in cfl growth. More...
 
int m_abortSteps
 Number of steps between checks for abort conditions. More...
 
bool m_explicitDiffusion
 Indicates if explicit or implicit treatment of diffusion is used. More...
 
bool m_explicitAdvection
 Indicates if explicit or implicit treatment of advection is used. More...
 
bool m_explicitReaction
 Indicates if explicit or implicit treatment of reaction is used. More...
 
int m_steadyStateSteps
 Check for steady state at step interval. More...
 
NekDouble m_steadyStateTol
 Tolerance to which steady state should be evaluated at. More...
 
int m_filtersInfosteps
 Number of time steps between outputting filters information. More...
 
std::vector< std::pair< std::string, FilterSharedPtr > > m_filters
 
bool m_homoInitialFwd
 Flag to determine if simulation should start in homogeneous forward transformed state. More...
 
std::ofstream m_errFile
 
NekDouble m_epsilon
 Diffusion coefficient. More...
 
- Protected Attributes inherited from Nektar::SolverUtils::EquationSystem
LibUtilities::CommSharedPtr m_comm
 Communicator. More...
 
bool m_verbose
 
LibUtilities::SessionReaderSharedPtr m_session
 The session reader. More...
 
std::map< std::string, SolverUtils::SessionFunctionSharedPtrm_sessionFunctions
 Map of known SessionFunctions. More...
 
LibUtilities::FieldIOSharedPtr m_fld
 Field input/output. More...
 
Array< OneD, MultiRegions::ExpListSharedPtrm_fields
 Array holding all dependent variables. More...
 
SpatialDomains::BoundaryConditionsSharedPtr m_boundaryConditions
 Pointer to boundary conditions object. More...
 
SpatialDomains::MeshGraphSharedPtr m_graph
 Pointer to graph defining mesh. More...
 
std::string m_sessionName
 Name of the session. More...
 
NekDouble m_time
 Current time of simulation. More...
 
int m_initialStep
 Number of the step where the simulation should begin. More...
 
NekDouble m_fintime
 Finish time of the simulation. More...
 
NekDouble m_timestep
 Time step size. More...
 
NekDouble m_lambda
 Lambda constant in real system if one required. More...
 
NekDouble m_checktime
 Time between checkpoints. More...
 
NekDouble m_lastCheckTime
 
NekDouble m_TimeIncrementFactor
 
int m_nchk
 Number of checkpoints written so far. More...
 
int m_steps
 Number of steps to take. More...
 
int m_checksteps
 Number of steps between checkpoints. More...
 
int m_infosteps
 Number of time steps between outputting status information. More...
 
int m_iterPIT = 0
 Number of parallel-in-time time iteration. More...
 
int m_windowPIT = 0
 Index of windows for parallel-in-time time iteration. More...
 
int m_spacedim
 Spatial dimension (>= expansion dim). More...
 
int m_expdim
 Expansion dimension. More...
 
bool m_singleMode
 Flag to determine if single homogeneous mode is used. More...
 
bool m_halfMode
 Flag to determine if half homogeneous mode is used. More...
 
bool m_multipleModes
 Flag to determine if use multiple homogenenous modes are used. More...
 
bool m_useFFT
 Flag to determine if FFT is used for homogeneous transform. More...
 
bool m_homogen_dealiasing
 Flag to determine if dealiasing is used for homogeneous simulations. More...
 
bool m_specHP_dealiasing
 Flag to determine if dealisising is usde for the Spectral/hp element discretisation. More...
 
enum MultiRegions::ProjectionType m_projectionType
 Type of projection; e.g continuous or discontinuous. More...
 
Array< OneD, Array< OneD, NekDouble > > m_traceNormals
 Array holding trace normals for DG simulations in the forwards direction. More...
 
Array< OneD, bool > m_checkIfSystemSingular
 Flag to indicate if the fields should be checked for singularity. More...
 
LibUtilities::FieldMetaDataMap m_fieldMetaDataMap
 Map to identify relevant solver info to dump in output fields. More...
 
Array< OneD, NekDoublem_movingFrameData
 Moving reference frame status in the inertial frame X, Y, Z, Theta_x, Theta_y, Theta_z, U, V, W, Omega_x, Omega_y, Omega_z, A_x, A_y, A_z, DOmega_x, DOmega_y, DOmega_z, pivot_x, pivot_y, pivot_z. More...
 
std::vector< std::string > m_strFrameData
 variable name in m_movingFrameData More...
 
int m_NumQuadPointsError
 Number of Quadrature points used to work out the error. More...
 
enum HomogeneousType m_HomogeneousType
 
NekDouble m_LhomX
 physical length in X direction (if homogeneous) More...
 
NekDouble m_LhomY
 physical length in Y direction (if homogeneous) More...
 
NekDouble m_LhomZ
 physical length in Z direction (if homogeneous) More...
 
int m_npointsX
 number of points in X direction (if homogeneous) More...
 
int m_npointsY
 number of points in Y direction (if homogeneous) More...
 
int m_npointsZ
 number of points in Z direction (if homogeneous) More...
 
int m_HomoDirec
 number of homogenous directions More...
 
- Protected Attributes inherited from Nektar::SolverUtils::ALEHelper
Array< OneD, MultiRegions::ExpListSharedPtrm_fieldsALE
 
Array< OneD, Array< OneD, NekDouble > > m_gridVelocity
 
Array< OneD, Array< OneD, NekDouble > > m_gridVelocityTrace
 
std::vector< ALEBaseShPtrm_ALEs
 
bool m_ALESolver = false
 
bool m_ImplicitALESolver = false
 
NekDouble m_prevStageTime = 0.0
 
int m_spaceDim
 

Friends

class MemoryManager< PulseWavePropagation >
 

Additional Inherited Members

- Protected Types inherited from Nektar::SolverUtils::EquationSystem
enum  HomogeneousType { eHomogeneous1D , eHomogeneous2D , eHomogeneous3D , eNotHomogeneous }
 Parameter for homogeneous expansions. More...
 
- Static Protected Attributes inherited from Nektar::SolverUtils::EquationSystem
static std::string equationSystemTypeLookupIds []
 
static std::string projectionTypeLookupIds []
 

Detailed Description

Set up the routines based on the weak formulation from "Computational Modelling of 1D blood flow with variable mechanical properties" by S. J. Sherwin et al. The weak formulation (1) reads: \( \sum_{e=1}^{N_{el}} \left[ \left( \frac{\partial \mathbf{U}^{\delta} }{\partial t} , \mathbf{\psi}^{\delta} \right)_{\Omega_e} - \left( \frac{\partial \mathbf{F(\mathbf{U})}^{\delta} } {\partial x}, \mathbf{\psi}^{\delta} \right)_{\Omega_e} + \left[ \mathbf{\psi}^{\delta} \cdot \{ \mathbf{F}^u - \mathbf{F}(\mathbf{U}^{\delta}) \} \right]_{x_e^l}^{x_eû} \right] = 0 \)

Definition at line 50 of file PulseWavePropagation.h.

Constructor & Destructor Documentation

◆ PulseWavePropagation()

Nektar::PulseWavePropagation::PulseWavePropagation ( const LibUtilities::SessionReaderSharedPtr pSession,
const SpatialDomains::MeshGraphSharedPtr pGraph 
)
protected

Definition at line 63 of file PulseWavePropagation.cpp.

66 : PulseWaveSystem(pSession, pGraph)
67{
68}
PulseWaveSystem(const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
Initialises PulseWaveSystem class members.

◆ ~PulseWavePropagation()

Nektar::PulseWavePropagation::~PulseWavePropagation ( )
overrideprotecteddefault

Member Function Documentation

◆ create()

static EquationSystemSharedPtr Nektar::PulseWavePropagation::create ( const LibUtilities::SessionReaderSharedPtr pSession,
const SpatialDomains::MeshGraphSharedPtr pGraph 
)
inlinestatic

Creates an instance of this class.

Definition at line 56 of file PulseWavePropagation.h.

59 {
62 pGraph);
63 p->InitObject();
64 return p;
65 }
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
std::shared_ptr< EquationSystem > EquationSystemSharedPtr
A shared pointer to an EquationSystem object.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and CellMLToNektar.cellml_metadata::p.

◆ DoOdeProjection()

void Nektar::PulseWavePropagation::DoOdeProjection ( const Array< OneD, const Array< OneD, NekDouble > > &  inarray,
Array< OneD, Array< OneD, NekDouble > > &  outarray,
const NekDouble  time 
)
protected

Definition at line 195 of file PulseWavePropagation.cpp.

199{
200 // Just copy over array
201 if (inarray != outarray)
202 {
203 for (size_t i = 0; i < m_nVariables; ++i)
204 {
205 Vmath::Vcopy(inarray[i].size(), inarray[i], 1, outarray[i], 1);
206 }
207 }
208}
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.hpp:825

References Nektar::PulseWaveSystem::m_nVariables, and Vmath::Vcopy().

Referenced by v_InitObject().

◆ DoOdeRhs()

void Nektar::PulseWavePropagation::DoOdeRhs ( const Array< OneD, const Array< OneD, NekDouble > > &  inarray,
Array< OneD, Array< OneD, NekDouble > > &  outarray,
const NekDouble  time 
)
protected

Computes the right hand side of (1). The RHS is everything except the term that contains the time derivative \(\frac{\partial \mathbf{U}}{\partial t}\). In case of a Discontinuous Galerkin projection, m_advObject->Advect will be called

Definition at line 142 of file PulseWavePropagation.cpp.

145{
146 size_t i;
147
148 Array<OneD, Array<OneD, NekDouble>> physarray(m_nVariables);
149
150 // Dummy array for WeakDG advection
151 Array<OneD, Array<OneD, NekDouble>> advVel(m_spacedim);
152
153 // Output array for advection
154 Array<OneD, Array<OneD, NekDouble>> out(m_nVariables);
155
156 size_t cnt = 0;
157
158 // Set up Inflow and Outflow boundary conditions.
159 SetPulseWaveBoundaryConditions(inarray, outarray, time);
160
161 // Set up any interface conditions and write into boundary condition
163
164 // do advection evaluation in all domains
165 for (size_t omega = 0; omega < m_nDomains; ++omega)
166 {
167 LibUtilities::Timer timer;
168 m_currentDomain = omega;
169 size_t nq = m_vessels[omega * m_nVariables]->GetTotPoints();
170
171 timer.Start();
172 for (i = 0; i < m_nVariables; ++i)
173 {
174 physarray[i] = inarray[i] + cnt; // note this is doing a hidden copy
175 out[i] = outarray[i] + cnt;
176 }
177
178 for (i = 0; i < m_nVariables; ++i)
179 {
180 m_fields[i] = m_vessels[omega * m_nVariables + i];
181 }
182
183 m_advObject->Advect(m_nVariables, m_fields, advVel, physarray, out,
184 time);
185 for (i = 0; i < m_nVariables; ++i)
186 {
187 Vmath::Neg(nq, out[i], 1);
188 }
189 timer.Stop();
190 timer.AccumulateRegion("PulseWavePropagation:_DoOdeRHS", 1);
191 cnt += nq;
192 }
193}
SolverUtils::AdvectionSharedPtr m_advObject
void SetPulseWaveBoundaryConditions(const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
Array< OneD, MultiRegions::ExpListSharedPtr > m_vessels
void EnforceInterfaceConditions(const Array< OneD, const Array< OneD, NekDouble > > &fields)
int m_spacedim
Spatial dimension (>= expansion dim).
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array holding all dependent variables.
void Neg(int n, T *x, const int incx)
Negate x = -x.
Definition: Vmath.hpp:292

References Nektar::LibUtilities::Timer::AccumulateRegion(), Nektar::PulseWaveSystem::EnforceInterfaceConditions(), m_advObject, Nektar::PulseWaveSystem::m_currentDomain, Nektar::SolverUtils::EquationSystem::m_fields, Nektar::PulseWaveSystem::m_nDomains, Nektar::PulseWaveSystem::m_nVariables, Nektar::SolverUtils::EquationSystem::m_spacedim, Nektar::PulseWaveSystem::m_vessels, Vmath::Neg(), SetPulseWaveBoundaryConditions(), Nektar::LibUtilities::Timer::Start(), and Nektar::LibUtilities::Timer::Stop().

Referenced by v_InitObject().

◆ GetA0()

Array< OneD, NekDouble > & Nektar::PulseWavePropagation::GetA0 ( )

Definition at line 380 of file PulseWavePropagation.cpp.

381{
383}
Array< OneD, Array< OneD, NekDouble > > m_A_0_trace

References Nektar::PulseWaveSystem::m_A_0_trace, and Nektar::PulseWaveSystem::m_currentDomain.

Referenced by v_InitObject().

◆ GetAlpha()

Array< OneD, NekDouble > & Nektar::PulseWavePropagation::GetAlpha ( )

Definition at line 390 of file PulseWavePropagation.cpp.

391{
393}
Array< OneD, Array< OneD, NekDouble > > m_alpha_trace

References Nektar::PulseWaveSystem::m_alpha_trace, and Nektar::PulseWaveSystem::m_currentDomain.

Referenced by v_InitObject().

◆ GetBeta()

Array< OneD, NekDouble > & Nektar::PulseWavePropagation::GetBeta ( )

Definition at line 385 of file PulseWavePropagation.cpp.

386{
388}
Array< OneD, Array< OneD, NekDouble > > m_beta_trace

References Nektar::PulseWaveSystem::m_beta_trace, and Nektar::PulseWaveSystem::m_currentDomain.

Referenced by v_InitObject().

◆ GetDomains()

NekDouble Nektar::PulseWavePropagation::GetDomains ( )

Definition at line 405 of file PulseWavePropagation.cpp.

406{
407 return m_nDomains;
408}

References Nektar::PulseWaveSystem::m_nDomains.

Referenced by v_InitObject().

◆ GetFluxVector()

void Nektar::PulseWavePropagation::GetFluxVector ( const Array< OneD, Array< OneD, NekDouble > > &  physfield,
Array< OneD, Array< OneD, Array< OneD, NekDouble > > > &  flux 
)
protected

DG Pulse Wave Propagation routines:

Calculates the second term of the weak form (1): \( \left( \frac{\partial \mathbf{F(\mathbf{U})}^{\delta} }{\partial x}, \mathbf{\psi}^{\delta} \right)_{\Omega_e} \) The variables of the system are $\mathbf{U} = [A,u]^T$ physfield[0] = A physfield[1] = u flux[0] = F[0] = A*u flux[1] = F[1] = u^2/2 + p/rho

Definition at line 304 of file PulseWavePropagation.cpp.

307{
308 size_t nq = m_vessels[m_currentDomain * m_nVariables]->GetTotPoints();
309 NekDouble domain = m_currentDomain;
310 m_pressure[domain] = Array<OneD, NekDouble>(nq);
311 Array<OneD, NekDouble> dAUdx(nq);
312 NekDouble viscoelasticGradient = 0.0;
313
314 LibUtilities::Timer timer;
315
316 timer.Start();
317 for (size_t j = 0; j < nq; ++j)
318 {
319 flux[0][0][j] = physfield[0][j] * physfield[1][j];
320 }
321 timer.Stop();
322 timer.AccumulateRegion("PulseWavePropagation:GetFluxVector-flux", 3);
323
324 // d/dx of AU, for the viscoelastic tube law and extra fields
325 m_fields[0]->PhysDeriv(flux[0][0], dAUdx);
326
327 for (size_t j = 0; j < nq; ++j)
328 {
329 if ((j == 0) || (j == nq - 1))
330 {
331 viscoelasticGradient = dAUdx[j];
332 }
333 else
334 {
335 viscoelasticGradient = (dAUdx[j] + dAUdx[j + 1]) / 2;
336 }
337
338 m_pressureArea->GetPressure(m_pressure[domain][j], m_beta[domain][j],
339 physfield[0][j], m_A_0[domain][j],
340 viscoelasticGradient, m_gamma[domain][j],
341 m_alpha[domain][j]);
342
343 flux[1][0][j] = physfield[1][j] * physfield[1][j] / 2 +
344 m_pressure[domain][j] / m_rho;
345 }
346
347 m_session->MatchSolverInfo("OutputExtraFields", "True", extraFields, true);
348
349 if (extraFields)
350 {
351 /*
352 Calculates wave speed and characteristic variables.
353
354 Ideally this should be moved to PulseWaveSystem, but it's easiest to
355 implement here.
356 */
357 size_t counter = 0;
358
359 m_PWV[domain] = Array<OneD, NekDouble>(nq);
360 m_W1[domain] = Array<OneD, NekDouble>(nq);
361 m_W2[domain] = Array<OneD, NekDouble>(nq);
362
363 for (size_t j = 0; j < nq; ++j)
364 {
365 m_pressureArea->GetC(m_PWV[domain][j], m_beta[domain][j],
366 physfield[0][counter + j], m_A_0[domain][j],
367 m_alpha[domain][j]);
368 m_pressureArea->GetW1(m_W1[domain][j], physfield[1][counter + j],
369 m_beta[domain][j], physfield[0][counter + j],
370 m_A_0[domain][j], m_alpha[domain][j]);
371 m_pressureArea->GetW2(m_W2[domain][j], physfield[1][counter + j],
372 m_beta[domain][j], physfield[0][counter + j],
373 m_A_0[domain][j], m_alpha[domain][j]);
374 }
375
376 counter += nq;
377 }
378}
Array< OneD, Array< OneD, NekDouble > > m_A_0
Array< OneD, Array< OneD, NekDouble > > m_W2
PulseWavePressureAreaSharedPtr m_pressureArea
Array< OneD, Array< OneD, NekDouble > > m_W1
Array< OneD, Array< OneD, NekDouble > > m_pressure
Array< OneD, Array< OneD, NekDouble > > m_gamma
Array< OneD, Array< OneD, NekDouble > > m_alpha
Array< OneD, Array< OneD, NekDouble > > m_PWV
Array< OneD, Array< OneD, NekDouble > > m_beta
LibUtilities::SessionReaderSharedPtr m_session
The session reader.
double NekDouble

References Nektar::LibUtilities::Timer::AccumulateRegion(), Nektar::PulseWaveSystem::extraFields, Nektar::PulseWaveSystem::m_A_0, Nektar::PulseWaveSystem::m_alpha, Nektar::PulseWaveSystem::m_beta, Nektar::PulseWaveSystem::m_currentDomain, Nektar::SolverUtils::EquationSystem::m_fields, Nektar::PulseWaveSystem::m_gamma, Nektar::PulseWaveSystem::m_nVariables, Nektar::PulseWaveSystem::m_pressure, Nektar::PulseWaveSystem::m_pressureArea, Nektar::PulseWaveSystem::m_PWV, Nektar::PulseWaveSystem::m_rho, Nektar::SolverUtils::EquationSystem::m_session, Nektar::PulseWaveSystem::m_vessels, Nektar::PulseWaveSystem::m_W1, Nektar::PulseWaveSystem::m_W2, Nektar::LibUtilities::Timer::Start(), and Nektar::LibUtilities::Timer::Stop().

Referenced by v_InitObject().

◆ GetN()

Array< OneD, NekDouble > & Nektar::PulseWavePropagation::GetN ( )

Definition at line 395 of file PulseWavePropagation.cpp.

396{
398}
Array< OneD, Array< OneD, NekDouble > > m_trace_fwd_normal

References Nektar::PulseWaveSystem::m_currentDomain, and Nektar::PulseWaveSystem::m_trace_fwd_normal.

Referenced by v_InitObject().

◆ GetRho()

NekDouble Nektar::PulseWavePropagation::GetRho ( )

Definition at line 400 of file PulseWavePropagation.cpp.

401{
402 return m_rho;
403}

References Nektar::PulseWaveSystem::m_rho.

Referenced by v_InitObject().

◆ SetPulseWaveBoundaryConditions()

void Nektar::PulseWavePropagation::SetPulseWaveBoundaryConditions ( const Array< OneD, const Array< OneD, NekDouble > > &  inarray,
Array< OneD, Array< OneD, NekDouble > > &  outarray,
const NekDouble  time 
)
protected
Does the projection between ... space and the ... space. Also checks for

Q-inflow boundary conditions at the inflow of the current arterial segment and applies the Q-inflow if specified

Definition at line 215 of file PulseWavePropagation.cpp.

220{
221 size_t omega;
222
223 Array<OneD, MultiRegions::ExpListSharedPtr> vessel(2);
224
225 size_t offset = 0;
226
227 // This will be moved to the RCR boundary condition once factory is setup
228 if (time == 0)
229 {
230 m_Boundary = Array<OneD, PulseWaveBoundarySharedPtr>(2 * m_nDomains);
231
232 for (omega = 0; omega < m_nDomains; ++omega)
233 {
234 vessel[0] = m_vessels[2 * omega];
235 vessel[1] = m_vessels[2 * omega + 1];
236
237 for (size_t j = 0; j < 2; ++j)
238 {
239 std::string BCType;
240
241 if (j < vessel[0]->GetBndConditions().size())
242 {
243 BCType = vessel[0]->GetBndConditions()[j]->GetUserDefined();
244 }
245
246 // if no condition given define it to be NoUserDefined
247 if (BCType.empty() || BCType == "Interface")
248 {
249 BCType = "NoUserDefined";
250 }
251
254
255 // turn on time dependent BCs
256 if (BCType == "Q-inflow")
257 {
258 vessel[0]->GetBndConditions()[j]->SetIsTimeDependent(true);
259 }
260 else if (BCType == "A-inflow")
261 {
262 vessel[0]->GetBndConditions()[j]->SetIsTimeDependent(true);
263 }
264 else if (BCType == "U-inflow")
265 {
266 vessel[1]->GetBndConditions()[j]->SetIsTimeDependent(true);
267 }
268 else if (BCType == "RCR-terminal")
269 {
270 vessel[0]->GetBndConditions()[j]->SetIsTimeDependent(true);
271 }
272 }
273 }
274 }
275
277
278 // Loop over all vessels and set boundary conditions
279 LibUtilities::Timer timer;
280 for (omega = 0; omega < m_nDomains; ++omega)
281 {
282 timer.Start();
283 for (size_t n = 0; n < 2; ++n)
284 {
285 m_Boundary[2 * omega + n]->DoBoundary(
286 inarray, m_A_0, m_beta, m_alpha, time, omega, offset, n);
287 }
288
289 offset += m_vessels[2 * omega]->GetTotPoints();
290 timer.Stop();
291 timer.AccumulateRegion("PulseWavePropagation:_SetBCs", 1);
292 }
293}
tBaseSharedPtr CreateInstance(tKey idKey, tParam... args)
Create an instance of the class referred to by idKey.
Array< OneD, PulseWaveBoundarySharedPtr > m_Boundary
SOLVER_UTILS_EXPORT void SetBoundaryConditions(NekDouble time)
Evaluates the boundary conditions at the given time.
BoundaryFactory & GetBoundaryFactory()

References Nektar::LibUtilities::Timer::AccumulateRegion(), Nektar::LibUtilities::NekFactory< tKey, tBase, tParam >::CreateInstance(), Nektar::GetBoundaryFactory(), Nektar::PulseWaveSystem::m_A_0, Nektar::PulseWaveSystem::m_alpha, Nektar::PulseWaveSystem::m_beta, m_Boundary, Nektar::PulseWaveSystem::m_nDomains, Nektar::PulseWaveSystem::m_pressureArea, Nektar::SolverUtils::EquationSystem::m_session, Nektar::PulseWaveSystem::m_vessels, Nektar::SolverUtils::EquationSystem::SetBoundaryConditions(), Nektar::LibUtilities::Timer::Start(), and Nektar::LibUtilities::Timer::Stop().

Referenced by DoOdeRhs().

◆ v_GenerateSummary()

void Nektar::PulseWavePropagation::v_GenerateSummary ( SolverUtils::SummaryList s)
overrideprotectedvirtual

Print summary routine, calls virtual routine reimplemented in UnsteadySystem

Reimplemented from Nektar::SolverUtils::UnsteadySystem.

Definition at line 414 of file PulseWavePropagation.cpp.

415{
417}
SOLVER_UTILS_EXPORT void v_GenerateSummary(SummaryList &s) override
Print a summary of time stepping parameters.

References Nektar::SolverUtils::UnsteadySystem::v_GenerateSummary().

◆ v_InitObject()

void Nektar::PulseWavePropagation::v_InitObject ( bool  DeclareField = false)
overrideprotectedvirtual

Init object for UnsteadySystem class.

Initialization object for UnsteadySystem class.

Reimplemented from Nektar::SolverUtils::UnsteadySystem.

Definition at line 70 of file PulseWavePropagation.cpp.

71{
72 // Will set up an array of vessels/fields in PulseWaveSystem::v_InitObject
73 // so set DeclareField to false so that the fields are not set up in
74 // EquationSystem unnecessarily. Note the number of fields in Equation
75 // system is related to the number of variables. The number of vessels is
76 // therefore held in PulwWaveSystem.
78
79 if (m_session->DefinesSolverInfo("PressureArea"))
80 {
82 m_session->GetSolverInfo("PressureArea"), m_vessels, m_session);
83 }
84 else
85 {
87 "Beta", m_vessels, m_session);
88 }
89
91 {
94 }
95 else
96 {
97 ASSERTL0(false, "Implicit Pulse Wave Propagation not set up.");
98 }
99
100 // Create advection object
101 std::string advName;
102 std::string riemName;
103 switch (m_upwindTypePulse)
104 {
105 case eUpwindPulse:
106 {
107 advName = "WeakDG";
108 riemName = "UpwindPulse";
109 }
110 break;
111 default:
112 {
113 ASSERTL0(false, "populate switch statement for upwind flux");
114 }
115 break;
116 }
121 riemName, m_session);
122 m_riemannSolver->SetScalar("A0", &PulseWavePropagation::GetA0, this);
123 m_riemannSolver->SetScalar("beta", &PulseWavePropagation::GetBeta, this);
124 m_riemannSolver->SetScalar("alpha", &PulseWavePropagation::GetAlpha, this);
125 m_riemannSolver->SetScalar("N", &PulseWavePropagation::GetN, this);
126 m_riemannSolver->SetParam("rho", &PulseWavePropagation::GetRho, this);
128 this);
129
130 m_advObject->SetRiemannSolver(m_riemannSolver);
131 m_advObject->InitObject(m_session, m_fields);
132}
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:208
void DefineProjection(FuncPointerT func, ObjectPointerT obj)
void DefineOdeRhs(FuncPointerT func, ObjectPointerT obj)
Array< OneD, NekDouble > & GetAlpha()
SolverUtils::RiemannSolverSharedPtr m_riemannSolver
void DoOdeProjection(const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
Array< OneD, NekDouble > & GetN()
Array< OneD, NekDouble > & GetA0()
void GetFluxVector(const Array< OneD, Array< OneD, NekDouble > > &physfield, Array< OneD, Array< OneD, Array< OneD, NekDouble > > > &flux)
DG Pulse Wave Propagation routines:
void DoOdeRhs(const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time)
Array< OneD, NekDouble > & GetBeta()
void v_InitObject(bool DeclareField=false) override
UpwindTypePulse m_upwindTypePulse
LibUtilities::TimeIntegrationSchemeOperators m_ode
The time integration scheme operators to use.
bool m_explicitAdvection
Indicates if explicit or implicit treatment of advection is used.
AdvectionFactory & GetAdvectionFactory()
Gets the factory for initialising advection objects.
Definition: Advection.cpp:43
RiemannSolverFactory & GetRiemannSolverFactory()
PressureAreaFactory & GetPressureAreaFactory()
@ eUpwindPulse
simple upwinding scheme

References ASSERTL0, Nektar::LibUtilities::NekFactory< tKey, tBase, tParam >::CreateInstance(), Nektar::LibUtilities::TimeIntegrationSchemeOperators::DefineOdeRhs(), Nektar::LibUtilities::TimeIntegrationSchemeOperators::DefineProjection(), DoOdeProjection(), DoOdeRhs(), Nektar::eUpwindPulse, GetA0(), Nektar::SolverUtils::GetAdvectionFactory(), GetAlpha(), GetBeta(), GetDomains(), GetFluxVector(), GetN(), Nektar::GetPressureAreaFactory(), GetRho(), Nektar::SolverUtils::GetRiemannSolverFactory(), m_advObject, Nektar::SolverUtils::UnsteadySystem::m_explicitAdvection, Nektar::SolverUtils::EquationSystem::m_fields, Nektar::SolverUtils::UnsteadySystem::m_ode, Nektar::PulseWaveSystem::m_pressureArea, m_riemannSolver, Nektar::SolverUtils::EquationSystem::m_session, Nektar::PulseWaveSystem::m_upwindTypePulse, Nektar::PulseWaveSystem::m_vessels, and Nektar::PulseWaveSystem::v_InitObject().

Friends And Related Function Documentation

◆ MemoryManager< PulseWavePropagation >

friend class MemoryManager< PulseWavePropagation >
friend

Definition at line 1 of file PulseWavePropagation.h.

Member Data Documentation

◆ className

std::string Nektar::PulseWavePropagation::className
static
Initial value:
=
"PulseWavePropagation", PulseWavePropagation::create,
"Pulse Wave Propagation equation.")
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
static EquationSystemSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
Creates an instance of this class.
EquationSystemFactory & GetEquationSystemFactory()

Name of class.

Definition at line 68 of file PulseWavePropagation.h.

◆ m_advObject

SolverUtils::AdvectionSharedPtr Nektar::PulseWavePropagation::m_advObject
protected

Definition at line 80 of file PulseWavePropagation.h.

Referenced by DoOdeRhs(), and v_InitObject().

◆ m_Boundary

Array<OneD, PulseWaveBoundarySharedPtr> Nektar::PulseWavePropagation::m_Boundary
protected

Definition at line 81 of file PulseWavePropagation.h.

Referenced by SetPulseWaveBoundaryConditions().

◆ m_riemannSolver

SolverUtils::RiemannSolverSharedPtr Nektar::PulseWavePropagation::m_riemannSolver
protected

Definition at line 79 of file PulseWavePropagation.h.

Referenced by v_InitObject().