IBAMR  IBAMR version 0.19.
Public Member Functions | Protected Member Functions | Static Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | List of all members
IBAMR::ConstraintIBMethod Class Referenceabstract

Class ConstraintIBMethod implements the rigidity constraint for rigid and deforming bodies using the constraint based IB method. More...

#include <ibamr/ConstraintIBMethod.h>

Inheritance diagram for IBAMR::ConstraintIBMethod:
Inheritance graph
[legend]

Public Member Functions

 ConstraintIBMethod (std::string object_name, SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > input_db, const int no_structures, bool register_for_restart=true)
 Constructor. More...
 
 ~ConstraintIBMethod ()
 Destructor. More...
 
void initializeHierarchyOperatorsandData ()
 Initialize Hierarchy operators and data at initial time. More...
 
virtual void registerEulerianVariables () override
 Register Eulerian variables with base IBStrategy class. More...
 
virtual void preprocessIntegrateData (double current_time, double new_time, int num_cycles) override
 Create Lagrangian workspace. More...
 
virtual void postprocessIntegrateData (double current_time, double new_time, int num_cycles) override
 Destroy Lagrangian workspace. More...
 
void registerConstraintIBKinematics (const std::vector< SAMRAI::tbox::Pointer< IBAMR::ConstraintIBKinematics > > &ib_kinematics_op)
 Register kinematics of the immersed structure(s) with this class. More...
 
void registerPreProcessSolveFluidEquationsCallBackFunction (void(*ptr_preprocess_callbackfnc)(const double, const double, const int, void *), void *ctx)
 Register any preprocess fluid solve callback functions. More...
 
virtual void preprocessSolveFluidEquations (double current_time, double new_time, int cycle_num) override
 Calculate any body forces for INS solver over here. More...
 
void registerPostProcessSolveFluidEquationsCallBackFunction (void(*ptr_postprocess_callbackfnc)(const double, const double, const int, void *), void *ctx)
 Register any postprocess fluid solve callback functions. More...
 
virtual void postprocessSolveFluidEquations (double current_time, double new_time, int cycle_num) override
 Apply the FuRMoRP algorithm in the postprocessSolveFluidEquations method. More...
 
virtual void forwardEulerStep (double current_time, double new_time) override
 Override the forwardEulerStep method of the base IBMethod class. More...
 
virtual void midpointStep (double current_time, double new_time) override
 Override the midpointStep method of the base IBMethod class. More...
 
virtual void putToDatabase (SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > db) override
 Override the putToDatabase method of the base Serializable class. More...
 
const std::vector< double > & getVolumeElement ()
 Get the volume element associated with material points of all structures. More...
 
const std::vector< std::vector< double > > & getCurrentCOMVelocity ()
 Get the current COM velocity associated with Lagrangian structures. More...
 
const std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > & getLagrangeMultiplier ()
 Get LData associated with Lagrange multiplier force field. More...
 
const std::vector< std::vector< double > > & getCurrentStructureCOM ()
 Get the current center of mass for all Lagrangian structures. More...
 
void setVelocityPhysBdryOp (IBTK::RobinPhysBdryPatchStrategy *u_phys_bdry_op)
 
void setVolumeElement (double vol_element, int struct_no)
 
void setVolumeElement (std::vector< double > vol_element)
 
const std::vector< double > & getStructureVolume ()
 
const std::vector< std::vector< double > > & getStructureMomentum ()
 Get the total linear momentum for all the Lagrangian structures. More...
 
const std::vector< std::vector< double > > & getStructureRotationalMomentum ()
 Get the total rotational momentum for all the Lagrangian structures with respect to their COM. More...
 
void registerIBLagrangianForceFunction (SAMRAI::tbox::Pointer< IBLagrangianForceStrategy > ib_force_fcn)
 
void registerIBLagrangianSourceFunction (SAMRAI::tbox::Pointer< IBLagrangianSourceStrategy > ib_source_fcn)
 
void registerLInitStrategy (SAMRAI::tbox::Pointer< IBTK::LInitStrategy > l_initializer)
 
void freeLInitStrategy ()
 
void registerIBMethodPostProcessor (SAMRAI::tbox::Pointer< IBMethodPostProcessStrategy > post_processor)
 
IBTK::LDataManagergetLDataManager () const
 
SAMRAI::tbox::Pointer< IBInstrumentPanelgetIBInstrumentPanel () const
 
void registerLSiloDataWriter (SAMRAI::tbox::Pointer< IBTK::LSiloDataWriter > silo_writer)
 
const SAMRAI::hier::IntVector< NDIM > & getMinimumGhostCellWidth () const override
 
void setupTagBuffer (SAMRAI::tbox::Array< int > &tag_buffer, SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > > gridding_alg) const override
 
virtual void inactivateLagrangianStructure (int structure_number=0, int level_number=std::numeric_limits< int >::max()) override
 
virtual void activateLagrangianStructure (int structure_number=0, int level_number=std::numeric_limits< int >::max()) override
 
virtual bool getLagrangianStructureIsActivated (int structure_number=0, int level_number=std::numeric_limits< int >::max()) const override
 
void createSolverVecs (Vec *X_vec, Vec *F_vec) override
 
void setupSolverVecs (Vec *X_vec, Vec *F_vec) override
 
void setUpdatedPosition (Vec &X_new_vec) override
 
void setLinearizedPosition (Vec &X_vec, double data_time) override
 
void computeResidual (Vec &R_vec) override
 
void computeLinearizedResidual (Vec &X_vec, Vec &R_vec) override
 
void updateFixedLEOperators () override
 
void interpolateVelocity (int u_data_idx, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > &u_synch_scheds, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &u_ghost_fill_scheds, double data_time) override
 
void interpolateLinearizedVelocity (int u_data_idx, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > &u_synch_scheds, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &u_ghost_fill_scheds, double data_time) override
 
void backwardEulerStep (double current_time, double new_time) override
 
void trapezoidalStep (double current_time, double new_time) override
 
void computeLagrangianForce (double data_time) override
 
void computeLinearizedLagrangianForce (Vec &X_vec, double data_time) override
 
void constructLagrangianForceJacobian (Mat &A, MatType mat_type, double data_time) override
 
void spreadForce (int f_data_idx, IBTK::RobinPhysBdryPatchStrategy *f_phys_bdry_op, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &f_prolongation_scheds, double data_time) override
 
void spreadLinearizedForce (int f_data_idx, IBTK::RobinPhysBdryPatchStrategy *f_phys_bdry_op, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &f_prolongation_scheds, double data_time) override
 
void constructInterpOp (Mat &J, void(*spread_fnc)(const double, double *), int stencil_width, const std::vector< int > &num_dofs_per_proc, int dof_index_idx, double data_time) override
 
bool hasFluidSources () const override
 
void computeLagrangianFluidSource (double data_time) override
 
void spreadFluidSource (int q_data_idx, IBTK::RobinPhysBdryPatchStrategy *q_phys_bdry_op, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &q_prolongation_scheds, double data_time) override
 
void interpolatePressure (int p_data_idx, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > &p_synch_scheds, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &p_ghost_fill_scheds, double data_time) override
 
void postprocessData () override
 
void initializePatchHierarchy (SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > > hierarchy, SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > > gridding_alg, int u_data_idx, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > &u_synch_scheds, const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &u_ghost_fill_scheds, int integrator_step, double init_data_time, bool initial_time) override
 
void registerLoadBalancer (SAMRAI::tbox::Pointer< SAMRAI::mesh::LoadBalancer< NDIM > > load_balancer, int workload_data_idx) override
 
void addWorkloadEstimate (SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > > hierarchy, const int workload_data_idx) override
 
void beginDataRedistribution (SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > > hierarchy, SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > > gridding_alg) override
 
void endDataRedistribution (SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > > hierarchy, SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > > gridding_alg) override
 
void initializeLevelData (SAMRAI::tbox::Pointer< SAMRAI::hier::BasePatchHierarchy< NDIM > > hierarchy, int level_number, double init_data_time, bool can_be_refined, bool initial_time, SAMRAI::tbox::Pointer< SAMRAI::hier::BasePatchLevel< NDIM > > old_level, bool allocate_data) override
 
virtual void initializeLevelData (const tbox::Pointer< hier::BasePatchHierarchy< DIM > > hierarchy, const int level_number, const double init_data_time, const bool can_be_refined, const bool initial_time, const tbox::Pointer< hier::BasePatchLevel< DIM > > old_level=tbox::Pointer< hier::BasePatchLevel< DIM > >(NULL), const bool allocate_data=true)=0
 
void resetHierarchyConfiguration (SAMRAI::tbox::Pointer< SAMRAI::hier::BasePatchHierarchy< NDIM > > hierarchy, int coarsest_level, int finest_level) override
 
virtual void resetHierarchyConfiguration (const tbox::Pointer< hier::BasePatchHierarchy< DIM > > hierarchy, const int coarsest_level, const int finest_level)=0
 
void applyGradientDetector (SAMRAI::tbox::Pointer< SAMRAI::hier::BasePatchHierarchy< NDIM > > hierarchy, int level_number, double error_data_time, int tag_index, bool initial_time, bool uses_richardson_extrapolation_too) override
 
virtual void applyGradientDetector (const tbox::Pointer< hier::BasePatchHierarchy< DIM > > hierarchy, const int level_number, const double error_data_time, const int tag_index, const bool initial_time, const bool uses_richardson_extrapolation_too)
 
double convertTimeEnumToDouble (IBTK::TimePoint time_pt)
 Convert the enum TimePoint to it respective value. More...
 
void getPositionData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **X_data, bool **X_needs_ghost_fill, IBTK::TimePoint time_pt)
 Get the structure position data at the specified time point. More...
 
void getVelocityData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **U_data, IBTK::TimePoint time_pt)
 Get the current structure velocity data at the specified time point. More...
 
void getForceData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **F_data, bool **F_needs_ghost_fill, IBTK::TimePoint time_pt)
 Get the current structure force data at the specified time point. More...
 
virtual void registerIBHierarchyIntegrator (IBHierarchyIntegrator *ib_solver)
 
virtual void registerEulerianCommunicationAlgorithms ()
 
virtual double getMaxPointDisplacement () const
 
void setUseFixedLEOperators (bool use_fixed_coupling_ops=true)
 
virtual void setUseMultistepTimeStepping (unsigned int n_previous_steps=1)
 
virtual void AB2Step (double current_time, double new_time)
 
virtual double getLevelDt (const tbox::Pointer< hier::BasePatchLevel< DIM > > level, const double dt_time, const bool initial_time)
 
virtual double advanceLevel (const tbox::Pointer< hier::BasePatchLevel< DIM > > level, const tbox::Pointer< hier::BasePatchHierarchy< DIM > > hierarchy, const double current_time, const double new_time, const bool first_step, const bool last_step, const bool regrid_advance=false)
 
virtual void resetTimeDependentData (const tbox::Pointer< hier::BasePatchLevel< DIM > > level, const double new_time, const bool can_be_refined)
 
virtual void resetDataToPreadvanceState (const tbox::Pointer< hier::BasePatchLevel< DIM > > level)
 
virtual void applyRichardsonExtrapolation (const tbox::Pointer< hier::PatchLevel< DIM > > level, const double error_data_time, const int tag_index, const double deltat, const int error_coarsen_ratio, const bool initial_time, const bool uses_gradient_detector_too)
 
virtual void coarsenDataForRichardsonExtrapolation (const tbox::Pointer< hier::PatchHierarchy< DIM > > hierarchy, const int level_number, const tbox::Pointer< hier::PatchLevel< DIM > > coarser_level, const double coarsen_data_time, const bool before_advance)
 

Protected Member Functions

void getPositionData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **X_data, bool **X_needs_ghost_fill, double data_time)
 
void getVelocityData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **U_data, double data_time)
 
void getForceData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **F_data, bool **F_needs_ghost_fill, double data_time)
 
void getLinearizedPositionData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **X_data, bool **X_needs_ghost_fill)
 
void getLECouplingPositionData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **X_LE_data, bool **X_LE_needs_ghost_fill, double data_time)
 
