proteus.mprans.RANS2P_IB module

class proteus.mprans.RANS2P_IB.Coefficients(epsFact=1.5, sigma=72.8, rho_0=998.2, nu_0=1.004e-06, rho_1=1.205, nu_1=1.5e-05, g=[0.0, 0.0, - 9.8], nd=3, LS_model=None, VF_model=None, KN_model=None, Closure_0_model=None, Closure_1_model=None, epsFact_density=None, stokes=False, sd=True, movingDomain=False, useVF=0.0, useRBLES=0.0, useMetrics=0.0, useConstant_he=False, dragAlpha=0.01, dragBeta=0.0, setParamsFunc=None, dragAlphaTypes=None, dragBetaTypes=None, porosityTypes=None, killNonlinearDrag=False, epsFact_source=1.0, epsFact_solid=None, eb_adjoint_sigma=1.0, eb_penalty_constant=10.0, forceStrongDirichlet=False, turbulenceClosureModel=0, smagorinskyConstant=0.1, barycenters=None, beamLocation=[], beamLength=[], beamRadius=[], EI=[], GJ=[], nBeamElements=4, beam_quadOrder=3, beamFilename='Beams', beam_useSparse=False, beam_Cd=1.2, beam_nlTol=1e-05, beamRigid=True)

Bases: proteus.mprans.RANS2P.Coefficients

Set the number of components (equations) of the PDE and initialize the dicitionaries describing the form of the coefficients. Strings naming each component (used for viewing and archiving) and a structure defining the sparsity pattern of diffusion tensors may also be provided.

attachModels(modelList)

Give the TC object access to other models in a loosely coupled split operator formulation (e.g. a transport equation for concentration might get velocity from a flow equation)

initializeElementQuadrature(t, cq)

Give the TC object access to the element quadrature storage

initializeGlobalExteriorElementBoundaryQuadrature(t, cebqe)

Give the TC object access to the exterior element boundary quadrature storage

initializeMesh(mesh)

Give the TC object access to the mesh for any mesh-dependent information.

postStep(t, firstStep=False)

Give the TC object an opportunity to modify itself after the time step.

delta_h(r)

updateBeamLoad()

updateBeams(t)

initializeBeams()

class proteus.mprans.RANS2P_IB.LevelModel(uDict, phiDict, testSpaceDict, matType, dofBoundaryConditionsDict, dofBoundaryConditionsSetterDict, coefficients, elementQuadrature, elementBoundaryQuadrature, fluxBoundaryConditionsDict=None, advectiveFluxBoundaryConditionsSetterDict=None, diffusiveFluxBoundaryConditionsSetterDictDict=None, stressTraceBoundaryConditionsSetterDictDict=None, stabilization=None, shockCapturing=None, conservativeFluxDict=None, numericalFluxType=None, TimeIntegrationClass=None, massLumping=False, reactionLumping=False, options=None, name='RANS2P', reuse_trial_and_test_quadrature=True, sd=True, movingDomain=False)

Bases: proteus.mprans.RANS2P.LevelModel

Allocate storage and initialize some variables.

Parameters

• uDict (dict) – Dictionary of proteus.FemTools.FiniteElementFunction objects.

• phiDict (dict) – Dictionary of proteus.FemTools.FiniteElementFunction objects.

• testSpaceDict (dict) – Dictionary of FiniteElementSpace objects

• dofBoundaryConditionsDict (dict) – Dictionary of DOFBoundaryConditions objects for the Dirichlet conditions.

• coefficients (proteus.TransportCoefficients.TC_base) – Problem’s Transport Coefficients class.

• elementQuadratureDict (dict) – Dictionary of dictionaries of quadrature rules for each element integral in each component equation.

• elementBoundaryQuadratureDict (dict) – Dictionary of dictionaries of quadrature rules for each element boundary integral in each component equation

• stabilization (bool) –

• shockCapturing (bool) –

• numericalFlux (bool) –

• bdyNullSpace (bool) – Indicates whether the boundary conditions create a global null space.

Notes

The constructor sets the input arguments, calculates dimensions, and allocates storage. The meanings of variable suffixes are

• global – per physical domain

• element – per element

• elementBoundary – per element boundary

The prefix n means ‘number of’.

Storage is divided into quantities required at different sets of points or geometric entities. Each type of storage has a dictionary for all the quantities of that type. The names and dimensions of the storage dictionaries are

• e – at element

• q – at element quadrature, unique to elements

• ebq – at element boundary quadrature, unique to elements

• ebq_global – at element boundary quadrature, unique to element boundary

• ebqe – at element boundary quadrature, unique to global, exterior element boundary

• phi_ip – at the generalized interpolation points required to build a nonlinear phi

getResidual(u, r)

Calculate the element residuals and add in to the global residual

getJacobian(jacobian)

calculateAuxiliaryQuantitiesAfterStep()

beamStep()