QuarterCircleDrivenCavityProblem< ELEMENT > Class Template Reference
+ Inheritance diagram for QuarterCircleDrivenCavityProblem< ELEMENT >:

Public Member Functions

 QuarterCircleDrivenCavityProblem (NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt body_force_fct_pt)
 Constructor. More...
 
 ~QuarterCircleDrivenCavityProblem ()
 Destructor: Empty – all memory gets cleaned up in base destructor. More...
 
void actions_after_newton_solve ()
 Update the after solve (empty) More...
 
void actions_before_newton_solve ()
 
void actions_after_adapt ()
 After adaptation: Unpin pressure and pin redudant pressure dofs. More...
 
void doc_solution (DocInfo &doc_info)
 Doc the solution. More...
 
 QuarterCircleDrivenCavityProblem (NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt body_force_fct_pt)
 Constructor. More...
 
 ~QuarterCircleDrivenCavityProblem ()
 Destructor: Empty. More...
 
void actions_after_newton_solve ()
 Update the after solve (empty) More...
 
void actions_before_newton_solve ()
 
void actions_after_adapt ()
 After adaptation: Unpin pressure and pin redudant pressure dofs. More...
 
void doc_solution (DocInfo &doc_info)
 Doc the solution. More...
 
 QuarterCircleDrivenCavityProblem (NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt body_force_fct_pt, const unsigned &driven_bdry)
 Constructor. More...
 
 ~QuarterCircleDrivenCavityProblem ()
 Destructor: Empty. More...
 
void actions_after_newton_solve ()
 Update the after solve (empty) More...
 
void actions_before_newton_solve ()
 
void actions_after_adapt ()
 After adaptation: Unpin pressure and pin redudant pressure dofs. More...
 
void doc_solution (DocInfo &doc_info)
 Doc the solution. More...
 
- Public Member Functions inherited from oomph::Problem
virtual void debug_hook_fct (const unsigned &i)
 
void set_analytic_dparameter (double *const &parameter_pt)
 
void unset_analytic_dparameter (double *const &parameter_pt)
 
bool is_dparameter_calculated_analytically (double *const &parameter_pt)
 
void set_analytic_hessian_products ()
 
void unset_analytic_hessian_products ()
 
bool are_hessian_products_calculated_analytically ()
 
void set_pinned_values_to_zero ()
 
bool distributed () const
 
virtual void actions_before_adapt ()
 
OomphCommunicatorcommunicator_pt ()
 access function to the oomph-lib communicator More...
 
const OomphCommunicatorcommunicator_pt () const
 access function to the oomph-lib communicator, const version More...
 
 Problem ()
 
 Problem (const Problem &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const Problem &)=delete
 Broken assignment operator. More...
 
virtual ~Problem ()
 Virtual destructor to clean up memory. More...
 
Mesh *& mesh_pt ()
 Return a pointer to the global mesh. More...
 
Mesh *const & mesh_pt () const
 Return a pointer to the global mesh (const version) More...
 
Mesh *& mesh_pt (const unsigned &imesh)
 
Mesh *const & mesh_pt (const unsigned &imesh) const
 Return a pointer to the i-th submesh (const version) More...
 
unsigned nsub_mesh () const
 Return number of submeshes. More...
 
unsigned add_sub_mesh (Mesh *const &mesh_pt)
 
void flush_sub_meshes ()
 
void build_global_mesh ()
 
void rebuild_global_mesh ()
 
LinearSolver *& linear_solver_pt ()
 Return a pointer to the linear solver object. More...
 
LinearSolver *const & linear_solver_pt () const
 Return a pointer to the linear solver object (const version) More...
 
LinearSolver *& mass_matrix_solver_for_explicit_timestepper_pt ()
 
LinearSolvermass_matrix_solver_for_explicit_timestepper_pt () const
 
EigenSolver *& eigen_solver_pt ()
 Return a pointer to the eigen solver object. More...
 
EigenSolver *const & eigen_solver_pt () const
 Return a pointer to the eigen solver object (const version) More...
 
Time *& time_pt ()
 Return a pointer to the global time object. More...
 
Timetime_pt () const
 Return a pointer to the global time object (const version). More...
 
doubletime ()
 Return the current value of continuous time. More...
 
double time () const
 Return the current value of continuous time (const version) More...
 
TimeStepper *& time_stepper_pt ()
 
const TimeSteppertime_stepper_pt () const
 
TimeStepper *& time_stepper_pt (const unsigned &i)
 Return a pointer to the i-th timestepper. More...
 
ExplicitTimeStepper *& explicit_time_stepper_pt ()
 Return a pointer to the explicit timestepper. More...
 
unsigned long set_timestepper_for_all_data (TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data=false)
 
virtual void shift_time_values ()
 Shift all values along to prepare for next timestep. More...
 
AssemblyHandler *& assembly_handler_pt ()
 Return a pointer to the assembly handler object. More...
 
AssemblyHandler *const & assembly_handler_pt () const
 Return a pointer to the assembly handler object (const version) More...
 
doubleminimum_dt ()
 Access function to min timestep in adaptive timestepping. More...
 
doublemaximum_dt ()
 Access function to max timestep in adaptive timestepping. More...
 
unsignedmax_newton_iterations ()
 Access function to max Newton iterations before giving up. More...
 
void problem_is_nonlinear (const bool &prob_lin)
 Access function to Problem_is_nonlinear. More...
 
doublemax_residuals ()
 
booltime_adaptive_newton_crash_on_solve_fail ()
 Access function for Time_adaptive_newton_crash_on_solve_fail. More...
 
doublenewton_solver_tolerance ()
 
void add_time_stepper_pt (TimeStepper *const &time_stepper_pt)
 
void set_explicit_time_stepper_pt (ExplicitTimeStepper *const &explicit_time_stepper_pt)
 
void initialise_dt (const double &dt)
 
void initialise_dt (const Vector< double > &dt)
 
Data *& global_data_pt (const unsigned &i)
 Return a pointer to the the i-th global data object. More...
 
void add_global_data (Data *const &global_data_pt)
 
void flush_global_data ()
 
LinearAlgebraDistribution *const & dof_distribution_pt () const
 Return the pointer to the dof distribution (read-only) More...
 