void getLinearizedVelocityData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **U_data)
 
void getLinearizedForceData (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **F_data, bool **F_needs_ghost_fill)
 
void reinitMidpointData (const std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > &current_data, const std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > &new_data, const std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > &half_data)
 
void resetAnchorPointValues (std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > U_data, int coarsest_ln, int finest_ln)
 
PetscErrorCode computeForce (Vec X, Vec F)
 
INSHierarchyIntegratorgetINSHierarchyIntegrator () const
 
SAMRAI::tbox::Pointer< SAMRAI::math::HierarchyDataOpsReal< NDIM, double > > getVelocityHierarchyDataOps () const
 
SAMRAI::tbox::Pointer< SAMRAI::math::HierarchyDataOpsReal< NDIM, double > > getPressureHierarchyDataOps () const
 
SAMRAI::tbox::Pointer< IBTK::HierarchyMathOpsgetHierarchyMathOps () const
 
void registerVariable (int &current_idx, int &new_idx, int &scratch_idx, SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > > variable, const SAMRAI::hier::IntVector< NDIM > &scratch_ghosts=SAMRAI::hier::IntVector< NDIM >(0), const std::string &coarsen_name="NO_COARSEN", const std::string &refine_name="NO_REFINE", SAMRAI::tbox::Pointer< IBTK::CartGridFunction > init_fcn=nullptr, const bool register_for_restart=true)
 
void registerVariable (int &idx, SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > > variable, const SAMRAI::hier::IntVector< NDIM > &ghosts=SAMRAI::hier::IntVector< NDIM >(0), SAMRAI::tbox::Pointer< SAMRAI::hier::VariableContext > ctx=SAMRAI::tbox::Pointer< SAMRAI::hier::VariableContext >(nullptr), const bool register_for_restart=true)
 
void registerGhostfillRefineAlgorithm (const std::string &name, SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineAlgorithm< NDIM > > ghostfill_alg, std::unique_ptr< SAMRAI::xfer::RefinePatchStrategy< NDIM > > ghostfill_patch_strategy=nullptr)
 
void registerProlongRefineAlgorithm (const std::string &name, SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineAlgorithm< NDIM > > prolong_alg, std::unique_ptr< SAMRAI::xfer::RefinePatchStrategy< NDIM > > prolong_patch_strategy=nullptr)
 
void registerCoarsenAlgorithm (const std::string &name, SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenAlgorithm< NDIM > > coarsen_alg, std::unique_ptr< SAMRAI::xfer::CoarsenPatchStrategy< NDIM > > coarsen_patch_strategy=nullptr)
 
SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineAlgorithm< NDIM > > getGhostfillRefineAlgorithm (const std::string &name) const
 
SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineAlgorithm< NDIM > > getProlongRefineAlgorithm (const std::string &name) const
 
SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenAlgorithm< NDIM > > getCoarsenAlgorithm (const std::string &name) const
 
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > & getGhostfillRefineSchedules (const std::string &name) const
 
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > & getProlongRefineSchedules (const std::string &name) const
 
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > & getCoarsenSchedules (const std::string &name) const
 

Static Protected Member Functions

static PetscErrorCode computeForce_SAMRAI (void *ctx, Vec X, Vec F)
 

Protected Attributes

bool d_do_log = false
 
SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > > d_hierarchy
 
SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > > d_gridding_alg
 
double d_current_time = std::numeric_limits<double>::quiet_NaN()
 
double d_new_time = std::numeric_limits<double>::quiet_NaN()
 
double d_half_time = std::numeric_limits<double>::quiet_NaN()
 
bool d_X_current_needs_ghost_fill = true
 
bool d_X_new_needs_ghost_fill = true
 
bool d_X_half_needs_ghost_fill = true
 
bool d_X_jac_needs_ghost_fill = true
 
bool d_X_LE_new_needs_ghost_fill = true
 
bool d_X_LE_half_needs_ghost_fill = true
 
bool d_F_current_needs_ghost_fill = true
 
bool d_F_new_needs_ghost_fill = true
 
bool d_F_half_needs_ghost_fill = true
 
bool d_F_jac_needs_ghost_fill = true
 
IBTK::LDataManagerd_l_data_manager
 
std::string d_interp_kernel_fcn = "IB_4"
 
std::string d_spread_kernel_fcn = "IB_4"
 
bool d_error_if_points_leave_domain = false
 
SAMRAI::hier::IntVector< NDIM > d_ghosts
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_X_current_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_X_new_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_X_half_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_X_jac_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_X_LE_new_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_X_LE_half_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_U_current_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_U_new_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_U_half_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_U_jac_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_F_current_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_F_new_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_F_half_data
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_F_jac_data
 
std::vector< std::set< int > > d_anchor_point_local_idxs
 
SAMRAI::tbox::Pointer< IBInstrumentPaneld_instrument_panel
 
std::vector< doubled_total_flow_volume
 
SAMRAI::tbox::Pointer< IBTK::LInitStrategyd_l_initializer
 
SAMRAI::tbox::Pointer< IBLagrangianForceStrategyd_ib_force_fcn
 
bool d_ib_force_fcn_needs_init = true
 
SAMRAI::tbox::Pointer< IBLagrangianSourceStrategyd_ib_source_fcn
 
bool d_ib_source_fcn_needs_init = true
 
std::vector< std::vector< IBTK::Point > > d_X_src
 
std::vector< std::vector< double > > d_r_src
 
std::vector< std::vector< double > > d_P_src
 
std::vector< std::vector< double > > d_Q_src
 
std::vector< intd_n_src
 
bool d_normalize_source_strength = false
 
SAMRAI::tbox::Pointer< IBMethodPostProcessStrategyd_post_processor
 
SAMRAI::tbox::Pointer< IBTK::LSiloDataWriterd_silo_writer
 
SAMRAI::tbox::Pointer< SAMRAI::mesh::LoadBalancer< NDIM > > d_load_balancer
 
int d_workload_idx = IBTK::invalid_index
 
std::string d_object_name
 
bool d_registered_for_restart
 
IBHierarchyIntegratord_ib_solver = nullptr
 
bool d_use_fixed_coupling_ops = false
 

Private Member Functions

 ConstraintIBMethod ()=delete
 Default constructor. More...
 
 ConstraintIBMethod (const ConstraintIBMethod &from)=delete
 Default copy constructor. More...
 
ConstraintIBMethodoperator= (const ConstraintIBMethod &that)=delete
 Default assignment operator. More...
 
void getFromInput (SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > input_db, const bool from_restart)
 Get values from input file. More...
 
void getFromRestart ()
 Get values from restart file. More...
 
void setInitialLagrangianVelocity ()
 Set initial Lagrangian velocity on material points. More...
 
void calculateCOMandMOIOfStructures ()
 Calculate center of mass and moment of inertia of immersed structures. More...
 
void calculateKinematicsVelocity ()
 Calculate the kinematics velocity for all structures handled by this class. More...
 
void calculateMomentumOfKinematicsVelocity (const int position_handle)
 Calculate momentum of kinematics velocity. This is extraneous momentum that needs to be subtracted from the kinematics velocity. More...
 
void calculateVolumeElement ()
 Calculate volume element associated with material points. More...
 
void setCounter ()
 Set the counter for this method. More...
 
void setFuRMoRPTime (const double current_time, const double new_time)
 Set the time at which FuRMoRP is applied. More...
 
void copyFluidVariable (int copy_from_idx, int copy_to)
 Copy vector. More...
 
void copyDensityVariable (int copy_from_idx, int copy_to)
 Copy density patch data. More...
 
void interpolateFluidSolveVelocity ()
 Interpolate fluid solve velocity from Eulerian grid onto the Lagrangian mesh. More...
 
void calculateRigidTranslationalMomentum ()
 Calculate the rigid translational velocity. More...
 
void calculateRigidRotationalMomentum ()
 Calculate the rigid rotational velocity. More...
 
void calculateCurrentLagrangianVelocity ()
 Calculate current velocity on the material points. More...
 
void correctVelocityOnLagrangianMesh ()
 Correct velocity on Lagrangian mesh. Set the velocity on Lagrangian mesh as U_lag_corr = U_trans + Omega X r + U_def - U_interpolated. More...
 
void spreadCorrectedLagrangianVelocity ()
 Spread the corrected velocity at the Lagrangian mesh to the Eulerian Grid. More...
 
void applyProjection ()
 The correction on Eulerian grid can lead to non-divergence free velocity field. We project the corrected velocity field onto a divergence free field. More...
 
void updateStructurePositionEulerStep ()
 Predict the position of structures according to forward Euler step method. More...
 
void updateStructurePositionMidPointStep ()
 Update the position of structures according to mid point step method. More...
 
void calculateMidPointVelocity ()
 Compute U_half = 0.5(U_current + U_new);. More...
 
void calculateDrag ()
 Calculate hydrodynamic drag on the immersed structures. More...
 
void calculateTorque ()
 Calculate hydrodynamic torque on the immersed structures. More...
 
void calculatePower ()
 Calculate power spent during swimming. More...
 
void calculateEulerianMomentum ()
 Calculate Eulerian Momentum. More...
 
void calculateStructureMomentum ()
 Calculate total translational momentum of all Lagrangian structures. More...
 
void calculateStructureRotationalMomentum ()
 Calculate the total rotational momentum of all Lagrangian structures with respect to their COM. More...
 
void resetLagrangianForceFunction (double init_data_time, bool initial_time)
 
void resetLagrangianSourceFunction (double init_data_time, bool initial_time)
 
void updateIBInstrumentationData (int timestep_num, double data_time)
 

Private Attributes

const int d_no_structures
 
std::vector< SAMRAI::tbox::Pointer< IBAMR::ConstraintIBKinematics > > d_ib_kinematics
 
double d_FuRMoRP_current_time = 0.0
 
double d_FuRMoRP_new_time = 0.0
 
std::vector< doubled_vol_element
 
std::vector< boold_vol_element_is_set
 
std::vector< doubled_structure_vol
 
std::vector< std::vector< double > > d_structure_mom
 
std::vector< std::vector< double > > d_structure_rotational_mom
 
bool d_needs_div_free_projection = false
 
std::vector< std::vector< double > > d_rigid_trans_vel_current
 
std::vector< std::vector< double > > d_rigid_trans_vel_new
 
std::vector< std::vector< double > > d_rigid_rot_vel_current
 
std::vector< std::vector< double > > d_rigid_rot_vel_new
 
std::vector< std::vector< double > > d_incremented_angle_from_reference_axis
 
std::vector< std::vector< double > > d_vel_com_def_current
 
std::vector< std::vector< double > > d_vel_com_def_new
 
std::vector< std::vector< double > > d_omega_com_def_current
 
std::vector< std::vector< double > > d_omega_com_def_new
 
std::vector< std::vector< double > > d_center_of_mass_current
 
std::vector< std::vector< double > > d_center_of_mass_new
 
std::vector< std::vector< double > > d_center_of_mass_unshifted_current
 
std::vector< std::vector< double > > d_center_of_mass_unshifted_new
 
std::vector< Eigen::Matrix3d > d_moment_of_inertia_current
 
std::vector< Eigen::Matrix3d > d_moment_of_inertia_new
 
std::vector< intd_tagged_pt_lag_idx
 
std::vector< std::vector< double > > d_tagged_pt_position
 
std::vector< doubled_rho_solid
 
double d_rho_fluid = std::numeric_limits<double>::quiet_NaN()
 
bool d_rho_is_const = true
 
bool d_calculate_structure_linear_mom = false
 
bool d_calculate_structure_rotational_mom = false
 
int d_timestep_counter = 0
 
int d_output_interval = 1
 
bool d_print_output = false
 
bool d_output_drag = false
 
bool d_output_torque = false
 
bool d_output_power = false
 
bool d_output_trans_vel = false
 
bool d_output_rot_vel = false
 
bool d_output_COM_coordinates = false
 
bool d_output_MOI = false
 
bool d_output_eul_mom = false
 
