FlowModelVelocityBC.h File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Macros
#define	_FLOWMODELVELOCITYBC_H_

Functions
T	fixedFluxContinuityBC (const StorageSite &faces, const Mesh &mesh, MultiFieldMatrix &matrix, MultiField &xField, MultiField &rField, const FlowBC< T > &bc)

Macro Definition Documentation

#define _FLOWMODELVELOCITYBC_H_

Definition at line 7 of file FlowModel_impl.h.

Function Documentation

T fixedFluxContinuityBC	(	const StorageSite &	faces,
		const Mesh &	mesh,
		MultiFieldMatrix &	matrix,
		MultiField &	xField,
		MultiField &	rField,
		const FlowBC< T > &	bc
	)

Definition at line 11 of file FlowModelVelocityBC.h.

References FlowBC< T >::bcType, dot(), Mesh::getCells(), StorageSite::getCount(), Mesh::getFaceCells(), MultiFieldMatrix::getMatrix(), FloatVarDict< T >::getVal(), and MultiFieldMatrix::hasMatrix().

Referenced by FlowModel< T >::Impl::linearizeContinuity().

  {
    const StorageSite& cells = mesh.getCells();

    const CRConnectivity& faceCells = mesh.getFaceCells(faces);

    MultiField::ArrayIndex mfIndex(&_flowFields.massFlux,&faces);
    MultiField::ArrayIndex vIndex(&_flowFields.velocity,&cells);
    MultiField::ArrayIndex pIndex(&_flowFields.pressure,&cells);

    PPMatrix& ppMatrix =
      dynamic_cast<PPMatrix&>(matrix.getMatrix(pIndex,pIndex));

    FMatrix& dFluxdP = dynamic_cast<FMatrix&>(matrix.getMatrix(mfIndex,pIndex));

    PPAssembler& ppAssembler = ppMatrix.getPairWiseAssembler(faceCells);
    PPDiagArray& ppDiag = ppMatrix.getDiag();

    const VectorT3Array& faceArea =
      dynamic_cast<const VectorT3Array&>(_geomFields.area[faces]);

    TArray& rCell = dynamic_cast<TArray&>(rField[pIndex]);
    TArray& massFlux = dynamic_cast<TArray&>(_flowFields.massFlux[faces]);
    TArray& pCell = dynamic_cast<TArray&>(_flowFields.pressure[cells]);
    const TArray& density = dynamic_cast<const TArray&>(_flowFields.density[cells]);

    FloatValEvaluator<VectorT3>
      bVelocity(bc.getVal("specifiedXVelocity"),
                bc.getVal("specifiedYVelocity"),
                bc.getVal("specifiedZVelocity"),
                faces);
    
    FloatValEvaluator<T>  bp(bc.getVal("specifiedPressure"),  faces);

    const bool fixPressure = (bc.bcType == "FixedBoundary");
    
    const int nFaces = faces.getCount();

    T netFlux(0.);
    
    for(int f=0; f<nFaces; f++)
    {
        const int c0 = faceCells(f,0);
        const int c1 = faceCells(f,1);

        massFlux[f] = density[c0]*dot(bVelocity[f],faceArea[f]);

        rCell[c0] -= massFlux[f];

        netFlux += massFlux[f];
        ppAssembler.getCoeff01(f) =0;
        ppDiag[c1] = -1;
        if (fixPressure)
        {
            rCell[c1] = bp[f]-pCell[c1];
            ppAssembler.getCoeff10(f) =0;
        }
        else
        {
          rCell[c1] = 0.;
          ppAssembler.getCoeff10(f) =1;
        }
        
        ppMatrix.setBoundary(c1);

        dFluxdP.setCoeffL(f,T(0.));
        dFluxdP.setCoeffR(f,T(0.));
        
    }
#ifdef PV_COUPLED
    if (matrix.hasMatrix(vIndex,pIndex))
    {
        VPMatrix& vpMatrix =
          dynamic_cast<VPMatrix&>(matrix.getMatrix(vIndex,pIndex));
        
        VPAssembler& vpAssembler = vpMatrix.getPairWiseAssembler(faceCells);
        VPDiagArray& vpDiag = vpMatrix.getDiag();

        for(int f=0; f<nFaces; f++)
        {
            const int c0 = faceCells(f,0);
            vpDiag[c0] += vpAssembler.getCoeff01(f);
            vpAssembler.getCoeff01(f) = NumTypeTraits<T>::getZero();
        }
    }
#endif
    return netFlux;
  }