unsigned long ndof () const
 Return the number of dofs. More...
 
unsigned ntime_stepper () const
 Return the number of time steppers. More...
 
unsigned nglobal_data () const
 Return the number of global data values. More...
 
unsigned self_test ()
 Self-test: Check meshes and global data. Return 0 for OK. More...
 
void enable_store_local_dof_pt_in_elements ()
 
void disable_store_local_dof_pt_in_elements ()
 
unsigned long assign_eqn_numbers (const bool &assign_local_eqn_numbers=true)
 
void describe_dofs (std::ostream &out= *(oomph_info.stream_pt())) const
 
void enable_discontinuous_formulation ()
 
void disable_discontinuous_formulation ()
 
void get_dofs (DoubleVector &dofs) const
 
void get_dofs (const unsigned &t, DoubleVector &dofs) const
 Return vector of the t'th history value of all dofs. More...
 
void set_dofs (const DoubleVector &dofs)
 Set the values of the dofs. More...
 
void set_dofs (const unsigned &t, DoubleVector &dofs)
 Set the history values of the dofs. More...
 
void set_dofs (const unsigned &t, Vector< double * > &dof_pt)
 
void add_to_dofs (const double &lambda, const DoubleVector &increment_dofs)
 Add lambda x incremenet_dofs[l] to the l-th dof. More...
 
doubleglobal_dof_pt (const unsigned &i)
 
doubledof (const unsigned &i)
 i-th dof in the problem More...
 
double dof (const unsigned &i) const
 i-th dof in the problem (const version) More...
 
double *& dof_pt (const unsigned &i)
 Pointer to i-th dof in the problem. More...
 
doubledof_pt (const unsigned &i) const
 Pointer to i-th dof in the problem (const version) More...
 
virtual void get_inverse_mass_matrix_times_residuals (DoubleVector &Mres)
 
virtual void get_dvaluesdt (DoubleVector &f)
 
virtual void get_residuals (DoubleVector &residuals)
 Get the total residuals Vector for the problem. More...
 
virtual void get_jacobian (DoubleVector &residuals, DenseDoubleMatrix &jacobian)
 
virtual void get_jacobian (DoubleVector &residuals, CRDoubleMatrix &jacobian)
 
virtual void get_jacobian (DoubleVector &residuals, CCDoubleMatrix &jacobian)
 
virtual void get_jacobian (DoubleVector &residuals, SumOfMatrices &jacobian)
 
void get_fd_jacobian (DoubleVector &residuals, DenseMatrix< double > &jacobian)
 Get the full Jacobian by finite differencing. More...
 
void get_derivative_wrt_global_parameter (double *const &parameter_pt, DoubleVector &result)
 
void get_hessian_vector_products (DoubleVectorWithHaloEntries const &Y, Vector< DoubleVectorWithHaloEntries > const &C, Vector< DoubleVectorWithHaloEntries > &product)
 
void solve_eigenproblem (const unsigned &n_eval, Vector< std::complex< double >> &eigenvalue, Vector< DoubleVector > &eigenvector, const bool &steady=true)
 Solve the eigenproblem. More...
 
void solve_eigenproblem (const unsigned &n_eval, Vector< std::complex< double >> &eigenvalue, const bool &steady=true)
 
virtual void get_eigenproblem_matrices (CRDoubleMatrix &mass_matrix, CRDoubleMatrix &main_matrix, const double &shift=0.0)
 
void assign_eigenvector_to_dofs (DoubleVector &eigenvector)
 Assign the eigenvector passed to the function to the dofs. More...
 
void add_eigenvector_to_dofs (const double &epsilon, const DoubleVector &eigenvector)
 
void store_current_dof_values ()
 Store the current values of the degrees of freedom. More...
 
void restore_dof_values ()
 Restore the stored values of the degrees of freedom. More...
 
void enable_jacobian_reuse ()
 
void disable_jacobian_reuse ()
 Disable recycling of Jacobian in Newton iteration. More...
 
bool jacobian_reuse_is_enabled ()
 Is recycling of Jacobian in Newton iteration enabled? More...
 
booluse_predictor_values_as_initial_guess ()
 
void newton_solve ()
 Use Newton method to solve the problem. More...
 
void enable_globally_convergent_newton_method ()
 enable globally convergent Newton method More...
 
void disable_globally_convergent_newton_method ()
 disable globally convergent Newton method More...
 
void newton_solve (unsigned const &max_adapt)
 
void steady_newton_solve (unsigned const &max_adapt=0)
 
void copy (Problem *orig_problem_pt)
 
virtual Problemmake_copy ()
 
virtual void read (std::ifstream &restart_file, bool &unsteady_restart)
 
virtual void read (std::ifstream &restart_file)
 
virtual void dump (std::ofstream &dump_file) const
 
void dump (const std::string &dump_file_name) const
 
void delete_all_external_storage ()
 
virtual void symmetrise_eigenfunction_for_adaptive_pitchfork_tracking ()
 
doublebifurcation_parameter_pt () const
 
void get_bifurcation_eigenfunction (Vector< DoubleVector > &eigenfunction)
 
void activate_fold_tracking (double *const &parameter_pt, const bool &block_solve=true)
 
void activate_bifurcation_tracking (double *const &parameter_pt, const DoubleVector &eigenvector, const bool &block_solve=true)
 
void activate_bifurcation_tracking (double *const &parameter_pt, const DoubleVector &eigenvector, const DoubleVector &normalisation, const bool &block_solve=true)
 
void activate_pitchfork_tracking (double *const &parameter_pt, const DoubleVector &symmetry_vector, const bool &block_solve=true)
 
void activate_hopf_tracking (double *const &parameter_pt, const bool &block_solve=true)
 
void activate_hopf_tracking (double *const &parameter_pt, const double &omega, const DoubleVector &null_real, const DoubleVector &null_imag, const bool &block_solve=true)
 
void deactivate_bifurcation_tracking ()
 
void reset_assembly_handler_to_default ()
 Reset the system to the standard non-augemented state. More...
 
double arc_length_step_solve (double *const &parameter_pt, const double &ds, const unsigned &max_adapt=0)
 
double arc_length_step_solve (Data *const &data_pt, const unsigned &data_index, const double &ds, const unsigned &max_adapt=0)
 