std::string d_dir_name = "./ConstraintIBMethodDump"
 
std::string d_base_output_filename = "ImmersedStructure"
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_l_data_U_interp
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_l_data_U_correction
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_l_data_U_new
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_l_data_U_current
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_l_data_U_half
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_l_data_X_half_Euler
 
std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > d_l_data_X_new_MidPoint
 
SAMRAI::tbox::Pointer< SAMRAI::math::HierarchySideDataOpsReal< NDIM, double > > d_hier_sc_data_ops
 
SAMRAI::tbox::Pointer< SAMRAI::math::HierarchyCellDataOpsReal< NDIM, double > > d_hier_cc_data_ops
 
SAMRAI::tbox::Pointer< IBTK::HierarchyGhostCellInterpolationd_no_fill_op
 
int d_wgt_cc_idx
 
int d_wgt_sc_idx
 
double d_volume
 
SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > > d_u_var
 
SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > > d_u_fluidSolve_var
 
SAMRAI::tbox::Pointer< SAMRAI::pdat::CellVariable< NDIM, double > > d_phi_var
 
SAMRAI::tbox::Pointer< SAMRAI::pdat::CellVariable< NDIM, double > > d_Div_u_var
 
SAMRAI::tbox::Pointer< SAMRAI::hier::VariableContextd_scratch_context
 
int d_u_fluidSolve_idx
 
int d_u_fluidSolve_cib_idx
 
int d_phi_idx
 
int d_Div_u_scratch_idx
 
SAMRAI::tbox::Pointer< SAMRAI::pdat::SideVariable< NDIM, double > > d_rho_var
 
int d_rho_ins_idx = IBTK::invalid_index
 
int d_rho_scratch_idx = IBTK::invalid_index
 
SAMRAI::solv::LocationIndexRobinBcCoefs< NDIM > d_velcorrection_projection_bc_coef
 
std::unique_ptr< SAMRAI::solv::PoissonSpecificationsd_velcorrection_projection_spec
 
SAMRAI::tbox::Pointer< IBTK::CCLaplaceOperatord_velcorrection_projection_op
 
SAMRAI::tbox::Pointer< IBTK::PETScKrylovPoissonSolverd_velcorrection_projection_solver
 
SAMRAI::tbox::Pointer< IBTK::CCPoissonPointRelaxationFACOperatord_velcorrection_projection_fac_op
 
SAMRAI::tbox::Pointer< SAMRAI::tbox::Databased_velcorrection_projection_fac_pc_db
 
SAMRAI::tbox::Pointer< IBTK::FACPreconditionerd_velcorrection_projection_fac_pc
 
std::vector< std::unique_ptr< std::ofstream > > d_trans_vel_stream
 
std::vector< std::unique_ptr< std::ofstream > > d_rot_vel_stream
 
std::vector< std::unique_ptr< std::ofstream > > d_drag_force_stream
 
std::vector< std::unique_ptr< std::ofstream > > d_moment_of_inertia_stream
 
std::vector< std::unique_ptr< std::ofstream > > d_torque_stream
 
std::vector< std::unique_ptr< std::ofstream > > d_position_COM_stream
 
std::vector< std::unique_ptr< std::ofstream > > d_power_spent_stream
 
std::fstream d_eulerian_mom_stream
 
std::vector< void(*)(const double, const double, const int, void *)> d_prefluidsolve_callback_fns
 
std::vector< void(*)(const double, const double, const int, void *)> d_postfluidsolve_callback_fns
 
std::vector< void * > d_prefluidsolve_callback_fns_ctx
 
std::vector< void * > d_postfluidsolve_callback_fns_ctx
 
IBTK::RobinPhysBdryPatchStrategyd_u_phys_bdry_op = nullptr
 
bool d_force_jac_mffd = false
 
Mat d_force_jac = nullptr
 
double d_force_jac_data_time
 

Detailed Description

References Bhalla et al. A unified mathematical framework and an adaptive numerical method for fluid-structure interaction with rigid, deforming, and elastic bodies. J Comput Phys, 250:446-476 (2013).

Constructor & Destructor Documentation

◆ ConstraintIBMethod() [1/3]

IBAMR::ConstraintIBMethod::ConstraintIBMethod ( std::string  object_name,
SAMRAI::tbox::Pointer< SAMRAI::tbox::Database input_db,
const int  no_structures,
bool  register_for_restart = true 
)

◆ ~ConstraintIBMethod()

IBAMR::ConstraintIBMethod::~ConstraintIBMethod ( )

◆ ConstraintIBMethod() [2/3]

IBAMR::ConstraintIBMethod::ConstraintIBMethod ( )
privatedelete
Note
This constructor is not implemented and should not be used.

◆ ConstraintIBMethod() [3/3]

IBAMR::ConstraintIBMethod::ConstraintIBMethod ( const ConstraintIBMethod from)
privatedelete
Note
This copy constructor is not implemented and should not be used.

Member Function Documentation

◆ initializeHierarchyOperatorsandData()

void IBAMR::ConstraintIBMethod::initializeHierarchyOperatorsandData ( )

◆ registerEulerianVariables()

virtual void IBAMR::ConstraintIBMethod::registerEulerianVariables ( )
overridevirtual

Reimplemented from IBAMR::IBStrategy.

◆ preprocessIntegrateData()

virtual void IBAMR::ConstraintIBMethod::preprocessIntegrateData ( double  current_time,
double  new_time,
int  num_cycles 
)
overridevirtual

Reimplemented from IBAMR::IBStrategy.

◆ postprocessIntegrateData()

virtual void IBAMR::ConstraintIBMethod::postprocessIntegrateData ( double  current_time,
double  new_time,
int  num_cycles 
)
overridevirtual

Reimplemented from IBAMR::IBStrategy.

◆ registerConstraintIBKinematics()

void IBAMR::ConstraintIBMethod::registerConstraintIBKinematics ( const std::vector< SAMRAI::tbox::Pointer< IBAMR::ConstraintIBKinematics > > &  ib_kinematics_op)

◆ registerPreProcessSolveFluidEquationsCallBackFunction()

void IBAMR::ConstraintIBMethod::registerPreProcessSolveFluidEquationsCallBackFunction ( void(*)(const double, const double, const int, void *)  ptr_preprocess_callbackfnc,
void *  ctx 
)
inline

◆ preprocessSolveFluidEquations()

virtual void IBAMR::ConstraintIBMethod::preprocessSolveFluidEquations ( double  current_time,
double  new_time,
int  cycle_num 
)
overridevirtual

Reimplemented from IBAMR::IBStrategy.

◆ registerPostProcessSolveFluidEquationsCallBackFunction()

void IBAMR::ConstraintIBMethod::registerPostProcessSolveFluidEquationsCallBackFunction ( void(*)(const double, const double, const int, void *)  ptr_postprocess_callbackfnc,
void *  ctx 
)
inline

◆ postprocessSolveFluidEquations()

virtual void IBAMR::ConstraintIBMethod::postprocessSolveFluidEquations ( double  current_time,
double  new_time,
int  cycle_num 
)
overridevirtual

Reimplemented from IBAMR::IBStrategy.

◆ forwardEulerStep()

virtual void IBAMR::ConstraintIBMethod::forwardEulerStep ( double  current_time,
double  new_time 
)
overridevirtual

Implements IBAMR::IBStrategy.

◆ midpointStep()

virtual void IBAMR::ConstraintIBMethod::midpointStep ( double  current_time,
double  new_time 
)
overridevirtual

Implements IBAMR::IBStrategy.

◆ putToDatabase()

virtual void IBAMR::ConstraintIBMethod::putToDatabase ( SAMRAI::tbox::Pointer< SAMRAI::tbox::Database db)
overridevirtual

◆ getVolumeElement()

const std::vector<double>& IBAMR::ConstraintIBMethod::getVolumeElement ( )
inline

◆ getCurrentCOMVelocity()

const std::vector<std::vector<double> >& IBAMR::ConstraintIBMethod::getCurrentCOMVelocity ( )
inline

◆ getLagrangeMultiplier()

const std::vector<SAMRAI::tbox::Pointer<IBTK::LData> >& IBAMR::ConstraintIBMethod::getLagrangeMultiplier ( )
inline

◆ getCurrentStructureCOM()

const std::vector<std::vector<double> >& IBAMR::ConstraintIBMethod::getCurrentStructureCOM ( )
inline

◆ setVelocityPhysBdryOp()

void IBAMR::ConstraintIBMethod::setVelocityPhysBdryOp ( IBTK::RobinPhysBdryPatchStrategy u_phys_bdry_op)
inline

◆ setVolumeElement() [1/2]

void IBAMR::ConstraintIBMethod::setVolumeElement ( double  vol_element,
int  struct_no 
)
inline

◆ setVolumeElement() [2/2]

void IBAMR::ConstraintIBMethod::setVolumeElement ( std::vector< double vol_element)
inline

◆ getStructureVolume()

const std::vector<double>& IBAMR::ConstraintIBMethod::getStructureVolume ( )
inline

◆ getStructureMomentum()

const std::vector<std::vector<double> >& IBAMR::ConstraintIBMethod::getStructureMomentum ( )
inline

◆ getStructureRotationalMomentum()

const std::vector<std::vector<double> >& IBAMR::ConstraintIBMethod::getStructureRotationalMomentum ( )
inline

◆ operator=()

ConstraintIBMethod& IBAMR::ConstraintIBMethod::operator= ( const ConstraintIBMethod that)
privatedelete
Note
This assignment operator is not implemented and should not be used.

◆ getFromInput()

void IBAMR::ConstraintIBMethod::getFromInput ( SAMRAI::tbox::Pointer< SAMRAI::tbox::Database input_db,
const bool  from_restart 
)
private

◆ getFromRestart()

void IBAMR::ConstraintIBMethod::getFromRestart ( )
private

◆ setInitialLagrangianVelocity()

void IBAMR::ConstraintIBMethod::setInitialLagrangianVelocity ( )
private

◆ calculateCOMandMOIOfStructures()

void IBAMR::ConstraintIBMethod::calculateCOMandMOIOfStructures ( )
private

◆ calculateKinematicsVelocity()

void IBAMR::ConstraintIBMethod::calculateKinematicsVelocity ( )
private

◆ calculateMomentumOfKinematicsVelocity()

void IBAMR::ConstraintIBMethod::calculateMomentumOfKinematicsVelocity ( const int  position_handle)
private

◆ calculateVolumeElement()

void IBAMR::ConstraintIBMethod::calculateVolumeElement ( )
private

◆ setCounter()

void IBAMR::ConstraintIBMethod::setCounter ( )
inlineprivate

◆ setFuRMoRPTime()

void IBAMR::ConstraintIBMethod::setFuRMoRPTime ( const double  current_time,
const double  new_time 
)
inlineprivate

◆ copyFluidVariable()

void IBAMR::ConstraintIBMethod::copyFluidVariable ( int  copy_from_idx,
int  copy_to 
)
private

◆ copyDensityVariable()

void IBAMR::ConstraintIBMethod::copyDensityVariable ( int  copy_from_idx,
int  copy_to 
)
private

◆ interpolateFluidSolveVelocity()

void IBAMR::ConstraintIBMethod::interpolateFluidSolveVelocity ( )
private

◆ calculateRigidTranslationalMomentum()

void IBAMR::ConstraintIBMethod::calculateRigidTranslationalMomentum ( )
private

◆ calculateRigidRotationalMomentum()

void IBAMR::ConstraintIBMethod::calculateRigidRotationalMomentum ( )
private

◆ calculateCurrentLagrangianVelocity()

void IBAMR::ConstraintIBMethod::calculateCurrentLagrangianVelocity ( )
private

◆ correctVelocityOnLagrangianMesh()

void IBAMR::ConstraintIBMethod::correctVelocityOnLagrangianMesh ( )
private

◆ spreadCorrectedLagrangianVelocity()

void IBAMR::ConstraintIBMethod::spreadCorrectedLagrangianVelocity ( )
private

◆ applyProjection()

