proteus.mprans.AddedMass module

class proteus.mprans.AddedMass.NumericalFlux(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)[source]

Bases: proteus.NumericalFlux.ConstantAdvection_Diffusion_SIPG_exterior

class proteus.mprans.AddedMass.Coefficients(nd=2, V_model=None, barycenters=None, flags_rigidbody=None)[source]

Bases: proteus.TransportCoefficients.TC_base

Parameters:
  • nd (int) – Number of space dimensions.
  • V_model (int) – Index of Navier-Stokes model.
  • barycenters (array_like) – List of domain barycenters.
  • flags_rigidbody (array_like) – Array of integers of length at least as long as maximum domain flag, with 0 at indices corresponding to flags for fixed walls, and 1 at indices corresponding to flags for rigid bodies. e.g. if only flag 2 is rigid body, and the max value of flag is 5, the array will look like np.array([0,0,1,0,0,0])

TODO

attachModels(modelList)[source]

Attach the model for velocity and density to PresureIncrement model

initializeMesh(mesh)[source]

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

initializeElementQuadrature(t, cq)[source]

Give the TC object access to the element quadrature storage

initializeElementBoundaryQuadrature(t, cebq, cebq_global)[source]

Give the TC object access to the element boundary quadrature storage

initializeGlobalExteriorElementBoundaryQuadrature(t, cebqe)[source]

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

initializeGeneralizedInterpolationPointQuadrature(t, cip)[source]

Give the TC object access to the generalized interpolation point storage. These points are used to project nonlinear potentials (phi).

preStep(t, firstStep=False)[source]

Move the current values to values_last to keep cached set of values for bdf1 algorithm

postStep(t, firstStep=False)[source]

Update the fluid velocities

evaluate(t, c)[source]
class proteus.mprans.AddedMass.LevelModel(uDict, phiDict, testSpaceDict, matType, dofBoundaryConditionsDict, dofBoundaryConditionsSetterDict, coefficients, elementQuadrature, elementBoundaryQuadrature, fluxBoundaryConditionsDict=None, advectiveFluxBoundaryConditionsSetterDict=None, diffusiveFluxBoundaryConditionsSetterDictDict=None, stressTraceBoundaryConditionsSetterDict=None, stabilization=None, shockCapturing=None, conservativeFluxDict=None, numericalFluxType=None, TimeIntegrationClass=None, massLumping=False, reactionLumping=False, options=None, name='defaultName', reuse_trial_and_test_quadrature=True, sd=True, movingDomain=False, bdyNullSpace=True)[source]

Bases: proteus.Transport.OneLevelTransport

nCalls = 0[source]
calculateCoefficients()[source]
calculateElementResidual()[source]
getResidual(u, r)[source]
getJacobian(jacobian)[source]
calculateElementQuadrature()[source]

Calculate the physical location and weights of the quadrature rules and the shape information at the quadrature points.

This function should be called only when the mesh changes.

calculateElementBoundaryQuadrature()[source]
calculateExteriorElementBoundaryQuadrature()[source]
estimate_mt()[source]
calculateAuxiliaryQuantitiesAfterStep()[source]
calculateSolutionAtQuadrature()[source]
updateAfterMeshMotion()[source]