FluidPropagateSurface Class Reference

#include <basic_actions.h>

Inheritance diagram for FluidPropagateSurface:

Collaboration diagram for FluidPropagateSurface:

Public Member Functions
	FluidPropagateSurface (FluidAgent *ag, BurnAgent *bag, const std::string b_rb, const std::string a_vm, double z)
void	init (double t)
void	run (double t, double dt, double alpha)
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
FluidAgent *	fagent
BurnAgent *	bagent
double	zoom
int	p_rb_hdl
int	b_rb_hdl
int	fb_rb_hdl
int	p_cnstr_type
int	p_pmesh_hdl
int	p_vm_hdl
int	a_vm_hdl
int	p_cflag_hdl
int	p_pos_hdl
int	p_bflag_hdl
int	MAP_reduce_maxabs
int	MAP_reduce_minabs
int	PROP_set_cnstr
int	PROP_propagate
int	PROPCON_find_intersections
int	PROPCON_constrain_displacements
int	PROPCON_burnout
int	PROPCON_burnout_filter

Additional Inherited Members
Protected Types inherited from Action
enum	{ IN =1, OUT =2, INOUT =3 }
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 Action
char *	action_name
char **	attr
int *	idx
int	count
void *	usr_ptr
std::vector< int >	inout

Detailed Description

Definition at line 171 of file basic_actions.h.

Constructor & Destructor Documentation

FluidPropagateSurface	(	FluidAgent *	ag,
		BurnAgent *	bag,
		const std::string	b_rb,
		const std::string	a_vm,
		double	z
	)

Definition at line 456 of file basic_actions.C.

References COM_DOUBLE, fagent, Action::IN, Action::OUT, FluidAgent::propBufAll, Agent::register_new_attribute(), Action::set_attr(), and Action::set_io().

                      : 
  Action(0, (const char**)NULL, NULL, NULL, (char *)"FluidPropagateSurface"), fagent(fag), bagent(bag), zoom(z)
{
  int io[] = {IN, OUT};
  set_io( 2, io); 
  
  std::string atts[2];
  atts[0] = b_rb;
  atts[1] = a_vm;
  set_attr(2, atts);

  fagent->register_new_attribute( fagent->propBufAll, ".vm", 'n', COM_DOUBLE, 3, "m/s");
}

Here is the call graph for this function:

Member Function Documentation

void init ( double  t )

virtual

Reimplemented from Action.

Definition at line 473 of file basic_actions.C.

References a_vm_hdl, b_rb_hdl, COM_get_attribute_handle(), COM_get_function_handle(), COM_get_window_handle(), fagent, fb_rb_hdl, FluidAgent::fluidBufB, Action::get_attribute_handle(), Agent::get_communicator(), Agent::get_coupling(), Coupling::get_rocmancontrol_param(), load_rocmap(), load_rocon(), load_rocprop(), MAP_reduce_maxabs, MAP_reduce_minabs, p_bflag_hdl, p_cflag_hdl, p_cnstr_type, p_pmesh_hdl, p_pos_hdl, p_rb_hdl, p_vm_hdl, RocmanControl_parameters::PROP_fom, PROP_propagate, PROP_set_cnstr, FluidAgent::propBuf, FluidAgent::propBufAll, PROPCON_burnout, PROPCON_burnout_filter, PROPCON_constrain_displacements, and PROPCON_find_intersections.

                                          {
  b_rb_hdl = get_attribute_handle( 0);
  a_vm_hdl = get_attribute_handle( 1);

  std::string propBuf = fagent->propBufAll;
  int PROP_fom = fagent->get_coupling()->get_rocmancontrol_param()->PROP_fom;
  if (!PROP_fom) propBuf = fagent->propBuf;
  p_rb_hdl     = COM_get_attribute_handle( propBuf+".rb");
  p_pmesh_hdl  = COM_get_attribute_handle( propBuf+".pmesh");
  p_vm_hdl     = COM_get_attribute_handle( propBuf+".vm");
  p_bflag_hdl  = COM_get_attribute_handle( propBuf+".bflag");
  p_cnstr_type = COM_get_attribute_handle( propBuf+".cnstr_type");
  p_cflag_hdl  = COM_get_attribute_handle( propBuf+".cflag");
  p_pos_hdl    = COM_get_attribute_handle( propBuf+".positions");
  fb_rb_hdl    = COM_get_attribute_handle( fagent->fluidBufB+".rb");

  load_rocprop(fagent->get_coupling()->get_rocmancontrol_param(), fagent->get_communicator());
  load_rocon(fagent->get_coupling()->get_rocmancontrol_param(),fagent->get_communicator());
  PROP_set_cnstr = COM_get_function_handle("PROP.set_constraints");
  PROP_propagate = COM_get_function_handle("PROP.propagate");
  if(COM_get_window_handle("PROPCON") >= 0) {
    PROPCON_find_intersections = COM_get_function_handle("PROPCON.find_intersections");
    PROPCON_constrain_displacements = COM_get_function_handle("PROPCON.constrain_displacements");
    PROPCON_burnout = COM_get_function_handle("PROPCON.burnout");
    PROPCON_burnout_filter = COM_get_function_handle("PROPCON.burnout_filter");
  }
  else {
    PROPCON_find_intersections = -1;
    PROPCON_constrain_displacements = -1;
    PROPCON_burnout = -1;
  }
  load_rocmap();
  MAP_reduce_maxabs = COM_get_function_handle("MAP.reduce_maxabs_on_shared_nodes");
  MAP_reduce_minabs = COM_get_function_handle("MAP.reduce_minabs_on_shared_nodes");
}

Here is the call graph for this function:

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

virtual

Reimplemented from Action.

Definition at line 509 of file basic_actions.C.

References a_vm_hdl, b_rb_hdl, COM_call_function(), RocBlas::copy, RocBlas::copy_scalar, RocBlas::div_scalar, fagent, fb_rb_hdl, Agent::get_communicator(), Agent::get_coupling(), Coupling::get_rocmancontrol_param(), MAN_DEBUG, man_verbose, MAP_reduce_maxabs, p_bflag_hdl, p_cflag_hdl, p_cnstr_type, p_pmesh_hdl, p_pos_hdl, p_rb_hdl, p_vm_hdl, RocmanControl_parameters::PROP_fom, PROP_propagate, PROP_set_cnstr, PROPCON_burnout, PROPCON_burnout_filter, PROPCON_constrain_displacements, PROPCON_find_intersections, rank, and zoom.

{
  if (zoom <= 0.0 ||  fagent->get_coupling()->initial_start() ) return;

  int rank = 0;
  MPI_Comm mycomm = fagent->get_communicator();
  MPI_Comm_rank(mycomm,&rank);
  if(rank == 0)
    MAN_DEBUG(2, ("Rocstar: FluidPropagateSurface::run() with t:%e dt:%e zoom:%f p_cnstr_type:%d.\n", t, dt, zoom, p_cnstr_type));

  double zero = 0.0;
  COM_call_function( RocBlas::copy_scalar, &zero, &p_rb_hdl);

  COM_call_function( RocBlas::copy, &b_rb_hdl, &fb_rb_hdl);

  // Burnout strategy:
  //
  // integer array (bflag) lives on interacting boundary 
  // faces.  bflag indicates whether the boundary face
  // is burning/reacting/injecting, etc.  In the current
  // implementation, bflag can be 1 or 0 for burning or
  // not-burning.
  //
  // integer array (cflag) lives on boundary nodes. 
  // cflag indicates the status of one or more geometrical
  // or other imposed constraints on the deformation or
  // solution space of the domain.
  // In the current implementation cflag = 1 for nodes
  // which have been constrained from further motion by
  // the physical boundary.  This happens when the domain
  // boundary and the physical boundary come in contact.
  //
  // Based on the value of bflag, certain models may
  // need to be used to determine physical quantities
  // on the boundary.   In the current implementation,
  // the burnrate solver (Rocburn) can determine whether
  // the propellant material should ignite, and if so
  // it sets bflag to 1 and sets a burning rate
  // for each burning face based on (aP^n), where
  // [P = fluid pressure], and the temperature of the
  // reacting surface.
  //
  //
  // This function sets bflag to 0 for every burned out
  // face.  Burned out faces are determined by faces
  // in which every node has cflag=1
  // BurnOut(cflag,bflag);
  if(PROPCON_burnout >= 0){
    //    std::cout << "Rocstar> Burning out surface elements with Rocon." 
    //        << std::endl;
    double data2 = 20000*dt;
    COM_call_function( PROP_propagate, &p_pmesh_hdl, &p_rb_hdl, &data2, &p_vm_hdl);
    MPI_Barrier(mycomm);
    COM_call_function(PROPCON_find_intersections,&p_pmesh_hdl,&p_vm_hdl,&p_pos_hdl,&p_cflag_hdl);
    COM_call_function( MAP_reduce_maxabs, &p_cflag_hdl);
    COM_call_function(PROPCON_burnout,&p_pmesh_hdl,&p_cflag_hdl,&p_bflag_hdl);
    COM_call_function( MAP_reduce_maxabs, &p_cflag_hdl); 
    COM_call_function(PROPCON_burnout,&p_pmesh_hdl,&p_cflag_hdl,&p_bflag_hdl);
    //    COM_call_function( MAP_reduce_maxabs,
  }

  //
  // This function sets the burning rate (rb) to 0.
  // The propagation code (Rocprop) sees rb as the 
  // (speed function).  If the speed function is 0
  // then Rocprop will not propagate the face.  This will
  // help Rocprop and Rocon to work together in keeping
  // the propagation constrained to the constraint surface.
  // p_rb_hdl = burning_rate
  if(PROPCON_burnout_filter >= 0)
    COM_call_function(PROPCON_burnout_filter,&p_bflag_hdl,&p_rb_hdl);
  //
  
  // Now original code to propagate the mesh:
  if ( p_cnstr_type > 0) 
    COM_call_function( PROP_set_cnstr, &p_cnstr_type);

  double data = dt * zoom;

  MPI_Barrier(mycomm);
  if(!rank && man_verbose > 2)
    std::cout << "Rocstar: Calling Rocprop"  << std::endl;
  // p_rb_hdl is the "speed function"
  COM_call_function( PROP_propagate, &p_pmesh_hdl, &p_rb_hdl, &data, &p_vm_hdl);
  MPI_Barrier(mycomm);
  if(!rank && man_verbose > 2)
    std::cout << "Rocstar: Rocprop done." << std::endl;
  // Rocon determines if the domain boundary and the physical boundary have
  // contacted and sets cflag (p_cflag_hdl) if so, and also constrains the 
  // nodes to not breach the physical boundary.
  if(PROPCON_constrain_displacements >= 0){
    MPI_Barrier(mycomm);
    if(!rank && man_verbose > 2)
      std::cout << "Rocstar: Finding constraint intersections." << std::endl;
    //    std::cout << "Constraining displacements with Rocon" << std::endl;
    //    COM_call_function(PROPCON_find_intersections,&p_pmesh_hdl,&p_vm_hdl,&p_pos_hdl,&p_cflag_hdl);
    MPI_Barrier(mycomm);
    if(!rank && man_verbose > 2)
      std::cout << "Rocstar: Reducing intersections." << std::endl;
    //    COM_call_function( MAP_reduce_maxabs, &p_cflag_hdl);
    //    COM_call_function( MAP_reduce_minabs, &p_cflag_hdl);
    MPI_Barrier(mycomm);
    if(!rank && man_verbose > 2)
      std::cout << "Rocstar: Constraining displacements." << std::endl;
    //    COM_call_function(PROPCON_constrain_displacements,&p_pmesh_hdl,&p_vm_hdl,&p_pos_hdl,&p_cflag_hdl);
    MPI_Barrier(mycomm);
    if(!rank && man_verbose > 2)
      std::cout << "Rocstar: Displacements constrained." << std::endl;
  }

  COM_call_function( RocBlas::div_scalar, &p_vm_hdl, &dt, &p_vm_hdl);

  const RocmanControl_parameters *param = fagent->get_coupling()->get_rocmancontrol_param();

  // Let's try this to help fix Rocflo symmetry plane motion
  COM_call_function( MAP_reduce_maxabs, &a_vm_hdl);
  //  COM_call_function(MAP_reduce_minabs,&a_vm_hdl);

  if ( ! param->PROP_fom) {
     // Move shared nodes between burning and nonburning panes.
     COM_call_function( MAP_reduce_maxabs, &a_vm_hdl);
  }

  //
  // Other functions that go elsewhere that are associated 
  // with burnout that I can think of are:
  // 
  // Sets mass flux to 0.0 on burned out faces
  // BurnOutFilter(bflag,0.0,mass_flux);
  // or BurnOutFilter(bflag,mass_flux_insulation,mass_flux);
  //
  // Sets Tflame to Tgas for a perfect insulator, or
  // something else if maybe the insulation chars or
  // something.
  // BurnOutFilter(bflag,gas_temp,tflame);
  // or maybe BurOutFilter(bflag,ins_flame_temp,tflame);
  //
  // 

//debug_print(fagent->propBufAll+".vm", 102, 0);
}

Here is the call graph for this function:

Member Data Documentation

int a_vm_hdl

private

Definition at line 180 of file basic_actions.h.

Referenced by init(), and run().

int b_rb_hdl

private

Definition at line 180 of file basic_actions.h.

Referenced by init(), and run().

BurnAgent* bagent

private

Definition at line 178 of file basic_actions.h.

FluidAgent* fagent

private

Definition at line 177 of file basic_actions.h.

Referenced by FluidPropagateSurface(), init(), and run().

int fb_rb_hdl

private

Definition at line 180 of file basic_actions.h.

Referenced by init(), and run().

int MAP_reduce_maxabs

private

Definition at line 182 of file basic_actions.h.

Referenced by init(), and run().

int MAP_reduce_minabs

private

Definition at line 183 of file basic_actions.h.

Referenced by init().

int p_bflag_hdl

private

Definition at line 181 of file basic_actions.h.

Referenced by init(), and run().

int p_cflag_hdl

private

Definition at line 181 of file basic_actions.h.

Referenced by init(), and run().

int p_cnstr_type

private

Definition at line 180 of file basic_actions.h.

Referenced by init(), and run().

int p_pmesh_hdl

private

Definition at line 180 of file basic_actions.h.

Referenced by init(), and run().

int p_pos_hdl

private

Definition at line 181 of file basic_actions.h.

Referenced by init(), and run().

int p_rb_hdl

private

Definition at line 180 of file basic_actions.h.

Referenced by init(), and run().

int p_vm_hdl

private

Definition at line 180 of file basic_actions.h.

Referenced by init(), and run().

int PROP_propagate

private

Definition at line 184 of file basic_actions.h.

Referenced by init(), and run().

int PROP_set_cnstr

private

Definition at line 184 of file basic_actions.h.

Referenced by init(), and run().

int PROPCON_burnout

private

Definition at line 186 of file basic_actions.h.

Referenced by init(), and run().

int PROPCON_burnout_filter

private

Definition at line 186 of file basic_actions.h.

Referenced by init(), and run().

int PROPCON_constrain_displacements

private

Definition at line 185 of file basic_actions.h.

Referenced by init(), and run().

int PROPCON_find_intersections

private

Definition at line 185 of file basic_actions.h.

Referenced by init(), and run().

double zoom

private

Definition at line 179 of file basic_actions.h.

Referenced by run().

---

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

• basic_actions.h
• basic_actions.C