Inheritance diagram for OpCalculateFStab< DIM_0, DIM_1 >:

Collaboration diagram for OpCalculateFStab< DIM_0, DIM_1 >:

Public Member Functions
	OpCalculateFStab (boost::shared_ptr< MatrixDouble > def_grad_ptr, boost::shared_ptr< MatrixDouble > def_grad_stab_ptr, boost::shared_ptr< MatrixDouble > def_grad_dot_ptr, double tau_F_ptr, double xi_F_ptr, boost::shared_ptr< MatrixDouble > grad_x_ptr, boost::shared_ptr< MatrixDouble > grad_vel_ptr)

MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Private Attributes
double	tauFPtr

double	xiF

boost::shared_ptr< MatrixDouble >	defGradPtr

boost::shared_ptr< MatrixDouble >	defGradStabPtr

boost::shared_ptr< MatrixDouble >	defGradDotPtr

boost::shared_ptr< MatrixDouble >	gradxPtr

boost::shared_ptr< MatrixDouble >	gradVelPtr

Detailed Description

template<int DIM_0, int DIM_1>
struct OpCalculateFStab< DIM_0, DIM_1 >

Examples: mofem/tutorials/adv-4_dynamic_first_order_con_law/dynamic_first_order_con_law.cpp.

Definition at line 104 of file dynamic_first_order_con_law.cpp.

Constructor & Destructor Documentation

◆ OpCalculateFStab()

template<int DIM_0, int DIM_1>

OpCalculateFStab< DIM_0, DIM_1 >::OpCalculateFStab	(	boost::shared_ptr< MatrixDouble >	def_grad_ptr,
		boost::shared_ptr< MatrixDouble >	def_grad_stab_ptr,
		boost::shared_ptr< MatrixDouble >	def_grad_dot_ptr,
		double	tau_F_ptr,
		double	xi_F_ptr,
		boost::shared_ptr< MatrixDouble >	grad_x_ptr,
		boost::shared_ptr< MatrixDouble >	grad_vel_ptr
	)

inline

Definition at line 105 of file dynamic_first_order_con_law.cpp.

      : ForcesAndSourcesCore::UserDataOperator(NOSPACE, OPLAST),
        defGradPtr(def_grad_ptr), defGradStabPtr(def_grad_stab_ptr),
        defGradDotPtr(def_grad_dot_ptr), tauFPtr(tau_F_ptr), xiF(xi_F_ptr),
        gradxPtr(grad_x_ptr), gradVelPtr(grad_vel_ptr) {}

Member Function Documentation

◆ doWork()

template<int DIM_0, int DIM_1>

MoFEMErrorCode OpCalculateFStab< DIM_0, DIM_1 >::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

inline

Examples: mofem/tutorials/adv-4_dynamic_first_order_con_law/dynamic_first_order_con_law.cpp.

Definition at line 116 of file dynamic_first_order_con_law.cpp.

                                                               {
    MoFEMFunctionBegin;
    // Define Indicies
    FTensor::Index<'i', SPACE_DIM> i;
    FTensor::Index<'j', SPACE_DIM> j;
 
    // Number of Gauss points
    const size_t nb_gauss_pts = getGaussPts().size2();
 
    defGradStabPtr->resize(nb_gauss_pts, DIM_0 * DIM_1, false);
    defGradStabPtr->clear();
 
    // Extract matrix from data matrix
    auto t_F = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*defGradPtr);
    auto t_Fstab = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*defGradStabPtr);
    auto t_F_dot = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*defGradDotPtr);
 
    // tau_F = alpha deltaT
    auto tau_F = tauFPtr;
    double xi_F = xiF;
    auto t_gradx = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*gradxPtr);
    auto t_gradVel = getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*gradVelPtr);
 
    for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
      // Stabilised Deformation Gradient
      t_Fstab(i, j) = t_F(i, j) + tau_F * (t_gradVel(i, j) - t_F_dot(i, j)) +
                      xi_F * (t_gradx(i, j) - t_F(i, j));
 
      ++t_F;
      ++t_Fstab;
      ++t_gradVel;
      ++t_F_dot;
 
      ++t_gradx;
    }
 
    MoFEMFunctionReturn(0);
  }