GetDeformedMesh_ALE Class Reference

#include <transfer_actions.h>

Inheritance diagram for GetDeformedMesh_ALE:

Collaboration diagram for GetDeformedMesh_ALE:

Public Member Functions
	GetDeformedMesh_ALE (FluidAgent *fag, SolidAgent *sag, const std::string s_x, const std::string uhat, const std::string s_y, double z)
void	init (double t)
void	run (double t, double dt, double alpha)
Public Member Functions inherited from InterMeshTransfer
	InterMeshTransfer (FluidAgent *fag, SolidAgent *sag, char *name=NULL)
Public Member Functions inherited from Action
	Action (void *p=0, char *name=NULL)
	Action (int n, const char *at[], int *i=NULL, void *p=0, char *name=NULL)
	Action (int n, const std::string at[], int *i=NULL, void *p=0, char *name=NULL)
virtual	~Action ()
virtual void	declare (Scheduler &)
virtual void	finalize ()
virtual char *	name ()
void	set_name (const char *name)
virtual void	print (FILE *f)
virtual void	print_toposort (FILE *f)
virtual void	schedule ()

Private Attributes
std::string	s_x_str
std::string	s_uhat_str
std::string	s_y_str
int	s_x_hdl
int	s_uhat_hdl
int	s_y_hdl
int	s_rb_hdl
int	s_mdot_hdl
int	s_areas_hdl
int	s_rhos_hdl
int	withALE
double	zoom
int	SURF_compute_bounded_volumes
int	SURF_compute_face_areas

Additional Inherited Members
Protected Types inherited from Action
enum	{ IN =1, OUT =2, INOUT =3 }
Protected Member Functions inherited from InterMeshTransfer
void	load_rocface (const RocmanControl_parameters *param)
Protected Member Functions inherited from Action
int	get_attribute_handle (int i)
int	get_attribute_handle_const (int i)
int	get_attribute_handle (const std::string str)
void	set_attr (int n, const std::string at[], int *id=NULL)
void	set_attr (int n, const char *at[], int *id=NULL)
void	set_io (int n, const int *io)
void	set_io (const char *io)
int	get_io (int i)
Protected Attributes inherited from InterMeshTransfer
FluidAgent *	fagent
SolidAgent *	sagent
Protected Attributes inherited from Action
char *	action_name
char **	attr
int *	idx
int	count
void *	usr_ptr
std::vector< int >	inout

Detailed Description

Definition at line 107 of file transfer_actions.h.

Constructor & Destructor Documentation

GetDeformedMesh_ALE ( FluidAgent *  fag, SolidAgent *  sag, const std::string  s_x, const std::string  uhat, const std::string  s_y, double  z  )

explicit

Definition at line 521 of file transfer_actions.C.

References Action::IN, Action::OUT, Action::set_attr(), and Action::set_io().

                                                                                                                                                    :
  InterMeshTransfer(fag, sag, (char *)"GetDeformedMesh_ALE"),
    s_x_str(s_x), s_uhat_str(s_uhat), s_y_str(s_y), zoom(z)
{
  int io[] = {IN, IN, OUT};
  set_io( 3, io); 
  
  std::string atts[4];
  atts[0] = s_x;
  atts[1] = s_uhat;
  atts[2] = s_y;
  set_attr(3, atts);

/*
    // register attributes
  std::string::size_type pos = s_x_str.find( ".");
  COM_assertion_msg(pos!=std::string::npos, "GetDeformedMesh_ALE::create_attribute failed!");
  std::string s = s_x_str.substr( 0, pos);
  std::string x = s_x_str.substr( pos, s_x_str.size());
  sagent->register_use_attribute( s, x, sagent->solidBufBase, ".nc");

    // s_y = solidBuf+".nc"
  pos = s_y_str.find( ".");
  COM_assertion_msg(pos!=std::string::npos, "GetDeformedMesh_ALE::create_attribute failed!");
  s = s_y_str.substr( 0, pos);
  std::string y = s_y_str.substr( pos, s_y_str.size());
  sagent->register_clone_attribute( 0, s, y, sagent->get_surface_window(), y);
*/
}

Here is the call graph for this function:

Member Function Documentation

void init ( double  t )

virtual

Reimplemented from Action.

Definition at line 551 of file transfer_actions.C.

References RocBlas::add, COM_call_function(), COM_get_attribute_handle(), COM_get_function_handle(), Action::get_attribute_handle(), load_rocsurf(), s_areas_hdl, s_mdot_hdl, s_rb_hdl, s_rhos_hdl, s_uhat_hdl, s_x_hdl, s_y_hdl, InterMeshTransfer::sagent, SolidAgent::solidBuf, SURF_compute_bounded_volumes, SURF_compute_face_areas, SolidAgent::withALE, and withALE.

                                        {
  s_x_hdl = get_attribute_handle(0);
  s_uhat_hdl = get_attribute_handle(1);
  s_y_hdl = get_attribute_handle(2);

  withALE = sagent->withALE;

  s_rb_hdl = COM_get_attribute_handle( sagent->solidBuf+".rb");
  s_mdot_hdl = COM_get_attribute_handle( sagent->solidBuf+".mdot");
  s_areas_hdl = COM_get_attribute_handle( sagent->solidBuf+".areas");
  s_rhos_hdl = COM_get_attribute_handle( sagent->solidBuf+".rhos");

  load_rocsurf();
  SURF_compute_bounded_volumes = COM_get_function_handle( "SURF.compute_bounded_volumes");
  SURF_compute_face_areas = COM_get_function_handle( "SURF.compute_element_areas");

#if 0
  // POST_INIT_SOLID
  COM_call_function( RocBlas::add, &s_x_hdl, &s_uhat_hdl, &s_y_hdl);
    // initial_start  ?????
  if (t == 0.0) {
//    double alpha = 0.0;
//     ????????????????? strange bug called in bcinitscheduler::init()
//    sagent->obtain_bc(&alpha);
  }
#endif
}