void reset_arc_length_parameters ()
 
intsign_of_jacobian ()
 
void explicit_timestep (const double &dt, const bool &shift_values=true)
 Take an explicit timestep of size dt. More...
 
void unsteady_newton_solve (const double &dt)
 
void unsteady_newton_solve (const double &dt, const bool &shift_values)
 
void unsteady_newton_solve (const double &dt, const unsigned &max_adapt, const bool &first, const bool &shift=true)
 
double doubly_adaptive_unsteady_newton_solve (const double &dt, const double &epsilon, const unsigned &max_adapt, const bool &first, const bool &shift=true)
 
double doubly_adaptive_unsteady_newton_solve (const double &dt, const double &epsilon, const unsigned &max_adapt, const unsigned &suppress_resolve_after_spatial_adapt_flag, const bool &first, const bool &shift=true)
 
double adaptive_unsteady_newton_solve (const double &dt_desired, const double &epsilon)
 
double adaptive_unsteady_newton_solve (const double &dt_desired, const double &epsilon, const bool &shift_values)
 
void assign_initial_values_impulsive ()
 
void assign_initial_values_impulsive (const double &dt)
 
void calculate_predictions ()
 Calculate predictions. More...
 
void enable_mass_matrix_reuse ()
 
void disable_mass_matrix_reuse ()
 
bool mass_matrix_reuse_is_enabled ()
 Return whether the mass matrix is being reused. More...
 
void refine_uniformly (const Vector< unsigned > &nrefine_for_mesh)
 
void refine_uniformly (const Vector< unsigned > &nrefine_for_mesh, DocInfo &doc_info)
 
void refine_uniformly_and_prune (const Vector< unsigned > &nrefine_for_mesh)
 
void refine_uniformly_and_prune (const Vector< unsigned > &nrefine_for_mesh, DocInfo &doc_info)
 
void refine_uniformly (DocInfo &doc_info)
 
void refine_uniformly_and_prune (DocInfo &doc_info)
 
void refine_uniformly ()
 
void refine_uniformly (const unsigned &i_mesh, DocInfo &doc_info)
 Do uniform refinement for submesh i_mesh with documentation. More...
 
void refine_uniformly (const unsigned &i_mesh)
 Do uniform refinement for submesh i_mesh without documentation. More...
 
void p_refine_uniformly (const Vector< unsigned > &nrefine_for_mesh)
 
void p_refine_uniformly (const Vector< unsigned > &nrefine_for_mesh, DocInfo &doc_info)
 
void p_refine_uniformly_and_prune (const Vector< unsigned > &nrefine_for_mesh)
 
void p_refine_uniformly_and_prune (const Vector< unsigned > &nrefine_for_mesh, DocInfo &doc_info)
 
void p_refine_uniformly (DocInfo &doc_info)
 
void p_refine_uniformly_and_prune (DocInfo &doc_info)
 
void p_refine_uniformly ()
 
void p_refine_uniformly (const unsigned &i_mesh, DocInfo &doc_info)
 Do uniform p-refinement for submesh i_mesh with documentation. More...
 
void p_refine_uniformly (const unsigned &i_mesh)
 Do uniform p-refinement for submesh i_mesh without documentation. More...
 
void refine_selected_elements (const Vector< unsigned > &elements_to_be_refined)
 
void refine_selected_elements (const Vector< RefineableElement * > &elements_to_be_refined_pt)
 
void refine_selected_elements (const unsigned &i_mesh, const Vector< unsigned > &elements_to_be_refined)
 
void refine_selected_elements (const unsigned &i_mesh, const Vector< RefineableElement * > &elements_to_be_refined_pt)
 
void refine_selected_elements (const Vector< Vector< unsigned >> &elements_to_be_refined)
 
void refine_selected_elements (const Vector< Vector< RefineableElement * >> &elements_to_be_refined_pt)
 
void p_refine_selected_elements (const Vector< unsigned > &elements_to_be_refined)
 
void p_refine_selected_elements (const Vector< PRefineableElement * > &elements_to_be_refined_pt)
 
void p_refine_selected_elements (const unsigned &i_mesh, const Vector< unsigned > &elements_to_be_refined)
 
void p_refine_selected_elements (const unsigned &i_mesh, const Vector< PRefineableElement * > &elements_to_be_refined_pt)
 
void p_refine_selected_elements (const Vector< Vector< unsigned >> &elements_to_be_refined)
 
void p_refine_selected_elements (const Vector< Vector< PRefineableElement * >> &elements_to_be_refined_pt)
 
unsigned unrefine_uniformly ()
 
unsigned unrefine_uniformly (const unsigned &i_mesh)
 
void p_unrefine_uniformly (DocInfo &doc_info)
 
void p_unrefine_uniformly (const unsigned &i_mesh, DocInfo &doc_info)
 Do uniform p-unrefinement for submesh i_mesh without documentation. More...
 
void adapt (unsigned &n_refined, unsigned &n_unrefined)
 
void adapt ()
 
void p_adapt (unsigned &n_refined, unsigned &n_unrefined)
 
void p_adapt ()
 
void adapt_based_on_error_estimates (unsigned &n_refined, unsigned &n_unrefined, Vector< Vector< double >> &elemental_error)
 
void adapt_based_on_error_estimates (Vector< Vector< double >> &elemental_error)
 
void get_all_error_estimates (Vector< Vector< double >> &elemental_error)
 
void doc_errors (DocInfo &doc_info)
 Get max and min error for all elements in submeshes. More...
 
void doc_errors ()
 Get max and min error for all elements in submeshes. More...
 
void enable_info_in_newton_solve ()
 
void disable_info_in_newton_solve ()
 Disable the output of information when in the newton solver. More...
 
- Public Member Functions inherited from oomph::ExplicitTimeSteppableObject
 ExplicitTimeSteppableObject ()
 Empty constructor. More...
 
 ExplicitTimeSteppableObject (const ExplicitTimeSteppableObject &)=delete
 Broken copy constructor. More...
 
void operator= (const ExplicitTimeSteppableObject &)=delete
 Broken assignment operator. More...
 
virtual ~ExplicitTimeSteppableObject ()
 Empty destructor. More...
 
virtual void actions_before_explicit_stage ()
 
virtual void actions_after_explicit_stage ()
 

Private Member Functions

void fix_pressure (const unsigned &e, const unsigned &pdof, const double &pvalue)
 Fix pressure in element e at pressure dof pdof and set to pvalue. More...
 
void fix_pressure (const unsigned &e, const unsigned &pdof, const double &pvalue)
 Fix pressure in element e at pressure dof pdof and set to pvalue. More...
 
void fix_pressure (const unsigned &e, const unsigned &pdof, const double &pvalue)
 Fix pressure in element e at pressure dof pdof and set to pvalue. More...
 

Private Attributes

NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt Body_force_fct_pt
 Pointer to body force function. More...
 
unsigned Driven_boundary
 Which boundary is driven? More...
 

Additional Inherited Members

- Public Types inherited from oomph::Problem
typedef void(* SpatialErrorEstimatorFctPt) (Mesh *&mesh_pt, Vector< double > &elemental_error)
 Function pointer for spatial error estimator. More...
 
typedef void(* SpatialErrorEstimatorWithDocFctPt) (Mesh *&mesh_pt, Vector< double > &elemental_error, DocInfo &doc_info)
 Function pointer for spatial error estimator with doc. More...
 
- Public Attributes inherited from oomph::Problem
bool Shut_up_in_newton_solve
 
- Static Public Attributes inherited from oomph::Problem
static bool Suppress_warning_about_actions_before_read_unstructured_meshes
 
- Protected Types inherited from oomph::Problem
enum  Assembly_method {
  Perform_assembly_using_vectors_of_pairs , Perform_assembly_using_two_vectors , Perform_assembly_using_maps , Perform_assembly_using_lists ,
  Perform_assembly_using_two_arrays
}
 Enumerated flags to determine which sparse assembly method is used. More...
 
- Protected Member Functions inherited from oomph::Problem
unsigned setup_element_count_per_dof ()
 
virtual void sparse_assemble_row_or_column_compressed (Vector< int * > &column_or_row_index, Vector< int * > &row_or_column_start, Vector< double * > &value, Vector< unsigned > &nnz, Vector< double * > &residual, bool compressed_row_flag)
 
virtual void actions_before_newton_convergence_check ()
 
virtual void actions_before_newton_step ()
 
virtual void actions_after_newton_step ()
 
virtual void actions_before_implicit_timestep ()
 
virtual void actions_after_implicit_timestep ()
 
virtual void actions_after_implicit_timestep_and_error_estimation ()
 
virtual void actions_before_explicit_timestep ()
 Actions that should be performed before each explicit time step. More...
 
virtual void actions_after_explicit_timestep ()
 Actions that should be performed after each explicit time step. More...
 
virtual void actions_before_read_unstructured_meshes ()
 
virtual void actions_after_read_unstructured_meshes ()
 
virtual void actions_after_change_in_global_parameter (double *const &parameter_pt)
 
virtual void actions_after_change_in_bifurcation_parameter ()
 
virtual void actions_after_parameter_increase (double *const &parameter_pt)
 
doubledof_derivative (const unsigned &i)
 
doubledof_current (const unsigned &i)
 
virtual void set_initial_condition ()
 
virtual double global_temporal_error_norm ()
 
unsigned newton_solve_continuation (double *const &parameter_pt)
 
unsigned newton_solve_continuation (double *const &parameter_pt, DoubleVector &z)
 
void calculate_continuation_derivatives (double *const &parameter_pt)
 
void calculate_continuation_derivatives (const DoubleVector &z)
 
void calculate_continuation_derivatives_fd (double *const &parameter_pt)
 
bool does_pointer_correspond_to_problem_data (double *const &parameter_pt)
 
void set_consistent_pinned_values_for_continuation ()
 
- Protected Attributes inherited from oomph::Problem
Vector< Problem * > Copy_of_problem_pt
 
std::map< double *, boolCalculate_dparameter_analytic
 
bool Calculate_hessian_products_analytic
 
LinearAlgebraDistributionDof_distribution_pt
 
Vector< double * > Dof_pt
 Vector of pointers to dofs. More...
 
DoubleVectorWithHaloEntries Element_count_per_dof
 
double Relaxation_factor
 
double Newton_solver_tolerance
 
unsigned Max_newton_iterations
 Maximum number of Newton iterations. More...
 
unsigned Nnewton_iter_taken
 
Vector< doubleMax_res
 Maximum residuals at start and after each newton iteration. More...
 
double Max_residuals
 
bool Time_adaptive_newton_crash_on_solve_fail
 
bool Jacobian_reuse_is_enabled
 Is re-use of Jacobian in Newton iteration enabled? Default: false. More...
 
bool Jacobian_has_been_computed
 
bool Problem_is_nonlinear
 
bool Pause_at_end_of_sparse_assembly
 
bool Doc_time_in_distribute
 
unsigned Sparse_assembly_method
 
unsigned Sparse_assemble_with_arrays_initial_allocation
 
unsigned Sparse_assemble_with_arrays_allocation_increment
 
Vector< Vector< unsigned > > Sparse_assemble_with_arrays_previous_allocation
 
double Numerical_zero_for_sparse_assembly
 
double FD_step_used_in_get_hessian_vector_products
 
bool Mass_matrix_reuse_is_enabled
 
bool Mass_matrix_has_been_computed
 
bool Discontinuous_element_formulation
 
double Minimum_dt
 Minimum desired dt: if dt falls below this value, exit. More...
 
double Maximum_dt
 Maximum desired dt. More...
 
double DTSF_max_increase
 
double DTSF_min_decrease
 
double Minimum_dt_but_still_proceed
 
bool Scale_arc_length
 Boolean to control whether arc-length should be scaled. More...
 
double Desired_proportion_of_arc_length
 Proportion of the arc-length to taken by the parameter. More...
 
double Theta_squared
 
int Sign_of_jacobian
 Storage for the sign of the global Jacobian. More...
 
double Continuation_direction
 
double Parameter_derivative
 Storage for the derivative of the global parameter wrt arc-length. More...
 
double Parameter_current
 Storage for the present value of the global parameter. More...
 
bool Use_continuation_timestepper
 Boolean to control original or new storage of dof stuff. More...
 
unsigned Dof_derivative_offset
 
unsigned Dof_current_offset
 
Vector< doubleDof_derivative
 Storage for the derivative of the problem variables wrt arc-length. More...
 
Vector< doubleDof_current
 Storage for the present values of the variables. More...
 
double Ds_current
 Storage for the current step value. More...
 
unsigned Desired_newton_iterations_ds
 
double Minimum_ds
 Minimum desired value of arc-length. More...
 
bool Bifurcation_detection
 Boolean to control bifurcation detection via determinant of Jacobian. More...
 
bool Bisect_to_find_bifurcation
 Boolean to control wheter bisection is used to located bifurcation. More...
 
bool First_jacobian_sign_change
 Boolean to indicate whether a sign change has occured in the Jacobian. More...
 
bool Arc_length_step_taken
 Boolean to indicate whether an arc-length step has been taken. More...
 
bool Use_finite_differences_for_continuation_derivatives
 
OomphCommunicatorCommunicator_pt
 The communicator for this problem. More...
 
bool Always_take_one_newton_step
 
double Timestep_reduction_factor_after_nonconvergence
 
bool Keep_temporal_error_below_tolerance
 
- Static Protected Attributes inherited from oomph::Problem
static ContinuationStorageScheme Continuation_time_stepper
 Storage for the single static continuation timestorage object. More...
 

Detailed Description

template<class ELEMENT>
class QuarterCircleDrivenCavityProblem< ELEMENT >

/////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////// Driven cavity problem in quarter circle domain, templated by element type.

Driven cavity problem in quarter circle domain, templated by element type.

Constructor & Destructor Documentation

◆ QuarterCircleDrivenCavityProblem() [1/3]

template<class ELEMENT >
QuarterCircleDrivenCavityProblem< ELEMENT >::QuarterCircleDrivenCavityProblem ( NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt  body_force_fct_pt)

Constructor.

Constructor for driven cavity problem in quarter circle domain.

189  :
190  Body_force_fct_pt(body_force_fct_pt)
191 {
192 
193  // Build geometric object that parametrises the curved boundary
194  // of the domain
195 
196  // Half axes for ellipse
197  double a_ellipse=1.0;
198  double b_ellipse=1.0;
199 
200  // Setup elliptical ring
201  GeomObject* Wall_pt=new Ellipse(a_ellipse,b_ellipse);
202 
203  // End points for wall
204  double xi_lo=0.0;
205  double xi_hi=2.0*atan(1.0);
206 
207  //Now create the mesh
208  double fract_mid=0.5;
209  Problem::mesh_pt() = new
211  Wall_pt,xi_lo,fract_mid,xi_hi);
212 
213  // Set error estimator
216  mesh_pt())->spatial_error_estimator_pt()=error_estimator_pt;
217 
218  // Set the boundary conditions for this problem: All nodes are
219  // free by default -- just pin the ones that have Dirichlet conditions
220  // here: All boundaries are Dirichlet boundaries.
221  unsigned num_bound = mesh_pt()->nboundary();
222  for(unsigned ibound=0;ibound<num_bound;ibound++)
223  {
224  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
225  for (unsigned inod=0;inod<num_nod;inod++)
226  {
227  // Loop over values (u and v velocities)
228  for (unsigned i=0;i<2;i++)
229  {
230  mesh_pt()->boundary_node_pt(ibound,inod)->pin(i);
231  }
232  }
233  } // end loop over boundaries
234 
235  //Find number of elements in mesh
236  unsigned n_element = mesh_pt()->nelement();
237 
238  // Loop over the elements to set up element-specific
239  // things that cannot be handled by constructor: Pass pointer to Reynolds
240  // number
241  for(unsigned e=0;e<n_element;e++)
242  {
243  // Upcast from GeneralisedElement to the present element
244  ELEMENT* el_pt = dynamic_cast<ELEMENT*>(mesh_pt()->element_pt(e));
245 
246  //Set the Reynolds number, etc
247  el_pt->re_pt() = &Global_Physical_Variables::Re;
248  //Set the Re/Fr
249  el_pt->re_invfr_pt() = &Global_Physical_Variables::Re_invFr;
250  //Set Gravity vector
251  el_pt->g_pt() = &Global_Physical_Variables::Gravity;
252  //set body force function
253  el_pt->body_force_fct_pt() = Body_force_fct_pt;
254 
255  } // end loop over elements
256 
257  // Initial refinement level
260 
261  // Pin redudant pressure dofs
263  pin_redundant_nodal_pressures(mesh_pt()->element_pt());
264 
265  // Now pin the first pressure dof in the first element and set it to 0.0
266  fix_pressure(0,0,0.0);
267 
268  // Setup equation numbering scheme
269  cout <<"Number of equations: " << assign_eqn_numbers() << std::endl;
270 
271 } // end_of_constructor
int i
Definition: BiCGSTAB_step_by_step.cpp:9
Array< double, 1, 3 > e(1./3., 0.5, 2.)
void fix_pressure(const unsigned &e, const unsigned &pdof, const double &pvalue)
Fix pressure in element e at pressure dof pdof and set to pvalue.
Definition: mpi/distribution/circular_driven_cavity/circular_driven_cavity.cc:172
NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt Body_force_fct_pt
Pointer to body force function.
Definition: mpi/distribution/circular_driven_cavity/circular_driven_cavity.cc:169
void pin(const unsigned &i)
Pin the i-th stored variable.
Definition: nodes.h:385
Definition: geom_objects.h:644
Definition: geom_objects.h:101
unsigned long nboundary_node(const unsigned &ibound) const
Return number of nodes on a particular boundary.
Definition: mesh.h:833
unsigned nboundary() const
Return number of boundaries.
Definition: mesh.h:827
GeneralisedElement *& element_pt(const unsigned long &e)
Return pointer to element e.
Definition: mesh.h:448
Node *& boundary_node_pt(const unsigned &b, const unsigned &n)
Return pointer to node n on boundary b.
Definition: mesh.h:493
unsigned long nelement() const
Return number of elements in the mesh.
Definition: mesh.h:590
Mesh *& mesh_pt()
Return a pointer to the global mesh.
Definition: problem.h:1280
unsigned long assign_eqn_numbers(const bool &assign_local_eqn_numbers=true)
Definition: problem.cc:1989
void refine_uniformly()
Definition: problem.h:2640
Definition: refineable_navier_stokes_elements.h:322
Definition: quarter_circle_sector_mesh.template.h:188
Definition: error_estimator.h:266
EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC bfloat16 atan(const bfloat16 &a)
Definition: BFloat16.h:636
double Gravity
Non-dim gravity.
Definition: meshing/quad_from_triangle_mesh/unstructured_two_d_solid.cc:214
double Re_invFr
Reynolds/Froude number.
Definition: mpi/distribution/circular_driven_cavity/circular_driven_cavity.cc:52
double Re
Reynolds number.
Definition: fibre.cc:55
Z2ErrorEstimator * error_estimator_pt
Definition: MortaringCantileverCompareToNonMortaring.cpp:190

References oomph::Problem::assign_eqn_numbers(), Eigen::bfloat16_impl::atan(), QuarterCircleDrivenCavityProblem< ELEMENT >::Body_force_fct_pt, oomph::Mesh::boundary_node_pt(), e(), oomph::Mesh::element_pt(), MeshRefinement::error_estimator_pt, QuarterCircleDrivenCavityProblem< ELEMENT >::fix_pressure(), Global_Physical_Variables::Gravity, i, oomph::Problem::mesh_pt(), oomph::Mesh::nboundary(), oomph::Mesh::nboundary_node(), oomph::Mesh::nelement(), oomph::Data::pin(), Global_Physical_Variables::Re, Global_Physical_Variables::Re_invFr, and oomph::Problem::refine_uniformly().

◆ ~QuarterCircleDrivenCavityProblem() [1/3]

template<class ELEMENT >
QuarterCircleDrivenCavityProblem< ELEMENT >::~QuarterCircleDrivenCavityProblem ( )
inline

Destructor: Empty – all memory gets cleaned up in base destructor.

94 {}

◆ QuarterCircleDrivenCavityProblem() [2/3]

template<class ELEMENT >
QuarterCircleDrivenCavityProblem< ELEMENT >::QuarterCircleDrivenCavityProblem ( NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt  body_force_fct_pt)

Constructor.

◆ ~QuarterCircleDrivenCavityProblem() [2/3]

template<class ELEMENT >
QuarterCircleDrivenCavityProblem< ELEMENT >::~QuarterCircleDrivenCavityProblem ( )
inline

Destructor: Empty.

94 {}

◆ QuarterCircleDrivenCavityProblem() [3/3]

template<class ELEMENT >
QuarterCircleDrivenCavityProblem< ELEMENT >::QuarterCircleDrivenCavityProblem ( NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt  body_force_fct_pt,
const unsigned driven_bdry 
)

Constructor.

Constructor for driven cavity problem in quarter circle domain.

218  :
219  Body_force_fct_pt(body_force_fct_pt), Driven_boundary(driven_bdry)
220 {
221 
222  // Build geometric object that parametrises the curved boundary
223  // of the domain
224 
225  // Half axes for ellipse
226  double a_ellipse=1.0;
227  double b_ellipse=1.0;
228 
229  // Setup elliptical ring
230  GeomObject* Wall_pt=new Ellipse(a_ellipse,b_ellipse);
231 
232  // End points for wall
233  double xi_lo=0.0;
234  double xi_hi=2.0*atan(1.0);
235 
236  //Now create the mesh
237  double fract_mid=0.5;
238  Problem::mesh_pt() = new
240  Wall_pt,xi_lo,fract_mid,xi_hi);
241 
242  // Raise the max refinement level
243  dynamic_cast<TreeBasedRefineableMeshBase*>(mesh_pt())->max_refinement_level()=7;
244 
245  // Set error estimator
248  mesh_pt())->spatial_error_estimator_pt()=error_estimator_pt;
249 
250  // Set the boundary conditions for this problem: All nodes are
251  // free by default -- just pin the ones that have Dirichlet conditions
252  // here: All boundaries are Dirichlet boundaries.
253  unsigned num_bound = mesh_pt()->nboundary();
254  for(unsigned ibound=0;ibound<num_bound;ibound++)
255  {
256  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
257  for (unsigned inod=0;inod<num_nod;inod++)
258  {
259  // Loop over values (u and v velocities)
260  for (unsigned i=0;i<2;i++)
261  {
262  mesh_pt()->boundary_node_pt(ibound,inod)->pin(i);
263  }
264  }
265  } // end loop over boundaries
266 
267  //Find number of elements in mesh
268  unsigned n_element = mesh_pt()->nelement();
269 
270  // Loop over the elements to set up element-specific
271  // things that cannot be handled by constructor: Pass pointer to Reynolds
272  // number
273  for(unsigned e=0;e<n_element;e++)
274  {
275  // Upcast from GeneralisedElement to the present element
276  ELEMENT* el_pt = dynamic_cast<ELEMENT*>(mesh_pt()->element_pt(e));
277 
278  //Set the Reynolds number, etc
279  el_pt->re_pt() = &Global_Physical_Variables::Re;
280  //Set the Re/Fr
281  el_pt->re_invfr_pt() = &Global_Physical_Variables::Re_invFr;
282  //Set Gravity vector
283  el_pt->g_pt() = &Global_Physical_Variables::Gravity;
284  //set body force function
285  el_pt->body_force_fct_pt() = Body_force_fct_pt;
286 
287  } // end loop over elements
288 
289  // Pin redudant pressure dofs
291  pin_redundant_nodal_pressures(mesh_pt()->element_pt());
292 
293  // Now pin the first pressure dof in the first element and set it to 0.0
294  fix_pressure(0,0,0.0);
295 
296  // Setup equation numbering scheme
297  cout <<"Number of equations: " << assign_eqn_numbers() << std::endl;
298 
299 } // end_of_constructor
unsigned Driven_boundary
Which boundary is driven?
Definition: circular_driven_cavity_hp_adapt.cc:197
Base class for tree-based refineable meshes.
Definition: refineable_mesh.h:376

References oomph::Problem::assign_eqn_numbers(), Eigen::bfloat16_impl::atan(), QuarterCircleDrivenCavityProblem< ELEMENT >::Body_force_fct_pt, oomph::Mesh::boundary_node_pt(), e(), oomph::Mesh::element_pt(), MeshRefinement::error_estimator_pt, QuarterCircleDrivenCavityProblem< ELEMENT >::fix_pressure(), Global_Physical_Variables::Gravity, i, oomph::Problem::mesh_pt(), oomph::Mesh::nboundary(), oomph::Mesh::nboundary_node(), oomph::Mesh::nelement(), oomph::Data::pin(), Global_Physical_Variables::Re, and Global_Physical_Variables::Re_invFr.

◆ ~QuarterCircleDrivenCavityProblem() [3/3]

template<class ELEMENT >
QuarterCircleDrivenCavityProblem< ELEMENT >::~QuarterCircleDrivenCavityProblem ( )
inline

Destructor: Empty.

91 {}

Member Function Documentation

◆ actions_after_adapt() [1/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_after_adapt ( )
inlinevirtual

After adaptation: Unpin pressure and pin redudant pressure dofs.

Reimplemented from oomph::Problem.

134  {
135  // Unpin all pressure dofs
137  unpin_all_pressure_dofs(mesh_pt()->element_pt());
138 
139  // Pin redundant pressure dofs
141  pin_redundant_nodal_pressures(mesh_pt()->element_pt());
142 
143  // Now pin the first pressure dof in the first element and set it to 0.0
144 
145  // Loop over all elements
146  unsigned n_element=mesh_pt()->nelement();
147  for (unsigned e=0;e<n_element;e++)
148  {
149  // If the lower left node of this element is (0,0), then fix the
150  // pressure dof in this element to zero
151  if (mesh_pt()->finite_element_pt(e)->node_pt(0)->x(0)==0.0 &&
152  mesh_pt()->finite_element_pt(e)->node_pt(0)->x(1)==0.0) // 2d problem
153  {
154  oomph_info << "I'm fixing the pressure " << std::endl;
155  // Fix the pressure in element e at pdof=0 to 0.0
156  unsigned pdof=0;
157  fix_pressure(e,pdof,0.0);
158  }
159  }
160 
161  } // end_of_actions_after_adapt
OomphInfo oomph_info
Definition: oomph_definitions.cc:319
list x
Definition: plotDoE.py:28

References e(), oomph::oomph_info, and plotDoE::x.

◆ actions_after_adapt() [2/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_after_adapt ( )
inlinevirtual

After adaptation: Unpin pressure and pin redudant pressure dofs.

Reimplemented from oomph::Problem.

134  {
135  // Unpin all pressure dofs
137  unpin_all_pressure_dofs(mesh_pt()->element_pt());
138 
139  // Pin redundant pressure dofs
141  pin_redundant_nodal_pressures(mesh_pt()->element_pt());
142 
143  // Now pin the first pressure dof in the first element and set it to 0.0
144  fix_pressure(0,0,0.0);
145  } // end_of_actions_after_adapt

◆ actions_after_adapt() [3/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_after_adapt ( )
inlinevirtual

After adaptation: Unpin pressure and pin redudant pressure dofs.

Reimplemented from oomph::Problem.

175  {
176  // Unpin all pressure dofs
178  unpin_all_pressure_dofs(mesh_pt()->element_pt());
179 
180  // Pin redundant pressure dofs
182  pin_redundant_nodal_pressures(mesh_pt()->element_pt());
183 
184  // Now pin the first pressure dof in the first element and set it to 0.0
185  fix_pressure(0,0,0.0);
186  } // end_of_actions_after_adapt

◆ actions_after_newton_solve() [1/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_after_newton_solve ( )
inlinevirtual

Update the after solve (empty)

Reimplemented from oomph::Problem.

97 {}

◆ actions_after_newton_solve() [2/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_after_newton_solve ( )
inlinevirtual

Update the after solve (empty)

Reimplemented from oomph::Problem.

97 {}

◆ actions_after_newton_solve() [3/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_after_newton_solve ( )
inlinevirtual

Update the after solve (empty)

Reimplemented from oomph::Problem.

94 {}

◆ actions_before_newton_solve() [1/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_before_newton_solve ( )
inlinevirtual

Update the problem specs before solve. (Re-)set velocity boundary conditions just to be on the safe side...

Reimplemented from oomph::Problem.

102  {
103  // Setup tangential flow along boundary 0:
104  unsigned ibound=0;
105  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
106  for (unsigned inod=0;inod<num_nod;inod++)
107  {
108  // Tangential flow
109  unsigned i=0;
110  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,1.0);
111  // No penetration
112  i=1;
113  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,0.0);
114  }
115 
116  // Overwrite with no flow along all other boundaries
117  unsigned num_bound = mesh_pt()->nboundary();
118  for(unsigned ibound=1;ibound<num_bound;ibound++)
119  {
120  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
121  for (unsigned inod=0;inod<num_nod;inod++)
122  {
123  for (unsigned i=0;i<2;i++)
124  {
125  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,0.0);
126  }
127  }
128  }
129  } // end_of_actions_before_newton_solve
void set_value(const unsigned &i, const double &value_)
Definition: nodes.h:271

References i.

◆ actions_before_newton_solve() [2/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_before_newton_solve ( )
inlinevirtual

Update the problem specs before solve. (Re-)set velocity boundary conditions just to be on the safe side...

Reimplemented from oomph::Problem.

102  {
103  // Setup tangential flow along boundary 0:
104  unsigned ibound=0;
105  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
106  for (unsigned inod=0;inod<num_nod;inod++)
107  {
108  // Tangential flow
109  unsigned i=0;
110  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,1.0);
111  // No penetration
112  i=1;
113  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,0.0);
114  }
115 
116  // Overwrite with no flow along all other boundaries
117  unsigned num_bound = mesh_pt()->nboundary();
118  for(unsigned ibound=1;ibound<num_bound;ibound++)
119  {
120  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
121  for (unsigned inod=0;inod<num_nod;inod++)
122  {
123  for (unsigned i=0;i<2;i++)
124  {
125  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,0.0);
126  }
127  }
128  }
129  } // end_of_actions_before_newton_solve

References i.

◆ actions_before_newton_solve() [3/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::actions_before_newton_solve ( )
inlinevirtual

Update the problem specs before solve. (Re-)set velocity boundary conditions just to be on the safe side...

Reimplemented from oomph::Problem.

99  {
100  if(Driven_boundary == 0)
101  {
102  // Setup tangential flow along boundary 0:
103  unsigned ibound=0;
104  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
105  for (unsigned inod=0;inod<num_nod;inod++)
106  {
107  // Tangential flow
108  unsigned i=0;
109  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,1.0);
110  // No penetration
111  i=1;
112  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,0.0);
113  }
114 
115  // Overwrite with no flow along all other boundaries
116  unsigned num_bound = mesh_pt()->nboundary();
117  for(unsigned ibound=1;ibound<num_bound;ibound++)
118  {
119  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
120  for (unsigned inod=0;inod<num_nod;inod++)
121  {
122  for (unsigned i=0;i<2;i++)
123  {
124  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,0.0);
125  }
126  }
127  }
128  }
129  else
130  {
131  // Setup tangential flow along boundary 1:
132  unsigned ibound=1;
133  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
134  for (unsigned inod=0;inod<num_nod;inod++)
135  {
136  // get coordinates
137  double x=mesh_pt()->boundary_node_pt(ibound,inod)->x(0);
138  double y=mesh_pt()->boundary_node_pt(ibound,inod)->x(1);
139  // find Lagrangian coordinate (the angle)
140  double zeta=0.0;
141  if (x!=0.0)
142  {
143  zeta=atan(y/x);
144  }
145  // Tangential flow u0
146  unsigned i=0;
147  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,-sin(zeta));
148  // Tangential flow u1
149  i=1;
150  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,cos(zeta));
151  }
152 
153  // Overwrite with no flow along all boundaries
154  unsigned num_bound = mesh_pt()->nboundary();
155  for(unsigned ibound=0;ibound<num_bound;ibound++)
156  {
157  if (ibound!=1)
158  {
159  unsigned num_nod= mesh_pt()->nboundary_node(ibound);
160  for (unsigned inod=0;inod<num_nod;inod++)
161  {
162  for (unsigned i=0;i<2;i++)
163  {
164  mesh_pt()->boundary_node_pt(ibound,inod)->set_value(i,0.0);
165  }
166  }
167  }
168  }
169  }
170  } // end_of_actions_before_newton_solve
AnnoyingScalar cos(const AnnoyingScalar &x)
Definition: AnnoyingScalar.h:136
AnnoyingScalar sin(const AnnoyingScalar &x)
Definition: AnnoyingScalar.h:137
double & x(const unsigned &i)
Return the i-th nodal coordinate.
Definition: nodes.h:1060
Scalar * y
Definition: level1_cplx_impl.h:128
EIGEN_STRONG_INLINE const Eigen::CwiseBinaryOp< Eigen::internal::scalar_zeta_op< typename DerivedX::Scalar >, const DerivedX, const DerivedQ > zeta(const Eigen::ArrayBase< DerivedX > &x, const Eigen::ArrayBase< DerivedQ > &q)
Definition: SpecialFunctionsArrayAPI.h:152

References Eigen::bfloat16_impl::atan(), cos(), i, sin(), plotDoE::x, y, and Eigen::zeta().

◆ doc_solution() [1/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::doc_solution ( DocInfo doc_info)

Doc the solution.

280 {
281 
282  ofstream some_file;
283  char filename[100];
284 
285  // Number of plot points
286  unsigned npts=5;
287 
288 
289  // Output solution
290  sprintf(filename,"%s/soln%i_on_proc%i.dat",doc_info.directory().c_str(),
291  doc_info.number(),this->communicator_pt()->my_rank());
292  some_file.open(filename);
293  mesh_pt()->output(some_file,npts);
294  some_file.close();
295 
296 } // end_of_doc_solution
std::string directory() const
Output directory.
Definition: oomph_utilities.h:524
unsigned & number()
Number used (e.g.) for labeling output files.
Definition: oomph_utilities.h:554
void output(std::ostream &outfile)
Output for all elements.
Definition: mesh.cc:2027
string filename
Definition: MergeRestartFiles.py:39

References oomph::DocInfo::directory(), MergeRestartFiles::filename, and oomph::DocInfo::number().

◆ doc_solution() [2/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::doc_solution ( DocInfo doc_info)

Doc the solution.

◆ doc_solution() [3/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::doc_solution ( DocInfo doc_info)

Doc the solution.

◆ fix_pressure() [1/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::fix_pressure ( const unsigned e,
const unsigned pdof,
const double pvalue 
)
inlineprivate

Fix pressure in element e at pressure dof pdof and set to pvalue.

174  {
175  //Cast to proper element and fix pressure
176  dynamic_cast<ELEMENT*>(mesh_pt()->element_pt(e))->
177  fix_pressure(pdof,pvalue);
178  } // end_of_fix_pressure

References e().

Referenced by QuarterCircleDrivenCavityProblem< ELEMENT >::QuarterCircleDrivenCavityProblem().

◆ fix_pressure() [2/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::fix_pressure ( const unsigned e,
const unsigned pdof,
const double pvalue 
)
inlineprivate

Fix pressure in element e at pressure dof pdof and set to pvalue.

158  {
159  //Cast to proper element and fix pressure
160  dynamic_cast<ELEMENT*>(mesh_pt()->element_pt(e))->
161  fix_pressure(pdof,pvalue);
162  } // end_of_fix_pressure

References e().

◆ fix_pressure() [3/3]

template<class ELEMENT >
void QuarterCircleDrivenCavityProblem< ELEMENT >::fix_pressure ( const unsigned e,
const unsigned pdof,
const double pvalue 
)
inlineprivate

Fix pressure in element e at pressure dof pdof and set to pvalue.

202  {
203  //Cast to proper element and fix pressure
204  dynamic_cast<ELEMENT*>(mesh_pt()->element_pt(e))->
205  fix_pressure(pdof,pvalue);
206  } // end_of_fix_pressure

References e().

Member Data Documentation

◆ Body_force_fct_pt

template<class ELEMENT >
NavierStokesEquations< 2 >::NavierStokesBodyForceFctPt QuarterCircleDrivenCavityProblem< ELEMENT >::Body_force_fct_pt
private

◆ Driven_boundary

template<class ELEMENT >
unsigned QuarterCircleDrivenCavityProblem< ELEMENT >::Driven_boundary
private

Which boundary is driven?


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