proteus.NumericalFlux module

A class hierarchy for numerical flux (numerical trace) computations

class proteus.NumericalFlux.NF_base(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, parallelPeriodic=False)

Bases: object

useWeakDirichletConditions = True

useStrongDirichletConstraints = False

hasInterior = True

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

setFromOptions(nOptions)

allow classes to set various numerical parameters

class proteus.NumericalFlux.DoNothing(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

setDirichletValues(ebqe)

class proteus.NumericalFlux.NoFlux(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, parallelPeriodic=False)

Bases: proteus.NumericalFlux.NF_base

calculateInteriorNumericalFlux(q, ebq, ebq_global)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.StrongDirichlet(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NF_base

useWeakDirichletConditions = False

hasInterior = False

setFluxBoundaryConditions()

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

proteus.NumericalFlux.StrongDirichletFactory(fluxBoundaryConditionsDict={})

class proteus.NumericalFlux.Advection_DiagonalUpwind(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NF_base

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Advection_Diagonal_average(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_IIPG(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_SIPG(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_IIPG

class proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_NIPG(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_IIPG

class proteus.NumericalFlux.Advection_DiagonalUpwind_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.Advection_DiagonalUpwind

hasInterior = False

calculateInteriorNumericalFlux(q, ebq, ebq_global)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.ConstantAdvection_Diffusion_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.ConstantAdvection_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

useStrongDirichletConstraints = True

hasInterior = False

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.MixedDarcy_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_NIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_IIPG_exterior

class proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_IIPG_exterior

class proteus.NumericalFlux.ConstantAdvection_Diffusion_NIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.ConstantAdvection_Diffusion_IIPG_exterior

class proteus.NumericalFlux.ConstantAdvection_Diffusion_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.ConstantAdvection_Diffusion_IIPG_exterior

class proteus.NumericalFlux.Advection_DiagonalUpwind_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Curvature_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Stokes_Advection_DiagonalUpwind_Diffusion_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

To use with regular Stokes, takes advantage of existence of ‘advectiveFlux’ flag even when there is no advective term

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Stokes_Advection_DiagonalUpwind_Diffusion_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.Stokes_Advection_DiagonalUpwind_Diffusion_IIPG_exterior

class proteus.NumericalFlux.StokesP_Advection_DiagonalUpwind_Diffusion_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

Weak boundary fluxes for StokesP

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.NavierStokes_Advection_DiagonalUpwind_Diffusion_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.NavierStokes_Advection_DiagonalUpwind_Diffusion_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NavierStokes_Advection_DiagonalUpwind_Diffusion_IIPG_exterior

hasInterior = False

class proteus.NumericalFlux.Diffusion_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Diffusion_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.Diffusion_IIPG_exterior

class proteus.NumericalFlux.DarcySplitPressure_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

weak dirichlet boundary conditions for Twophase_split_pressure class

Diffusion_IIPG_exterior is ok except for need to switch between psi_w and psi_n bc types

hasInterior = False

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Diffusion_LDG(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

initial LDG numerical flux for pure diffusion

calculateDiffusionSplittings(q, ebq, ebq_global)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_LDG(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.Diffusion_LDG

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.RusanovNumericalFlux_Diagonal(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicConditions)

Bases: proteus.NumericalFlux.Advection_DiagonalUpwind

apply numerical flux $f_{num}(a,b)=1/2(f(a)+f(b)-\overline{\lambda }(b-a)$ where $\lambda >=max|f^{\prime }|fora<=u<=b$

this one applies flux to each component of flux separately

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

class proteus.NumericalFlux.RusanovNumericalFlux_Diagonal_Diffusion_IIPG(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.Advection_DiagonalUpwind_Diffusion_IIPG

apply numerical flus $f_{num}(a,b)=1/2(f(a)+f(b)-\overline{\lambda }(b-a)$ where $\lambda >=max|f^{\prime }|$ for $a<=u<=b$ this one applies flux to each component of flux separately

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

class proteus.NumericalFlux.ConvexOneSonicPointNumericalFlux(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions, sonicPoint=0.0, sonicFlux=0.0)

Bases: proteus.NumericalFlux.Advection_DiagonalUpwind

basic Godunov flux $f_{num}(a,b)=ma{x}_{b<=u<=a}f(u)$ if $a>=b=mi{n}_{a<=u<=b}f(u)$ otherwise where there is only one sonic point, $u_s$ with $f^{\prime }(u_s)=0$ and $f$ is convex so $f(u_s)$ is a minimum

This class typically has to be “wrapped” for a given problem to specify the correct sonic point and sonic flux

calculateInteriorNumericalFlux(q, ebq, ebq_global)

class proteus.NumericalFlux.HamiltonJacobi_DiagonalLesaintRaviart(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, speedEvaluationType=2)

Bases: proteus.NumericalFlux.NF_base

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.HamiltonJacobi_DiagonalLesaintRaviart_Diffusion_IIPG(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, speedEvaluationType=1)

Bases: proteus.NumericalFlux.NF_base

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.HamiltonJacobi_DiagonalLesaintRaviart_Diffusion_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, speedEvaluationType=1)

Bases: proteus.NumericalFlux.Diffusion_SIPG_exterior

setDirichletValues(ebqe)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.HamiltonJacobi_Pressure_DiagonalLesaintRaviart_Diffusion_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, speedEvaluationType=1)

Bases: proteus.NumericalFlux.HamiltonJacobi_DiagonalLesaintRaviart_Diffusion_SIPG_exterior

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

class proteus.NumericalFlux.DarcyFCFF_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

weak dirichlet boundary conditions for Twophase_fc_ff class

hasInterior = False

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.DarcyFC_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

weak dirichlet boundary conditions for Twophase_fc class

hasInterior = False

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.DarcyFCPP_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

weak dirichlet boundary conditions for Twophase_fc class

hasInterior = False

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.ShallowWater_1D(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, h_eps=1e-08, tol_u=1e-08)

Bases: proteus.NumericalFlux.NF_base

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.ShallowWaterHLL_1D(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, h_eps=1e-08, tol_u=1e-08)

Bases: proteus.NumericalFlux.NF_base

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.ShallowWater_2D(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, h_eps=1e-08, tol_u=1e-08)

Bases: proteus.NumericalFlux.NF_base

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.RusanovNumericalFlux(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.RusanovNumericalFlux_Diagonal

Base class for Rusanov scheme for generic systems, relies on a user-specified estimate for the maximum (magnitude) eigenvalue for the system.

Default is to take

$$\overline{\text{lamda}}=max|{\overrightarrow{f}}^{\prime }\cdot n\Vert$$

Then we just apply apply numerical flux

$$f_{num}(a,b)=1/2(f(a)+f(b)-\overline{\lambda }(b-a)$$

For now, we will try to piggy back on the cfl calculation, assuming the user has implemented this correctly for the system in question

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

class proteus.NumericalFlux.RusanovLDG(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.Diffusion_LDG

combine Rusanov flux for advection and LDG for diffusion right now advection jacobian is not calculated

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

class proteus.NumericalFlux.HamiltonJacobi_DiagonalChengShu(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None, speedEvaluationType=0)

Bases: proteus.NumericalFlux.NF_base

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Stress_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

setDirichletValues(ebqe)

calculateInteriorNumericalFlux(q, ebq, ebq_global)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateInteriorNumericalFluxJacobian(l2g, q, ebq, ebq_global, dphi, fluxJacobian, fluxJacobian_eb, fluxJacobian_hj)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Stress_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions, getPeriodicBoundaryConditions=None)

Bases: proteus.NumericalFlux.Stress_IIPG_exterior

hasInterior = False

class proteus.NumericalFlux.Richards_IIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.NF_base

hasInterior = False

setDirichletValues(ebqe)

calculateExteriorNumericalFlux(inflowFlag, q, ebqe)

updateExteriorNumericalFluxJacobian(l2g, inflowFlag, q, ebqe, dphi, fluxJacobian_exterior, fluxJacobian_eb, fluxJacobian_hj)

class proteus.NumericalFlux.Richards_SIPG_exterior(vt, getPointwiseBoundaryConditions, getAdvectiveFluxBoundaryConditions, getDiffusiveFluxBoundaryConditions)

Bases: proteus.NumericalFlux.Richards_IIPG_exterior

hasInterior = False