void IBAMR::ConstraintIBMethod::applyProjection ( )
private

◆ updateStructurePositionEulerStep()

void IBAMR::ConstraintIBMethod::updateStructurePositionEulerStep ( )
private

◆ updateStructurePositionMidPointStep()

void IBAMR::ConstraintIBMethod::updateStructurePositionMidPointStep ( )
private

◆ calculateMidPointVelocity()

void IBAMR::ConstraintIBMethod::calculateMidPointVelocity ( )
private

◆ calculateDrag()

void IBAMR::ConstraintIBMethod::calculateDrag ( )
private

◆ calculateTorque()

void IBAMR::ConstraintIBMethod::calculateTorque ( )
private

◆ calculatePower()

void IBAMR::ConstraintIBMethod::calculatePower ( )
private

◆ calculateEulerianMomentum()

void IBAMR::ConstraintIBMethod::calculateEulerianMomentum ( )
private

◆ calculateStructureMomentum()

void IBAMR::ConstraintIBMethod::calculateStructureMomentum ( )
private

◆ calculateStructureRotationalMomentum()

void IBAMR::ConstraintIBMethod::calculateStructureRotationalMomentum ( )
private

◆ registerIBLagrangianForceFunction()

void IBAMR::IBMethod::registerIBLagrangianForceFunction ( SAMRAI::tbox::Pointer< IBLagrangianForceStrategy ib_force_fcn)
inherited

Supply a Lagrangian force object.

◆ registerIBLagrangianSourceFunction()

void IBAMR::IBMethod::registerIBLagrangianSourceFunction ( SAMRAI::tbox::Pointer< IBLagrangianSourceStrategy ib_source_fcn)
inherited

Supply a Lagrangian source object.

◆ registerLInitStrategy()

void IBAMR::IBMethod::registerLInitStrategy ( SAMRAI::tbox::Pointer< IBTK::LInitStrategy l_initializer)
inherited

Supply a Lagrangian initialization object.

◆ freeLInitStrategy()

void IBAMR::IBMethod::freeLInitStrategy ( )
inherited

Free references to Lagrangian initialization objects.

◆ registerIBMethodPostProcessor()

void IBAMR::IBMethod::registerIBMethodPostProcessor ( SAMRAI::tbox::Pointer< IBMethodPostProcessStrategy post_processor)
inherited

Supply a post processor object.

◆ getLDataManager()

IBTK::LDataManager* IBAMR::IBMethod::getLDataManager ( ) const
inherited

Return a pointer to the Lagrangian data manager object.

◆ getIBInstrumentPanel()

SAMRAI::tbox::Pointer<IBInstrumentPanel> IBAMR::IBMethod::getIBInstrumentPanel ( ) const
inherited

Return a pointer to the instrumentation manager object.

◆ registerLSiloDataWriter()

void IBAMR::IBMethod::registerLSiloDataWriter ( SAMRAI::tbox::Pointer< IBTK::LSiloDataWriter silo_writer)
inherited

Register a Lagrangian Silo data writer so this class will write plot files that may be postprocessed with the VisIt visualization tool.

◆ getMinimumGhostCellWidth()

const SAMRAI::hier::IntVector<NDIM>& IBAMR::IBMethod::getMinimumGhostCellWidth ( ) const
overridevirtualinherited

Return the number of ghost cells required by the Lagrangian-Eulerian interaction routines.

Implements IBAMR::IBStrategy.

◆ setupTagBuffer()

void IBAMR::IBMethod::setupTagBuffer ( SAMRAI::tbox::Array< int > &  tag_buffer,
SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > >  gridding_alg 
) const
overridevirtualinherited

Setup the tag buffer.

Reimplemented from IBAMR::IBStrategy.

◆ inactivateLagrangianStructure()

virtual void IBAMR::IBMethod::inactivateLagrangianStructure ( int  structure_number = 0,
int  level_number = std::numeric_limits< int >::max() 
)
overridevirtualinherited

Inactivate a structure/part. See IBAMR::IBStrategy::inactivateLagrangianStructure().

Reimplemented from IBAMR::IBStrategy.

◆ activateLagrangianStructure()

virtual void IBAMR::IBMethod::activateLagrangianStructure ( int  structure_number = 0,
int  level_number = std::numeric_limits< int >::max() 
)
overridevirtualinherited

Activate a previously inactivated structure/part to be used again in FSI calculations. See IBAMR::IBStrategy::activateLagrangianStructure().

Reimplemented from IBAMR::IBStrategy.

◆ getLagrangianStructureIsActivated()

virtual bool IBAMR::IBMethod::getLagrangianStructureIsActivated ( int  structure_number = 0,
int  level_number = std::numeric_limits< int >::max() 
) const
overridevirtualinherited

Determine whether or not the given structure or part is currently activated. See IBAMR::IBStrategy::getLagrangianStructureIsActivated().

Reimplemented from IBAMR::IBStrategy.

◆ createSolverVecs()

void IBAMR::IBMethod::createSolverVecs ( Vec *  X_vec,
Vec *  F_vec 
)
overridevirtualinherited

Create solution and rhs data on the specified level of the patch hierarchy.

Implements IBAMR::IBImplicitStrategy.

◆ setupSolverVecs()

void IBAMR::IBMethod::setupSolverVecs ( Vec *  X_vec,
Vec *  F_vec 
)
overridevirtualinherited

Setup solution and rhs data on the specified level of the patch hierarchy.

Implements IBAMR::IBImplicitStrategy.

◆ setUpdatedPosition()

void IBAMR::IBMethod::setUpdatedPosition ( Vec &  X_new_vec)
overridevirtualinherited

Set the value of the updated position vector.

Implements IBAMR::IBImplicitStrategy.

◆ setLinearizedPosition()

void IBAMR::IBMethod::setLinearizedPosition ( Vec &  X_vec,
double  data_time 
)
overridevirtualinherited

Set the value of the intermediate position vector used in evaluating the linearized problem.

Implements IBAMR::IBImplicitStrategy.

◆ computeResidual()

void IBAMR::IBMethod::computeResidual ( Vec &  R_vec)
overridevirtualinherited

Compute the residual on the specified level of the patch hierarchy.

Implements IBAMR::IBImplicitStrategy.

◆ computeLinearizedResidual()

void IBAMR::IBMethod::computeLinearizedResidual ( Vec &  X_vec,
Vec &  R_vec 
)
overridevirtualinherited

Compute the linearized residual for the given intermediate position vector.

Implements IBAMR::IBImplicitStrategy.

◆ updateFixedLEOperators()

void IBAMR::IBMethod::updateFixedLEOperators ( )
overridevirtualinherited

Update the positions used for the "fixed" interpolation and spreading operators.

Reimplemented from IBAMR::IBStrategy.

◆ interpolateVelocity()

void IBAMR::IBMethod::interpolateVelocity ( int  u_data_idx,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > &  u_synch_scheds,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &  u_ghost_fill_scheds,
double  data_time 
)
overridevirtualinherited

Interpolate the Eulerian velocity to the curvilinear mesh at the specified time within the current time interval.

Implements IBAMR::IBStrategy.

◆ interpolateLinearizedVelocity()

void IBAMR::IBMethod::interpolateLinearizedVelocity ( int  u_data_idx,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > &  u_synch_scheds,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &  u_ghost_fill_scheds,
double  data_time 
)
overridevirtualinherited

Interpolate the Eulerian velocity to the curvilinear mesh at the specified time within the current time interval for use in evaluating the residual of the linearized problem.

Implements IBAMR::IBImplicitStrategy.

◆ backwardEulerStep()

void IBAMR::IBMethod::backwardEulerStep ( double  current_time,
double  new_time 
)
overridevirtualinherited

Advance the positions of the Lagrangian structure using the backward Euler method.

Reimplemented from IBAMR::IBStrategy.

◆ trapezoidalStep()

void IBAMR::IBMethod::trapezoidalStep ( double  current_time,
double  new_time 
)
overridevirtualinherited

Advance the positions of the Lagrangian structure using the trapezoidal rule.

Implements IBAMR::IBStrategy.

Reimplemented in IBAMR::PenaltyIBMethod.

◆ computeLagrangianForce()

void IBAMR::IBMethod::computeLagrangianForce ( double  data_time)
overridevirtualinherited

Compute the Lagrangian force at the specified time within the current time interval.

Implements IBAMR::IBStrategy.

Reimplemented in IBAMR::PenaltyIBMethod.

◆ computeLinearizedLagrangianForce()

void IBAMR::IBMethod::computeLinearizedLagrangianForce ( Vec &  X_vec,
double  data_time 
)
overridevirtualinherited

Compute the Lagrangian force of the linearized problem for the specified configuration of the updated position vector.

Implements IBAMR::IBImplicitStrategy.

◆ constructLagrangianForceJacobian()

void IBAMR::IBMethod::constructLagrangianForceJacobian ( Mat &  A,
MatType  mat_type,
double  data_time 
)
overridevirtualinherited

Construct the linearized Lagrangian force Jacobian.

Implements IBAMR::IBImplicitStrategy.

◆ spreadForce()

void IBAMR::IBMethod::spreadForce ( int  f_data_idx,
IBTK::RobinPhysBdryPatchStrategy f_phys_bdry_op,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &  f_prolongation_scheds,
double  data_time 
)
overridevirtualinherited

Spread the Lagrangian force to the Cartesian grid at the specified time within the current time interval.

Implements IBAMR::IBStrategy.

◆ spreadLinearizedForce()

void IBAMR::IBMethod::spreadLinearizedForce ( int  f_data_idx,
IBTK::RobinPhysBdryPatchStrategy f_phys_bdry_op,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &  f_prolongation_scheds,
double  data_time 
)
overridevirtualinherited

Spread the Lagrangian force of the linearized problem to the Cartesian grid at the specified time within the current time interval.

Implements IBAMR::IBImplicitStrategy.

◆ constructInterpOp()

void IBAMR::IBMethod::constructInterpOp ( Mat &  J,
void(*)(const double, double *)  spread_fnc,
int  stencil_width,
const std::vector< int > &  num_dofs_per_proc,
int  dof_index_idx,
double  data_time 
)
overridevirtualinherited

Construct the IB interpolation operator.

Implements IBAMR::IBImplicitStrategy.

◆ hasFluidSources()

bool IBAMR::IBMethod::hasFluidSources ( ) const
overridevirtualinherited

Indicate whether there are any internal fluid sources/sinks.

Reimplemented from IBAMR::IBStrategy.

◆ computeLagrangianFluidSource()

void IBAMR::IBMethod::computeLagrangianFluidSource ( double  data_time)
overridevirtualinherited

Compute the Lagrangian source/sink density at the specified time within the current time interval.

Reimplemented from IBAMR::IBStrategy.

◆ spreadFluidSource()

void IBAMR::IBMethod::spreadFluidSource ( int  q_data_idx,
IBTK::RobinPhysBdryPatchStrategy q_phys_bdry_op,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &  q_prolongation_scheds,
double  data_time 
)
overridevirtualinherited

Spread the Lagrangian source/sink density to the Cartesian grid at the specified time within the current time interval.

Reimplemented from IBAMR::IBStrategy.

◆ interpolatePressure()

void IBAMR::IBMethod::interpolatePressure ( int  p_data_idx,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > &  p_synch_scheds,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &  p_ghost_fill_scheds,
double  data_time 
)
overridevirtualinherited

Compute the pressures at the positions of any distributed internal fluid sources or sinks.

Reimplemented from IBAMR::IBStrategy.

◆ postprocessData()

void IBAMR::IBMethod::postprocessData ( )
overridevirtualinherited

Execute user-defined post-processing operations.

Reimplemented from IBAMR::IBStrategy.

◆ initializePatchHierarchy()

void IBAMR::IBMethod::initializePatchHierarchy ( SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > >  hierarchy,
SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > >  gridding_alg,
int  u_data_idx,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenSchedule< NDIM > > > &  u_synch_scheds,
const std::vector< SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineSchedule< NDIM > > > &  u_ghost_fill_scheds,
int  integrator_step,
double  init_data_time,
bool  initial_time 
)
overridevirtualinherited

Initialize Lagrangian data corresponding to the given AMR patch hierarchy at the start of a computation. If the computation is begun from a restart file, data may be read from the restart databases.

A patch data descriptor is provided for the Eulerian velocity in case initialization requires interpolating Eulerian data. Ghost cells for Eulerian data will be filled upon entry to this function.

Reimplemented from IBAMR::IBStrategy.

Reimplemented in IBAMR::PenaltyIBMethod.

◆ registerLoadBalancer()

void IBAMR::IBMethod::registerLoadBalancer ( SAMRAI::tbox::Pointer< SAMRAI::mesh::LoadBalancer< NDIM > >  load_balancer,
int  workload_data_idx 
)
overridevirtualinherited

Register a load balancer and work load patch data index with the IB strategy object.

Deprecated:
This method is no longer necessary with the current workload estimation scheme.

Reimplemented from IBAMR::IBStrategy.

◆ addWorkloadEstimate()

void IBAMR::IBMethod::addWorkloadEstimate ( SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > >  hierarchy,
const int  workload_data_idx 
)
overridevirtualinherited

Add the estimated computational work from the current object per cell into the specified workload_data_idx.

Reimplemented from IBAMR::IBStrategy.

◆ beginDataRedistribution()

void IBAMR::IBMethod::beginDataRedistribution ( SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > >  hierarchy,
SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > >  gridding_alg 
)
overridevirtualinherited

Begin redistributing Lagrangian data prior to regridding the patch hierarchy.

Reimplemented from IBAMR::IBStrategy.

◆ endDataRedistribution()

void IBAMR::IBMethod::endDataRedistribution ( SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > >  hierarchy,
SAMRAI::tbox::Pointer< SAMRAI::mesh::GriddingAlgorithm< NDIM > >  gridding_alg 
)
overridevirtualinherited

Complete redistributing Lagrangian data following regridding the patch hierarchy.

Reimplemented from IBAMR::IBStrategy.

◆ initializeLevelData() [1/2]

void IBAMR::IBMethod::initializeLevelData ( SAMRAI::tbox::Pointer< SAMRAI::hier::BasePatchHierarchy< NDIM > >  hierarchy,
int  level_number,
double  init_data_time,
bool  can_be_refined,
bool  initial_time,
SAMRAI::tbox::Pointer< SAMRAI::hier::BasePatchLevel< NDIM > >  old_level,
bool  allocate_data 
)
overrideinherited

Initialize data on a new level after it is inserted into an AMR patch hierarchy by the gridding algorithm.

See also
SAMRAI::mesh::StandardTagAndInitStrategy::initializeLevelData

◆ initializeLevelData() [2/2]

virtual void SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::initializeLevelData ( const tbox::Pointer< hier::BasePatchHierarchy< DIM > >  hierarchy,
const int  level_number,
const double  init_data_time,
const bool  can_be_refined,
const bool  initial_time,
const tbox::Pointer< hier::BasePatchLevel< DIM > >  old_level = tbox::Pointerhier::BasePatchLevel<DIM> >(NULL),
const bool  allocate_data = true 
)
pure virtualinherited

Initialize data on a new level after it is inserted into an AMR patch hierarchy by the gridding algorithm. The level number indicates that of the new level.

Generally, when data is set, it is interpolated from coarser levels in the hierarchy. If the old level pointer in the argument list is non-null, then data is copied from the old level to the new level on regions of intersection between those levels before interpolation occurs. In this case, the level number must match that of the old level. The specific operations that occur when initializing level data are determined by the particular solution methods in use; i.e., in the subclass of this abstract base class.

The boolean argument initial_time indicates whether the level is being introduced for the first time (i.e., at initialization time), or after some regrid process during the calculation beyond the initial hierarchy construction. This information is provided since the initialization of the data may be different in each of those circumstances. The can_be_refined boolean argument indicates whether the level is the finest allowable level in the hierarchy.

◆ resetHierarchyConfiguration() [1/2]

void IBAMR::IBMethod::resetHierarchyConfiguration ( SAMRAI::tbox::Pointer< SAMRAI::hier::BasePatchHierarchy< NDIM > >  hierarchy,
int  coarsest_level,
int  finest_level 
)
overrideinherited

Reset cached hierarchy dependent data.

See also
SAMRAI::mesh::StandardTagAndInitStrategy::resetHierarchyConfiguration

◆ resetHierarchyConfiguration() [2/2]

virtual void SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::resetHierarchyConfiguration ( const tbox::Pointer< hier::BasePatchHierarchy< DIM > >  hierarchy,
const int  coarsest_level,
const int  finest_level 
)
pure virtualinherited

After hierarchy levels have changed and data has been initialized on the new levels, this routine can be used to reset any information needed by the solution method that is particular to the hierarchy configuration. For example, the solution procedure may cache communication schedules to amortize the cost of data movement on the AMR patch hierarchy. This function will be called by the gridding algorithm after the initialization occurs so that the algorithm-specific subclass can reset such things. Also, if the solution method must make the solution consistent across multiple levels after the hierarchy is changed, this process may be invoked by this routine. Of course the details of these processes are determined by the particular solution methods in use.

The level number arguments indicate the coarsest and finest levels in the current hierarchy configuration that have changed. It should be assumed that all intermediate levels have changed as well.

◆ applyGradientDetector() [1/2]

void IBAMR::IBMethod::applyGradientDetector ( SAMRAI::tbox::Pointer< SAMRAI::hier::BasePatchHierarchy< NDIM > >  hierarchy,
int  level_number,
double  error_data_time,
int  tag_index,
bool  initial_time,
bool  uses_richardson_extrapolation_too 
)
overrideinherited

Set integer tags to "one" in cells where refinement of the given level should occur according to user-supplied feature detection criteria.

See also
SAMRAI::mesh::StandardTagAndInitStrategy::applyGradientDetector

◆ applyGradientDetector() [2/2]

virtual void SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::applyGradientDetector ( const tbox::Pointer< hier::BasePatchHierarchy< DIM > >  hierarchy,
const int  level_number,
const double  error_data_time,
const int  tag_index,
const bool  initial_time,
const bool  uses_richardson_extrapolation_too 
)
virtualinherited

Set integer tags to "one" in cells where refinement of the given level should occur according to some user-supplied gradient criteria. The double time argument is the regrid time. The integer "tag_index" argument is the patch descriptor index of the cell-centered integer tag array on each patch in the hierarchy. The boolean argument initial_time indicates whether the level is being subject to refinement at the initial simulation time. If it is false, then the error estimation process is being invoked at some later time after the AMR hierarchy was initially constructed. Typically, this information is passed to the user's patch tagging routines since the error estimator or gradient detector may be different in each case.

The boolean uses_richardson_extrapolation_too is true when Richardson extrapolation error estimation is used in addition to the gradient detector, and false otherwise. This argument helps the user to manage multiple regridding criteria.

This routine is only when gradient detector is being used. It is virtual with an empty implementation here (rather than pure virtual) so that users are not required to provide an implementation when the function is not needed.

◆ convertTimeEnumToDouble()

double IBAMR::IBMethod::convertTimeEnumToDouble ( IBTK::TimePoint  time_pt)
inherited

If TimePoint is not one of CURRENT_TIME, HALF_TIME, or NEW_TIME, this returns NaN.

◆ getPositionData() [1/2]

void IBAMR::IBMethod::getPositionData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  X_data,
bool **  X_needs_ghost_fill,
IBTK::TimePoint  time_pt 
)
inherited

The time point should be one of CURRENT_TIME, HALF_TIME, or NEW_TIME. If this condition is met, X_data is set to the data at that respective time, otherwise the X_data pointers are unchanged.

◆ getPositionData() [2/2]

void IBAMR::IBMethod::getPositionData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  X_data,
bool **  X_needs_ghost_fill,
double  data_time 
)
protectedinherited

Get the current structure position data.

data_time must be equal to one of current time, new time, or half time. If this condition is met, X_data is set to the data at that respective time, otherwise the X_data pointers are unchanged.

◆ getVelocityData() [1/2]

void IBAMR::IBMethod::getVelocityData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  U_data,
IBTK::TimePoint  time_pt 
)
inherited

The time point should be one of CURRENT_TIME, HALF_TIME, or NEW_TIME. If this condition is met, U_data is set to the data at that respective time, otherwise the U_data pointers are unchanged.

◆ getVelocityData() [2/2]

void IBAMR::IBMethod::getVelocityData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  U_data,
double  data_time 
)
protectedinherited

Get the current structure velocity data.

data_time must be equal to one of current time, new time, or half time. If this condition is met, U_data is set to the data at that respective time, otherwise the U_data pointers are unchanged.

◆ getForceData() [1/2]

void IBAMR::IBMethod::getForceData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  F_data,
bool **  F_needs_ghost_fill,
IBTK::TimePoint  time_pt 
)
inherited

The time point should be one of CURRENT_TIME, HALF_TIME, or NEW_TIME. If this condition is met, F_data is set to the data at that respective time, otherwise the F_data pointers are unchanged.

◆ getForceData() [2/2]

void IBAMR::IBMethod::getForceData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  F_data,
bool **  F_needs_ghost_fill,
double  data_time 
)
protectedinherited

Get the current structure force data.

data_time must be equal to one of current time, new time, or half time. If this condition is met, F_data is set to the data at that respective time, otherwise the F_data pointers are unchanged.

◆ getLinearizedPositionData()

void IBAMR::IBMethod::getLinearizedPositionData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  X_data,
bool **  X_needs_ghost_fill 
)
protectedinherited

Get the linearized structure position data.

If the linearized position data does not exist, it will be created.

◆ getLECouplingPositionData()

void IBAMR::IBMethod::getLECouplingPositionData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  X_LE_data,
bool **  X_LE_needs_ghost_fill,
double  data_time 
)
protectedinherited

Get the current interpolation/spreading position data.

data_time must be equal to one of current time, new time, or half time. If this condition is met, X_LE_data is set to the data at that respective time, otherwise the X_LE_data is unchanged.

If this class is not set up to use fixed coupling, this returns data from getPositionData().

◆ getLinearizedVelocityData()

void IBAMR::IBMethod::getLinearizedVelocityData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  U_data)
protectedinherited

Get the linearized structure velocity data.

If the linearized velocity data does not exist, it will be created.

◆ getLinearizedForceData()

void IBAMR::IBMethod::getLinearizedForceData ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > **  F_data,
bool **  F_needs_ghost_fill 
)
protectedinherited

Get the linearized structure force data.

If the linearized force data does not exist, it will be created.

◆ reinitMidpointData()

void IBAMR::IBMethod::reinitMidpointData ( const std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > &  current_data,
const std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > &  new_data,
const std::vector< SAMRAI::tbox::Pointer< IBTK::LData > > &  half_data 
)
protectedinherited

Interpolate the current and new data to obtain values at the midpoint of the time interval.

◆ resetAnchorPointValues()

void IBAMR::IBMethod::resetAnchorPointValues ( std::vector< SAMRAI::tbox::Pointer< IBTK::LData > >  U_data,
int  coarsest_ln,
int  finest_ln 
)
protectedinherited

Set the elements of the Lagrangian vector to zero at anchored nodes of the curvilinear mesh.

◆ computeForce_SAMRAI()

static PetscErrorCode IBAMR::IBMethod::computeForce_SAMRAI ( void *  ctx,
Vec  X,
Vec  F 
)
staticprotectedinherited

◆ computeForce()

PetscErrorCode IBAMR::IBMethod::computeForce ( Vec  X,
Vec  F 
)
protectedinherited

◆ resetLagrangianForceFunction()

void IBAMR::IBMethod::resetLagrangianForceFunction ( double  init_data_time,
bool  initial_time 
)
privateinherited

Reset the Lagrangian force function object.

◆ resetLagrangianSourceFunction()

void IBAMR::IBMethod::resetLagrangianSourceFunction ( double  init_data_time,
bool  initial_time 
)
privateinherited

Reset the Lagrangian source function object.

◆ updateIBInstrumentationData()

void IBAMR::IBMethod::updateIBInstrumentationData ( int  timestep_num,
double  data_time 
)
privateinherited

Compute the flow rates and pressures in the internal flow meters and pressure gauges.

◆ registerIBHierarchyIntegrator()

virtual void IBAMR::IBStrategy::registerIBHierarchyIntegrator ( IBHierarchyIntegrator ib_solver)
virtualinherited

Register the IBHierarchyIntegrator object that is using this strategy class.

Reimplemented in IBAMR::IBLevelSetMethod, and IBAMR::IBStrategySet.

◆ registerEulerianCommunicationAlgorithms()

virtual void IBAMR::IBStrategy::registerEulerianCommunicationAlgorithms ( )
virtualinherited

Register Eulerian refinement or coarsening algorithms with the parent IBHierarchyIntegrator using the two versions of the protected methods IBStrategy::registerGhostfillRefineAlgorithm(), IBStrategy::registerProlongRefineAlgorithm(), and IBStrategy::registerCoarsenAlgorithm().

An empty default implementation is provided.

Reimplemented in IBAMR::CIBMethod, IBAMR::IBLevelSetMethod, IBAMR::IBInterpolantMethod, IBAMR::IBStrategySet, and IBAMR::GeneralizedIBMethod.

◆ getMaxPointDisplacement()

virtual double IBAMR::IBStrategy::getMaxPointDisplacement ( ) const
virtualinherited

Get the ratio of the maximum point displacement of all the structures owned by the current class to the cell width of the grid level on which the structure is assigned. This value is useful for determining if the Eulerian patch hierarchy needs to be regridded.

Note
The process of regridding is distinct, for some IBStrategy objects (like IBFEMethod), from forming (or reforming) the association between Lagrangian structures and patches. In particular, this function computes the distance between the current position of the structure and the structure at the point of the last regrid, which may not be the same point at which we last rebuilt the structure-to-patch mappings. The reassociation check should be implemented in postprocessIntegrateData().

Reimplemented in IBAMR::IBFEMethod, and IBAMR::IBStrategySet.

◆ setUseFixedLEOperators()

void IBAMR::IBStrategy::setUseFixedLEOperators ( bool  use_fixed_coupling_ops = true)
inherited

Indicate whether "fixed" interpolation and spreading operators should be used during Lagrangian-Eulerian interaction.

◆ setUseMultistepTimeStepping()

virtual void IBAMR::IBStrategy::setUseMultistepTimeStepping ( unsigned int  n_previous_steps = 1)
virtualinherited

Indicate that multistep time stepping will be used.

A default implementation is provided that emits an unrecoverable exception.

Parameters
[in]n_previous_stepsNumber of previous solution values that can be used by the multistep scheme.

Reimplemented in IBAMR::IBFEMethod.

◆ AB2Step()

virtual void IBAMR::IBStrategy::AB2Step ( double  current_time,
double  new_time 
)
virtualinherited

Advance the positions of the Lagrangian structure using the standard 2nd-order Adams-Bashforth rule.

A default implementation is provided that emits an unrecoverable exception.

Reimplemented in IBAMR::IBFEMethod.

◆ getINSHierarchyIntegrator()

INSHierarchyIntegrator* IBAMR::IBStrategy::getINSHierarchyIntegrator ( ) const
protectedinherited

Return a pointer to the INSHierarchyIntegrator object being used with the IBHierarchyIntegrator class registered with this IBStrategy object.

◆ getVelocityHierarchyDataOps()

SAMRAI::tbox::Pointer<SAMRAI::math::HierarchyDataOpsReal<NDIM, double> > IBAMR::IBStrategy::getVelocityHierarchyDataOps ( ) const
protectedinherited

Return a pointer to the HierarchyDataOpsReal object associated with velocity-like variables.

◆ getPressureHierarchyDataOps()

SAMRAI::tbox::Pointer<SAMRAI::math::HierarchyDataOpsReal<NDIM, double> > IBAMR::IBStrategy::getPressureHierarchyDataOps ( ) const
protectedinherited

Return a pointer to the HierarchyDataOpsReal object associated with pressure-like variables.

◆ getHierarchyMathOps()

SAMRAI::tbox::Pointer<IBTK::HierarchyMathOps> IBAMR::IBStrategy::getHierarchyMathOps ( ) const
protectedinherited

Return a pointer to a HierarchyMathOps object.

◆ registerVariable() [1/2]

void IBAMR::IBStrategy::registerVariable ( int current_idx,
int new_idx,
int scratch_idx,
SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > >  variable,
const SAMRAI::hier::IntVector< NDIM > &  scratch_ghosts = SAMRAI::hier::IntVector< NDIM >(0),
const std::string &  coarsen_name = "NO_COARSEN",
const std::string &  refine_name = "NO_REFINE",
SAMRAI::tbox::Pointer< IBTK::CartGridFunction init_fcn = nullptr,
const bool  register_for_restart = true 
)
protectedinherited

Register a state variable with the integrator. When a refine operator is specified, the data for the variable are automatically maintained as the patch hierarchy evolves.

All state variables are registered with three contexts: current, new, and scratch. The current context of a state variable is maintained from time step to time step and, if the necessary coarsen and refine operators are specified, as the patch hierarchy evolves.

◆ registerVariable() [2/2]

void IBAMR::IBStrategy::registerVariable ( int idx,
SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > >  variable,
const SAMRAI::hier::IntVector< NDIM > &  ghosts = SAMRAI::hier::IntVector< NDIM >(0),
SAMRAI::tbox::Pointer< SAMRAI::hier::VariableContext ctx = SAMRAI::tbox::PointerSAMRAI::hier::VariableContext >(nullptr),
const bool  register_for_restart = true 
)
protectedinherited

Register a variable with the integrator that may not be maintained from time step to time step.

By default, variables are registered with the scratch context, which is deallocated after each time step.

◆ registerGhostfillRefineAlgorithm()

void IBAMR::IBStrategy::registerGhostfillRefineAlgorithm ( const std::string &  name,
SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineAlgorithm< NDIM > >  ghostfill_alg,
std::unique_ptr< SAMRAI::xfer::RefinePatchStrategy< NDIM > >  ghostfill_patch_strategy = nullptr 
)
protectedinherited

Register a ghost cell-filling refine algorithm.

◆ registerProlongRefineAlgorithm()

void IBAMR::IBStrategy::registerProlongRefineAlgorithm ( const std::string &  name,
SAMRAI::tbox::Pointer< SAMRAI::xfer::RefineAlgorithm< NDIM > >  prolong_alg,
std::unique_ptr< SAMRAI::xfer::RefinePatchStrategy< NDIM > >  prolong_patch_strategy = nullptr 
)
protectedinherited

Register a data-prolonging refine algorithm.

◆ registerCoarsenAlgorithm()

void IBAMR::IBStrategy::registerCoarsenAlgorithm ( const std::string &  name,
SAMRAI::tbox::Pointer< SAMRAI::xfer::CoarsenAlgorithm< NDIM > >  coarsen_alg,
std::unique_ptr< SAMRAI::xfer::CoarsenPatchStrategy< NDIM > >  coarsen_patch_strategy = nullptr 
)
protectedinherited

Register a coarsen algorithm.

◆ getGhostfillRefineAlgorithm()

SAMRAI::tbox::Pointer<SAMRAI::xfer::RefineAlgorithm<NDIM> > IBAMR::IBStrategy::getGhostfillRefineAlgorithm ( const std::string &  name) const
protectedinherited

Get ghost cell-filling refine algorithm.

◆ getProlongRefineAlgorithm()

SAMRAI::tbox::Pointer<SAMRAI::xfer::RefineAlgorithm<NDIM> > IBAMR::IBStrategy::getProlongRefineAlgorithm ( const std::string &  name) const
protectedinherited

Get data-prolonging refine algorithm.

◆ getCoarsenAlgorithm()

SAMRAI::tbox::Pointer<SAMRAI::xfer::CoarsenAlgorithm<NDIM> > IBAMR::IBStrategy::getCoarsenAlgorithm ( const std::string &  name) const
protectedinherited

Get coarsen algorithm.

◆ getGhostfillRefineSchedules()

const std::vector<SAMRAI::tbox::Pointer<SAMRAI::xfer::RefineSchedule<NDIM> > >& IBAMR::IBStrategy::getGhostfillRefineSchedules ( const std::string &  name) const
protectedinherited

Get ghost cell-filling refine schedules.

◆ getProlongRefineSchedules()

const std::vector<SAMRAI::tbox::Pointer<SAMRAI::xfer::RefineSchedule<NDIM> > >& IBAMR::IBStrategy::getProlongRefineSchedules ( const std::string &  name) const
protectedinherited

Get data-prolonging refine schedules.

Note
These schedules are allocated only for level numbers >= 1.

◆ getCoarsenSchedules()

const std::vector<SAMRAI::tbox::Pointer<SAMRAI::xfer::CoarsenSchedule<NDIM> > >& IBAMR::IBStrategy::getCoarsenSchedules ( const std::string &  name) const
protectedinherited

Get coarsen schedules.

Note
These schedules are allocated only for level numbers >= 1.

◆ getLevelDt()

virtual double SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::getLevelDt ( const tbox::Pointer< hier::BasePatchLevel< DIM > >  level,
const double  dt_time,
const bool  initial_time 
)
virtualinherited

Determine time increment to advance data on level. The recompute_dt option specifies whether to compute the timestep using the current level data or to return the value stored by the time integrator. The default true setting means the timestep will be computed if no value is supplied.

This routine is only when Richardson extrapolation is being used. It is virtual with an empty implementation here (rather than pure virtual) so that users are not required to provide an implementation when the function is not needed.

◆ advanceLevel()

virtual double SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::advanceLevel ( const tbox::Pointer< hier::BasePatchLevel< DIM > >  level,
const tbox::Pointer< hier::BasePatchHierarchy< DIM > >  hierarchy,
const double  current_time,
const double  new_time,
const bool  first_step,
const bool  last_step,
const bool  regrid_advance = false 
)
virtualinherited

Advance data on all patches on specified patch level from current time (current_time) to new time (new_time). This routine is called only during time-dependent regridding procedures, such as Richardson extrapolation. It is virtual with an empty implementation here (rather than pure virtual) so that users are not required to provide an implementation when the function is not needed. The boolean arguments are used to determine the state of the algorithm and the data when the advance routine is called. Note that this advance function is also used during normal time integration steps.

When this function is called, the level data required to begin the advance must be allocated and be defined appropriately. Typically, this is equivalent to what is needed to initialize a new level after regridding. Upon exiting this routine, both current and new data may exist on the level. This data is needed until level synchronization occurs, in general. Current and new data may be reset by calling the member function resetTimeDependentData().

This routine is called from two different points within the Richardson exptrapolation process: to advance a temporary level that is coarser than the hierarchy level on which error estimation is performed, and to advance the hierarchy level itself. In the first case, the values of the boolean flags are:

  • first_step = true.
  • last_step = true.
  • regrid_advance = true.

In the second case, the values of the boolean flags are:

  • first_step (when regridding during time integration sequence) = true when the level is not coarsest level to synchronize immediately before the regridding process; else, false. (when generating initial hierarchy construction) = true, even though there may be multiple advance steps.
  • last_step = true when the advance is the last in the Richardson extrapolation step sequence; else false.
  • regrid_advance = true.

◆ resetTimeDependentData()

virtual void SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::resetTimeDependentData ( const tbox::Pointer< hier::BasePatchLevel< DIM > >  level,
const double  new_time,
const bool  can_be_refined 
)
virtualinherited

Reset time-dependent data storage for the specified patch level.

This routine only applies when Richardson extrapolation is being used. It is virtual with an empty implementation here (rather than pure virtual) so that users are not required to provide an implementation when the function is not needed.

◆ resetDataToPreadvanceState()

virtual void SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::resetDataToPreadvanceState ( const tbox::Pointer< hier::BasePatchLevel< DIM > >  level)
virtualinherited

Reset data on the patch level by destroying all patch data other than that which is needed to initialize the solution on that level. In other words, this is the data needed to begin a time integration step on the level.

This routine is only when Richardson extrapolation is being used. It is virtual with an empty implementation here (rather than pure virtual) so that users are not required to provide an implementation when the function is not needed.

◆ applyRichardsonExtrapolation()

virtual void SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::applyRichardsonExtrapolation ( const tbox::Pointer< hier::PatchLevel< DIM > >  level,
const double  error_data_time,
const int  tag_index,
const double  deltat,
const int  error_coarsen_ratio,
const bool  initial_time,
const bool  uses_gradient_detector_too 
)
virtualinherited

Set integer tags to "one" in cells where refinement of the given level should occur according to some user-supplied Richardson extrapolation criteria. The "error_data_time" argument is the regrid time. The "deltat" argument is the time increment to advance the solution on the level to be refined. Note that that level is finer than the level in the argument list, in general. The ratio between the argument level and the actual hierarchy level is given by the integer "coarsen ratio".

The integer "tag_index" argument is the patch descriptor index of the cell-centered integer tag array on each patch in the hierarchy.

The boolean argument initial_time indicates whether the level is being subject to refinement at the initial simulation time. If it is false, then the error estimation process is being invoked at some later time after the AMR hierarchy was initially constructed. Typically, this information is passed to the user's patch tagging routines since the application of the Richardson extrapolation process may be different in each case.

The boolean uses_gradient_detector_too is true when a gradient detector procedure is used in addition to Richardson extrapolation, and false otherwise. This argument helps the user to manage multiple regridding criteria.

This routine is only when Richardson extrapolation is being used. It is virtual with an empty implementation here (rather than pure virtual) so that users are not required to provide an implementation when the function is not needed.

◆ coarsenDataForRichardsonExtrapolation()

virtual void SAMRAI::mesh::StandardTagAndInitStrategy< DIM >::coarsenDataForRichardsonExtrapolation ( const tbox::Pointer< hier::PatchHierarchy< DIM > >  hierarchy,
const int  level_number,
const tbox::Pointer< hier::PatchLevel< DIM > >  coarser_level,
const double  coarsen_data_time,
const bool  before_advance 
)
virtualinherited

Coarsen solution data from level to coarse_level for Richardson extrapolation. Note that this routine will be called twice during the Richardson extrapolation error estimation process, once to set data on the coarser level and once to coarsen data from after advancing the fine level. The init_coarse_level boolean argument indicates whether data is set on the coarse level by coarsening the "old" time level solution or by coarsening the "new" solution on the fine level (i.e., after it has been advanced).

This routine is only when Richardson extrapolation is being used. It is virtual with an empty implementation here (rather than pure virtual) so that users are not required to provide an implementation when the function is not needed.

Member Data Documentation

◆ d_no_structures

const int IBAMR::ConstraintIBMethod::d_no_structures
private

No of immersed structures.

◆ d_ib_kinematics

std::vector<SAMRAI::tbox::Pointer<IBAMR::ConstraintIBKinematics> > IBAMR::ConstraintIBMethod::d_ib_kinematics
private

Pointer to the kinematics of the immersed structures.

◆ d_FuRMoRP_current_time

double IBAMR::ConstraintIBMethod::d_FuRMoRP_current_time = 0.0
private

FuRMoRP apply time.

◆ d_FuRMoRP_new_time

double IBAMR::ConstraintIBMethod::d_FuRMoRP_new_time = 0.0
private

◆ d_vol_element

std::vector<double> IBAMR::ConstraintIBMethod::d_vol_element
private

Volume element associated with material points.

◆ d_vol_element_is_set

std::vector<bool> IBAMR::ConstraintIBMethod::d_vol_element_is_set
private

◆ d_structure_vol

std::vector<double> IBAMR::ConstraintIBMethod::d_structure_vol
private

Volume associated with each immersed structure

◆ d_structure_mom

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_structure_mom
private

Linear momentum associated with each immersed structure

◆ d_structure_rotational_mom

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_structure_rotational_mom
private

Rotational momentum associated with each immersed structure with respect to their COM

◆ d_needs_div_free_projection

bool IBAMR::ConstraintIBMethod::d_needs_div_free_projection = false
private

If divergence free projection is needed after FuRMoRP algorithm?

◆ d_rigid_trans_vel_current

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_rigid_trans_vel_current
private

Rigid translational velocity of the structures.

◆ d_rigid_trans_vel_new

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_rigid_trans_vel_new
private

◆ d_rigid_rot_vel_current

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_rigid_rot_vel_current
private

Rigid rotational velocity of the structures.

◆ d_rigid_rot_vel_new

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_rigid_rot_vel_new
private

◆ d_incremented_angle_from_reference_axis

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_incremented_angle_from_reference_axis
private

Incremented angle from x, y and z axis when the body is rotating.

◆ d_vel_com_def_current

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_vel_com_def_current
private

Translational velocity of the structures due to deformational kinematics.

◆ d_vel_com_def_new

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_vel_com_def_new
private

◆ d_omega_com_def_current

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_omega_com_def_current
private

Rotational velocity of the structures due to deformational kinematics.

◆ d_omega_com_def_new

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_omega_com_def_new
private

◆ d_center_of_mass_current

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_center_of_mass_current
private

Center of mass of the immersed structures.

◆ d_center_of_mass_new

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_center_of_mass_new
private

◆ d_center_of_mass_unshifted_current

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_center_of_mass_unshifted_current
private

unshifted center of mass

◆ d_center_of_mass_unshifted_new

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_center_of_mass_unshifted_new
private

◆ d_moment_of_inertia_current

std::vector<Eigen::Matrix3d> IBAMR::ConstraintIBMethod::d_moment_of_inertia_current
private

Moment of inertia of the structures.

◆ d_moment_of_inertia_new

std::vector<Eigen::Matrix3d> IBAMR::ConstraintIBMethod::d_moment_of_inertia_new
private

◆ d_tagged_pt_lag_idx

std::vector<int> IBAMR::ConstraintIBMethod::d_tagged_pt_lag_idx
private

Tag a Lagrangian point (generally eye of the fish) of the immersed structures.

◆ d_tagged_pt_position

std::vector<std::vector<double> > IBAMR::ConstraintIBMethod::d_tagged_pt_position
private

Coordinates of the tagged points of different structures.

◆ d_rho_solid

std::vector<double> IBAMR::ConstraintIBMethod::d_rho_solid
private

Density of the structures.

◆ d_rho_fluid

double IBAMR::ConstraintIBMethod::d_rho_fluid = std::numeric_limits<double>::quiet_NaN()
private

Density of the fluid in constant coefficient case.

◆ d_rho_is_const

bool IBAMR::ConstraintIBMethod::d_rho_is_const = true
private

Whether or not the density from the integrator is constant

◆ d_calculate_structure_linear_mom

bool IBAMR::ConstraintIBMethod::d_calculate_structure_linear_mom = false
private

Bools for computing linear and rotational momentums of the body

◆ d_calculate_structure_rotational_mom

bool IBAMR::ConstraintIBMethod::d_calculate_structure_rotational_mom = false
private

◆ d_timestep_counter

int IBAMR::ConstraintIBMethod::d_timestep_counter = 0
private

Iteration_counter for printing stuff.

◆ d_output_interval

int IBAMR::ConstraintIBMethod::d_output_interval = 1
private

◆ d_print_output

bool IBAMR::ConstraintIBMethod::d_print_output = false
private

Bools for outputing stuff which is calculated on the fly.

◆ d_output_drag

bool IBAMR::ConstraintIBMethod::d_output_drag = false
private

◆ d_output_torque

bool IBAMR::ConstraintIBMethod::d_output_torque = false
private

◆ d_output_power

bool IBAMR::ConstraintIBMethod::d_output_power = false
private

◆ d_output_trans_vel

bool IBAMR::ConstraintIBMethod::d_output_trans_vel = false
private

◆ d_output_rot_vel

bool IBAMR::ConstraintIBMethod::d_output_rot_vel = false
private

◆ d_output_COM_coordinates

bool IBAMR::ConstraintIBMethod::d_output_COM_coordinates = false
private

◆ d_output_MOI

bool IBAMR::ConstraintIBMethod::d_output_MOI = false
private

◆ d_output_eul_mom

bool IBAMR::ConstraintIBMethod::d_output_eul_mom = false
private

◆ d_dir_name

std::string IBAMR::ConstraintIBMethod::d_dir_name = "./ConstraintIBMethodDump"
private

output file name string.

◆ d_base_output_filename

std::string IBAMR::ConstraintIBMethod::d_base_output_filename = "ImmersedStructure"
private

◆ d_l_data_U_interp

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::ConstraintIBMethod::d_l_data_U_interp
private

Store LData for only those levels which contain immersed structures.

◆ d_l_data_U_correction

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::ConstraintIBMethod::d_l_data_U_correction
private

◆ d_l_data_U_new

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::ConstraintIBMethod::d_l_data_U_new
private

◆ d_l_data_U_current

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::ConstraintIBMethod::d_l_data_U_current
private

◆ d_l_data_U_half

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::ConstraintIBMethod::d_l_data_U_half
private

◆ d_l_data_X_half_Euler

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::ConstraintIBMethod::d_l_data_X_half_Euler
private

◆ d_l_data_X_new_MidPoint

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::ConstraintIBMethod::d_l_data_X_new_MidPoint
private

◆ d_hier_sc_data_ops

SAMRAI::tbox::Pointer<SAMRAI::math::HierarchySideDataOpsReal<NDIM, double> > IBAMR::ConstraintIBMethod::d_hier_sc_data_ops
private

Hierarchy operations object. Needed for projection step.

◆ d_hier_cc_data_ops

SAMRAI::tbox::Pointer<SAMRAI::math::HierarchyCellDataOpsReal<NDIM, double> > IBAMR::ConstraintIBMethod::d_hier_cc_data_ops
private

◆ d_no_fill_op

SAMRAI::tbox::Pointer<IBTK::HierarchyGhostCellInterpolation> IBAMR::ConstraintIBMethod::d_no_fill_op
private

◆ d_wgt_cc_idx

int IBAMR::ConstraintIBMethod::d_wgt_cc_idx
private

◆ d_wgt_sc_idx

int IBAMR::ConstraintIBMethod::d_wgt_sc_idx
private

◆ d_volume

double IBAMR::ConstraintIBMethod::d_volume
private

◆ d_u_var

SAMRAI::tbox::Pointer<SAMRAI::hier::Variable<NDIM> > IBAMR::ConstraintIBMethod::d_u_var
private

Variables and variable contexts associated with calculating divergence free projection.

◆ d_u_fluidSolve_var

SAMRAI::tbox::Pointer<SAMRAI::hier::Variable<NDIM> > IBAMR::ConstraintIBMethod::d_u_fluidSolve_var
private

◆ d_phi_var

SAMRAI::tbox::Pointer<SAMRAI::pdat::CellVariable<NDIM, double> > IBAMR::ConstraintIBMethod::d_phi_var
private

◆ d_Div_u_var

SAMRAI::tbox::Pointer<SAMRAI::pdat::CellVariable<NDIM, double> > IBAMR::ConstraintIBMethod::d_Div_u_var
private

◆ d_scratch_context

SAMRAI::tbox::Pointer<SAMRAI::hier::VariableContext> IBAMR::ConstraintIBMethod::d_scratch_context
private

◆ d_u_fluidSolve_idx

int IBAMR::ConstraintIBMethod::d_u_fluidSolve_idx
private

◆ d_u_fluidSolve_cib_idx

int IBAMR::ConstraintIBMethod::d_u_fluidSolve_cib_idx
private

◆ d_phi_idx

int IBAMR::ConstraintIBMethod::d_phi_idx
private

◆ d_Div_u_scratch_idx

int IBAMR::ConstraintIBMethod::d_Div_u_scratch_idx
private

◆ d_rho_var

SAMRAI::tbox::Pointer<SAMRAI::pdat::SideVariable<NDIM, double> > IBAMR::ConstraintIBMethod::d_rho_var
private

