|
IBAMR
IBAMR version 0.19.
|
Class ConstraintIBMethod implements the rigidity constraint for rigid and deforming bodies using the constraint based IB method. More...
#include <ibamr/ConstraintIBMethod.h>

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::LDataManager * | getLDataManager () const |
| SAMRAI::tbox::Pointer< IBInstrumentPanel > | getIBInstrumentPanel () 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) |
Static Protected Member Functions | |
| static PetscErrorCode | computeForce_SAMRAI (void *ctx, Vec X, Vec F) |
Private Member Functions | |
| ConstraintIBMethod ()=delete | |
| Default constructor. More... | |
| ConstraintIBMethod (const ConstraintIBMethod &from)=delete | |
| Default copy constructor. More... | |
| ConstraintIBMethod & | operator= (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) |
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).
| IBAMR::ConstraintIBMethod::ConstraintIBMethod | ( | std::string | object_name, |
| SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > | input_db, | ||
| const int | no_structures, | ||
| bool | register_for_restart = true |
||
| ) |
| IBAMR::ConstraintIBMethod::~ConstraintIBMethod | ( | ) |
|
privatedelete |
|
privatedelete |
| void IBAMR::ConstraintIBMethod::initializeHierarchyOperatorsandData | ( | ) |
|
overridevirtual |
Reimplemented from IBAMR::IBStrategy.
|
overridevirtual |
Reimplemented from IBAMR::IBStrategy.
|
overridevirtual |
Reimplemented from IBAMR::IBStrategy.
| void IBAMR::ConstraintIBMethod::registerConstraintIBKinematics | ( | const std::vector< SAMRAI::tbox::Pointer< IBAMR::ConstraintIBKinematics > > & | ib_kinematics_op | ) |
|
inline |
|
overridevirtual |
Reimplemented from IBAMR::IBStrategy.
|
inline |
|
overridevirtual |
Reimplemented from IBAMR::IBStrategy.
|
overridevirtual |
Implements IBAMR::IBStrategy.
|
overridevirtual |
Implements IBAMR::IBStrategy.
|
overridevirtual |
Implements SAMRAI::tbox::Serializable.
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
|
privatedelete |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
inlineprivate |
|
inlineprivate |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
inherited |
Supply a Lagrangian force object.
|
inherited |
Supply a Lagrangian source object.
|
inherited |
Supply a Lagrangian initialization object.
|
inherited |
Free references to Lagrangian initialization objects.
|
inherited |
Supply a post processor object.
|
inherited |
Return a pointer to the Lagrangian data manager object.
|
inherited |
Return a pointer to the instrumentation manager object.
|
inherited |
Register a Lagrangian Silo data writer so this class will write plot files that may be postprocessed with the VisIt visualization tool.
|
overridevirtualinherited |
Return the number of ghost cells required by the Lagrangian-Eulerian interaction routines.
Implements IBAMR::IBStrategy.
|
overridevirtualinherited |
Setup the tag buffer.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Inactivate a structure/part. See IBAMR::IBStrategy::inactivateLagrangianStructure().
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Activate a previously inactivated structure/part to be used again in FSI calculations. See IBAMR::IBStrategy::activateLagrangianStructure().
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Determine whether or not the given structure or part is currently activated. See IBAMR::IBStrategy::getLagrangianStructureIsActivated().
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Create solution and rhs data on the specified level of the patch hierarchy.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Setup solution and rhs data on the specified level of the patch hierarchy.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Set the value of the updated position vector.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Set the value of the intermediate position vector used in evaluating the linearized problem.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Compute the residual on the specified level of the patch hierarchy.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Compute the linearized residual for the given intermediate position vector.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Update the positions used for the "fixed" interpolation and spreading operators.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Interpolate the Eulerian velocity to the curvilinear mesh at the specified time within the current time interval.
Implements IBAMR::IBStrategy.
|
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.
|
overridevirtualinherited |
Advance the positions of the Lagrangian structure using the backward Euler method.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Advance the positions of the Lagrangian structure using the trapezoidal rule.
Implements IBAMR::IBStrategy.
Reimplemented in IBAMR::PenaltyIBMethod.
|
overridevirtualinherited |
Compute the Lagrangian force at the specified time within the current time interval.
Implements IBAMR::IBStrategy.
Reimplemented in IBAMR::PenaltyIBMethod.
|
overridevirtualinherited |
Compute the Lagrangian force of the linearized problem for the specified configuration of the updated position vector.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Construct the linearized Lagrangian force Jacobian.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Spread the Lagrangian force to the Cartesian grid at the specified time within the current time interval.
Implements IBAMR::IBStrategy.
|
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.
|
overridevirtualinherited |
Construct the IB interpolation operator.
Implements IBAMR::IBImplicitStrategy.
|
overridevirtualinherited |
Indicate whether there are any internal fluid sources/sinks.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Compute the Lagrangian source/sink density at the specified time within the current time interval.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Spread the Lagrangian source/sink density to the Cartesian grid at the specified time within the current time interval.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Compute the pressures at the positions of any distributed internal fluid sources or sinks.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Execute user-defined post-processing operations.
Reimplemented from IBAMR::IBStrategy.
|
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.
|
overridevirtualinherited |
Register a load balancer and work load patch data index with the IB strategy object.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Add the estimated computational work from the current object per cell into the specified workload_data_idx.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Begin redistributing Lagrangian data prior to regridding the patch hierarchy.
Reimplemented from IBAMR::IBStrategy.
|
overridevirtualinherited |
Complete redistributing Lagrangian data following regridding the patch hierarchy.
Reimplemented from IBAMR::IBStrategy.
|
overrideinherited |
Initialize data on a new level after it is inserted into an AMR patch hierarchy by the gridding algorithm.
|
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.
|
overrideinherited |
Reset cached hierarchy dependent data.
|
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.
|
overrideinherited |
Set integer tags to "one" in cells where refinement of the given level should occur according to user-supplied feature detection criteria.
|
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.
|
inherited |
If TimePoint is not one of CURRENT_TIME, HALF_TIME, or NEW_TIME, this returns NaN.
|
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.
|
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.
|
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.
|
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.
|
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.
|
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.
|
protectedinherited |
Get the linearized structure position data.
If the linearized position data does not exist, it will be created.
|
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().
|
protectedinherited |
Get the linearized structure velocity data.
If the linearized velocity data does not exist, it will be created.
|
protectedinherited |
Get the linearized structure force data.
If the linearized force data does not exist, it will be created.
|
protectedinherited |
Interpolate the current and new data to obtain values at the midpoint of the time interval.
|
protectedinherited |
Set the elements of the Lagrangian vector to zero at anchored nodes of the curvilinear mesh.
|
staticprotectedinherited |
|
protectedinherited |
|
privateinherited |
Reset the Lagrangian force function object.
|
privateinherited |
Reset the Lagrangian source function object.
|
privateinherited |
Compute the flow rates and pressures in the internal flow meters and pressure gauges.
|
virtualinherited |
Register the IBHierarchyIntegrator object that is using this strategy class.
Reimplemented in IBAMR::IBLevelSetMethod, and IBAMR::IBStrategySet.
|
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.
|
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.
Reimplemented in IBAMR::IBFEMethod, and IBAMR::IBStrategySet.
|
inherited |
Indicate whether "fixed" interpolation and spreading operators should be used during Lagrangian-Eulerian interaction.
|
virtualinherited |
Indicate that multistep time stepping will be used.
A default implementation is provided that emits an unrecoverable exception.
| [in] | n_previous_steps | Number of previous solution values that can be used by the multistep scheme. |
Reimplemented in IBAMR::IBFEMethod.
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.
|
protectedinherited |
Return a pointer to the INSHierarchyIntegrator object being used with the IBHierarchyIntegrator class registered with this IBStrategy object.
|
protectedinherited |
Return a pointer to the HierarchyDataOpsReal object associated with velocity-like variables.
|
protectedinherited |
Return a pointer to the HierarchyDataOpsReal object associated with pressure-like variables.
|
protectedinherited |
Return a pointer to a HierarchyMathOps object.
|
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.
|
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.
|
protectedinherited |
Register a ghost cell-filling refine algorithm.
|
protectedinherited |
Register a data-prolonging refine algorithm.
|
protectedinherited |
Register a coarsen algorithm.
|
protectedinherited |
Get ghost cell-filling refine algorithm.
|
protectedinherited |
Get data-prolonging refine algorithm.
|
protectedinherited |
Get coarsen algorithm.
|
protectedinherited |
Get ghost cell-filling refine schedules.
|
protectedinherited |
Get data-prolonging refine schedules.
|
protectedinherited |
Get coarsen schedules.
|
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.
|
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:
In the second case, the values of the boolean flags are:
|
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.
|
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.
|
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.
|
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.
|
private |
No of immersed structures.
|
private |
Pointer to the kinematics of the immersed structures.
|
private |
FuRMoRP apply time.
|
private |
|
private |
Volume element associated with material points.
|
private |
|
private |
Volume associated with each immersed structure
|
private |
Linear momentum associated with each immersed structure
|
private |
Rotational momentum associated with each immersed structure with respect to their COM
|
private |
If divergence free projection is needed after FuRMoRP algorithm?
|
private |
Rigid translational velocity of the structures.
|
private |
|
private |
Rigid rotational velocity of the structures.
|
private |
|
private |
Incremented angle from x, y and z axis when the body is rotating.
|
private |
Translational velocity of the structures due to deformational kinematics.
|
private |
|
private |
Rotational velocity of the structures due to deformational kinematics.
|
private |
|
private |
Center of mass of the immersed structures.
|
private |
|
private |
unshifted center of mass
|
private |
|
private |
Moment of inertia of the structures.
|
private |
|
private |
Tag a Lagrangian point (generally eye of the fish) of the immersed structures.
|
private |
Coordinates of the tagged points of different structures.
|
private |
Density of the structures.
Density of the fluid in constant coefficient case.
|
private |
Whether or not the density from the integrator is constant
|
private |
Bools for computing linear and rotational momentums of the body
|
private |
|
private |
Iteration_counter for printing stuff.
|
private |
|
private |
Bools for outputing stuff which is calculated on the fly.
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
output file name string.
|
private |
|
private |
Store LData for only those levels which contain immersed structures.
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
Hierarchy operations object. Needed for projection step.
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
Variables and variable contexts associated with calculating divergence free projection.
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
Variables associated with the spatially varying density field, which is maintained by an integrator.
|
private |
|
private |
|
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.
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
File streams associated for the output.
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
|
private |
Stream for calculating Eulerian momentum.
|
private |
Pre and post fluid solve call back functions and contexts.
|
private |
|
private |
|
private |
|
private |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
privateinherited |
Jacobian data.
|
privateinherited |
|
privateinherited |
|
protectedinherited |
The IBHierarchyIntegrator object that is using this strategy class.
|
protectedinherited |
Whether to use "fixed" Lagrangian-Eulerian coupling operators.
1.8.17