DynamicPipe Class Reference

#include <pipe.h>

Inheritance diagram for DynamicPipe:

Collaboration diagram for DynamicPipe:

Public Member Functions
	DynamicPipe (unsigned int uL, const DofOwner *pD, HydraulicFluid *hf, const PressureNode *p1, const PressureNode *p2, doublereal Dh, doublereal A, doublereal L, flag transition, doublereal q0, flag fOut)
	~DynamicPipe (void)
virtual HydraulicElem::Type	GetHydraulicType (void) const
virtual std::ostream &	Restart (std::ostream &out) const
virtual unsigned int	iGetNumDof (void) const
virtual DofOrder::Order	GetDofType (unsigned int i) const
virtual DofOrder::Order	GetEqType (unsigned int i) const
virtual void	WorkSpaceDim (integer *piNumRows, integer *piNumCols) const
VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP)
virtual void	Output (OutputHandler &OH) const
virtual void	SetValue (DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *ph=0)
virtual void	GetConnectedNodes (std::vector< const Node * > &connectedNodes) const
Public Member Functions inherited from Elem
	Elem (unsigned int uL, flag fOut)
virtual	~Elem (void)
virtual std::ostream &	DescribeDof (std::ostream &out, const char *prefix="", bool bInitial=false) const
virtual void	DescribeDof (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const
virtual std::ostream &	DescribeEq (std::ostream &out, const char *prefix="", bool bInitial=false) const
virtual void	DescribeEq (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const
virtual void	AssMats (VariableSubMatrixHandler &WorkMatA, VariableSubMatrixHandler &WorkMatB, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual bool	bInverseDynamics (void) const
void	SetInverseDynamicsFlags (unsigned uIDF)
unsigned	GetInverseDynamicsFlags (void) const
bool	bIsErgonomy (void) const
bool	bIsRightHandSide (void) const
virtual VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)
virtual SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, const VectorHandler &XPrimePrimeCurr, InverseDynamics::Order iOrder=InverseDynamics::INVERSE_DYNAMICS)
virtual int	GetNumConnectedNodes (void) const
Public Member Functions inherited from WithLabel
	WithLabel (unsigned int uL=0, const std::string &sN="")
virtual	~WithLabel (void)
void	PutLabel (unsigned int uL)
void	PutName (const std::string &sN)
unsigned int	GetLabel (void) const
const std::string &	GetName (void) const
Public Member Functions inherited from SimulationEntity
	SimulationEntity (void)
virtual	~SimulationEntity (void)
virtual bool	bIsValidIndex (unsigned int i) const
virtual Hint *	ParseHint (DataManager *pDM, const char *s) const
virtual void	BeforePredict (VectorHandler &, VectorHandler &, VectorHandler &, VectorHandler &) const
virtual void	AfterPredict (VectorHandler &X, VectorHandler &XP)
virtual void	Update (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	DerivativesUpdate (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	Update (const VectorHandler &XCurr, InverseDynamics::Order iOrder)
virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP, const VectorHandler &XPP)
virtual unsigned int	iGetNumPrivData (void) const
virtual unsigned int	iGetPrivDataIdx (const char *s) const
virtual doublereal	dGetPrivData (unsigned int i) const
virtual std::ostream &	OutputAppend (std::ostream &out) const
virtual void	ReadInitialState (MBDynParser &HP)
Public Member Functions inherited from ToBeOutput
	ToBeOutput (flag fOut=fDefaultOut)
virtual	~ToBeOutput (void)
virtual void	OutputPrepare (OutputHandler &OH)
virtual void	Output (OutputHandler &OH, const VectorHandler &X, const VectorHandler &XP) const
virtual flag	fToBeOutput (void) const
virtual bool	bToBeOutput (void) const
virtual void	SetOutputFlag (flag f=flag(1))
Public Member Functions inherited from HydraulicElem
	HydraulicElem (unsigned int uL, const DofOwner *pDO, HydraulicFluid *hf, flag fOut)
virtual	~HydraulicElem (void)
virtual Elem::Type	GetElemType (void) const
Public Member Functions inherited from ElemWithDofs
	ElemWithDofs (unsigned int uL, const DofOwner *pDO, flag fOut)
virtual	~ElemWithDofs (void)
Public Member Functions inherited from DofOwnerOwner
	DofOwnerOwner (const DofOwner *pDO)
virtual	~DofOwnerOwner ()
virtual const DofOwner *	pGetDofOwner (void) const
virtual integer	iGetFirstIndex (void) const
virtual void	SetInitialValue (VectorHandler &X)

Private Attributes
const PressureNode *	pNode1
const PressureNode *	pNode2
doublereal	diameter
doublereal	area
doublereal	length
flag	turbulent
doublereal	q0
doublereal	Re
doublereal	p1
doublereal	p2
doublereal	p1p
doublereal	p2p
doublereal	q1
doublereal	q2
doublereal	q1p
doublereal	q2p
doublereal	density1
doublereal	density0
doublereal	density2
doublereal	densityDPres1
doublereal	densityDPres2
doublereal	viscosity
doublereal	dKlam
doublereal	dKtrb
doublereal	dKtra

Additional Inherited Members
Public Types inherited from Elem
enum	Type {   UNKNOWN = -1, AIRPROPERTIES = 0, INDUCEDVELOCITY, AUTOMATICSTRUCTURAL,   GRAVITY, BODY, JOINT, JOINT_REGULARIZATION,   BEAM, PLATE, FORCE, INERTIA,   ELECTRICBULK, ELECTRIC, THERMAL, HYDRAULIC,   BULK, LOADABLE, DRIVEN, EXTERNAL,   AEROMODAL, AERODYNAMIC, GENEL, SOCKETSTREAM_OUTPUT,   RTAI_OUTPUT = SOCKETSTREAM_OUTPUT, LASTELEMTYPE }
Public Types inherited from SimulationEntity
typedef std::vector< Hint * >	Hints
Public Types inherited from ToBeOutput
enum	{ OUTPUT = 0x1U, OUTPUT_MASK = 0xFU, OUTPUT_PRIVATE = 0x10U, OUTPUT_PRIVATE_MASK = ~OUTPUT_MASK }
Public Types inherited from HydraulicElem
enum	Type {   UNKNOWN = -1, MINOR_LOSS = 0, THREEWAYMINORLOSS, CONTROL_VALVE,   DYNAMIC_CONTROL_VALVE, PRESSURE_FLOW_CONTROL_VALVE, PRESSURE_VALVE, FLOW_VALVE,   ORIFICE, ACCUMULATOR, TANK, PIPE,   DYNAMIC_PIPE, HYDRAULIC_ACTUATOR, ACTUATOR, LASTHYDRAULICTYPE }
Protected Attributes inherited from WithLabel
unsigned int	uLabel
std::string	sName
Protected Attributes inherited from ToBeOutput
flag	fOutput
Protected Attributes inherited from HydraulicElem
HydraulicFluid *	HF

Detailed Description

Definition at line 201 of file pipe.h.

Constructor & Destructor Documentation

DynamicPipe::DynamicPipe	(	unsigned int	uL,
		const DofOwner *	pD,
		HydraulicFluid *	hf,
		const PressureNode *	p1,
		const PressureNode *	p2,
		doublereal	Dh,
		doublereal	A,
		doublereal	L,
		flag	transition,
		doublereal	q0,
		flag	fOut
	)

Definition at line 1006 of file pipe.cc.

References area, ASSERT, diameter, dKlam, dKtra, dKtrb, PressureNode::GetNodeType(), Node::HYDRAULIC, length, pNode1, pNode2, and grad::pow().

: Elem(uL, fOut), 
HydraulicElem(uL, pDO, hf, fOut),
pNode1(p1), 
pNode2(p2),
diameter(Dh), 
area(A),
length(L), 
turbulent(transition), 
q0(q0)
{
   ASSERT(pNode1 != NULL);
   ASSERT(pNode1->GetNodeType() == Node::HYDRAULIC);
   ASSERT(pNode2 != NULL);
   ASSERT(pNode2->GetNodeType() == Node::HYDRAULIC);
   
   ASSERT(Dh > std::numeric_limits<doublereal>::epsilon());
   ASSERT(A > std::numeric_limits<doublereal>::epsilon());
   ASSERT(L > std::numeric_limits<doublereal>::epsilon());
   
   dKlam = 16.*length/(diameter*diameter);
   dKtra = .25*length/(diameter*area);
   dKtrb = .5*length*.1582/(pow(diameter, 1.25)*pow(area, .75));
}

Here is the call graph for this function:

DynamicPipe::~DynamicPipe

(

void

)

Definition at line 1041 of file pipe.cc.

References NO_OP.

{
   NO_OP;
}

Member Function Documentation

void DynamicPipe::AfterConvergence ( const VectorHandler &  X, const VectorHandler &  XP  )

virtual

Reimplemented from SimulationEntity.

Definition at line 1353 of file pipe.cc.

References DEBUGCOUT, HydraulicFluid::dGetRe(), WithLabel::GetLabel(), HydraulicElem::HF, HydraulicFluid::LOWER, Re, turbulent, and HydraulicFluid::UPPER.

{
   if (Re < HF->dGetRe(HydraulicFluid::LOWER)) {
      turbulent = 0;
   } else if (Re > HF->dGetRe(HydraulicFluid::UPPER)) {
      turbulent = 1;
   }
   
#ifdef HYDR_DEVEL
   DEBUGCOUT("DynamicPipe(" << GetLabel() << "): turbulent mode =  " 
                   << turbulent << std::endl);
#endif /* HYDR_DEVEL */
}

Here is the call graph for this function:

VariableSubMatrixHandler & DynamicPipe::AssJac ( VariableSubMatrixHandler &  WorkMat, doublereal  dCoef, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr  )

virtual

Implements Elem.

Definition at line 1084 of file pipe.cc.

References a, area, c, density0, density1, density2, densityDPres1, densityDPres2, HydraulicFluid::dGetDensityDPres(), HydraulicFluid::dGetRe(), diameter, dKlam, dKtrb, grad::fabs(), HydraulicElem::HF, Node::iGetFirstColIndex(), DofOwnerOwner::iGetFirstIndex(), Node::iGetFirstRowIndex(), length, HydraulicFluid::LOWER, p1, p2, pNode1, pNode2, grad::pow(), FullSubMatrixHandler::PutCoef(), FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), q1, q2, Re, FullSubMatrixHandler::ResizeReset(), VariableSubMatrixHandler::SetFull(), turbulent, HydraulicFluid::UPPER, and viscosity.

{
   FullSubMatrixHandler& WM = WorkMat.SetFull();
   WM.ResizeReset(6, 6);
  
   integer iNode1RowIndex = pNode1->iGetFirstRowIndex()+1;
   integer iNode2RowIndex = pNode2->iGetFirstRowIndex()+1;
   integer iNode1ColIndex = pNode1->iGetFirstColIndex()+1;
   integer iNode2ColIndex = pNode2->iGetFirstColIndex()+1;
   
   integer iFirstIndex = iGetFirstIndex();  
   
   WM.PutRowIndex(1, iNode1RowIndex);
   WM.PutRowIndex(2, iNode2RowIndex);
   WM.PutRowIndex(3, iFirstIndex+1);
   WM.PutRowIndex(4, iFirstIndex+2);
   WM.PutRowIndex(5, iFirstIndex+3);
   WM.PutRowIndex(6, iFirstIndex+4);
   
   WM.PutColIndex(1, iNode1ColIndex);
   WM.PutColIndex(2, iNode2ColIndex);
   WM.PutColIndex(3, iFirstIndex+1);
   WM.PutColIndex(4, iFirstIndex+2);
   WM.PutColIndex(5, iFirstIndex+3);
   WM.PutColIndex(6, iFirstIndex+4);

   doublereal dRDP1 = HF->dGetDensityDPres(p1);
   doublereal dRDP2 = HF->dGetDensityDPres(p2);
   doublereal dRDP0 = HF->dGetDensityDPres(.5*(p1+p2));

   doublereal qq1 = .75*q1+.25*q2;
   doublereal qq2 = .25*q1+.75*q2;
   
   doublereal dd1 = .5*(density1+density0);
   doublereal dd2 = .5*(density0+density2);      
   
   doublereal dLoss11 = 0.;
   doublereal dLoss12 = 0.;
   doublereal dLoss21 = 0.;
   doublereal dLoss22 = 0.;
   
   if ((Re < HF->dGetRe(HydraulicFluid::LOWER)) || ((turbulent == 0) 
                           && (Re < 1.25*HF->dGetRe(HydraulicFluid::LOWER)))) {
      /* laminar, or ascending transition ascendente (continuation) */
      doublereal dk = dKlam*viscosity*dCoef;

      dLoss11 = .75*dk/dd1;
      dLoss12 = .25*dk/dd1;
      dLoss21 = .25*dk/dd2;
      dLoss22 = .75*dk/dd2;
   } else if ((Re > HF->dGetRe(HydraulicFluid::UPPER)) || ((turbulent == 1) 
                           && (Re > .775*HF->dGetRe(HydraulicFluid::UPPER)))) {
      /* turbulent, or descending transition (continuation) */
      
      doublereal dk = dKtrb*1.75*pow(viscosity, .25)*dCoef;
      doublereal dl1 = dk*pow(fabs(qq1), .75)/dd1;
      doublereal dl2 = dk*pow(fabs(qq2), .75)/dd2;
      
      dLoss11 = .75*dl1;
      dLoss12 = .25*dl1;
      dLoss21 = .25*dl2;
      dLoss22 = .75*dl2;
   } else {
      /* transition */
      doublereal dva = diameter/(viscosity*area);
      doublereal dva2 = dva*dva;
      doublereal dva3 = dva2*dva;
      
      doublereal a, b, c, d;
      
      if (turbulent == 0) {
         /* ascending */
         a = 7.e-13*dva3;
         b = 8.e-10*dva2;
         c = -1.8e-5*dva;
         d = .0542;
      } else /* if (turbulent == 1) */ {
         /* descending */
         a = 9.e-13*dva3;
         b = 2.e-9*dva2;
         c = -1.8e-5*dva;
         d = .0528;
      }
         
      doublereal fqq1 = fabs(qq1);
      doublereal fqq2 = fabs(qq2);
      
      doublereal fax1 = (((5.*a*fqq1+4.*b)*fqq1+3.*c)*fqq1+2.*d)*fqq1;
      doublereal fax2 = (((5.*a*fqq2+4.*b)*fqq2+3.*c)*fqq2+2.*d)*fqq2;
      
      doublereal dk = dKtrb*dCoef;
      
      doublereal dl1 = dk*fax1/dd1;
      doublereal dl2 = dk*fax2/dd2;
      
      dLoss11 = .75*dl1;
      dLoss12 = .25*dl1;
      dLoss21 = .25*dl2;
      dLoss22 = .75*dl2;
   }
   
   /* primo blocco: conservazione massa */
   doublereal dr = .5*dCoef;
   WM.PutCoef(1, 5, -dr);
   WM.PutCoef(1, 6, -dr);
   WM.PutCoef(2, 5, dr);
   WM.PutCoef(2, 6, dr);
 
   dr = area*length/8.;
   WM.PutCoef(1, 3, -dr*3.*densityDPres1);
   WM.PutCoef(1, 4, -dr*densityDPres1);
   WM.PutCoef(2, 3, -dr*densityDPres2);
   WM.PutCoef(2, 4, -dr*3.*densityDPres2);
   
   
   /* secondo blocco: equazione quantita' di moto */
   dr = .5*dCoef*area;
   doublereal q12 = .5*(q1+q2);
   doublereal ddq12 = .5*(q12*q12)/(density0*density0*area)*dRDP0*dCoef;
   doublereal ddq1 = q1*q1/(density1*density1*area)*dRDP1*dCoef;
   doublereal ddq2 = q2*q2/(density2*density2*area)*dRDP2*dCoef;
   WM.PutCoef(3, 3, -dr + ddq1 - ddq12);
   WM.PutCoef(3, 4, dr - ddq12);
   WM.PutCoef(4, 3, -dr + ddq12);
   WM.PutCoef(4, 4, dr + ddq12 - ddq2);
   
   /* manca la viscosita' */
   dr = length/8.;
   doublereal dq12 = .5*(q1+q2)/density0;   
   WM.PutCoef(3, 5, 3.*dr+(dq12-2.*q1/density1)/area*dCoef  + dLoss11);
   WM.PutCoef(3, 6, dr+dq12/area*dCoef + dLoss12);
   WM.PutCoef(4, 5, dr-dq12/area*dCoef + dLoss21);
   WM.PutCoef(4, 6, 3.*dr+(2.*q2/density2-dq12)/area*dCoef + dLoss22);
   
   /* terzo blocco: definizione delle pressioni nodali interne */   
   WM.PutCoef(5, 1, -1.);
   WM.PutCoef(6, 2, -1.);
   
   WM.PutCoef(5, 3, dCoef);
   WM.PutCoef(6, 4, dCoef);
   
   return WorkMat;
}

Here is the call graph for this function:

SubVectorHandler & DynamicPipe::AssRes ( SubVectorHandler &  WorkVec, doublereal  dCoef, const VectorHandler &  XCurr, const VectorHandler &  XPrimeCurr  )

virtual

Implements Elem.

Definition at line 1233 of file pipe.cc.

References a, area, c, DEBUGCOUT, density0, density1, density2, densityDPres1, densityDPres2, HydraulicFluid::dGetDensity(), HydraulicFluid::dGetDensityDPres(), HydraulicFluid::dGetRe(), HydraulicFluid::dGetViscosity(), ScalarAlgebraicNode::dGetX(), diameter, dKlam, dKtra, dKtrb, grad::fabs(), HydraulicElem::HF, DofOwnerOwner::iGetFirstIndex(), Node::iGetFirstRowIndex(), length, HydraulicFluid::LOWER, p0, p1, p1p, p2, p2p, pNode1, pNode2, grad::pow(), SubVectorHandler::PutItem(), q1, q1p, q2, q2p, Re, VectorHandler::Resize(), turbulent, HydraulicFluid::UPPER, and viscosity.

{
   DEBUGCOUT("Entering DynamicPipe::AssRes()" << std::endl);
   WorkVec.Resize(6);
   
   integer iNode1RowIndex = pNode1->iGetFirstRowIndex()+1;
   integer iNode2RowIndex = pNode2->iGetFirstRowIndex()+1;
   
   doublereal pn1 = pNode1->dGetX();
   doublereal pn2 = pNode2->dGetX();
   
   integer iFirstIndex = iGetFirstIndex();
   
   p1 = XCurr(iFirstIndex+1);        /* pressione */
   p2 = XCurr(iFirstIndex+2);        /* pressione */
   p1p = XPrimeCurr(iFirstIndex+1);  /* derivata pressione */
   p2p = XPrimeCurr(iFirstIndex+2);  /* derivata pressione */
   
   q1 = XCurr(iFirstIndex+3);        /* portata nodo 1 */
   q2 = XCurr(iFirstIndex+4);        /* portata nodo 2 */
   q1p = XPrimeCurr(iFirstIndex+3);  /* derivata portata nodo 1 */
   q2p = XPrimeCurr(iFirstIndex+4);  /* derivata portata nodo 2 */
   
   doublereal p0 = .5*(p1+p2);
   
   density1 = HF->dGetDensity(p1);            /* densita' all'inizio del tubo */
   density2 = HF->dGetDensity(p2);            /* densita' alla fine del nodo */
   density0 = HF->dGetDensity(p0);            /* densita' a meta' tubo */
      
   /* viscosity = HF->dGetViscosity(); */
   densityDPres1 = HF->dGetDensityDPres(.75*p1+.25*p2);
   densityDPres2 = HF->dGetDensityDPres(.25*p1+.75*p2);
   
   viscosity = HF->dGetViscosity(p0);
   
   doublereal qq1 = .75*q1+.25*q2;
   doublereal qq2 = .25*q1+.75*q2;
   
   doublereal dd1 = .5*(density1+density0);
   doublereal dd2 = .5*(density0+density2);
   
   doublereal q12 = .5*(q1+q2);

   /* determines whether the flow will be turbulent or laminar */
   Re = fabs(q12)/area*diameter/viscosity;

   doublereal dLoss1 = 0.;
   doublereal dLoss2 = 0.;

   if ((Re < HF->dGetRe(HydraulicFluid::LOWER)) || ((turbulent == 0) 
                           && (Re < 1.25*HF->dGetRe(HydraulicFluid::LOWER)))) {
      /* laminar, or ascending transition (continuation) */
      doublereal dk = dKlam*viscosity;
      
      dLoss1 = dk*qq1/dd1;
      dLoss2 = dk*qq2/dd2;
   } else if ((Re > HF->dGetRe(HydraulicFluid::UPPER)) || ((turbulent == 1) 
                           && (Re > .775*HF->dGetRe(HydraulicFluid::UPPER)))) {
      /* turbulent, or descending transition (continuation) */
      doublereal dk = dKtrb*pow(viscosity, .25);
      
      dLoss1 = dk*qq1*pow(fabs(qq1), .75)/dd1;
      dLoss2 = dk*qq2*pow(fabs(qq2), .75)/dd2;
   } else {
      /* transition */
      doublereal dva = diameter/(viscosity*area);
      doublereal dva2 = dva*dva;
      doublereal dva3 = dva2*dva;
      
      doublereal a, b, c, d;
      
      if (turbulent == 0) {
         /* ascending */
         a = 7.e-13*dva3;
         b = 8.e-10*dva2;
         c = -1.8e-5*dva;
         d = .0542;
      } else /* if (turbulent == 1) */ {
         /* descending */
         a = 9.e-13*dva3;
         b = 2.e-9*dva2;
         c = -1.8e-5*dva;
         d = .0528;
      }
         
      doublereal fqq1 = fabs(qq1);
      doublereal fqq2 = fabs(qq2);
      
      doublereal fax1 = qq1*(((a*fqq1+b)*fqq1+c)*fqq1+d)*fqq1;
      doublereal fax2 = qq2*(((a*fqq2+b)*fqq2+c)*fqq2+d)*fqq2;
      
      dLoss1 = dKtra*fax1/dd1;
      dLoss2 = dKtra*fax2/dd2;
   }

   /* mass conservation */
   doublereal dr = area*length/8.;
   WorkVec.PutItem(1, iNode1RowIndex , q12+dr*densityDPres1*(3.*p1p+p2p));
   WorkVec.PutItem(2, iNode2RowIndex , -q12+dr*densityDPres2*(p1p+3.*p2p));
   
   /* momentum balance */
   doublereal dp = .5*area*(p2-p1);
   doublereal dq12 = q12*q12/(density0*area);
   dr = length/8.;
   WorkVec.PutItem(3, iFirstIndex+1,
                    -dr*(3.*q1p+q2p)-dq12+q1*q1/(density1*area)-dp-dLoss1);
   WorkVec.PutItem(4, iFirstIndex+2,
                    -dr*(q1p+3.*q2p)-q2*q2/(density2*area)+dq12-dp-dLoss2);
   
   /* differential pressure definition */
   WorkVec.PutItem(5, iFirstIndex+3, pn1-p1);
   WorkVec.PutItem(6, iFirstIndex+4, pn2-p2);
      
   return WorkVec;
}

Here is the call graph for this function:

virtual void DynamicPipe::GetConnectedNodes ( std::vector< const Node * > &  connectedNodes ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 279 of file pipe.h.

References pNode1, and pNode2.

                                                                                 {
      connectedNodes.resize(2);
      connectedNodes[0] = pNode1;
      connectedNodes[1] = pNode2;
   };

DofOrder::Order DynamicPipe::GetDofType ( unsigned int  i ) const

virtual

Reimplemented from Elem.

Definition at line 1063 of file pipe.cc.

References ASSERT, and DofOrder::DIFFERENTIAL.

{
   ASSERT(i >= 0 && i < 4);
   return DofOrder::DIFFERENTIAL;
}

DofOrder::Order DynamicPipe::GetEqType ( unsigned int  i ) const

virtual

Reimplemented from SimulationEntity.

Definition at line 1070 of file pipe.cc.

References ASSERT, and DofOrder::DIFFERENTIAL.

{
   ASSERT(i >= 0 && i < 4);
   return DofOrder::DIFFERENTIAL;
}

HydraulicElem::Type DynamicPipe::GetHydraulicType ( void  ) const

virtual

Implements HydraulicElem.

Definition at line 1047 of file pipe.cc.

References HydraulicElem::DYNAMIC_PIPE.

                                                          {
   return HydraulicElem::DYNAMIC_PIPE;
}

unsigned int DynamicPipe::iGetNumDof ( void  ) const

virtual

Reimplemented from Elem.

Definition at line 1057 of file pipe.cc.

{
   return 4;
}

void DynamicPipe::Output ( OutputHandler &  OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 1367 of file pipe.cc.

References ToBeOutput::bToBeOutput(), density0, density1, density2, densityDPres1, densityDPres2, WithLabel::GetLabel(), OutputHandler::Hydraulic(), p1, p1p, p2, p2p, q1, q1p, q2, q2p, Re, and turbulent.

{
   if (bToBeOutput()) {
      std::ostream& out = OH.Hydraulic();
      out 
        << std::setw(8) << GetLabel()   /*  1 */
        << " " << p1                    /*  2 */
        << " " << p2                    /*  3 */
        << " " << p1p                   /*  4 */
        << " " << p2p                   /*  5 */
        << " " << q1                    /*  6 */
        << " " << q2                    /*  7 */
        << " " << q1p                   /*  8 */
        << " " << q2p                   /*  9 */
        << " " << density1              /* 10 */
        << " " << density0              /* 11 */
        << " " << density2              /* 12 */
        << " " << densityDPres1         /* 13 */
        << " " << densityDPres2         /* 14 */
        << " " << Re                    /* 15 */
        << " " << turbulent             /* 16 */
        << std::endl;
   }
}

Here is the call graph for this function:

std::ostream & DynamicPipe::Restart ( std::ostream &  out ) const

virtual

Implements Elem.

Definition at line 1052 of file pipe.cc.

{
   return out << "Pipe not implemented yet!" << std::endl;
}

void DynamicPipe::SetValue ( DataManager *  pDM, VectorHandler &  X, VectorHandler &  XP, SimulationEntity::Hints *  ph = 0  )

virtual

Reimplemented from SimulationEntity.

Definition at line 1393 of file pipe.cc.

References ScalarAlgebraicNode::dGetX(), DofOwnerOwner::iGetFirstIndex(), p1, p2, pNode1, pNode2, VectorHandler::PutCoef(), and q0.

{
   integer i = iGetFirstIndex();
   
   doublereal p1 = pNode1->dGetX();
   doublereal p2 = pNode2->dGetX();
   
   X.PutCoef(i+1, p1);
   X.PutCoef(i+2, p2);
   X.PutCoef(i+3, q0);
   X.PutCoef(i+4, -q0);
   
   XP.PutCoef(i+1, 0.);
   XP.PutCoef(i+2, 0.);
   XP.PutCoef(i+3, 0.);
   XP.PutCoef(i+4, 0.);
}

Here is the call graph for this function:

void DynamicPipe::WorkSpaceDim ( integer *  piNumRows, integer *  piNumCols  ) const

virtual

Implements Elem.

Definition at line 1077 of file pipe.cc.

{
   *piNumRows = 6;
   *piNumCols = 6; 
}

Member Data Documentation

doublereal DynamicPipe::area

private

Definition at line 207 of file pipe.h.

Referenced by AssJac(), AssRes(), and DynamicPipe().

doublereal DynamicPipe::density0

private

Definition at line 226 of file pipe.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal DynamicPipe::density1

private

Definition at line 225 of file pipe.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal DynamicPipe::density2

private

Definition at line 227 of file pipe.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal DynamicPipe::densityDPres1

private

Definition at line 228 of file pipe.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal DynamicPipe::densityDPres2

private

Definition at line 229 of file pipe.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal DynamicPipe::diameter

private

Definition at line 206 of file pipe.h.

Referenced by AssJac(), AssRes(), and DynamicPipe().

doublereal DynamicPipe::dKlam

private

Definition at line 233 of file pipe.h.

Referenced by AssJac(), AssRes(), and DynamicPipe().

doublereal DynamicPipe::dKtra

private

Definition at line 235 of file pipe.h.

Referenced by AssRes(), and DynamicPipe().

doublereal DynamicPipe::dKtrb

private

Definition at line 234 of file pipe.h.

Referenced by AssJac(), AssRes(), and DynamicPipe().

doublereal DynamicPipe::length

private

Definition at line 208 of file pipe.h.

Referenced by AssJac(), AssRes(), and DynamicPipe().

doublereal DynamicPipe::p1

private

Definition at line 215 of file pipe.h.

Referenced by AssJac(), AssRes(), Output(), and SetValue().

doublereal DynamicPipe::p1p

private

Definition at line 217 of file pipe.h.

Referenced by AssRes(), and Output().

doublereal DynamicPipe::p2

private

Definition at line 216 of file pipe.h.

Referenced by AssJac(), AssRes(), Output(), and SetValue().

doublereal DynamicPipe::p2p

private

Definition at line 218 of file pipe.h.

Referenced by AssRes(), and Output().

const PressureNode* DynamicPipe::pNode1

private

Definition at line 203 of file pipe.h.

Referenced by AssJac(), AssRes(), DynamicPipe(), GetConnectedNodes(), and SetValue().

const PressureNode* DynamicPipe::pNode2

private

Definition at line 204 of file pipe.h.

Referenced by AssJac(), AssRes(), DynamicPipe(), GetConnectedNodes(), and SetValue().

doublereal DynamicPipe::q0

private

Definition at line 211 of file pipe.h.

Referenced by SetValue().

doublereal DynamicPipe::q1

private

Definition at line 220 of file pipe.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal DynamicPipe::q1p

private

Definition at line 222 of file pipe.h.

Referenced by AssRes(), and Output().

doublereal DynamicPipe::q2

private

Definition at line 221 of file pipe.h.

Referenced by AssJac(), AssRes(), and Output().

doublereal DynamicPipe::q2p

private

Definition at line 223 of file pipe.h.

Referenced by AssRes(), and Output().

doublereal DynamicPipe::Re

private

Definition at line 213 of file pipe.h.

Referenced by AfterConvergence(), AssJac(), AssRes(), and Output().

flag DynamicPipe::turbulent

private

Definition at line 210 of file pipe.h.

Referenced by AfterConvergence(), AssJac(), AssRes(), and Output().

doublereal DynamicPipe::viscosity

private

Definition at line 231 of file pipe.h.

Referenced by AssJac(), and AssRes().

---

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

• pipe.h
• pipe.cc