#include <womersley_elements.h>
Inheritance diagram for oomph::WomersleyImpedanceTubeBase< ELEMENT, DIM >:
Public Types | |
| typedef double(* | PrescribedVolumeFluxFctPt) (const double &time) |
| Function pointer to fct that prescribes volume flux q=fct(t) – mainly used for validation purposes. More... | |
Public Member Functions | |
| WomersleyImpedanceTubeBase (const double &length, PrescribedVolumeFluxFctPt prescribed_volume_flux_fct_pt) | |
| Constructor: Specify length of tube and pointer to function that specifies the prescribed volume flux. Outlet pressure is set to zero. More... | |
| WomersleyImpedanceTubeBase (const double &length, Mesh *navier_stokes_outflow_mesh_pt) | |
| Constructor: Specify length of tube and the pointer to the mesh of either NavierStokesImpedanceTractionElements or NavierStokesFluxControlElements that are attached to the outflow cross-section of a (higher-dimensional) Navier Stokes mesh and provide the inflow into the ImpedanceTube. Outlet pressure is set to zero. More... | |
| double & | p_out () |
| Access fct to outlet pressure. More... | |
| virtual Mesh * | build_mesh_and_apply_boundary_conditions (TimeStepper *time_stepper_pt)=0 |
| Pure virtual function in which the user of a derived class must create the mesh of WomersleyElements (of the type specified by the class's template argument) and apply the boundary conditions. The Womersley elements use the timestepper specified as the input argument. More... | |
| void | setup () |
| Set up the Womersley tubes so that a subsequent call to get_response(...) computes the inlet pressure for the currently prescribed instantaneous flow rate. Steady version! More... | |
| void | setup (double *re_st_pt, const double &dt, const double &q_initial, TimeStepper *time_stepper_pt=0) |
| Set up the Womersley tubes so that a subsequent call to get_response(...) computes the inlet pressure for the currently prescribed instantaneous flow rate, assuming that at all previous times the tube conveyed steady, fully-developed flow with flowrate q_initial. dt specifies the timestep for the subsequent time integration. Specify: Womersley number, (constant) timestep, the initial volume flux (from which the subsequent impulsive start is performed) and, optionally the pointer to the timestepper to be used in the Womersley elements (defaults to BDF<2>). More... | |
| WomersleyProblem< ELEMENT, DIM > * | womersley_problem_pt () |
| Access to underlying Womersley problem. More... | |
| void | shift_time_values (const double &dt) |
| Shift history values to allow coputation of next timestep. Note: When used with a full Navier-Stokes problem this function must be called in actions_before_implicit_timestep() More... | |
| double | total_volume_flux_into_impedance_tube () |
| Compute total current volume flux into the "impedance tube" that provides the flow resistance (flux is either obtained from the function that specifies it externally or by by adding up the flux through all NavierStokesImpedanceTractionElements in the mesh pointed to by the Navier_stokes_outflow_mesh_pt. More... | |
| void | get_response (double &p_in, double &dp_in_dq) |
| Compute inlet pressure, p_in, required to achieve the currently imposed, instantaneous volume flux q prescribed by total_volume_flux_into_impedance_tube(), and its derivative, dp_in/dq. More... | |
| Public Member Functions inherited from oomph::TemplateFreeWomersleyImpedanceTubeBase | |
| TemplateFreeWomersleyImpedanceTubeBase () | |
| Empty constructor. More... | |
| virtual | ~TemplateFreeWomersleyImpedanceTubeBase () |
| Empty virtual destructor. More... | |
Protected Member Functions | |
| void | precompute_aux_integrals () |
| Precompute auxiliary integrals required for the computation of the Jacobian in the NavierStokesImpedanceTractionElement. Also pass the pointer to the pre-computed integrals to the elements in the Navier_stokes_outflow_mesh_pt so they can refer to it. More... | |
Protected Attributes | |
| double | Length |
| Length of the tube. More... | |
| double | Dp_in_dq |
| Derivative of inflow pressure w.r.t. instantaenous volume flux (Note: Can be pre-computed) More... | |
| double * | Current_volume_flux_pt |
| Pointer to double that specifies the currently imposed instantaneous volume flux into the impedance tube. This is used to communicate with the Womersley elements which require access to the flux via a pointer to a double. More... | |
| WomersleyProblem< ELEMENT, DIM > * | Womersley_problem_pt |
| Pointer to Womersley problem that determines the pressure gradient along the tube. More... | |
| double | P_out |
| Outlet pressure. More... | |
| PrescribedVolumeFluxFctPt | Prescribed_volume_flux_fct_pt |
| Pointer to function that specifies the prescribed volume flux. More... | |
| Mesh * | Navier_stokes_outflow_mesh_pt |
| Pointer to the mesh of NavierStokesImpedanceTractionElements that are attached to the outflow cross-section of the higher-dimensional Navier Stokes mesh and provide the inflow into the Impedance tube. More... | |
| std::map< unsigned, double > * | Aux_integral_pt |
| Pointer to auxiliary integral, containing the derivative of the total volume flux through the outflow boundary of the (higher-dimensional) Navier-Stokes mesh w.r.t. to the discrete (global) (velocity) degrees of freedom. More... | |
Private Attributes | |
| bool | Using_flux_control_elements |
Additional Inherited Members | |
| Static Public Attributes inherited from oomph::TemplateFreeWomersleyImpedanceTubeBase | |
| static double | Zero =0.0 |
| Zero! More... | |
Detailed Description
template<class ELEMENT, unsigned DIM>
class oomph::WomersleyImpedanceTubeBase< ELEMENT, DIM >
Base class for Womersley impedance tube. Allows the computation of the inlet pressure p_in into a uniform tube of specified length that is assumed to convey fully-developed, but time-dependent flow with a presribed instantaneous flow rate, q. Also computes the derivative dp_in/dq required when this is used to determine impedance-type outlet boundary conditions in a Navier-Stokes computation.
Definition at line 1120 of file womersley_elements.h.
Member Typedef Documentation
◆ PrescribedVolumeFluxFctPt
| typedef double(* oomph::WomersleyImpedanceTubeBase< ELEMENT, DIM >::PrescribedVolumeFluxFctPt) (const double &time) |
Function pointer to fct that prescribes volume flux q=fct(t) – mainly used for validation purposes.
Definition at line 1128 of file womersley_elements.h.
Constructor & Destructor Documentation
◆ WomersleyImpedanceTubeBase() [1/2]
|
inline |
Constructor: Specify length of tube and pointer to function that specifies the prescribed volume flux. Outlet pressure is set to zero.
Definition at line 1133 of file womersley_elements.h.
◆ WomersleyImpedanceTubeBase() [2/2]
|
inline |
Constructor: Specify length of tube and the pointer to the mesh of either NavierStokesImpedanceTractionElements or NavierStokesFluxControlElements that are attached to the outflow cross-section of a (higher-dimensional) Navier Stokes mesh and provide the inflow into the ImpedanceTube. Outlet pressure is set to zero.
Definition at line 1156 of file womersley_elements.h.
Member Function Documentation
◆ build_mesh_and_apply_boundary_conditions()
|
pure virtual |
Pure virtual function in which the user of a derived class must create the mesh of WomersleyElements (of the type specified by the class's template argument) and apply the boundary conditions. The Womersley elements use the timestepper specified as the input argument.
Implemented in oomph::WomersleyOutflowImpedanceTube< ELEMENT, DIM >.
◆ get_response()
|
inlinevirtual |
Compute inlet pressure, p_in, required to achieve the currently imposed, instantaneous volume flux q prescribed by total_volume_flux_into_impedance_tube(), and its derivative, dp_in/dq.
Implements oomph::TemplateFreeWomersleyImpedanceTubeBase.
Definition at line 1428 of file womersley_elements.h.
◆ p_out()
|
inline |
Access fct to outlet pressure.
Definition at line 1224 of file womersley_elements.h.
◆ precompute_aux_integrals()
|
inlineprotected |
Precompute auxiliary integrals required for the computation of the Jacobian in the NavierStokesImpedanceTractionElement. Also pass the pointer to the pre-computed integrals to the elements in the Navier_stokes_outflow_mesh_pt so they can refer to it.
Definition at line 1454 of file womersley_elements.h.
◆ setup() [1/2]
|
inline |
Set up the Womersley tubes so that a subsequent call to get_response(...) computes the inlet pressure for the currently prescribed instantaneous flow rate. Steady version!
Definition at line 1238 of file womersley_elements.h.
◆ setup() [2/2]
|
inline |
Set up the Womersley tubes so that a subsequent call to get_response(...) computes the inlet pressure for the currently prescribed instantaneous flow rate, assuming that at all previous times the tube conveyed steady, fully-developed flow with flowrate q_initial. dt specifies the timestep for the subsequent time integration. Specify: Womersley number, (constant) timestep, the initial volume flux (from which the subsequent impulsive start is performed) and, optionally the pointer to the timestepper to be used in the Womersley elements (defaults to BDF<2>).
By default, we do want to suppress the output from the Newton solver
Definition at line 1258 of file womersley_elements.h.
◆ shift_time_values()
|
inline |
Shift history values to allow coputation of next timestep. Note: When used with a full Navier-Stokes problem this function must be called in actions_before_implicit_timestep()
Definition at line 1364 of file womersley_elements.h.
◆ total_volume_flux_into_impedance_tube()
|
inline |
Compute total current volume flux into the "impedance tube" that provides the flow resistance (flux is either obtained from the function that specifies it externally or by by adding up the flux through all NavierStokesImpedanceTractionElements in the mesh pointed to by the Navier_stokes_outflow_mesh_pt.
Definition at line 1390 of file womersley_elements.h.
◆ womersley_problem_pt()
|
inline |
Access to underlying Womersley problem.
Definition at line 1354 of file womersley_elements.h.
Member Data Documentation
◆ Aux_integral_pt
|
protected |
Pointer to auxiliary integral, containing the derivative of the total volume flux through the outflow boundary of the (higher-dimensional) Navier-Stokes mesh w.r.t. to the discrete (global) (velocity) degrees of freedom.
Definition at line 1520 of file womersley_elements.h.
◆ Current_volume_flux_pt
|
protected |
Pointer to double that specifies the currently imposed instantaneous volume flux into the impedance tube. This is used to communicate with the Womersley elements which require access to the flux via a pointer to a double.
Definition at line 1499 of file womersley_elements.h.
◆ Dp_in_dq
|
protected |
Derivative of inflow pressure w.r.t. instantaenous volume flux (Note: Can be pre-computed)
Definition at line 1493 of file womersley_elements.h.
◆ Length
|
protected |
Length of the tube.
Definition at line 1489 of file womersley_elements.h.
◆ Navier_stokes_outflow_mesh_pt
|
protected |
Pointer to the mesh of NavierStokesImpedanceTractionElements that are attached to the outflow cross-section of the higher-dimensional Navier Stokes mesh and provide the inflow into the Impedance tube.
Definition at line 1514 of file womersley_elements.h.
◆ P_out
|
protected |
Outlet pressure.
Definition at line 1506 of file womersley_elements.h.
◆ Prescribed_volume_flux_fct_pt
|
protected |
Pointer to function that specifies the prescribed volume flux.
Definition at line 1509 of file womersley_elements.h.
◆ Using_flux_control_elements
|
private |
Definition at line 1526 of file womersley_elements.h.
◆ Womersley_problem_pt
|
protected |
Pointer to Womersley problem that determines the pressure gradient along the tube.
Definition at line 1503 of file womersley_elements.h.
The documentation for this class was generated from the following file: