proteus.mprans.RDLS module¶

class proteus.mprans.RDLS.SubgridError(coefficients, nd)[source]¶

Bases: proteus.SubgridError.SGE_base

initializeElementQuadrature(mesh, t, cq)[source]¶

calculateSubgridError(q)[source]¶

updateSubgridErrorHistory(initializationPhase=False)[source]¶

class proteus.mprans.RDLS.ShockCapturing(coefficients, nd, shockCapturingFactor=0.25, lag=True, nStepsToDelay=None)[source]¶

Bases: proteus.ShockCapturing.ShockCapturing_base

initializeElementQuadrature(mesh, t, cq)[source]¶

updateShockCapturingHistory()[source]¶

class proteus.mprans.RDLS.PsiTC(model, nOptions)[source]¶

Bases: proteus.StepControl.SC_base

stepExact_model(tOut)[source]¶

initialize_dt_model(t0, tOut)[source]¶

updateSubstep()[source]¶

choose_dt_model()[source]¶

class proteus.mprans.RDLS.Coefficients(applyRedistancing=True, epsFact=2.0, nModelId=None, u0=None, rdModelId=0, penaltyParameter=0.0, useMetrics=0.0, useConstantH=False, weakDirichletFactor=10.0, backgroundDiffusionFactor=0.01, computeMetrics=False, ELLIPTIC_REDISTANCING=0, alpha=1000000000.0, backgroundDissipationEllipticRedist=1.0, outputQuantDOFs=False)[source]¶

Bases: proteus.TransportCoefficients.TC_base

attachModels(modelList)[source]¶

initializeMesh(mesh)[source]¶

initializeElementQuadrature(t, cq)[source]¶

initializeElementBoundaryQuadrature(t, cebq, cebq_global)[source]¶

initializeGlobalExteriorElementBoundaryQuadrature(t, cebqe)[source]¶

preStep(t, firstStep=False)[source]¶

postStep(t, firstStep=False)[source]¶

getICDofs(cj)[source]¶

updateToMovingDomain(t, c)[source]¶

evaluate(t, c)[source]¶

static setZeroLSweakDirichletBCs(vt)[source]¶

static setZeroLSweakDirichletBCs2(vt)[source]¶

static setZeroLSweakDirichletBCs3(vt)[source]¶

redistanceLevelSetCoefficientsEvaluate()[source]¶

evaluate the coefficients for level set redistancing equation

class proteus.mprans.RDLS.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=False)[source]¶

Bases: proteus.Transport.OneLevelTransport

nCalls = 0[source]¶

runAtEOS()[source]¶

getMetricsAtEOS(u_exact)[source]¶

calculateCoefficients()[source]¶

calculateElementResidual()[source]¶

getNormalReconstruction()[source]¶

getResidual(u, r)[source]¶

getJacobian(jacobian)[source]¶

getJacobianNew(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]¶

Calculate the physical location and weights of the quadrature rules and the shape information at the quadrature points on global element boundaries.

This function should be called only when the mesh changes.

estimate_mt()[source]¶

calculateSolutionAtQuadrature()[source]¶

calculateAuxiliaryQuantitiesAfterStep()[source]¶

updateAfterMeshMotion()[source]¶

proteus.mprans.RDLS.setZeroLSweakDirichletBCs(RDLSvt)[source]¶

proteus.mprans.RDLS.setZeroLSweakDirichletBCsSimple(RDLSvt)[source]¶