Variables associated with the spatially varying density field, which is maintained by an integrator.

◆ d_rho_ins_idx

int IBAMR::ConstraintIBMethod::d_rho_ins_idx = IBTK::invalid_index
private

◆ d_rho_scratch_idx

int IBAMR::ConstraintIBMethod::d_rho_scratch_idx = IBTK::invalid_index
private

◆ d_velcorrection_projection_bc_coef

SAMRAI::solv::LocationIndexRobinBcCoefs<NDIM> IBAMR::ConstraintIBMethod::d_velcorrection_projection_bc_coef
private

The following variables are needed to solve cell centered poison equation for \( \phi \) ,which is used to project the corrected background fluid velocity on divergence free field to remove kinkiness introduced via FuRMoRP algorithm.

◆ d_velcorrection_projection_spec

std::unique_ptr<SAMRAI::solv::PoissonSpecifications> IBAMR::ConstraintIBMethod::d_velcorrection_projection_spec
private

◆ d_velcorrection_projection_op

SAMRAI::tbox::Pointer<IBTK::CCLaplaceOperator> IBAMR::ConstraintIBMethod::d_velcorrection_projection_op
private

◆ d_velcorrection_projection_solver

SAMRAI::tbox::Pointer<IBTK::PETScKrylovPoissonSolver> IBAMR::ConstraintIBMethod::d_velcorrection_projection_solver
private

◆ d_velcorrection_projection_fac_op

SAMRAI::tbox::Pointer<IBTK::CCPoissonPointRelaxationFACOperator> IBAMR::ConstraintIBMethod::d_velcorrection_projection_fac_op
private

◆ d_velcorrection_projection_fac_pc_db

SAMRAI::tbox::Pointer<SAMRAI::tbox::Database> IBAMR::ConstraintIBMethod::d_velcorrection_projection_fac_pc_db
private

◆ d_velcorrection_projection_fac_pc

SAMRAI::tbox::Pointer<IBTK::FACPreconditioner> IBAMR::ConstraintIBMethod::d_velcorrection_projection_fac_pc
private

◆ d_trans_vel_stream

std::vector<std::unique_ptr<std::ofstream> > IBAMR::ConstraintIBMethod::d_trans_vel_stream
private

File streams associated for the output.

◆ d_rot_vel_stream

std::vector<std::unique_ptr<std::ofstream> > IBAMR::ConstraintIBMethod::d_rot_vel_stream
private

◆ d_drag_force_stream

std::vector<std::unique_ptr<std::ofstream> > IBAMR::ConstraintIBMethod::d_drag_force_stream
private

◆ d_moment_of_inertia_stream

std::vector<std::unique_ptr<std::ofstream> > IBAMR::ConstraintIBMethod::d_moment_of_inertia_stream
private

◆ d_torque_stream

std::vector<std::unique_ptr<std::ofstream> > IBAMR::ConstraintIBMethod::d_torque_stream
private

◆ d_position_COM_stream

std::vector<std::unique_ptr<std::ofstream> > IBAMR::ConstraintIBMethod::d_position_COM_stream
private

◆ d_power_spent_stream

std::vector<std::unique_ptr<std::ofstream> > IBAMR::ConstraintIBMethod::d_power_spent_stream
private

◆ d_eulerian_mom_stream

std::fstream IBAMR::ConstraintIBMethod::d_eulerian_mom_stream
private

Stream for calculating Eulerian momentum.

◆ d_prefluidsolve_callback_fns

std::vector<void (*)(const double, const double, const int, void*)> IBAMR::ConstraintIBMethod::d_prefluidsolve_callback_fns
private

Pre and post fluid solve call back functions and contexts.

◆ d_postfluidsolve_callback_fns

std::vector<void (*)(const double, const double, const int, void*)> IBAMR::ConstraintIBMethod::d_postfluidsolve_callback_fns
private

◆ d_prefluidsolve_callback_fns_ctx

std::vector<void*> IBAMR::ConstraintIBMethod::d_prefluidsolve_callback_fns_ctx
private

◆ d_postfluidsolve_callback_fns_ctx

std::vector<void*> IBAMR::ConstraintIBMethod::d_postfluidsolve_callback_fns_ctx
private

◆ d_u_phys_bdry_op

IBTK::RobinPhysBdryPatchStrategy* IBAMR::ConstraintIBMethod::d_u_phys_bdry_op = nullptr
private

◆ d_do_log

bool IBAMR::IBMethod::d_do_log = false
protectedinherited

◆ d_hierarchy

SAMRAI::tbox::Pointer<SAMRAI::hier::PatchHierarchy<NDIM> > IBAMR::IBMethod::d_hierarchy
protectedinherited

◆ d_gridding_alg

SAMRAI::tbox::Pointer<SAMRAI::mesh::GriddingAlgorithm<NDIM> > IBAMR::IBMethod::d_gridding_alg
protectedinherited

◆ d_current_time

double IBAMR::IBMethod::d_current_time = std::numeric_limits<double>::quiet_NaN()
protectedinherited

◆ d_new_time

double IBAMR::IBMethod::d_new_time = std::numeric_limits<double>::quiet_NaN()
protectedinherited

◆ d_half_time

double IBAMR::IBMethod::d_half_time = std::numeric_limits<double>::quiet_NaN()
protectedinherited

◆ d_X_current_needs_ghost_fill

bool IBAMR::IBMethod::d_X_current_needs_ghost_fill = true
protectedinherited

◆ d_X_new_needs_ghost_fill

bool IBAMR::IBMethod::d_X_new_needs_ghost_fill = true
protectedinherited

◆ d_X_half_needs_ghost_fill

bool IBAMR::IBMethod::d_X_half_needs_ghost_fill = true
protectedinherited

◆ d_X_jac_needs_ghost_fill

bool IBAMR::IBMethod::d_X_jac_needs_ghost_fill = true
protectedinherited

◆ d_X_LE_new_needs_ghost_fill

bool IBAMR::IBMethod::d_X_LE_new_needs_ghost_fill = true
protectedinherited

◆ d_X_LE_half_needs_ghost_fill

bool IBAMR::IBMethod::d_X_LE_half_needs_ghost_fill = true
protectedinherited

◆ d_F_current_needs_ghost_fill

bool IBAMR::IBMethod::d_F_current_needs_ghost_fill = true
protectedinherited

◆ d_F_new_needs_ghost_fill

bool IBAMR::IBMethod::d_F_new_needs_ghost_fill = true
protectedinherited

◆ d_F_half_needs_ghost_fill

bool IBAMR::IBMethod::d_F_half_needs_ghost_fill = true
protectedinherited

◆ d_F_jac_needs_ghost_fill

bool IBAMR::IBMethod::d_F_jac_needs_ghost_fill = true
protectedinherited

◆ d_l_data_manager

IBTK::LDataManager* IBAMR::IBMethod::d_l_data_manager
protectedinherited

◆ d_interp_kernel_fcn

std::string IBAMR::IBMethod::d_interp_kernel_fcn = "IB_4"
protectedinherited

◆ d_spread_kernel_fcn

std::string IBAMR::IBMethod::d_spread_kernel_fcn = "IB_4"
protectedinherited

◆ d_error_if_points_leave_domain

bool IBAMR::IBMethod::d_error_if_points_leave_domain = false
protectedinherited

◆ d_ghosts

SAMRAI::hier::IntVector<NDIM> IBAMR::IBMethod::d_ghosts
protectedinherited

◆ d_X_current_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_X_current_data
protectedinherited

◆ d_X_new_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_X_new_data
protectedinherited

◆ d_X_half_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_X_half_data
protectedinherited

◆ d_X_jac_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_X_jac_data
protectedinherited

◆ d_X_LE_new_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_X_LE_new_data
protectedinherited

◆ d_X_LE_half_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_X_LE_half_data
protectedinherited

◆ d_U_current_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_U_current_data
protectedinherited

◆ d_U_new_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_U_new_data
protectedinherited

◆ d_U_half_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_U_half_data
protectedinherited

◆ d_U_jac_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_U_jac_data
protectedinherited

◆ d_F_current_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_F_current_data
protectedinherited

◆ d_F_new_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_F_new_data
protectedinherited

◆ d_F_half_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_F_half_data
protectedinherited

◆ d_F_jac_data

std::vector<SAMRAI::tbox::Pointer<IBTK::LData> > IBAMR::IBMethod::d_F_jac_data
protectedinherited

◆ d_anchor_point_local_idxs

std::vector<std::set<int> > IBAMR::IBMethod::d_anchor_point_local_idxs
protectedinherited

◆ d_instrument_panel

SAMRAI::tbox::Pointer<IBInstrumentPanel> IBAMR::IBMethod::d_instrument_panel
protectedinherited

◆ d_total_flow_volume

std::vector<double> IBAMR::IBMethod::d_total_flow_volume
protectedinherited

◆ d_l_initializer

SAMRAI::tbox::Pointer<IBTK::LInitStrategy> IBAMR::IBMethod::d_l_initializer
protectedinherited

◆ d_ib_force_fcn

SAMRAI::tbox::Pointer<IBLagrangianForceStrategy> IBAMR::IBMethod::d_ib_force_fcn
protectedinherited

◆ d_ib_force_fcn_needs_init

bool IBAMR::IBMethod::d_ib_force_fcn_needs_init = true
protectedinherited

◆ d_ib_source_fcn

SAMRAI::tbox::Pointer<IBLagrangianSourceStrategy> IBAMR::IBMethod::d_ib_source_fcn
protectedinherited

◆ d_ib_source_fcn_needs_init

bool IBAMR::IBMethod::d_ib_source_fcn_needs_init = true
protectedinherited

◆ d_X_src

std::vector<std::vector<IBTK::Point> > IBAMR::IBMethod::d_X_src
protectedinherited

◆ d_r_src

std::vector<std::vector<double> > IBAMR::IBMethod::d_r_src
protectedinherited

◆ d_P_src

std::vector<std::vector<double> > IBAMR::IBMethod::d_P_src
protectedinherited

◆ d_Q_src

std::vector<std::vector<double> > IBAMR::IBMethod::d_Q_src
protectedinherited

◆ d_n_src

std::vector<int> IBAMR::IBMethod::d_n_src
protectedinherited

◆ d_normalize_source_strength

bool IBAMR::IBMethod::d_normalize_source_strength = false
protectedinherited

◆ d_post_processor

SAMRAI::tbox::Pointer<IBMethodPostProcessStrategy> IBAMR::IBMethod::d_post_processor
protectedinherited

◆ d_silo_writer

SAMRAI::tbox::Pointer<IBTK::LSiloDataWriter> IBAMR::IBMethod::d_silo_writer
protectedinherited

◆ d_load_balancer

SAMRAI::tbox::Pointer<SAMRAI::mesh::LoadBalancer<NDIM> > IBAMR::IBMethod::d_load_balancer
protectedinherited

◆ d_workload_idx

int IBAMR::IBMethod::d_workload_idx = IBTK::invalid_index
protectedinherited

◆ d_object_name

std::string IBAMR::IBMethod::d_object_name
protectedinherited

◆ d_registered_for_restart

bool IBAMR::IBMethod::d_registered_for_restart
protectedinherited

◆ d_force_jac_mffd

bool IBAMR::IBMethod::d_force_jac_mffd = false
privateinherited

Jacobian data.

◆ d_force_jac

Mat IBAMR::IBMethod::d_force_jac = nullptr
privateinherited

◆ d_force_jac_data_time

double IBAMR::IBMethod::d_force_jac_data_time
privateinherited

◆ d_ib_solver

IBHierarchyIntegrator* IBAMR::IBStrategy::d_ib_solver = nullptr
protectedinherited

The IBHierarchyIntegrator object that is using this strategy class.

◆ d_use_fixed_coupling_ops

bool IBAMR::IBStrategy::d_use_fixed_coupling_ops = false
protectedinherited

Whether to use "fixed" Lagrangian-Eulerian coupling operators.


The documentation for this class was generated from the following file: