Orifice Class Reference

#include <hminor.h>

Inheritance diagram for Orifice:

Collaboration diagram for Orifice:

Public Member Functions
	Orifice (unsigned int uL, const DofOwner *pD, HydraulicFluid *hf, const PressureNode *p1, const PressureNode *p2, doublereal Dh, doublereal A_diaf, doublereal A_pipe, doublereal ReCR, flag fOut)
	~Orifice (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 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	Output (OutputHandler &OH) const
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 DofOrder::Order	GetEqType (unsigned int i) const
virtual void	SetValue (DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *h=0)
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)
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 *	m_pNode1
const PressureNode *	m_pNode2
doublereal	diameter
doublereal	viscosity
doublereal	m_Area_diaf
doublereal	m_Area_pipe
doublereal	ReCr
doublereal	CriticJump
doublereal	delta
doublereal	Cd
doublereal	flow
doublereal	vel
doublereal	Re
bool	turbulent

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 170 of file hminor.h.

Constructor & Destructor Documentation

Orifice::Orifice	(	unsigned int	uL,
		const DofOwner *	pD,
		HydraulicFluid *	hf,
		const PressureNode *	p1,
		const PressureNode *	p2,
		doublereal	Dh,
		doublereal	A_diaf,
		doublereal	A_pipe,
		doublereal	ReCR,
		flag	fOut
	)

Definition at line 420 of file hminor.cc.

References ASSERT, Cd, CriticJump, delta, HydraulicFluid::dGetDensity(), HydraulicFluid::dGetViscosity(), ScalarAlgebraicNode::dGetX(), diameter, WithLabel::GetLabel(), PressureNode::GetNodeType(), HydraulicElem::HF, Node::HYDRAULIC, m_Area_diaf, m_Area_pipe, m_pNode1, m_pNode2, MBDYN_EXCEPT_ARGS, grad::pow(), grad::sqrt(), and viscosity.

: Elem(uL, fOut),
HydraulicElem(uL, pDO, hf, fOut),
m_pNode1(p1), m_pNode2(p2),
diameter(Dh), m_Area_diaf(A_diaf),
m_Area_pipe(A_pipe), ReCr(ReCr),
flow(0.),
vel(0.),
Re(0.),
turbulent(false)
{
        ASSERT(m_pNode1 != NULL);
        ASSERT(m_pNode1->GetNodeType() == Node::HYDRAULIC);
        ASSERT(m_pNode2 != NULL);
        ASSERT(m_pNode2->GetNodeType() == Node::HYDRAULIC);
        ASSERT(Dh > std::numeric_limits<doublereal>::epsilon());
        ASSERT(A_diaf > std::numeric_limits<doublereal>::epsilon());
        ASSERT(A_pipe > std::numeric_limits<doublereal>::epsilon());

        /* se |p1-p2| < CriticJump avrò sicuramente moto laminare se no turbolento */

        viscosity = HF->dGetViscosity();
        /*
         * Merritt, pp. 43-45:
         *
         * delta = 0.2 for a sharp-edged round orifice
         * delta = 0.157 for a sharp-edged slit orifice
         */
        delta = 0.2;
        // ReCr = pow(0.611/delta, 2);
        doublereal density = HF->dGetDensity((m_pNode2->dGetX() + m_pNode1->dGetX())/2.);
        CriticJump = m_Area_pipe*ReCr*pow(viscosity/(diameter*delta), 2.)/(2.*density*m_Area_diaf);

        /* calcolo del Cd */

        /* Cc = funzione delle due aree; */
        doublereal rapp = m_Area_diaf/m_Area_pipe;
        /* parametro adimensionale in funzione del rapporto m_Area_diaf/m_Area_pipe */
        doublereal Cc = (((.4855*rapp - .4971)*rapp + .158)*rapp + .1707)*rapp + .6005;

        doublereal base = Cc*rapp;
        doublereal rad = 1. - base*base;
        if (rad < 1.e3*std::numeric_limits<doublereal>::epsilon()) {
                silent_cerr("Orifice(" << GetLabel() << "): error computing Cd" << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }

        doublereal Cv = .98; /* costante adimensionale */
        Cd = Cv*Cc/sqrt(rad);
}

Here is the call graph for this function:

Orifice::~Orifice

(

void

)

Definition at line 475 of file hminor.cc.

References NO_OP.

{
        NO_OP;
}

Member Function Documentation

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

virtual

Implements Elem.

Definition at line 515 of file hminor.cc.

References Cd, CriticJump, DEBUGCOUT, delta, HydraulicFluid::dGetDensity(), ScalarAlgebraicNode::dGetX(), diameter, grad::fabs(), HydraulicElem::HF, Node::iGetFirstColIndex(), Node::iGetFirstRowIndex(), m_Area_diaf, m_pNode1, m_pNode2, FullSubMatrixHandler::PutCoef(), FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::Resize(), VariableSubMatrixHandler::SetFull(), grad::sqrt(), and viscosity.

{
        DEBUGCOUT("Entering Orifice::AssJac()" << std::endl);

        FullSubMatrixHandler& WM = WorkMat.SetFull();
        WM.Resize(2, 2);

        integer iNode1RowIndex = m_pNode1->iGetFirstRowIndex() + 1;
        integer iNode2RowIndex = m_pNode2->iGetFirstRowIndex() + 1;
        integer iNode1ColIndex = m_pNode1->iGetFirstColIndex() + 1;
        integer iNode2ColIndex = m_pNode2->iGetFirstColIndex() + 1;

        WM.PutRowIndex(1, iNode1RowIndex);
        WM.PutRowIndex(2, iNode2RowIndex);
        WM.PutColIndex(1, iNode1ColIndex);
        WM.PutColIndex(2, iNode2ColIndex);

        doublereal p1 = m_pNode1->dGetX();
        doublereal p2 = m_pNode2->dGetX();
        doublereal jumpPres = fabs(p1-p2);

        doublereal Jac = 0.;
        doublereal density = HF->dGetDensity((p1 + p2)/2.);

        if (jumpPres < CriticJump) {
                /*  moto sicuramente laminare (jacobiano) */
                Jac = -density*2.*(delta*delta)*diameter*m_Area_diaf/viscosity;

        } else {
                /*  moto sicuramente turbolento  (jacobiano) */
                Jac = -Cd*m_Area_diaf/sqrt(2.*jumpPres/density);
        }

        /* Jac *= dCoef; */

        WM.PutCoef(1, 1, Jac);
        WM.PutCoef(1, 2, -Jac);
        WM.PutCoef(2, 1, -Jac);
        WM.PutCoef(2, 2, Jac);

        return WorkMat;
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 562 of file hminor.cc.

References Cd, copysign(), CriticJump, DEBUGCOUT, delta, HydraulicFluid::dGetDensity(), ScalarAlgebraicNode::dGetX(), diameter, grad::fabs(), flow, HydraulicElem::HF, Node::iGetFirstRowIndex(), m_Area_diaf, m_Area_pipe, m_pNode1, m_pNode2, SubVectorHandler::PutItem(), Re, VectorHandler::Resize(), grad::sqrt(), turbulent, vel, and viscosity.

{
        DEBUGCOUT("Entering Orifice::AssRes()" << std::endl);

        WorkVec.Resize(2);

        integer iNode1RowIndex = m_pNode1->iGetFirstRowIndex() + 1;
        integer iNode2RowIndex = m_pNode2->iGetFirstRowIndex() + 1;

        doublereal p1 = m_pNode1->dGetX();
        doublereal p2 = m_pNode2->dGetX();

        doublereal jumpPres = fabs(p1 - p2);

        doublereal density = HF->dGetDensity((p1 + p2)/2.);

        if (jumpPres < CriticJump) {
                /*  moto sicuramente laminare (residuo) */
#ifdef HYDR_DEVEL
                DEBUGCOUT("we are in orifice laminar" << std::endl);
#endif
                flow = density*2.*(delta*delta)*diameter*m_Area_diaf*(p1 - p2)/viscosity;
                turbulent = false;

        } else {
                /*  moto sicuramente turbolento  (residuo) */
#ifdef HYDR_DEVEL
                DEBUGCOUT("we are in orifice turbulent:" << std::endl);
#endif
                flow = density*Cd*m_Area_diaf*copysign(sqrt(2.*jumpPres/density), p1 - p2);
                turbulent = true;
        }

        vel = flow/(density*m_Area_pipe);
        Re = fabs(density*vel*diameter/viscosity);

#ifdef HYDR_DEVEL
        DEBUGCOUT("jumpPres:       " << jumpPres << std::endl);
        DEBUGCOUT("density:        " << density << std::endl);
        DEBUGCOUT("Cd:             " << Cd << std::endl);
        DEBUGCOUT("p1:             " << p1 << std::endl);
        DEBUGCOUT("p2:             " << p2 << std::endl);
        DEBUGCOUT("Cc:             " << Cc << std::endl);
        DEBUGCOUT("CriticJump:     " << CriticJump << std::endl);
        DEBUGCOUT("delta:          " << delta << std::endl);
        DEBUGCOUT("viscosity:      " << viscosity << std::endl);
        DEBUGCOUT("rad:            " << rad << std::endl);
        DEBUGCOUT("area_pipe:      " << m_Area_pipe << std::endl);
        DEBUGCOUT("area_diaf:      " << m_Area_diaf << std::endl);
        DEBUGCOUT("RES area_pipe : " << m_Area_pipe << std::endl);
        DEBUGCOUT("RES flow:       " << flow << std::endl);
        DEBUGCOUT("RES Reynolds:   " << Re << std::endl);
        DEBUGCOUT("******************************************" << std::endl);
        DEBUGCOUT("RES velocita': " << vel << std::endl);
        DEBUGCOUT("    se positiva il fluido va dal nodo 1 al nodo 2" << std::endl);
        DEBUGCOUT("*********************************************" << std::endl);
#endif

        WorkVec.PutItem(1, iNode1RowIndex, flow);
        WorkVec.PutItem(2, iNode2RowIndex, -flow);

        return WorkVec;
}

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from Elem.

Definition at line 223 of file hminor.h.

References m_pNode1, and m_pNode2.

                                                                                 {
     connectedNodes.resize(2);
     connectedNodes[0] = m_pNode1;
     connectedNodes[1] = m_pNode2;
   };

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

virtual

Reimplemented from Elem.

Definition at line 501 of file hminor.cc.

References MBDYN_EXCEPT_ARGS.

{
        silent_cerr("Orifice has no dofs!" << std::endl);
        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
}

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

virtual

Implements HydraulicElem.

Definition at line 482 of file hminor.cc.

References HydraulicElem::ORIFICE.

{
        return HydraulicElem::ORIFICE;
}

unsigned int Orifice::iGetNumDof ( void  ) const

virtual

Reimplemented from Elem.

Definition at line 495 of file hminor.cc.

{
        return 0;
}

void Orifice::Output ( OutputHandler &  OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 630 of file hminor.cc.

References ToBeOutput::bToBeOutput(), flow, WithLabel::GetLabel(), OutputHandler::Hydraulic(), Re, turbulent, and vel.

{
        if (bToBeOutput()) {
                std::ostream& out = OH.Hydraulic();
                out << std::setw(8) << GetLabel()       /*  1 */
                        << " " << vel                   /*  2 */
                        << " " << flow                  /*  3 */
                        << " " << Re                    /*  4 */
                        << " " << turbulent             /*  5 */
                        << std::endl;
        }
}

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 489 of file hminor.cc.

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

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

virtual

Implements Elem.

Definition at line 508 of file hminor.cc.

{
        *piNumRows = 2;
        *piNumCols = 2;
}

Member Data Documentation

doublereal Orifice::Cd

private

Definition at line 182 of file hminor.h.

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

doublereal Orifice::CriticJump

private

Definition at line 180 of file hminor.h.

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

doublereal Orifice::delta

private

Definition at line 181 of file hminor.h.

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

doublereal Orifice::diameter

private

Definition at line 174 of file hminor.h.

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

doublereal Orifice::flow

private

Definition at line 183 of file hminor.h.

Referenced by AssRes(), and Output().

doublereal Orifice::m_Area_diaf

private

Definition at line 176 of file hminor.h.

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

doublereal Orifice::m_Area_pipe

private

Definition at line 177 of file hminor.h.

Referenced by AssRes(), and Orifice().

const PressureNode* Orifice::m_pNode1

private

Definition at line 172 of file hminor.h.

Referenced by AssJac(), AssRes(), GetConnectedNodes(), and Orifice().

const PressureNode* Orifice::m_pNode2

private

Definition at line 173 of file hminor.h.

Referenced by AssJac(), AssRes(), GetConnectedNodes(), and Orifice().

doublereal Orifice::Re

private

Definition at line 185 of file hminor.h.

Referenced by AssRes(), and Output().

doublereal Orifice::ReCr

private

Definition at line 178 of file hminor.h.

bool Orifice::turbulent

private

Definition at line 186 of file hminor.h.

Referenced by AssRes(), and Output().

doublereal Orifice::vel

private

Definition at line 184 of file hminor.h.

Referenced by AssRes(), and Output().

doublereal Orifice::viscosity

private

Definition at line 175 of file hminor.h.

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

---

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

• hminor.h
• hminor.cc