Here is the call graph for this function:

void run ( double  t, double  dt, double  alpha  )

virtual

Reimplemented from Action.

Definition at line 580 of file transfer_actions.C.

References RocBlas::add, COM_call_function(), RocBlas::copy, RocBlas::div, RocBlas::div_scalar, InterMeshTransfer::fagent, Agent::get_coupling(), Coupling::get_rocmancontrol_param(), MAN_DEBUG, RocBlas::maxof_scalar, RocBlas::min_scalar_MPI, MPI_COMM_SELF, MPI_COMM_WORLD, RocBlas::mul, RocmanControl_parameters::PROP_fom, s_areas_hdl, s_mdot_hdl, s_rb_hdl, s_rhos_hdl, s_uhat_hdl, s_x_hdl, s_y_hdl, SURF_compute_bounded_volumes, SURF_compute_face_areas, withALE, and zoom.

                                                                {
  MAN_DEBUG(3, ("Rocstar: calling GetDeformedMesh_ALE::run() with t:%e dt:%e alpha:%e withALE=%d zoom=%e.\n", t, dt, alpha, withALE, zoom));

  double mdot_min;

  if ( withALE && zoom>0 && !fagent->get_coupling()->get_rocmancontrol_param()->PROP_fom) {
      // compute volumes only for faces with nonzero burning rate
    COM_call_function( RocBlas::copy, &s_rb_hdl, &s_mdot_hdl);

    //debug_print(sagent->solidBuf+".nc", 202, 1);
    //debug_print(sagent->solidBuf+".x", 202, 1);

      // Compute mass injection using s_x and s_y as buffers.
      //  So far s_y stores the undeformed configuration s_x.
    int one = 1;
    COM_call_function( SURF_compute_bounded_volumes, &s_x_hdl, &s_y_hdl, &s_mdot_hdl, &one);

    //debug_print(sagent->solidBuf+".nc", 202, 1);
    //debug_print(attr[2], 202, 1);

#if 0
      //  Check whether the burning rates are nonnegative.
    int comm = MPI_COMM_SELF;
    COM_call_function( RocBlas::min_scalar_MPI, &s_mdot_hdl, &mdot_min, &comm);

    if ( mdot_min < 0) {
          printf("Rocstar ERROR: Negative mdot found %e. Aborting...\n", mdot_min);
          MPI_Abort( MPI_COMM_WORLD, -1);
    }
#else
    double zero = 0.0;
    COM_call_function( RocBlas::maxof_scalar, &s_mdot_hdl, &zero, &s_mdot_hdl);
#endif

    // Compute deformed configuration
    COM_call_function( RocBlas::add, &s_x_hdl, &s_uhat_hdl, &s_y_hdl);
    COM_call_function( SURF_compute_face_areas, &s_areas_hdl);

    // Compute mdot=dV*rhos/area/dt
    COM_call_function( RocBlas::mul, &s_mdot_hdl, &s_rhos_hdl, &s_mdot_hdl);
    COM_call_function( RocBlas::div, &s_mdot_hdl, &s_areas_hdl, &s_mdot_hdl);
    COM_call_function( RocBlas::div_scalar, &s_mdot_hdl, &dt, &s_mdot_hdl);
  }
  else {
      // Compute deformed configuration
    COM_call_function( RocBlas::add, &s_x_hdl, &s_uhat_hdl, &s_y_hdl);
      // Compute mdot as rb*rhos
    if ( withALE)
        COM_call_function( RocBlas::mul, &s_rb_hdl, &s_rhos_hdl, &s_mdot_hdl);
  }
}

Here is the call graph for this function:

Member Data Documentation

int s_areas_hdl

private

Definition at line 115 of file transfer_actions.h.

Referenced by init(), and run().

int s_mdot_hdl

private

Definition at line 115 of file transfer_actions.h.

Referenced by init(), and run().

int s_rb_hdl

private

Definition at line 115 of file transfer_actions.h.

Referenced by init(), and run().

int s_rhos_hdl

private

Definition at line 115 of file transfer_actions.h.

Referenced by init(), and run().

int s_uhat_hdl

private

Definition at line 114 of file transfer_actions.h.

Referenced by init(), and run().

std::string s_uhat_str

private

Definition at line 113 of file transfer_actions.h.

int s_x_hdl

private

Definition at line 114 of file transfer_actions.h.

Referenced by init(), and run().

std::string s_x_str

private

Definition at line 113 of file transfer_actions.h.

int s_y_hdl

private

Definition at line 114 of file transfer_actions.h.

Referenced by init(), and run().

std::string s_y_str

private

Definition at line 113 of file transfer_actions.h.

int SURF_compute_bounded_volumes

private

Definition at line 118 of file transfer_actions.h.

Referenced by init(), and run().

int SURF_compute_face_areas

private

Definition at line 118 of file transfer_actions.h.

Referenced by init(), and run().

int withALE

private

Definition at line 116 of file transfer_actions.h.

Referenced by init(), and run().

double zoom

private

Definition at line 117 of file transfer_actions.h.

Referenced by run().

---

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

• transfer_actions.h
• transfer_actions.C