m numerical modeling / IICSM

THE NUMERICAL MODELING WORKFLOW.

From equations to reliable results.

Browse all techniques ↓
DISCRETIZE · SOLVE · VERIFY
4 −1 0−1 4 −10 −1 4MESHAx = bERROR → REFINEMENT
Schematic workflow, not a computed result

01 / DISCRETIZE

Discretize the problem

Turn continuous equations into finite-dimensional approximations.

Workflow stages are a navigation guide; numerical techniques often span several stages. Select a technique to open its description; use “Go back to the slider” to return here.

HOW THE TECHNIQUES CONNECT

A tree of numerical techniques.

Explore workflow stage → method family → technique. Branches organize all 68 entries; select a technique to see its description and relationships. Connections between families appear in each technique’s “Relationships to other techniques” section.

This is a browsing hierarchy, not a universal dependency tree. “Same discipline” means shared subject area; specific links distinguish approximations, extensions, closures, numerical methods, and coupling.

68 techniques in the tree

Numerical modeling techniques 68 entries
Discretize 8
Spatial discretization 8
Solve algebra 20
Linear algebra 12
Nonlinear equations 8
Advance in time 8
Time integration 8
Optimize 8
Optimization & inverse problems 8
Approximate & sample 16
Approximation & reduction 8
Quadrature & sampling 8
Verify & assess 8
Error analysis & verification 8
↑ Go back to the slider

THE PHYSICAL WORLD, REPRESENTED

Small particles.
Big systems.
A world of models.

Explore the ideas we use to understand, predict, and engineer reality — from molecular interactions to planetary systems.

Explore the catalog ↓

68 entries · 8 disciplines · References for every entry

84 cumulative daily visits · View visitor map ↘

CHANGE YOUR PERSPECTIVE

Explore by workflow

THE REFERENCE COLLECTION

Find your numerical technique.

Download catalog JSON ↓

A curated, expandable atlas of numerical techniques. It is not exhaustive. Workflow labels indicate a typical role, and techniques can be combined across stages.

Showing 68 of 68 entries

Spatial discretization001

Finite difference method

Approximates derivatives with weighted values on a grid.

DiscretizeNumerical technique
Formulation & short derivation

Representative numerical formulation

u′′(xi)≈ui+1−2ui+ui−1h2u''(x_i)\approx\frac{u_{i+1}-2u_i+u_{i-1}}{h^2}

Derivation / construction sketch

  1. Expand neighboring values in Taylor series.
  2. Add the expansions to cancel odd derivatives.
  3. Divide by h squared to obtain the centered second derivative.

Symbols & assumptions

Second-order accuracy requires a smooth solution and a uniform grid; boundaries need compatible stencils.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Spatial discretization

Practical use

Temperature fields in heat sinks and plates.

Engineering application examples

Temperature fields in heat sinks and plates.

Software product or implementation route

Custom finite-difference solver ↗

Implementation route: assemble stencils in Python/NumPy or a compiled solver; this reference is a textbook, not a packaged application.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Temperature fields in heat sinks and plates.

Method limitations

Second-order accuracy requires a smooth solution and a uniform grid; boundaries need compatible stencils.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • DiscretizationQuantum harmonic oscillator ↗

    Describes a quantum degree of freedom in a quadratic potential.

    For fields on structured grids; design boundary stencils and check mesh convergence.

  • DiscretizationParticle-in-a-box model ↗

    Confines a quantum particle within idealized boundaries.

    For fields on structured grids; design boundary stencils and check mesh convergence.

Relationships to other techniques

Discretize → Spatial discretization

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Spatial discretization002

Finite volume method

Balances conserved quantities over control volumes.

DiscretizeNumerical technique
Formulation & short derivation

Representative numerical formulation

VidUidt+∑fF^f⋅nfAf=ViSiV_i\frac{dU_i}{dt}+\sum_f\widehat F_f\cdot n_f A_f=V_iS_i

Derivation / construction sketch

  1. Integrate the conservation law over a cell.
  2. Use the divergence theorem to convert volume flux divergence to surface flux.
  3. Approximate each face flux consistently with its neighbor.

Symbols & assumptions

Accuracy depends on reconstruction, numerical flux, mesh quality, and time integration.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Spatial discretization

Practical use

Pipe-flow and aerodynamic CFD workflows.

Engineering application examples

Pipe-flow and aerodynamic CFD workflows.

Software product or implementation route

OpenFOAM ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Pipe-flow and aerodynamic CFD workflows.

Method limitations

Accuracy depends on reconstruction, numerical flux, mesh quality, and time integration.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Conservative discretizationPopulation balance model ↗

    Tracks the distribution of particle sizes or other internal properties.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationMass-action reaction kinetics ↗

    Relates reaction rates to species concentrations and reaction orders.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationFickian diffusion ↗

    Relates diffusive flux to concentration gradients.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationMaxwell–Stefan diffusion ↗

    Represents multicomponent diffusion through interspecies friction.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationAdvection–diffusion–reaction model ↗

    Combines bulk transport, diffusion and reaction sources.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationNavier–Stokes model ↗

    Conserves mass and momentum for a viscous continuum fluid.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationEuler flow model ↗

    Neglects viscous stresses in compressible or incompressible flow.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationStokes creeping-flow model ↗

    Neglects inertial terms relative to viscosity.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationBoundary-layer model ↗

    Resolves thin near-wall regions with scale-based simplifications.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationLubrication approximation ↗

    Simplifies viscous flow in thin gaps.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationOldroyd-B model ↗

    Combines solvent viscosity with an elastic polymer stress.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationReynolds-averaged Navier–Stokes (RANS) ↗

    Models mean flow with closure for unresolved turbulent stresses.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationSpalart–Allmaras model ↗

    Uses a transported turbulence variable to obtain eddy viscosity.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationk–epsilon model ↗

    Uses turbulent kinetic energy and dissipation rate to close mean flow.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationk–omega model ↗

    Uses turbulent kinetic energy and specific dissipation rate.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationSST k–omega model ↗

    Blends near-wall and outer-flow behavior with a shear-stress limiter.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationReynolds-stress transport model ↗

    Transports individual turbulent stress components.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationLarge-eddy simulation (LES) ↗

    Resolves larger turbulent motions and models subgrid effects.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationSmagorinsky subgrid model ↗

    Relates subgrid eddy viscosity to resolved strain and filter scale.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationDetached-eddy simulation (DES) ↗

    Combines RANS near walls with LES-like treatment away from them.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationVolume-of-fluid (VOF) representation ↗

    Tracks phase volume fractions to represent an interface.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationEuler–Euler two-fluid model ↗

    Treats phases as interpenetrating continua with exchange terms.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationLagrangian particle tracking ↗

    Tracks discrete particles through a carrier flow.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationRadiative transfer equation ↗

    Tracks radiation intensity through emission, absorption and scattering.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationPoisson–Nernst–Planck model ↗

    Couples electrostatics to diffusion and migration of ions.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationDoyle–Fuller–Newman (DFN/P2D) model ↗

    Combines porous-electrode transport and particle diffusion.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationSingle-particle battery model (SPM) ↗

    Represents each electrode by a representative active-material particle.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationSingle-particle model with electrolyte (SPMe) ↗

    Adds electrolyte concentration effects to a single-particle approximation.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationSaint-Venant shallow-water model ↗

    Depth-averages mass and momentum in free-surface flow.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationKinematic-wave routing ↗

    Simplifies flow routing by approximating dominant slope and friction balance.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationGroundwater flow model ↗

    Combines water conservation with porous-flow relations.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationAdvection–dispersion groundwater model ↗

    Represents contaminant transport and spreading through an aquifer.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationNumerical weather prediction ↗

    Evolves atmospheric dynamics and thermodynamics from an analyzed initial state.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationGeneral circulation model (GCM) ↗

    Represents large-scale atmospheric or oceanic circulation.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationEarth system model (ESM) ↗

    Couples atmosphere, ocean, land, ice and biogeochemical processes.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationOcean circulation model ↗

    Evolves ocean momentum, temperature and salinity.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationSea-ice thermodynamic-dynamic model ↗

    Couples freezing, melting and ice motion.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationElastic seismic-wave model ↗

    Propagates elastic disturbances through Earth materials.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationBoltzmann kinetic equation ↗

    Evolves a particle distribution under transport and collisions.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationVlasov–Poisson model ↗

    Couples collisionless distribution dynamics to electrostatic fields.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationVlasov–Maxwell model ↗

    Couples collisionless kinetic distributions to electromagnetic fields.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationMagnetohydrodynamics (MHD) ↗

    Treats a conducting fluid coupled to a magnetic field.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationNeutron diffusion approximation ↗

    Simplifies neutron transport to a diffusion description.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

  • Conservative discretizationGeneral relativity model ↗

    Relates spacetime curvature to matter and energy.

    For transport/balance-law formulations; use consistent face fluxes and appropriate reconstruction.

Relationships to other techniques

Discretize → Spatial discretization

Specific connections

  • Can be time-integrated by Embedded Runge-Kutta RK45

    A method-of-lines discretization can be advanced by a suitable ODE integrator; stability must be checked.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Spatial discretization003

Finite element method

Uses piecewise basis functions and a weak formulation on a mesh.

DiscretizeNumerical technique
Formulation & short derivation

Representative numerical formulation

uh=∑jUjNju_h=\sum_j U_jN_jKij=∫Ω∇Ni⋅k∇Nj dΩK_{ij}=\int_\Omega\nabla N_i\cdot k\nabla N_j\,d\OmegaKU=fKU=f

Derivation / construction sketch

  1. Multiply a diffusion equation by a test function.
  2. Integrate by parts and impose boundary conditions.
  3. Expand in basis functions and assemble the element contributions.

Symbols & assumptions

Diffusion example shown; elements and function spaces must suit the PDE and boundary conditions.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Spatial discretization

Practical use

Structural brackets, heat exchangers, and electromagnetic components.

Engineering application examples

Structural brackets, heat exchangers, and electromagnetic components.

Software product or implementation route

MFEM ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Structural brackets, heat exchangers, and electromagnetic components.

Method limitations

Diffusion example shown; elements and function spaces must suit the PDE and boundary conditions.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • DiscretizationSchrödinger model ↗

    Evolves a nonrelativistic quantum state using a Hamiltonian.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationDirac model ↗

    Describes relativistic spin-half particles with a spinor wave equation.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationBorn–Oppenheimer approximation ↗

    Separates electronic motion from slower nuclear motion.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationPotential-flow model ↗

    Represents irrotational velocity using a scalar potential.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationTransient heat equation ↗

    Balances thermal storage, conduction and heat sources.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationStefan phase-change problem ↗

    Couples heat transport to a moving melting or freezing boundary.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationEnthalpy–porosity model ↗

    Represents melting using enthalpy and a porous resistance in the mushy zone.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationLinear elasticity (Hooke model) ↗

    Relates stress linearly to small elastic strain.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationOrthotropic elasticity ↗

    Uses direction-dependent elastic properties along material axes.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationNeo-Hookean hyperelasticity ↗

    Models large elastic deformation with a strain-energy function.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationMooney–Rivlin hyperelasticity ↗

    Uses multiple strain invariants to fit rubber-like response.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationOgden hyperelasticity ↗

    Uses powers of principal stretches to represent nonlinear elasticity.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationEuler–Bernoulli beam model ↗

    Describes slender-beam bending while neglecting transverse shear deformation.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationTimoshenko beam model ↗

    Includes transverse shear deformation and rotational effects.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationKirchhoff–Love plate model ↗

    Describes thin-plate bending with normals remaining normal.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationMindlin–Reissner plate model ↗

    Includes transverse shear deformation in plate bending.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationShell model ↗

    Combines membrane and bending behavior on a curved surface.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationTruss model ↗

    Represents a structure with axial-force members joined at idealized nodes.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationCable and membrane models ↗

    Represent slender or thin structures dominated by tension.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • Discretizationvon Mises J2 plasticity ↗

    Uses deviatoric stress to define yielding in an isotropic ductile material.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationTresca yield model ↗

    Defines yield using maximum shear stress.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationDrucker–Prager plasticity ↗

    Uses a smooth pressure-dependent yield surface.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationMohr–Coulomb model ↗

    Relates frictional shear strength to normal stress and cohesion.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationJohnson–Cook model ↗

    Uses empirical strain, strain-rate and temperature factors.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationCrystal plasticity ↗

    Represents plastic flow through crystallographic slip systems.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationPhase-field fracture model ↗

    Represents cracks with a continuous damage-like field.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationLinear acoustic wave model ↗

    Describes small pressure perturbations about an equilibrium state.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationHelmholtz acoustic model ↗

    Represents harmonic acoustic fields at one frequency.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationTransmission-line acoustic model ↗

    Uses distributed wave propagation in a narrow duct or tube.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationMaxwell electromagnetic model ↗

    Couples electric and magnetic fields with charges and currents.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationElectrostatic Poisson model ↗

    Relates electric potential to charge density.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationMagnetostatic model ↗

    Represents steady magnetic fields driven by currents and magnetization.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationEddy-current model ↗

    Models induced conducting currents in a time-varying magnetic field.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationDarcy porous-flow model ↗

    Relates averaged fluid flux to hydraulic gradient.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationBrinkman porous-flow model ↗

    Adds a viscous shear term to a Darcy-like resistance model.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationRichards equation ↗

    Describes variably saturated water movement in porous media.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationBiot poroelasticity ↗

    Couples solid deformation and pore-fluid pressure.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationTerzaghi consolidation model ↗

    Describes time-dependent settlement from pore-pressure dissipation.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationModified Cam-Clay model ↗

    Uses critical-state plasticity for idealized clay behavior.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationPennes bioheat model ↗

    Adds perfusion and metabolic heat to tissue heat transfer.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationReaction–diffusion morphogenesis model ↗

    Couples reacting substances with diffusion.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationHomogenization ↗

    Derives effective properties or equations from smaller-scale structure.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationRepresentative volume element (RVE) ↗

    Uses a finite microstructural sample to estimate bulk response.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationFE² computational homogenization ↗

    Solves microscale problems within a macroscale finite-element calculation.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationQM/MM coupling ↗

    Combines quantum mechanics in a selected region with molecular mechanics around it.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationAtomistic–continuum coupling ↗

    Connects particle-level and continuum descriptions.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationFluid–structure interaction (FSI) ↗

    Couples fluid loads with structural motion or deformation.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

  • DiscretizationThermomechanical coupling ↗

    Couples temperature evolution and mechanical response.

    For a suitable weak-form spatial problem; choose function spaces and boundary conditions for the operator.

Relationships to other techniques

Discretize → Spatial discretization

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Spatial discretization004

Discontinuous Galerkin method

Combines element-local trial functions with numerical fluxes across interfaces.

DiscretizeNumerical technique
Formulation & short derivation

Representative numerical formulation

∫Kv∂tuh−∫K∇v⋅F(uh)+∫∂KvF^⋅n=0\int_K v\partial_tu_h-\int_K\nabla v\cdot F(u_h)+\int_{\partial K}v\widehat F\cdot n=0

Derivation / construction sketch

  1. Use independent polynomial spaces on each element.
  2. Integrate the conservation law against local test functions.
  3. Couple neighbors through a stable numerical flux.

Symbols & assumptions

Flux choice and stabilization matter; explicit high-order schemes can need small time steps.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Spatial discretization

Practical use

Compressible-flow simulation for nozzles and wings.

Engineering application examples

Compressible-flow simulation for nozzles and wings.

Software product or implementation route

MFEM ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Compressible-flow simulation for nozzles and wings.

Method limitations

Flux choice and stabilization matter; explicit high-order schemes can need small time steps.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • DiscretizationNavier–Stokes model ↗

    Conserves mass and momentum for a viscous continuum fluid.

    For element-local high-order transport or wave formulations; stable interface fluxes and time steps are essential.

  • DiscretizationLinear acoustic wave model ↗

    Describes small pressure perturbations about an equilibrium state.

    For element-local high-order transport or wave formulations; stable interface fluxes and time steps are essential.

  • DiscretizationHelmholtz acoustic model ↗

    Represents harmonic acoustic fields at one frequency.

    For element-local high-order transport or wave formulations; stable interface fluxes and time steps are essential.

  • DiscretizationTransmission-line acoustic model ↗

    Uses distributed wave propagation in a narrow duct or tube.

    For element-local high-order transport or wave formulations; stable interface fluxes and time steps are essential.

Relationships to other techniques

Discretize → Spatial discretization

Specific connections

  • Discontinuous finite-element family of Finite element method

    Uses element-local spaces with interface flux coupling.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Spatial discretization005

Boundary element method

Transfers suitable linear PDE problems to boundary integral equations.

DiscretizeNumerical technique
Formulation & short derivation

Representative numerical formulation

c(x)u(x)+∫Γu∂nG dΓ=∫ΓG∂nu dΓc(x)u(x)+\int_\Gamma u\partial_nG\,d\Gamma=\int_\Gamma G\partial_nu\,d\Gamma

Derivation / construction sketch

  1. Apply Green's identity with a fundamental solution.
  2. Take the observation point to the boundary.
  3. Discretize the boundary fields and integrals.

Symbols & assumptions

Representative Laplace formulation; singular quadrature and dense matrices require care.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Spatial discretization

Practical use

Exterior acoustics, electrostatics, and scattering models.

Engineering application examples

Exterior acoustics, electrostatics, and scattering models.

Software product or implementation route

Bempp ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Exterior acoustics, electrostatics, and scattering models.

Method limitations

Representative Laplace formulation; singular quadrature and dense matrices require care.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Boundary formulationPotential-flow model ↗

    Represents irrotational velocity using a scalar potential.

    For a linear homogeneous-domain formulation with a known fundamental solution; general nonlinear/inhomogeneous problems need extensions.

  • Boundary formulationLinear acoustic wave model ↗

    Describes small pressure perturbations about an equilibrium state.

    For a linear homogeneous-domain formulation with a known fundamental solution; general nonlinear/inhomogeneous problems need extensions.

  • Boundary formulationHelmholtz acoustic model ↗

    Represents harmonic acoustic fields at one frequency.

    For a linear homogeneous-domain formulation with a known fundamental solution; general nonlinear/inhomogeneous problems need extensions.

  • Boundary formulationTransmission-line acoustic model ↗

    Uses distributed wave propagation in a narrow duct or tube.

    For a linear homogeneous-domain formulation with a known fundamental solution; general nonlinear/inhomogeneous problems need extensions.

  • Boundary formulationElectrostatic Poisson model ↗

    Relates electric potential to charge density.

    For a linear homogeneous-domain formulation with a known fundamental solution; general nonlinear/inhomogeneous problems need extensions.

  • Boundary formulationMagnetostatic model ↗

    Represents steady magnetic fields driven by currents and magnetization.

    For a linear homogeneous-domain formulation with a known fundamental solution; general nonlinear/inhomogeneous problems need extensions.

Relationships to other techniques

Discretize → Spatial discretization

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Spatial discretization006

Spectral collocation

Approximates smooth fields globally and enforces the equation at selected nodes.

DiscretizeNumerical technique
Formulation & short derivation

Representative numerical formulation

uN(x)=∑j=0NUjℓj(x)u_N(x)=\sum_{j=0}^NU_j\ell_j(x)uN′(xi)=∑jDijUju_N'(x_i)=\sum_jD_{ij}U_j

Derivation / construction sketch

  1. Interpolate a field with a global polynomial or Fourier basis.
  2. Differentiate the basis analytically.
  3. Enforce the PDE at collocation nodes.

Symbols & assumptions

Rapid convergence requires sufficient smoothness; discontinuities cause oscillations and aliasing.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Spatial discretization

Practical use

Smooth fluid instabilities and wave propagation.

Engineering application examples

Smooth fluid instabilities and wave propagation.

Software product or implementation route

Dedalus ↗

Dedalus provides spectral methods; the paper discusses its formulation and implementation. Collocation details depend on basis and problem.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Smooth fluid instabilities and wave propagation.

Method limitations

Rapid convergence requires sufficient smoothness; discontinuities cause oscillations and aliasing.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Smooth-field discretizationSchrödinger model ↗

    Evolves a nonrelativistic quantum state using a Hamiltonian.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationDirac model ↗

    Describes relativistic spin-half particles with a spinor wave equation.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationBorn–Oppenheimer approximation ↗

    Separates electronic motion from slower nuclear motion.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationTime-dependent DFT (TDDFT) ↗

    Evolves electron density to approximate excited-state response.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationQuantum harmonic oscillator ↗

    Describes a quantum degree of freedom in a quadratic potential.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationParticle-in-a-box model ↗

    Confines a quantum particle within idealized boundaries.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationCahn–Hilliard model ↗

    Evolves a conserved composition field through chemical-potential gradients.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationAllen–Cahn model ↗

    Evolves a nonconserved order parameter toward lower free energy.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationPhase-field crystal model ↗

    Uses a periodic density-like field to represent crystalline ordering.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

  • Smooth-field discretizationScalar diffraction model ↗

    Uses a scalar wave approximation for light diffraction.

    For sufficiently smooth fields in compatible geometries; use suitable bases, dealiasing, and boundary treatment.

Relationships to other techniques

Discretize → Spatial discretization

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Spatial discretization007

Radial basis function discretization

Builds meshfree interpolants or derivative stencils from radial kernels.

DiscretizeNumerical technique
Formulation & short derivation

Representative numerical formulation

uh(x)=∑jcjϕ(∥x−xj∥)u_h(x)=\sum_jc_j\phi(\lVert x-x_j\rVert)(LΦ)c=f(L\Phi)c=f

Derivation / construction sketch

  1. Choose kernel centers and a radial basis.
  2. Apply the differential operator to the basis functions.
  3. Collocate interior and boundary equations to solve for coefficients.

Symbols & assumptions

Polynomial augmentation and shape parameters affect solvability and conditioning.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Spatial discretization

Practical use

Scattered geometry interpolation and meshfree field approximation.

Engineering application examples

Scattered geometry interpolation and meshfree field approximation.

Software product or implementation route

SciPy RBFInterpolator ↗

Interpolation implementation; a PDE collocation solver requires applying operators and boundary conditions separately.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Scattered geometry interpolation and meshfree field approximation.

Method limitations

Polynomial augmentation and shape parameters affect solvability and conditioning.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Scattered-field approximationFickian diffusion ↗

    Relates diffusive flux to concentration gradients.

    For scattered samples or a suitable meshfree collocation formulation; test conditioning and boundary accuracy.

Relationships to other techniques

Discretize → Spatial discretization

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Spatial discretization008

Adaptive mesh refinement

Concentrates degrees of freedom where a numerical error indicator is large.

DiscretizeNumerical technique
Formulation & short derivation

Representative numerical formulation

η2=∑KηK2\eta^2=\sum_K\eta_K^2∑K∈MηK2≥θη2\sum_{K\in\mathcal M}\eta_K^2\ge\theta\eta^2

Derivation / construction sketch

  1. Solve on the current mesh.
  2. Estimate local error and mark cells carrying a chosen error fraction.
  3. Refine, transfer the solution, and repeat.

Symbols & assumptions

Bulk-marking criterion shown; an indicator is not automatically a rigorous error bound.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Spatial discretization

Practical use

Resolving stress concentrations and localized heat sources.

Engineering application examples

Resolving stress concentrations and localized heat sources.

Software product or implementation route

MFEM ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Resolving stress concentrations and localized heat sources.

Method limitations

Bulk-marking criterion shown; an indicator is not automatically a rigorous error bound.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Spatial error controlDiscrete dislocation dynamics ↗

    Tracks line defects and their interactions.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlNeo-Hookean hyperelasticity ↗

    Models large elastic deformation with a strain-energy function.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlMooney–Rivlin hyperelasticity ↗

    Uses multiple strain invariants to fit rubber-like response.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlOgden hyperelasticity ↗

    Uses powers of principal stretches to represent nonlinear elasticity.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlvon Mises J2 plasticity ↗

    Uses deviatoric stress to define yielding in an isotropic ductile material.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlTresca yield model ↗

    Defines yield using maximum shear stress.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlDrucker–Prager plasticity ↗

    Uses a smooth pressure-dependent yield surface.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlMohr–Coulomb model ↗

    Relates frictional shear strength to normal stress and cohesion.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlJohnson–Cook model ↗

    Uses empirical strain, strain-rate and temperature factors.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlCrystal plasticity ↗

    Represents plastic flow through crystallographic slip systems.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlPhase-field fracture model ↗

    Represents cracks with a continuous damage-like field.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlMaxwell electromagnetic model ↗

    Couples electric and magnetic fields with charges and currents.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlEddy-current model ↗

    Models induced conducting currents in a time-varying magnetic field.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlBoltzmann kinetic equation ↗

    Evolves a particle distribution under transport and collisions.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlVlasov–Poisson model ↗

    Couples collisionless distribution dynamics to electrostatic fields.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlVlasov–Maxwell model ↗

    Couples collisionless kinetic distributions to electromagnetic fields.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlMagnetohydrodynamics (MHD) ↗

    Treats a conducting fluid coupled to a magnetic field.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlNeutron diffusion approximation ↗

    Simplifies neutron transport to a diffusion description.

    Refine localized gradients or error indicators after choosing the PDE discretization.

  • Spatial error controlGeneral relativity model ↗

    Relates spacetime curvature to matter and energy.

    Refine localized gradients or error indicators after choosing the PDE discretization.

Relationships to other techniques

Discretize → Spatial discretization

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra009

LU factorization

Solves a linear system through triangular factors with pivoting.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

PA=LUPA=LULy=PbLy=PbUx=yUx=y

Derivation / construction sketch

  1. Eliminate entries below the diagonal.
  2. Record elimination multipliers in a lower triangular factor.
  3. Apply the row permutation and solve two triangular systems.

Symbols & assumptions

Pivoting is important; sparse fill-in and ill-conditioning can dominate cost and accuracy.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Circuit matrices and repeated small-to-medium engineering solves.

Engineering application examples

Circuit matrices and repeated small-to-medium engineering solves.

Software product or implementation route

SciPy linalg.lu ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Circuit matrices and repeated small-to-medium engineering solves.

Method limitations

Pivoting is important; sparse fill-in and ill-conditioning can dominate cost and accuracy.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Linear solveLumped-capacitance thermal model ↗

    Represents a body with one spatially uniform temperature.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveThermal resistance-capacitance network ↗

    Represents heat paths and storage with connected lumped elements.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveNewton cooling model ↗

    Uses a heat-transfer coefficient between a surface and a fluid.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveSurface-to-surface radiosity model ↗

    Balances diffuse radiation exchange between surfaces.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveNeo-Hookean hyperelasticity ↗

    Models large elastic deformation with a strain-energy function.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveMooney–Rivlin hyperelasticity ↗

    Uses multiple strain invariants to fit rubber-like response.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveOgden hyperelasticity ↗

    Uses powers of principal stretches to represent nonlinear elasticity.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveLumped RLC circuit model ↗

    Uses resistors, capacitors and inductors connected by Kirchhoff laws.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveTransmission-line electrical model ↗

    Represents distributed inductance, capacitance and losses.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveDrift–diffusion semiconductor model ↗

    Combines electrostatics with carrier drift, diffusion and continuity.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveHydrodynamic carrier model ↗

    Adds carrier-energy or momentum information to transport.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveEquivalent-circuit battery model ↗

    Uses fitted electrical elements to approximate terminal behavior.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveAC power-flow model ↗

    Balances complex power on an electrical network.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveDC power-flow approximation ↗

    Linearizes active-power flow under restrictive grid assumptions.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveState-space model ↗

    Represents system evolution with internal states, inputs and outputs.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveTransfer-function model ↗

    Relates linear time-invariant input and output in the transform domain.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveBond-graph model ↗

    Represents energy exchange across mechanical, electrical and other domains.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveSystem-dynamics stock-flow model ↗

    Represents accumulated quantities and their rates of change.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveMarkov state model ↗

    Represents probabilistic transitions between a finite set of states.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveProper orthogonal decomposition (POD) ↗

    Builds a compact basis from representative field snapshots.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

  • Linear solveReduced basis model ↗

    Projects a parameterized governing model onto a small approximation space.

    For assembled nonsingular linear systems; pivoting, conditioning, and sparse fill-in determine reliability and cost.

Relationships to other techniques

Solve algebra → Linear algebra

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra010

Cholesky factorization

Factors a symmetric positive-definite matrix efficiently.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

A=LLTA=LL^{\mathsf T}Ly=bLy=bLTx=yL^{\mathsf T}x=y

Derivation / construction sketch

  1. Match entries in the product of a lower triangular factor and its transpose.
  2. Compute each positive diagonal square root.
  3. Use forward and backward substitution.

Symbols & assumptions

Requires positive definiteness; the complex analogue uses a conjugate transpose.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Elasticity systems with sufficient constraints and covariance calculations.

Engineering application examples

Elasticity systems with sufficient constraints and covariance calculations.

Software product or implementation route

SciPy linalg.cholesky ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Elasticity systems with sufficient constraints and covariance calculations.

Method limitations

Requires positive definiteness; the complex analogue uses a conjugate transpose.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Linear solveLinear elasticity (Hooke model) ↗

    Relates stress linearly to small elastic strain.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveOrthotropic elasticity ↗

    Uses direction-dependent elastic properties along material axes.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveEuler–Bernoulli beam model ↗

    Describes slender-beam bending while neglecting transverse shear deformation.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveTimoshenko beam model ↗

    Includes transverse shear deformation and rotational effects.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveKirchhoff–Love plate model ↗

    Describes thin-plate bending with normals remaining normal.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveMindlin–Reissner plate model ↗

    Includes transverse shear deformation in plate bending.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveShell model ↗

    Combines membrane and bending behavior on a curved surface.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveTruss model ↗

    Represents a structure with axial-force members joined at idealized nodes.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveCable and membrane models ↗

    Represent slender or thin structures dominated by tension.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveKalman state estimator ↗

    Combines a dynamical model with noisy observations using covariance updates.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

  • Linear solveGaussian-process surrogate ↗

    Predicts responses with a probabilistic function model fitted to samples.

    Only for symmetric positive-definite assembled systems after constraints are handled; not for general coupled saddle-point systems.

Relationships to other techniques

Solve algebra → Linear algebra

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra011

QR factorization

Uses an orthogonal factorization to solve least-squares systems.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

A=QRA=QRRx=QTbR x=Q^{\mathsf T}b

Derivation / construction sketch

  1. Apply orthogonal transformations to eliminate subdiagonal entries.
  2. Orthogonality preserves the residual norm.
  3. Solve the triangular least-squares problem.

Symbols & assumptions

Thin full-rank form shown; rank-deficient problems need pivoting or SVD.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Sensor calibration and polynomial parameter fitting.

Engineering application examples

Sensor calibration and polynomial parameter fitting.

Software product or implementation route

SciPy linalg.qr ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Sensor calibration and polynomial parameter fitting.

Method limitations

Thin full-rank form shown; rank-deficient problems need pivoting or SVD.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Stable fitting / linear solveKalman state estimator ↗

    Combines a dynamical model with noisy observations using covariance updates.

    For a linearized least-squares or calibration problem; use pivoting or SVD when rank is uncertain.

  • Stable fitting / linear solveProper orthogonal decomposition (POD) ↗

    Builds a compact basis from representative field snapshots.

    For a linearized least-squares or calibration problem; use pivoting or SVD when rank is uncertain.

  • Stable fitting / linear solveReduced basis model ↗

    Projects a parameterized governing model onto a small approximation space.

    For a linearized least-squares or calibration problem; use pivoting or SVD when rank is uncertain.

Relationships to other techniques

Solve algebra → Linear algebra

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra012

Singular value decomposition

Separates matrix directions by their amplification strengths.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

A=UΣVTA=U\Sigma V^{\mathsf T}x=VΣ+UTbx=V\Sigma^+U^{\mathsf T}b

Derivation / construction sketch

  1. Diagonalize the action of the matrix into orthogonal input and output directions.
  2. Invert retained nonzero singular values.
  3. Discard or regularize poorly determined directions if required.

Symbols & assumptions

Threshold choice affects effective rank; SVD can be expensive for large dense systems.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Inverse imaging, reduced-order models, and signal denoising.

Engineering application examples

Inverse imaging, reduced-order models, and signal denoising.

Software product or implementation route

SciPy linalg.svd ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Inverse imaging, reduced-order models, and signal denoising.

Method limitations

Threshold choice affects effective rank; SVD can be expensive for large dense systems.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Rank / inverse analysisTight-binding model ↗

    Represents electronic states with localized orbitals and hopping parameters.

    For reduced bases, rank diagnosis, or regularized inverse fitting; select truncation using the data and error budget.

  • Rank / inverse analysisHubbard model ↗

    Models competition between particle hopping and local electron interactions.

    For reduced bases, rank diagnosis, or regularized inverse fitting; select truncation using the data and error budget.

  • Rank / inverse analysisProper orthogonal decomposition (POD) ↗

    Builds a compact basis from representative field snapshots.

    For reduced bases, rank diagnosis, or regularized inverse fitting; select truncation using the data and error budget.

  • Rank / inverse analysisReduced basis model ↗

    Projects a parameterized governing model onto a small approximation space.

    For reduced bases, rank diagnosis, or regularized inverse fitting; select truncation using the data and error budget.

  • Rank / inverse analysisTight-binding electronic model ↗

    Builds crystal electronic bands from localized orbitals and intersite hopping.

    Use SVD on H(k)-E I to check null states, condition analysis to assess sensitivity, and quadrature for Brillouin-zone averages. General bands require a Hermitian eigensolver, not SVD singular values interpreted as signed energies. For reduced bases, rank diagnosis, or regularized inverse fitting; select truncation using the data and error budget.

  • Rank / inverse analysisNearly-free-electron model ↗

    Predicts band gaps by perturbing free electrons with a weak periodic potential.

    Use SVD for null-state diagnostics, sensitivity checks for small gaps, and quadrature for band averages. Diagonalize the Hermitian plane-wave Hamiltonian for actual energies. For reduced bases, rank diagnosis, or regularized inverse fitting; select truncation using the data and error budget.

  • Rank / inverse analysisHarmonic lattice dynamics ↗

    Computes phonon modes from a quadratic expansion of crystal potential energy.

    Fit force constants with SVD, cross-check real-time harmonic motion with Verlet, and test displacement/time-step refinement. Obtain phonon frequencies with a Hermitian dynamical-matrix eigensolver. For reduced bases, rank diagnosis, or regularized inverse fitting; select truncation using the data and error budget.

Relationships to other techniques

Solve algebra → Linear algebra

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra013

Jacobi iteration

Updates each unknown from the previous iterate using the diagonal.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

A=D+L+UA=D+L+Uxk+1=D−1[b−(L+U)xk]x^{k+1}=D^{-1}[b-(L+U)x^k]

Derivation / construction sketch

  1. Split the matrix into diagonal and off-diagonal parts.
  2. Move the off-diagonal contribution to the right-hand side.
  3. Evaluate all new components from the old iterate.

Symbols & assumptions

Convergence requires the iteration matrix spectral radius below one; diagonal dominance is a sufficient condition.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Parallel smoothing and simple stationary iterative solvers.

Engineering application examples

Parallel smoothing and simple stationary iterative solvers.

Software product or implementation route

PETSc PCJACOBI ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Parallel smoothing and simple stationary iterative solvers.

Method limitations

Convergence requires the iteration matrix spectral radius below one; diagonal dominance is a sufficient condition.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Relaxation / preconditioningFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    Use as a diagonal preconditioner or suitable smoother; standalone iteration requires a convergence check.

Relationships to other techniques

Solve algebra → Linear algebra

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra014

Gauss-Seidel iteration

Uses newly updated values immediately within each sweep.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

(D+L)xk+1=b−Uxk(D+L)x^{k+1}=b-Ux^k

Derivation / construction sketch

  1. Split the matrix into lower triangular and upper parts.
  2. Sweep through the unknowns in a fixed order.
  3. Reuse new components as soon as they are available.

Symbols & assumptions

Ordering influences convergence and parallelism; convergence is not guaranteed for arbitrary matrices.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Smoothers inside multigrid and structured-grid elliptic solvers.

Engineering application examples

Smoothers inside multigrid and structured-grid elliptic solvers.

Software product or implementation route

PETSc PCSOR ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Smoothers inside multigrid and structured-grid elliptic solvers.

Method limitations

Ordering influences convergence and parallelism; convergence is not guaranteed for arbitrary matrices.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Relaxation / smoothingFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    Useful for suitable diffusion-like matrices or multigrid smoothing; ordering and parallelism matter.

Relationships to other techniques

Solve algebra → Linear algebra

Specific connections

  • Generalized by relaxation in Successive over-relaxation

    Setting omega to one recovers the Gauss-Seidel iteration.

  • Can smooth within Geometric multigrid

    Relaxation damps components of the error before and after coarse correction.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra015

Successive over-relaxation

Relaxes a Gauss-Seidel correction with a tunable weight.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

xk+1=(D+ωL)−1[ωb−((ω−1)D+ωU)xk]x^{k+1}=(D+\omega L)^{-1}[\omega b-((\omega-1)D+\omega U)x^k]

Derivation / construction sketch

  1. Form the Gauss-Seidel component update.
  2. Blend its correction using relaxation weight omega.
  3. Choose the weight to improve convergence for the problem class.

Symbols & assumptions

For symmetric positive-definite systems, 0 < omega < 2 gives convergence; an optimal weight is problem dependent.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Elliptic potential and diffusion calculations.

Engineering application examples

Elliptic potential and diffusion calculations.

Software product or implementation route

PETSc PCSOR ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Elliptic potential and diffusion calculations.

Method limitations

For symmetric positive-definite systems, 0 < omega < 2 gives convergence; an optimal weight is problem dependent.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • RelaxationFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    For suitable elliptic systems; choose relaxation carefully and verify convergence rather than assuming acceleration.

Relationships to other techniques

Solve algebra → Linear algebra

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra016

Conjugate gradient

Solves symmetric positive-definite systems using conjugate search directions.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

αk=rkTrkpkTApk\alpha_k=\frac{r_k^{\mathsf T}r_k}{p_k^{\mathsf T}Ap_k}xk+1=xk+αkpkx_{k+1}=x_k+\alpha_kp_krk+1=rk−αkApkr_{k+1}=r_k-\alpha_kAp_k

Derivation / construction sketch

  1. Minimize the quadratic energy along a search direction.
  2. Choose successive directions to be A-conjugate.
  3. Update residuals and stop using a scaled tolerance.

Symbols & assumptions

Unpreconditioned step shown; A must be symmetric positive definite, with compatible preconditioning.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Large constrained elasticity and diffusion systems.

Engineering application examples

Large constrained elasticity and diffusion systems.

Software product or implementation route

PETSc KSPCG ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Large constrained elasticity and diffusion systems.

Method limitations

Unpreconditioned step shown; A must be symmetric positive definite, with compatible preconditioning.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Linear solvePotential-flow model ↗

    Represents irrotational velocity using a scalar potential.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveTransient heat equation ↗

    Balances thermal storage, conduction and heat sources.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveStefan phase-change problem ↗

    Couples heat transport to a moving melting or freezing boundary.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveEnthalpy–porosity model ↗

    Represents melting using enthalpy and a porous resistance in the mushy zone.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveLinear elasticity (Hooke model) ↗

    Relates stress linearly to small elastic strain.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveOrthotropic elasticity ↗

    Uses direction-dependent elastic properties along material axes.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveEuler–Bernoulli beam model ↗

    Describes slender-beam bending while neglecting transverse shear deformation.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveTimoshenko beam model ↗

    Includes transverse shear deformation and rotational effects.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveKirchhoff–Love plate model ↗

    Describes thin-plate bending with normals remaining normal.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveMindlin–Reissner plate model ↗

    Includes transverse shear deformation in plate bending.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveShell model ↗

    Combines membrane and bending behavior on a curved surface.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveTruss model ↗

    Represents a structure with axial-force members joined at idealized nodes.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveCable and membrane models ↗

    Represent slender or thin structures dominated by tension.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveDC power-flow approximation ↗

    Linearizes active-power flow under restrictive grid assumptions.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveHomogenization ↗

    Derives effective properties or equations from smaller-scale structure.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveRepresentative volume element (RVE) ↗

    Uses a finite microstructural sample to estimate bulk response.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

  • Linear solveFE² computational homogenization ↗

    Solves microscale problems within a macroscale finite-element calculation.

    Only when both the matrix and preconditioner satisfy the required symmetry and positive-definiteness conditions.

Relationships to other techniques

Solve algebra → Linear algebra

Specific connections

  • Can be preconditioned by Algebraic multigrid

    The multigrid cycle must preserve the symmetry and positivity required by CG.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra017

GMRES

Minimizes the residual over a Krylov subspace for nonsymmetric systems.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

xm=x0+Vmymx_m=x_0+V_my_mym=arg⁡min⁡y∥βe1−Hˉmy∥2y_m=\arg\min_y\lVert\beta e_1-\bar H_my\rVert_2

Derivation / construction sketch

  1. Build an orthonormal Krylov basis with Arnoldi iteration.
  2. Represent the matrix action by an upper Hessenberg matrix.
  3. Solve the small residual-minimization problem.

Symbols & assumptions

Restarting limits storage but may slow or stall convergence; preconditioning is often essential.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Advection-diffusion and coupled multiphysics linearizations.

Engineering application examples

Advection-diffusion and coupled multiphysics linearizations.

Software product or implementation route

PETSc KSPGMRES ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Advection-diffusion and coupled multiphysics linearizations.

Method limitations

Restarting limits storage but may slow or stall convergence; preconditioning is often essential.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Linear solveNavier–Stokes model ↗

    Conserves mass and momentum for a viscous continuum fluid.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveEuler flow model ↗

    Neglects viscous stresses in compressible or incompressible flow.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveStokes creeping-flow model ↗

    Neglects inertial terms relative to viscosity.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveBoundary-layer model ↗

    Resolves thin near-wall regions with scale-based simplifications.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveLubrication approximation ↗

    Simplifies viscous flow in thin gaps.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveOldroyd-B model ↗

    Combines solvent viscosity with an elastic polymer stress.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveReynolds-averaged Navier–Stokes (RANS) ↗

    Models mean flow with closure for unresolved turbulent stresses.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveSpalart–Allmaras model ↗

    Uses a transported turbulence variable to obtain eddy viscosity.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solvek–epsilon model ↗

    Uses turbulent kinetic energy and dissipation rate to close mean flow.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solvek–omega model ↗

    Uses turbulent kinetic energy and specific dissipation rate.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveSST k–omega model ↗

    Blends near-wall and outer-flow behavior with a shear-stress limiter.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveReynolds-stress transport model ↗

    Transports individual turbulent stress components.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveLarge-eddy simulation (LES) ↗

    Resolves larger turbulent motions and models subgrid effects.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveSmagorinsky subgrid model ↗

    Relates subgrid eddy viscosity to resolved strain and filter scale.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveDetached-eddy simulation (DES) ↗

    Combines RANS near walls with LES-like treatment away from them.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveVolume-of-fluid (VOF) representation ↗

    Tracks phase volume fractions to represent an interface.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveEuler–Euler two-fluid model ↗

    Treats phases as interpenetrating continua with exchange terms.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveLagrangian particle tracking ↗

    Tracks discrete particles through a carrier flow.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveSurface-to-surface radiosity model ↗

    Balances diffuse radiation exchange between surfaces.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveMaxwell electromagnetic model ↗

    Couples electric and magnetic fields with charges and currents.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

  • Linear solveEddy-current model ↗

    Models induced conducting currents in a time-varying magnetic field.

    For nonsymmetric linearized or discretized systems; plan preconditioning and restart/storage settings.

Relationships to other techniques

Solve algebra → Linear algebra

Specific connections

  • Can solve inner systems in Newton-Krylov method

    Matrix-free Jacobian-vector products can drive a preconditioned GMRES iteration.

  • Can be preconditioned by Algebraic multigrid

    A fixed suitable multigrid operator can accelerate a GMRES solve; variable operators call for flexible variants.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra018

BiCGSTAB

Uses short recurrences to stabilize a nonsymmetric Krylov iteration.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

rk=pk(A)r0r_k=p_k(A)r_0pk(0)=1p_k(0)=1

Derivation / construction sketch

  1. Construct a bi-conjugate residual polynomial.
  2. Combine it with local residual-smoothing factors.
  3. Use short vector recurrences to avoid storing a full Krylov basis.

Symbols & assumptions

Polynomial viewpoint only; breakdown and irregular convergence are possible.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Memory-limited nonsymmetric CFD matrix solves.

Engineering application examples

Memory-limited nonsymmetric CFD matrix solves.

Software product or implementation route

PETSc KSPBCGS ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Memory-limited nonsymmetric CFD matrix solves.

Method limitations

Polynomial viewpoint only; breakdown and irregular convergence are possible.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Linear solveNavier–Stokes model ↗

    Conserves mass and momentum for a viscous continuum fluid.

    A short-recurrence option for nonsymmetric systems; monitor breakdown and irregular residual convergence.

Relationships to other techniques

Solve algebra → Linear algebra

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra019

Geometric multigrid

Removes error at multiple mesh resolutions.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

rh=bh−Ahxhr_h=b_h-A_hx_hAHeH=RrhA_He_H=Rr_hxh←xh+PeHx_h\leftarrow x_h+Pe_H

Derivation / construction sketch

  1. Smooth high-frequency error on the fine grid.
  2. Restrict the residual and solve for coarse-grid error.
  3. Prolong the correction and apply additional smoothing.

Symbols & assumptions

Transfers, coarse operators, boundary conditions, and smoothers must work together.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Large Poisson and diffusion problems on mesh hierarchies.

Engineering application examples

Large Poisson and diffusion problems on mesh hierarchies.

Software product or implementation route

PETSc PCMG ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Large Poisson and diffusion problems on mesh hierarchies.

Method limitations

Transfers, coarse operators, boundary conditions, and smoothers must work together.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Linear accelerationFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    For suitable elliptic operators with a mesh hierarchy and compatible transfer operators and smoothers.

  • Linear accelerationElectrostatic Poisson model ↗

    Relates electric potential to charge density.

    For suitable elliptic operators with a mesh hierarchy and compatible transfer operators and smoothers.

  • Linear accelerationMagnetostatic model ↗

    Represents steady magnetic fields driven by currents and magnetization.

    For suitable elliptic operators with a mesh hierarchy and compatible transfer operators and smoothers.

Relationships to other techniques

Solve algebra → Linear algebra

Specific connections

  • Can use as smoother Gauss-Seidel iteration

    Relaxation damps components of the error before and after coarse correction.

  • Has algebraic alternative Algebraic multigrid

    Both combine smoothing and coarse corrections; the hierarchy is constructed differently.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Linear algebra020

Algebraic multigrid

Constructs coarse spaces from matrix structure rather than an explicit mesh hierarchy.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

Ac=RAPA_c=RAPx←x+PAc−1R(b−Ax)x\leftarrow x+PA_c^{-1}R(b-Ax)

Derivation / construction sketch

  1. Identify algebraically strong connections.
  2. Construct interpolation and a coarse operator.
  3. Combine coarse corrections with relaxation.

Symbols & assumptions

One correction shown; performance depends on operator structure and coarsening choices.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Linear algebra

Practical use

Large sparse elliptic systems on complex engineering meshes.

Engineering application examples

Large sparse elliptic systems on complex engineering meshes.

Software product or implementation route

PETSc PCGAMG ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Large sparse elliptic systems on complex engineering meshes.

Method limitations

One correction shown; performance depends on operator structure and coarsening choices.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Linear accelerationDarcy porous-flow model ↗

    Relates averaged fluid flux to hydraulic gradient.

    For suitable sparse elliptic blocks; coupled, indefinite, or strongly anisotropic operators need tailored treatment.

  • Linear accelerationBrinkman porous-flow model ↗

    Adds a viscous shear term to a Darcy-like resistance model.

    For suitable sparse elliptic blocks; coupled, indefinite, or strongly anisotropic operators need tailored treatment.

  • Linear accelerationRichards equation ↗

    Describes variably saturated water movement in porous media.

    For suitable sparse elliptic blocks; coupled, indefinite, or strongly anisotropic operators need tailored treatment.

  • Linear accelerationBiot poroelasticity ↗

    Couples solid deformation and pore-fluid pressure.

    For suitable sparse elliptic blocks; coupled, indefinite, or strongly anisotropic operators need tailored treatment.

  • Linear accelerationTerzaghi consolidation model ↗

    Describes time-dependent settlement from pore-pressure dissipation.

    For suitable sparse elliptic blocks; coupled, indefinite, or strongly anisotropic operators need tailored treatment.

  • Linear accelerationModified Cam-Clay model ↗

    Uses critical-state plasticity for idealized clay behavior.

    For suitable sparse elliptic blocks; coupled, indefinite, or strongly anisotropic operators need tailored treatment.

Relationships to other techniques

Solve algebra → Linear algebra

Specific connections

  • Can precondition GMRES

    A fixed suitable multigrid operator can accelerate a GMRES solve; variable operators call for flexible variants.

  • Can precondition Conjugate gradient

    The multigrid cycle must preserve the symmetry and positivity required by CG.

  • Algebraic alternative to Geometric multigrid

    Both combine smoothing and coarse corrections; the hierarchy is constructed differently.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Nonlinear equations021

Bisection

Reliably narrows a continuous scalar root bracket.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

ck=ak+bk2c_k=\frac{a_k+b_k}{2}bk−ak=b0−a02kb_k-a_k=\frac{b_0-a_0}{2^k}

Derivation / construction sketch

  1. Begin with opposite endpoint signs.
  2. Evaluate the midpoint.
  3. Retain the half interval that preserves a sign change.

Symbols & assumptions

Continuity and a valid bracket are required; even-multiplicity roots may not change sign.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Nonlinear equations

Practical use

Finding operating points from monotone balance equations.

Engineering application examples

Finding operating points from monotone balance equations.

Software product or implementation route

SciPy optimize.bisect ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Finding operating points from monotone balance equations.

Method limitations

Continuity and a valid bracket are required; even-multiplicity roots may not change sign.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Scalar rootNRTL activity model ↗

    Uses local-composition parameters to describe nonideal liquid mixtures.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootUNIQUAC activity model ↗

    Combines molecular size, shape and interaction contributions.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootDebye–Hückel model ↗

    Approximates ionic activity using screened electrostatic interactions.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootArrhenius rate model ↗

    Relates a rate coefficient to temperature through an activation energy.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootTransition-state theory ↗

    Estimates reaction rates from a free-energy barrier.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootMichaelis–Menten kinetics ↗

    Approximates enzyme reaction rates with substrate saturation.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootLangmuir adsorption isotherm ↗

    Models adsorption on equivalent sites with finite occupancy.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootLangmuir–Hinshelwood kinetics ↗

    Models surface reactions involving adsorbed reactants.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootHagen–Poiseuille model ↗

    Predicts fully developed laminar flow in a circular pipe.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootDarcy–Weisbach model ↗

    Relates pipe pressure loss to friction factor and flow speed.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootNon-Newtonian power-law fluid ↗

    Relates shear stress to a power of shear rate.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootBingham plastic model ↗

    Represents a material with a yield stress and post-yield viscosity.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootHerschel–Bulkley model ↗

    Combines yield stress with nonlinear post-yield flow.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootHybrid dynamical model ↗

    Combines continuous dynamics with discrete state changes.

    For a continuous scalar closure or balance with a known sign-changing bracket.

  • Scalar rootCarnahan-Starling hard-sphere equation of state ↗

    Approximates the compressibility factor of a monodisperse hard-sphere fluid from its packing fraction.

    A robust bracketed alternative for EOS inversion; restrict the bracket to the physically intended fluid branch.

Relationships to other techniques

Solve algebra → Nonlinear equations

Specific connections

  • Provides safeguard for Brent root finding

    Bracketing is retained while interpolation proposes faster steps.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Nonlinear equations022

Newton-Raphson method

Solves nonlinear equations by repeated local linearization.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

J(xk)sk=−F(xk)J(x_k)s_k=-F(x_k)xk+1=xk+skx_{k+1}=x_k+s_k

Derivation / construction sketch

  1. Expand the residual to first order at the current state.
  2. Set the linearized residual to zero.
  3. Solve for a correction and repeat, with damping if needed.

Symbols & assumptions

Fast local convergence requires a suitable initial guess, smoothness, and a nonsingular Jacobian.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Nonlinear equations

Practical use

Nonlinear material equilibria and implicit time-step equations.

Engineering application examples

Nonlinear material equilibria and implicit time-step equations.

Software product or implementation route

SciPy optimize.newton ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Nonlinear material equilibria and implicit time-step equations.

Method limitations

Fast local convergence requires a suitable initial guess, smoothness, and a nonsingular Jacobian.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Nonlinear solveIdeal gas equation of state ↗

    Relates pressure, volume and temperature for a dilute noninteracting gas.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveVan der Waals equation of state ↗

    Adds molecular attraction and excluded volume to an ideal gas model.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solvePeng–Robinson equation of state ↗

    Uses a cubic equation of state for real-fluid behavior.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveSoave–Redlich–Kwong equation of state ↗

    Uses a temperature-dependent attraction correction in a cubic fluid model.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveVirial equation of state ↗

    Represents nonideal behavior as a density or pressure expansion.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveGibbs-energy minimization ↗

    Finds equilibrium by minimizing free energy under conservation constraints.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveCALPHAD model ↗

    Combines assessed phase free energies to predict equilibria.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveMass-action reaction kinetics ↗

    Relates reaction rates to species concentrations and reaction orders.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveContinuous stirred-tank reactor (CSTR) ↗

    Assumes a well-mixed reactor with inlet and outlet flows.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solvePlug-flow reactor (PFR) ↗

    Approximates axial evolution without axial back-mixing.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveBatch reactor model ↗

    Evolves composition and energy in a closed reacting charge.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solvevon Mises J2 plasticity ↗

    Uses deviatoric stress to define yielding in an isotropic ductile material.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveTresca yield model ↗

    Defines yield using maximum shear stress.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveDrucker–Prager plasticity ↗

    Uses a smooth pressure-dependent yield surface.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveMohr–Coulomb model ↗

    Relates frictional shear strength to normal stress and cohesion.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveJohnson–Cook model ↗

    Uses empirical strain, strain-rate and temperature factors.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveCrystal plasticity ↗

    Represents plastic flow through crystallographic slip systems.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solvePhase-field fracture model ↗

    Represents cracks with a continuous damage-like field.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveLumped RLC circuit model ↗

    Uses resistors, capacitors and inductors connected by Kirchhoff laws.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveTransmission-line electrical model ↗

    Represents distributed inductance, capacitance and losses.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveShockley diode model ↗

    Approximates diode current with an exponential voltage relation.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveEbers–Moll transistor model ↗

    Models coupled junction currents in a bipolar transistor.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveMOSFET square-law model ↗

    Approximates long-channel transistor current from terminal voltages.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveBSIM compact-model family ↗

    Uses detailed parameterized MOS transistor relations.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveDrift–diffusion semiconductor model ↗

    Combines electrostatics with carrier drift, diffusion and continuity.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveHydrodynamic carrier model ↗

    Adds carrier-energy or momentum information to transport.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveNernst equilibrium potential ↗

    Relates electrochemical equilibrium potential to species activities.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveButler–Volmer kinetics ↗

    Relates interfacial current to electrochemical overpotential.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveTafel approximation ↗

    Approximates high-overpotential behavior of Butler–Volmer kinetics.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveSix-degree-of-freedom flight model ↗

    Evolves vehicle translation and rotation using aerodynamic and propulsion forces.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveBicycle vehicle model ↗

    Combines left and right wheels into a planar steering model.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveQuarter-car suspension model ↗

    Represents one wheel assembly and a fraction of vehicle body mass.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveAC power-flow model ↗

    Balances complex power on an electrical network.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveSwing-equation generator model ↗

    Represents rotor-angle dynamics from mechanical-electrical power imbalance.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveStellar structure model ↗

    Couples hydrostatic balance, energy transport and energy generation.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

  • Nonlinear solveFLRW cosmological model ↗

    Assumes a homogeneous and isotropic expanding spacetime.

    For differentiable residuals with a suitable initial guess; use globalization and consistent Jacobians.

Relationships to other techniques

Solve algebra → Nonlinear equations

Specific connections

  • Can solve implicit step in Backward Euler

    An implicit step is a nonlinear equation unless the dynamics are linear.

  • Has quasi-Newton alternative Broyden method

    Rank-one secant updates reduce repeated Jacobian evaluation.

  • Can use Krylov implementation Newton-Krylov method

    A Krylov method approximately solves each linearized Newton correction.

  • Has derivative-free scalar variant Secant method

    A divided difference replaces the scalar derivative.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Nonlinear equations023

Secant method

Approximates a scalar derivative from two previous points.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

xk+1=xk−f(xk)xk−xk−1f(xk)−f(xk−1)x_{k+1}=x_k-f(x_k)\frac{x_k-x_{k-1}}{f(x_k)-f(x_{k-1})}

Derivation / construction sketch

  1. Replace Newton's derivative with a divided difference.
  2. Intersect the resulting secant line with the horizontal axis.
  3. Advance the pair of iterates.

Symbols & assumptions

Does not preserve a root bracket; small denominator differences can cause large steps.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Nonlinear equations

Practical use

Scalar balance equations with expensive derivatives.

Engineering application examples

Scalar balance equations with expensive derivatives.

Software product or implementation route

SciPy optimize.newton ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Scalar balance equations with expensive derivatives.

Method limitations

Does not preserve a root bracket; small denominator differences can cause large steps.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Scalar rootStefan–Boltzmann surface model ↗

    Relates idealized surface radiant emission to the fourth power of temperature.

    For smooth scalar equations when derivatives are costly; it does not preserve a bracket and can fail.

Relationships to other techniques

Solve algebra → Nonlinear equations

Specific connections

  • Provides interpolation step for Brent root finding

    Secant and inverse-quadratic proposals are accepted only when suitable.

  • Derivative approximation to Newton-Raphson method

    A divided difference replaces the scalar derivative.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Nonlinear equations024

Brent root finding

Combines bracket reliability with interpolation-based acceleration.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

f(ak)f(bk)≤0f(a_k)f(b_k)\le0ck∈(ak,bk)c_k\in(a_k,b_k)

Derivation / construction sketch

  1. Maintain a sign-changing bracket.
  2. Try secant or inverse-quadratic interpolation when its step is acceptable.
  3. Fall back to bisection when interpolation is unsafe.

Symbols & assumptions

Invariant shown rather than a full algorithm; the function must be continuous on a valid bracket.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Nonlinear equations

Practical use

Robust engineering threshold and equilibrium calculations.

Engineering application examples

Robust engineering threshold and equilibrium calculations.

Software product or implementation route

SciPy optimize.brentq ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Robust engineering threshold and equilibrium calculations.

Method limitations

Invariant shown rather than a full algorithm; the function must be continuous on a valid bracket.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Scalar rootIdeal gas equation of state ↗

    Relates pressure, volume and temperature for a dilute noninteracting gas.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootVan der Waals equation of state ↗

    Adds molecular attraction and excluded volume to an ideal gas model.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootPeng–Robinson equation of state ↗

    Uses a cubic equation of state for real-fluid behavior.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootSoave–Redlich–Kwong equation of state ↗

    Uses a temperature-dependent attraction correction in a cubic fluid model.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootVirial equation of state ↗

    Represents nonideal behavior as a density or pressure expansion.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootNRTL activity model ↗

    Uses local-composition parameters to describe nonideal liquid mixtures.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootUNIQUAC activity model ↗

    Combines molecular size, shape and interaction contributions.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootDebye–Hückel model ↗

    Approximates ionic activity using screened electrostatic interactions.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootArrhenius rate model ↗

    Relates a rate coefficient to temperature through an activation energy.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootTransition-state theory ↗

    Estimates reaction rates from a free-energy barrier.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootMichaelis–Menten kinetics ↗

    Approximates enzyme reaction rates with substrate saturation.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootLangmuir adsorption isotherm ↗

    Models adsorption on equivalent sites with finite occupancy.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootLangmuir–Hinshelwood kinetics ↗

    Models surface reactions involving adsorbed reactants.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootHagen–Poiseuille model ↗

    Predicts fully developed laminar flow in a circular pipe.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootDarcy–Weisbach model ↗

    Relates pipe pressure loss to friction factor and flow speed.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootNon-Newtonian power-law fluid ↗

    Relates shear stress to a power of shear rate.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootBingham plastic model ↗

    Represents a material with a yield stress and post-yield viscosity.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootHerschel–Bulkley model ↗

    Combines yield stress with nonlinear post-yield flow.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootStefan–Boltzmann surface model ↗

    Relates idealized surface radiant emission to the fourth power of temperature.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootNorton creep law ↗

    Relates creep rate to a power of stress.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootLinear elastic fracture mechanics (LEFM) ↗

    Uses crack-tip intensity parameters in an elastic body.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootCohesive-zone model ↗

    Uses traction-separation relations across a fracture process zone.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootParis fatigue crack-growth law ↗

    Relates cyclic crack-growth rate to stress-intensity-factor range.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootMiner cumulative damage rule ↗

    Adds fractions of fatigue life consumed by load cycles.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootArchard wear model ↗

    Relates wear volume to load, sliding distance and hardness.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootMagnetic-circuit model ↗

    Uses reluctance and magnetomotive force in lumped magnetic paths.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootJiles–Atherton hysteresis model ↗

    Represents path-dependent magnetization with phenomenological parameters.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootGaussian beam model ↗

    Represents a paraxial beam with a Gaussian transverse profile.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootDrude–Lorentz optical model ↗

    Represents free-carrier and bound-charge contributions to permittivity.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootShockley diode model ↗

    Approximates diode current with an exponential voltage relation.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootEbers–Moll transistor model ↗

    Models coupled junction currents in a bipolar transistor.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootMOSFET square-law model ↗

    Approximates long-channel transistor current from terminal voltages.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootBSIM compact-model family ↗

    Uses detailed parameterized MOS transistor relations.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootNernst equilibrium potential ↗

    Relates electrochemical equilibrium potential to species activities.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootButler–Volmer kinetics ↗

    Relates interfacial current to electrochemical overpotential.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootTafel approximation ↗

    Approximates high-overpotential behavior of Butler–Volmer kinetics.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootForchheimer model ↗

    Adds inertial resistance to porous flow.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootvan Genuchten retention model ↗

    Relates water saturation to pressure head with fitted parameters.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootLifting-line model ↗

    Approximates finite-wing lift using a spanwise circulation distribution.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootBlade-element momentum model ↗

    Combines blade-section loads with momentum balances.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootHybrid dynamical model ↗

    Combines continuous dynamics with discrete state changes.

    For continuous scalar closures with a valid sign-changing bracket; verify the intended physical root.

  • Scalar rootCarnahan-Starling hard-sphere equation of state ↗

    Approximates the compressibility factor of a monodisperse hard-sphere fluid from its packing fraction.

    Invert the hard-sphere EOS for packing fraction at a specified pressure using a bracket within the fluid regime. Direct pressure evaluation requires no root solver.

Relationships to other techniques

Solve algebra → Nonlinear equations

Specific connections

  • Can use interpolation from Secant method

    Secant and inverse-quadratic proposals are accepted only when suitable.

  • Uses safeguarding from Bisection

    Bracketing is retained while interpolation proposes faster steps.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Nonlinear equations025

Fixed-point iteration

Iterates a rearranged equation until the state stops changing.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

xk+1=g(xk)x_{k+1}=g(x_k)∥g(x)−g(y)∥≤q∥x−y∥,q<1\lVert g(x)-g(y)\rVert\le q\lVert x-y\rVert,\quad q<1

Derivation / construction sketch

  1. Rewrite the equation as x = g(x).
  2. Use the latest state on the right-hand side.
  3. A contraction on an invariant complete set guarantees convergence.

Symbols & assumptions

A poor rearrangement can diverge even when the original equation has a root.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Nonlinear equations

Practical use

Partitioned coupling and simple nonlinear balance solvers.

Engineering application examples

Partitioned coupling and simple nonlinear balance solvers.

Software product or implementation route

SciPy optimize.fixed_point ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Partitioned coupling and simple nonlinear balance solvers.

Method limitations

A poor rearrangement can diverge even when the original equation has a root.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Self-consistency / couplingHartree–Fock model ↗

    Approximates a many-electron wavefunction by one self-consistent Slater determinant.

    For a contractive or suitably relaxed fixed-point formulation; monitor residuals and possible divergence.

  • Self-consistency / couplingDensity functional theory (DFT) ↗

    Uses electron density to determine ground-state properties with an approximate exchange-correlation functional.

    For a contractive or suitably relaxed fixed-point formulation; monitor residuals and possible divergence.

  • Self-consistency / couplingLifting-line model ↗

    Approximates finite-wing lift using a spanwise circulation distribution.

    For a contractive or suitably relaxed fixed-point formulation; monitor residuals and possible divergence.

  • Self-consistency / couplingBlade-element momentum model ↗

    Combines blade-section loads with momentum balances.

    For a contractive or suitably relaxed fixed-point formulation; monitor residuals and possible divergence.

  • Self-consistency / couplingOrnstein-Zernike equation ↗

    Relates total and direct pair correlations in a homogeneous liquid, linking microscopic structure to scattering.

    Iterate the coupled OZ/closure equations with damping; check residuals and grid convergence, especially at high density.

  • Self-consistency / couplingPercus-Yevick closure ↗

    Closes the liquid integral equation using an approximate relation between pair correlations and interactions.

    Iterate the coupled OZ/closure equations with damping; check residuals and grid convergence, especially at high density.

  • Self-consistency / couplingHypernetted-chain (HNC) closure ↗

    Approximates liquid pair structure by neglecting bridge diagrams in the exact closure.

    Iterate the coupled OZ/closure equations with damping; check residuals and grid convergence, especially at high density.

Relationships to other techniques

Solve algebra → Nonlinear equations

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Nonlinear equations026

Broyden method

Updates an approximate Jacobian from observed changes.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

Bk+1=Bk+(yk−Bksk)skTskTskB_{k+1}=B_k+\frac{(y_k-B_ks_k)s_k^{\mathsf T}}{s_k^{\mathsf T}s_k}yk=F(xk+1)−F(xk)y_k=F(x_{k+1})-F(x_k)

Derivation / construction sketch

  1. Take a step using the current Jacobian approximation.
  2. Measure the residual change.
  3. Apply a rank-one correction satisfying the new secant condition.

Symbols & assumptions

Good-Broyden Jacobian form shown; scaling, safeguards, and initial approximation matter.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Nonlinear equations

Practical use

Nonlinear systems when full Jacobians are costly.

Engineering application examples

Nonlinear systems when full Jacobians are costly.

Software product or implementation route

SciPy optimize.broyden1 ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Nonlinear systems when full Jacobians are costly.

Method limitations

Good-Broyden Jacobian form shown; scaling, safeguards, and initial approximation matter.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Nonlinear solveHartree–Fock model ↗

    Approximates a many-electron wavefunction by one self-consistent Slater determinant.

    For smooth nonlinear residuals when repeated full Jacobians are expensive; scale and safeguard the iteration.

  • Nonlinear solveDensity functional theory (DFT) ↗

    Uses electron density to determine ground-state properties with an approximate exchange-correlation functional.

    For smooth nonlinear residuals when repeated full Jacobians are expensive; scale and safeguard the iteration.

  • Nonlinear solveOrnstein-Zernike equation ↗

    Relates total and direct pair correlations in a homogeneous liquid, linking microscopic structure to scattering.

    Solve the discretized OZ/closure residual with quasi-Newton mixing when simple iteration is slow; enforce core conditions and inspect convergence.

  • Nonlinear solvePercus-Yevick closure ↗

    Closes the liquid integral equation using an approximate relation between pair correlations and interactions.

    Solve the discretized OZ/closure residual with quasi-Newton mixing when simple iteration is slow; enforce core conditions and inspect convergence.

  • Nonlinear solveHypernetted-chain (HNC) closure ↗

    Approximates liquid pair structure by neglecting bridge diagrams in the exact closure.

    Solve the discretized OZ/closure residual with quasi-Newton mixing when simple iteration is slow; enforce core conditions and inspect convergence.

Relationships to other techniques

Solve algebra → Nonlinear equations

Specific connections

  • Quasi-Newton alternative to Newton-Raphson method

    Rank-one secant updates reduce repeated Jacobian evaluation.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Nonlinear equations027

Newton-Krylov method

Solves each Newton correction approximately with a Krylov method.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

J(x)v≈F(x+εv)−F(x)εJ(x)v\approx\frac{F(x+\varepsilon v)-F(x)}{\varepsilon}Js=−F(x)Js=-F(x)

Derivation / construction sketch

  1. Linearize the nonlinear residual.
  2. Supply Jacobian-vector products without necessarily forming a matrix.
  3. Use a preconditioned Krylov solve and globalize the Newton step.

Symbols & assumptions

Finite-difference perturbations balance truncation and roundoff; preconditioning remains critical.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Nonlinear equations

Practical use

Large nonlinear PDE systems and implicit multiphysics.

Engineering application examples

Large nonlinear PDE systems and implicit multiphysics.

Software product or implementation route

SciPy optimize.newton_krylov ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Large nonlinear PDE systems and implicit multiphysics.

Method limitations

Finite-difference perturbations balance truncation and roundoff; preconditioning remains critical.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Large nonlinear solvePoisson–Nernst–Planck model ↗

    Couples electrostatics to diffusion and migration of ions.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveDoyle–Fuller–Newman (DFN/P2D) model ↗

    Combines porous-electrode transport and particle diffusion.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveSingle-particle battery model (SPM) ↗

    Represents each electrode by a representative active-material particle.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveSingle-particle model with electrolyte (SPMe) ↗

    Adds electrolyte concentration effects to a single-particle approximation.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveDarcy porous-flow model ↗

    Relates averaged fluid flux to hydraulic gradient.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveBrinkman porous-flow model ↗

    Adds a viscous shear term to a Darcy-like resistance model.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveRichards equation ↗

    Describes variably saturated water movement in porous media.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveBiot poroelasticity ↗

    Couples solid deformation and pore-fluid pressure.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveTerzaghi consolidation model ↗

    Describes time-dependent settlement from pore-pressure dissipation.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveModified Cam-Clay model ↗

    Uses critical-state plasticity for idealized clay behavior.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveFluid–structure interaction (FSI) ↗

    Couples fluid loads with structural motion or deformation.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

  • Large nonlinear solveThermomechanical coupling ↗

    Couples temperature evolution and mechanical response.

    For large smooth residual systems; matrix-free products still need effective preconditioning and globalization.

Relationships to other techniques

Solve algebra → Nonlinear equations

Specific connections

  • Can use inner linear solver GMRES

    Matrix-free Jacobian-vector products can drive a preconditioned GMRES iteration.

  • Krylov implementation of Newton-Raphson method

    A Krylov method approximately solves each linearized Newton correction.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Nonlinear equations028

Pseudo-arclength continuation

Tracks solution branches through turning points by augmenting the nonlinear system.

Solve algebraNumerical technique
Formulation & short derivation

Representative numerical formulation

F(x,λ)=0F(x,\lambda)=0txT(x−x0)+tλ(λ−λ0)=Δst_x^{\mathsf T}(x-x_0)+t_\lambda(\lambda-\lambda_0)=\Delta s

Derivation / construction sketch

  1. Compute a tangent to the known solution branch.
  2. Predict along that tangent.
  3. Correct using the original residual and an arclength constraint.

Symbols & assumptions

Step-size adaptation and branch switching need additional logic; the constraint is local.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Nonlinear equations

Practical use

Buckling, bifurcation, and nonlinear operating-envelope analysis.

Engineering application examples

Buckling, bifurcation, and nonlinear operating-envelope analysis.

Software product or implementation route

AUTO-07p ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Buckling, bifurcation, and nonlinear operating-envelope analysis.

Method limitations

Step-size adaptation and branch switching need additional logic; the constraint is local.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Branch followingDuffing oscillator ↗

    Adds nonlinear stiffness to an oscillator.

    For equilibrium branches or parameter sweeps near turning points; it is not a time integrator.

Relationships to other techniques

Solve algebra → Nonlinear equations

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Time integration029

Forward Euler

Advances an ODE using the current slope.

Advance in timeNumerical technique
Formulation & short derivation

Representative numerical formulation

yn+1=yn+hf(tn,yn)y_{n+1}=y_n+h f(t_n,y_n)

Derivation / construction sketch

  1. Integrate the ODE over one time interval.
  2. Approximate the integral with its left-endpoint slope.
  3. Repeat with a step size that satisfies stability and accuracy needs.

Symbols & assumptions

First-order global accuracy; the explicit stability region is limited.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Time integration

Practical use

Simple dynamics prototypes and teaching simulations.

Engineering application examples

Simple dynamics prototypes and teaching simulations.

Software product or implementation route

PETSc TSEULER ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Simple dynamics prototypes and teaching simulations.

Method limitations

First-order global accuracy; the explicit stability region is limited.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Time integrationTransient heat equation ↗

    Balances thermal storage, conduction and heat sources.

    For a nonstiff ODE or semidiscrete equation when the explicit stability bound and error budget permit; usually a baseline rather than the most efficient choice.

Relationships to other techniques

Advance in time → Time integration

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Time integration030

Backward Euler

Uses the next-step slope and solves an implicit equation.

Advance in timeNumerical technique
Formulation & short derivation

Representative numerical formulation

yn+1=yn+hf(tn+1,yn+1)y_{n+1}=y_n+h f(t_{n+1},y_{n+1})

Derivation / construction sketch

  1. Approximate the time integral with its right-endpoint slope.
  2. Rearrange as a nonlinear residual for the new state.
  3. Solve that residual each step.

Symbols & assumptions

First order and strongly damping; an implicit solve does not guarantee an accurate large step.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Time integration

Practical use

Stiff thermal and dissipative systems.

Engineering application examples

Stiff thermal and dissipative systems.

Software product or implementation route

PETSc TSBEULER ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Stiff thermal and dissipative systems.

Method limitations

First order and strongly damping; an implicit solve does not guarantee an accurate large step.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Time integrationFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationTransient heat equation ↗

    Balances thermal storage, conduction and heat sources.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationLumped-capacitance thermal model ↗

    Represents a body with one spatially uniform temperature.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationThermal resistance-capacitance network ↗

    Represents heat paths and storage with connected lumped elements.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationNewton cooling model ↗

    Uses a heat-transfer coefficient between a surface and a fluid.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationStefan phase-change problem ↗

    Couples heat transport to a moving melting or freezing boundary.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationEnthalpy–porosity model ↗

    Represents melting using enthalpy and a porous resistance in the mushy zone.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationRainfall–runoff model ↗

    Converts precipitation and catchment storage into streamflow.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationEnergy-balance climate model ↗

    Balances incoming and outgoing energy in a simplified climate system.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationBuilding thermal-zone model ↗

    Balances heat gains, losses and storage within building zones.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationPoint reactor kinetics ↗

    Approximates time-dependent neutron population with delayed-neutron groups.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

  • Time integrationBateman decay-chain model ↗

    Evolves coupled radioactive parent and daughter populations.

    For dissipative stiff evolution when first-order accuracy and damping are acceptable; solve each implicit step.

Relationships to other techniques

Advance in time → Time integration

Specific connections

  • Can supply implicit component of IMEX integration

    The representative IMEX Euler formula combines forward and backward Euler parts.

  • Is first-order member of Backward differentiation formulas

    Backward Euler is the first-order BDF formula.

  • May require nonlinear solve by Newton-Raphson method

    An implicit step is a nonlinear equation unless the dynamics are linear.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Time integration031

Crank-Nicolson

Averages endpoint slopes to obtain a second-order implicit step.

Advance in timeNumerical technique
Formulation & short derivation

Representative numerical formulation

yn+1=yn+h2[f(tn,yn)+f(tn+1,yn+1)]y_{n+1}=y_n+\frac h2[f(t_n,y_n)+f(t_{n+1},y_{n+1})]

Derivation / construction sketch

  1. Integrate over one time step.
  2. Use trapezoidal quadrature for the right-hand side.
  3. Solve the resulting implicit system.

Symbols & assumptions

Second order for smooth solutions; A-stable for the test equation but not L-stable, so stiff oscillations can persist.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Time integration

Practical use

Transient diffusion and parabolic PDE discretizations.

Engineering application examples

Transient diffusion and parabolic PDE discretizations.

Software product or implementation route

PETSc TSCN ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Transient diffusion and parabolic PDE discretizations.

Method limitations

Second order for smooth solutions; A-stable for the test equation but not L-stable, so stiff oscillations can persist.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Time integrationTime-dependent DFT (TDDFT) ↗

    Evolves electron density to approximate excited-state response.

    For smooth evolution where a second-order implicit scheme fits; stiff transients may ring without adequate resolution.

Relationships to other techniques

Advance in time → Time integration

Specific connections

  • Uses quadrature principle of Composite trapezoidal rule

    Trapezoidal integration of the time derivative produces the endpoint-average method.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Time integration032

Classical Runge-Kutta RK4

Combines four explicit slope evaluations in a fourth-order step.

Advance in timeNumerical technique
Formulation & short derivation

Representative numerical formulation

k1=f(tn,yn)k_1=f(t_n,y_n)k2=f(tn+h/2,yn+hk1/2)k_2=f(t_n+h/2,y_n+hk_1/2)k3=f(tn+h/2,yn+hk2/2)k_3=f(t_n+h/2,y_n+hk_2/2)k4=f(tn+h,yn+hk3)k_4=f(t_n+h,y_n+hk_3)yn+1=yn+h6(k1+2k2+2k3+k4)y_{n+1}=y_n+\frac h6(k_1+2k_2+2k_3+k_4)

Derivation / construction sketch

  1. Evaluate beginning, midpoint, and endpoint slopes.
  2. Choose weights to match the Taylor expansion through fourth order.
  3. Combine the stages to update the state.

Symbols & assumptions

Explicit and not suited to strongly stiff problems without very small steps; no embedded error estimate.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Time integration

Practical use

Nonstiff flight, vibration, and control simulations.

Engineering application examples

Nonstiff flight, vibration, and control simulations.

Software product or implementation route

PETSc TSRK ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Nonstiff flight, vibration, and control simulations.

Method limitations

Explicit and not suited to strongly stiff problems without very small steps; no embedded error estimate.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Time integrationNewton–Euler rigid-body model ↗

    Balances forces and moments on translating and rotating bodies.

    For smooth nonstiff ODEs with a carefully selected fixed step; perform a time-step convergence study.

  • Time integrationLagrangian mechanics ↗

    Derives motion from kinetic and potential energy with constraints.

    For smooth nonstiff ODEs with a carefully selected fixed step; perform a time-step convergence study.

  • Time integrationHamiltonian mechanics ↗

    Describes dynamics in generalized coordinates and momenta.

    For smooth nonstiff ODEs with a carefully selected fixed step; perform a time-step convergence study.

  • Time integrationMass–spring–damper model ↗

    Represents inertia, stiffness and dissipation with lumped elements.

    For smooth nonstiff ODEs with a carefully selected fixed step; perform a time-step convergence study.

  • Time integrationModal superposition model ↗

    Expands linear structural response into vibration modes.

    For smooth nonstiff ODEs with a carefully selected fixed step; perform a time-step convergence study.

  • Time integrationDuffing oscillator ↗

    Adds nonlinear stiffness to an oscillator.

    For smooth nonstiff ODEs with a carefully selected fixed step; perform a time-step convergence study.

  • Time integrationMultibody dynamics ↗

    Couples rigid or flexible bodies through joints and force elements.

    For smooth nonstiff ODEs with a carefully selected fixed step; perform a time-step convergence study.

Relationships to other techniques

Advance in time → Time integration

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Time integration033

Embedded Runge-Kutta RK45

Uses two related formulas to estimate local error and adapt the step.

Advance in timeNumerical technique
Formulation & short derivation

Representative numerical formulation

yn+1(p)=yn+h∑ibikiy_{n+1}^{(p)}=y_n+h\sum_i b_i k_ie=h∑i(bi−b^i)kie=h\sum_i(b_i-\widehat b_i)k_i

Derivation / construction sketch

  1. Reuse a shared set of stage slopes in two formulas of different order.
  2. Compare their updates to estimate local error.
  3. Accept or reject and adjust the step using scaled tolerances.

Symbols & assumptions

Representative embedded-pair form; RK45 commonly uses a Dormand-Prince 5(4) pair. Local control is not a global error guarantee.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Time integration

Practical use

Adaptive nonstiff ODE simulation.

Engineering application examples

Adaptive nonstiff ODE simulation.

Software product or implementation route

SciPy integrate.RK45 ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Adaptive nonstiff ODE simulation.

Method limitations

Representative embedded-pair form; RK45 commonly uses a Dormand-Prince 5(4) pair. Local control is not a global error guarantee.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Adaptive time integrationDiscrete dislocation dynamics ↗

    Tracks line defects and their interactions.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationContinuous stirred-tank reactor (CSTR) ↗

    Assumes a well-mixed reactor with inlet and outlet flows.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationPlug-flow reactor (PFR) ↗

    Approximates axial evolution without axial back-mixing.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationBatch reactor model ↗

    Evolves composition and energy in a closed reacting charge.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationLumped-capacitance thermal model ↗

    Represents a body with one spatially uniform temperature.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationThermal resistance-capacitance network ↗

    Represents heat paths and storage with connected lumped elements.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationNewton cooling model ↗

    Uses a heat-transfer coefficient between a surface and a fluid.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationMaxwell viscoelastic model ↗

    Combines an elastic spring and viscous dashpot in series.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationKelvin–Voigt model ↗

    Combines an elastic spring and viscous dashpot in parallel.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationStandard linear solid ↗

    Combines elastic and viscoelastic branches.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationNewton–Euler rigid-body model ↗

    Balances forces and moments on translating and rotating bodies.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationLagrangian mechanics ↗

    Derives motion from kinetic and potential energy with constraints.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationHamiltonian mechanics ↗

    Describes dynamics in generalized coordinates and momenta.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationMass–spring–damper model ↗

    Represents inertia, stiffness and dissipation with lumped elements.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationModal superposition model ↗

    Expands linear structural response into vibration modes.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationDuffing oscillator ↗

    Adds nonlinear stiffness to an oscillator.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationMultibody dynamics ↗

    Couples rigid or flexible bodies through joints and force elements.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationGeometrical optics ↗

    Approximates light propagation as rays.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationEquivalent-circuit battery model ↗

    Uses fitted electrical elements to approximate terminal behavior.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationRainfall–runoff model ↗

    Converts precipitation and catchment storage into streamflow.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationEnergy-balance climate model ↗

    Balances incoming and outgoing energy in a simplified climate system.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationSix-degree-of-freedom flight model ↗

    Evolves vehicle translation and rotation using aerodynamic and propulsion forces.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationBicycle vehicle model ↗

    Combines left and right wheels into a planar steering model.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationQuarter-car suspension model ↗

    Represents one wheel assembly and a fraction of vehicle body mass.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationSwing-equation generator model ↗

    Represents rotor-angle dynamics from mechanical-electrical power imbalance.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationBuilding thermal-zone model ↗

    Balances heat gains, losses and storage within building zones.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationState-space model ↗

    Represents system evolution with internal states, inputs and outputs.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationTransfer-function model ↗

    Relates linear time-invariant input and output in the transform domain.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationHybrid dynamical model ↗

    Combines continuous dynamics with discrete state changes.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationBond-graph model ↗

    Represents energy exchange across mechanical, electrical and other domains.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationSystem-dynamics stock-flow model ↗

    Represents accumulated quantities and their rates of change.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationHodgkin–Huxley membrane model ↗

    Uses voltage-dependent ion-channel conductances.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationFitzHugh–Nagumo model ↗

    Simplifies excitation and recovery into two dynamical variables.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationHill muscle model ↗

    Represents muscle mechanics with active and passive elements.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationWindkessel circulation model ↗

    Represents vascular resistance and compliance with lumped elements.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationMonod growth model ↗

    Relates microbial growth to a limiting substrate.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationPhysiologically based compartment model ↗

    Represents exchange between anatomically motivated compartments.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationNewtonian gravitational N-body model ↗

    Evolves masses under mutual inverse-square attraction.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationStellar structure model ↗

    Couples hydrostatic balance, energy transport and energy generation.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

  • Adaptive time integrationFLRW cosmological model ↗

    Assumes a homogeneous and isotropic expanding spacetime.

    For nonstiff ODEs; detect events and check whether constraints or fast scales require another solver.

Relationships to other techniques

Advance in time → Time integration

Specific connections

  • Can advance Finite volume method

    A method-of-lines discretization can be advanced by a suitable ODE integrator; stability must be checked.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Time integration034

Backward differentiation formulas

Use several past states to approximate the new-time derivative.

Advance in timeNumerical technique
Formulation & short derivation

Representative numerical formulation

3yn+1−4yn+yn−12h=f(tn+1,yn+1)\frac{3y_{n+1}-4y_n+y_{n-1}}{2h}=f(t_{n+1},y_{n+1})

Derivation / construction sketch

  1. Interpolate recent states with a polynomial.
  2. Differentiate that polynomial at the newest time.
  3. Solve the implicit residual for the new state.

Symbols & assumptions

Constant-step BDF2 shown; higher-order formulas have different stability limits and require startup/history handling.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Time integration

Practical use

Chemical kinetics, batteries, and other stiff systems.

Engineering application examples

Chemical kinetics, batteries, and other stiff systems.

Software product or implementation route

SciPy integrate.BDF ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Chemical kinetics, batteries, and other stiff systems.

Method limitations

Constant-step BDF2 shown; higher-order formulas have different stability limits and require startup/history handling.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Stiff time integrationMass-action reaction kinetics ↗

    Relates reaction rates to species concentrations and reaction orders.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationContinuous stirred-tank reactor (CSTR) ↗

    Assumes a well-mixed reactor with inlet and outlet flows.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationPlug-flow reactor (PFR) ↗

    Approximates axial evolution without axial back-mixing.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationBatch reactor model ↗

    Evolves composition and energy in a closed reacting charge.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationMaxwell viscoelastic model ↗

    Combines an elastic spring and viscous dashpot in series.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationKelvin–Voigt model ↗

    Combines an elastic spring and viscous dashpot in parallel.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationStandard linear solid ↗

    Combines elastic and viscoelastic branches.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationLumped RLC circuit model ↗

    Uses resistors, capacitors and inductors connected by Kirchhoff laws.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationTransmission-line electrical model ↗

    Represents distributed inductance, capacitance and losses.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationDrift–diffusion semiconductor model ↗

    Combines electrostatics with carrier drift, diffusion and continuity.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationHydrodynamic carrier model ↗

    Adds carrier-energy or momentum information to transport.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationPoisson–Nernst–Planck model ↗

    Couples electrostatics to diffusion and migration of ions.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationDoyle–Fuller–Newman (DFN/P2D) model ↗

    Combines porous-electrode transport and particle diffusion.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationSingle-particle battery model (SPM) ↗

    Represents each electrode by a representative active-material particle.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationSingle-particle model with electrolyte (SPMe) ↗

    Adds electrolyte concentration effects to a single-particle approximation.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationHodgkin–Huxley membrane model ↗

    Uses voltage-dependent ion-channel conductances.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationFitzHugh–Nagumo model ↗

    Simplifies excitation and recovery into two dynamical variables.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationHill muscle model ↗

    Represents muscle mechanics with active and passive elements.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationWindkessel circulation model ↗

    Represents vascular resistance and compliance with lumped elements.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationMonod growth model ↗

    Relates microbial growth to a limiting substrate.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationPhysiologically based compartment model ↗

    Represents exchange between anatomically motivated compartments.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationPoint reactor kinetics ↗

    Approximates time-dependent neutron population with delayed-neutron groups.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

  • Stiff time integrationBateman decay-chain model ↗

    Evolves coupled radioactive parent and daughter populations.

    For stiff ODEs; constrained or algebraic formulations require an appropriate DAE-capable implementation, not a plain ODE routine.

Relationships to other techniques

Advance in time → Time integration

Specific connections

  • Multistep generalization of Backward Euler

    Backward Euler is the first-order BDF formula.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Time integration035

Velocity Verlet

Advances positions and velocities with a symmetric force update.

Advance in timeNumerical technique
Formulation & short derivation

Representative numerical formulation

qn+1=qn+hvn+h22a(qn)q_{n+1}=q_n+h v_n+\frac{h^2}{2}a(q_n)vn+1=vn+h2[a(qn)+a(qn+1)]v_{n+1}=v_n+\frac h2[a(q_n)+a(q_{n+1})]

Derivation / construction sketch

  1. Apply a half velocity kick.
  2. Drift the position with the intermediate velocity.
  3. Apply the remaining half kick using the new force.

Symbols & assumptions

Second order and symplectic for suitable separable Hamiltonian systems at fixed step; arbitrary damping changes these properties.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Time integration

Practical use

Molecular dynamics and long-time particle mechanics.

Engineering application examples

Molecular dynamics and long-time particle mechanics.

Software product or implementation route

LAMMPS ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Molecular dynamics and long-time particle mechanics.

Method limitations

Second order and symplectic for suitable separable Hamiltonian systems at fixed step; arbitrary damping changes these properties.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Mechanical time integrationClassical molecular dynamics (MD) ↗

    Integrates atomic motion under specified interaction forces.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationAb initio molecular dynamics ↗

    Computes interatomic forces from electronic-structure calculations during motion.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationLennard–Jones potential ↗

    Combines short-range repulsion with an inverse-sixth-power attraction.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationMorse potential ↗

    Represents an anharmonic bond with a finite dissociation energy.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationEmbedded-atom method (EAM) ↗

    Combines pair interactions with an embedding energy dependent on local electron density.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationModified embedded-atom method (MEAM) ↗

    Extends embedding models with angular information.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationTersoff bond-order potential ↗

    Makes bond strength depend on the local bonding environment.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationStillinger–Weber potential ↗

    Uses two-body and three-body terms to favor local tetrahedral structure.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationReaxFF reactive force field ↗

    Uses variable bond orders and charge equilibration to represent chemical reactions.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationAMBER force-field family ↗

    Uses parameterized bonded and nonbonded interactions for biomolecules.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationCHARMM force-field family ↗

    Models biomolecular interactions with chemistry-specific parameter sets.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationOPLS force-field family ↗

    Uses parameterized molecular interactions developed for condensed phases.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationSPC/E water model ↗

    Approximates water using a rigid three-site classical model.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationTIP4P water-model family ↗

    Uses a four-site geometry with an off-oxygen charge site.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationDrude polarizable model ↗

    Uses auxiliary charged particles to represent induced polarization.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationMachine-learned interatomic potential ↗

    Fits atomic energies and forces from reference data using statistical learning.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationCoarse-grained molecular model ↗

    Groups atoms into effective interaction sites.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationMartini coarse-grained model ↗

    Uses mapped molecular beads and parameterized interactions.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationNewton–Euler rigid-body model ↗

    Balances forces and moments on translating and rotating bodies.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationLagrangian mechanics ↗

    Derives motion from kinetic and potential energy with constraints.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationHamiltonian mechanics ↗

    Describes dynamics in generalized coordinates and momenta.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationMass–spring–damper model ↗

    Represents inertia, stiffness and dissipation with lumped elements.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationModal superposition model ↗

    Expands linear structural response into vibration modes.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationDuffing oscillator ↗

    Adds nonlinear stiffness to an oscillator.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationMultibody dynamics ↗

    Couples rigid or flexible bodies through joints and force elements.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationNewtonian gravitational N-body model ↗

    Evolves masses under mutual inverse-square attraction.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationQM/MM coupling ↗

    Combines quantum mechanics in a selected region with molecular mechanics around it.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationAtomistic–continuum coupling ↗

    Connects particle-level and continuum descriptions.

    For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

  • Mechanical time integrationGreen-Kubo viscosity relation ↗

    Obtains equilibrium shear viscosity from the time integral of microscopic shear-stress fluctuations.

    Generate equilibrium MD trajectories and shear-pressure samples with a compatible ensemble and force field; this time integrator alone does not estimate viscosity.

  • Mechanical time integrationHarmonic lattice dynamics ↗

    Computes phonon modes from a quadratic expansion of crystal potential energy.

    Fit force constants with SVD, cross-check real-time harmonic motion with Verlet, and test displacement/time-step refinement. Obtain phonon frequencies with a Hermitian dynamical-matrix eigensolver. For compatible position-dependent conservative forces; constraints, thermostats, and stochastic forces require specialized extensions.

Relationships to other techniques

Advance in time → Time integration

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Time integration036

IMEX integration

Treats stiff terms implicitly and nonstiff terms explicitly.

Advance in timeNumerical technique
Formulation & short derivation

Representative numerical formulation

yn+1=yn+hfE(yn)+hfI(yn+1)y_{n+1}=y_n+h f_E(y_n)+h f_I(y_{n+1})

Derivation / construction sketch

  1. Split the dynamics by stiffness or computational structure.
  2. Use an explicit update for one part and an implicit update for the other.
  3. Combine with compatible order and stability conditions.

Symbols & assumptions

First-order IMEX Euler shown; splitting and stability constraints remain problem dependent.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Time integration

Practical use

Reaction-transport and advection-diffusion systems.

Engineering application examples

Reaction-transport and advection-diffusion systems.

Software product or implementation route

PETSc TSARKIMEX ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Reaction-transport and advection-diffusion systems.

Method limitations

First-order IMEX Euler shown; splitting and stability constraints remain problem dependent.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Split time integrationCahn–Hilliard model ↗

    Evolves a conserved composition field through chemical-potential gradients.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationAllen–Cahn model ↗

    Evolves a nonconserved order parameter toward lower free energy.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationPhase-field crystal model ↗

    Uses a periodic density-like field to represent crystalline ordering.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationPopulation balance model ↗

    Tracks the distribution of particle sizes or other internal properties.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationFickian diffusion ↗

    Relates diffusive flux to concentration gradients.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationMaxwell–Stefan diffusion ↗

    Represents multicomponent diffusion through interspecies friction.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationAdvection–diffusion–reaction model ↗

    Combines bulk transport, diffusion and reaction sources.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationReynolds-averaged Navier–Stokes (RANS) ↗

    Models mean flow with closure for unresolved turbulent stresses.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationSpalart–Allmaras model ↗

    Uses a transported turbulence variable to obtain eddy viscosity.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationk–epsilon model ↗

    Uses turbulent kinetic energy and dissipation rate to close mean flow.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationk–omega model ↗

    Uses turbulent kinetic energy and specific dissipation rate.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationSST k–omega model ↗

    Blends near-wall and outer-flow behavior with a shear-stress limiter.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationReynolds-stress transport model ↗

    Transports individual turbulent stress components.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationLarge-eddy simulation (LES) ↗

    Resolves larger turbulent motions and models subgrid effects.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationSmagorinsky subgrid model ↗

    Relates subgrid eddy viscosity to resolved strain and filter scale.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationDetached-eddy simulation (DES) ↗

    Combines RANS near walls with LES-like treatment away from them.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationVolume-of-fluid (VOF) representation ↗

    Tracks phase volume fractions to represent an interface.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationEuler–Euler two-fluid model ↗

    Treats phases as interpenetrating continua with exchange terms.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationLagrangian particle tracking ↗

    Tracks discrete particles through a carrier flow.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationSaint-Venant shallow-water model ↗

    Depth-averages mass and momentum in free-surface flow.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationKinematic-wave routing ↗

    Simplifies flow routing by approximating dominant slope and friction balance.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationGroundwater flow model ↗

    Combines water conservation with porous-flow relations.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationAdvection–dispersion groundwater model ↗

    Represents contaminant transport and spreading through an aquifer.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationNumerical weather prediction ↗

    Evolves atmospheric dynamics and thermodynamics from an analyzed initial state.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationGeneral circulation model (GCM) ↗

    Represents large-scale atmospheric or oceanic circulation.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationEarth system model (ESM) ↗

    Couples atmosphere, ocean, land, ice and biogeochemical processes.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationOcean circulation model ↗

    Evolves ocean momentum, temperature and salinity.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationSea-ice thermodynamic-dynamic model ↗

    Couples freezing, melting and ice motion.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationElastic seismic-wave model ↗

    Propagates elastic disturbances through Earth materials.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationPennes bioheat model ↗

    Adds perfusion and metabolic heat to tissue heat transfer.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationReaction–diffusion morphogenesis model ↗

    Couples reacting substances with diffusion.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationBoltzmann kinetic equation ↗

    Evolves a particle distribution under transport and collisions.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationVlasov–Poisson model ↗

    Couples collisionless distribution dynamics to electrostatic fields.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationVlasov–Maxwell model ↗

    Couples collisionless kinetic distributions to electromagnetic fields.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationMagnetohydrodynamics (MHD) ↗

    Treats a conducting fluid coupled to a magnetic field.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationNeutron diffusion approximation ↗

    Simplifies neutron transport to a diffusion description.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationGeneral relativity model ↗

    Relates spacetime curvature to matter and energy.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationFluid–structure interaction (FSI) ↗

    Couples fluid loads with structural motion or deformation.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

  • Split time integrationThermomechanical coupling ↗

    Couples temperature evolution and mechanical response.

    When a meaningful stiff/nonstiff split exists; check the stability and coupling error of both parts.

Relationships to other techniques

Advance in time → Time integration

Specific connections

  • Uses implicit step in first-order form Backward Euler

    The representative IMEX Euler formula combines forward and backward Euler parts.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Optimization & inverse problems037

Gradient descent

Moves downhill along the negative objective gradient.

OptimizeNumerical technique
Formulation & short derivation

Representative numerical formulation

xk+1=xk−αk∇f(xk)x_{k+1}=x_k-\alpha_k\nabla f(x_k)

Derivation / construction sketch

  1. Use a first-order local objective approximation.
  2. Choose a descent direction opposite to the gradient.
  3. Select a step using a line search or a justified fixed rule.

Symbols & assumptions

Convergence depends on smoothness and step size; nonconvex objectives may have local minima.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Optimization & inverse problems

Practical use

Parameter tuning and inverse-model fitting.

Engineering application examples

Parameter tuning and inverse-model fitting.

Software product or implementation route

SciPy optimize.line_search ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Parameter tuning and inverse-model fitting.

Method limitations

Convergence depends on smoothness and step size; nonconvex objectives may have local minima.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Parameter fittingPhysics-informed neural network (PINN) ↗

    Trains a neural approximation using data and governing-equation residuals.

    For a differentiable calibration objective with a justified step rule; slow convergence is possible under poor scaling.

  • Parameter fittingNeural operator ↗

    Learns a map between function-valued inputs and outputs.

    For a differentiable calibration objective with a justified step rule; slow convergence is possible under poor scaling.

Relationships to other techniques

Optimize → Optimization & inverse problems

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Optimization & inverse problems038

Newton optimization

Uses curvature to compute a local stationary-point correction.

OptimizeNumerical technique
Formulation & short derivation

Representative numerical formulation

∇2f(xk)pk=−∇f(xk)\nabla^2 f(x_k)p_k=-\nabla f(x_k)xk+1=xk+αkpkx_{k+1}=x_k+\alpha_kp_k

Derivation / construction sketch

  1. Build a second-order Taylor model of the objective.
  2. Set its gradient to zero.
  3. Globalize the step with line search or a trust region.

Symbols & assumptions

An indefinite Hessian may not produce descent; regularization or trust-region safeguards may be needed.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Optimization & inverse problems

Practical use

Smooth engineering design optimization.

Engineering application examples

Smooth engineering design optimization.

Software product or implementation route

SciPy Newton-CG ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Smooth engineering design optimization.

Method limitations

An indefinite Hessian may not produce descent; regularization or trust-region safeguards may be needed.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Design / calibrationMachine-learned interatomic potential ↗

    Fits atomic energies and forces from reference data using statistical learning.

    For a smooth objective with usable curvature; distinguish optimization from solving the physical state equations.

Relationships to other techniques

Optimize → Optimization & inverse problems

Specific connections

  • Can use gradients from Adjoint sensitivity analysis

    An adjoint supplies objective sensitivities; Hessian information requires additional work.

  • Has quasi-Newton alternative BFGS quasi-Newton

    Gradient differences replace direct Hessian evaluation.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Optimization & inverse problems039

BFGS quasi-Newton

Updates an inverse-Hessian approximation using gradient differences.

OptimizeNumerical technique
Formulation & short derivation

Representative numerical formulation

Hk+1=(I−ρsyT)Hk(I−ρysT)+ρssTH_{k+1}=(I-\rho s y^{\mathsf T})H_k(I-\rho y s^{\mathsf T})+\rho ss^{\mathsf T}ρ=1yTs\rho=\frac1{y^{\mathsf T}s}

Derivation / construction sketch

  1. Take a descent step and measure the gradient change.
  2. Enforce a secant relation while preserving symmetry.
  3. Use the BFGS update and a curvature-respecting line search.

Symbols & assumptions

Positive definiteness needs positive curvature y dot s; noisy gradients can be problematic.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Optimization & inverse problems

Practical use

Smooth unconstrained calibration and shape optimization.

Engineering application examples

Smooth unconstrained calibration and shape optimization.

Software product or implementation route

SciPy BFGS ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Smooth unconstrained calibration and shape optimization.

Method limitations

Positive definiteness needs positive curvature y dot s; noisy gradients can be problematic.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Design / calibrationMachine-learned interatomic potential ↗

    Fits atomic energies and forces from reference data using statistical learning.

    For smooth moderate-size parameter fitting without explicit Hessians; use accurate gradients and a line search.

Relationships to other techniques

Optimize → Optimization & inverse problems

Specific connections

  • Has limited-memory relative L-BFGS-B

    Stores curvature pairs and incorporates bound-aware steps.

  • Curvature approximation to Newton optimization

    Gradient differences replace direct Hessian evaluation.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Optimization & inverse problems040

L-BFGS-B

Uses limited curvature history with bound constraints.

OptimizeNumerical technique
Formulation & short derivation

Representative numerical formulation

min⁡xf(x)\min_x f(x)li≤xi≤uil_i\le x_i\le u_ipk≈−Hk∇f(xk)p_k\approx-H_k\nabla f(x_k)

Derivation / construction sketch

  1. Store a limited set of displacement and gradient-difference pairs.
  2. Identify bound-active variables and a feasible search direction.
  3. Apply a line search within bounds.

Symbols & assumptions

Schematic direction shown; the algorithm includes generalized Cauchy and subspace steps. Bounds do not make a nonconvex problem globally solvable.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Optimization & inverse problems

Practical use

Large parameter estimation with physical parameter bounds.

Engineering application examples

Large parameter estimation with physical parameter bounds.

Software product or implementation route

SciPy L-BFGS-B ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Large parameter estimation with physical parameter bounds.

Method limitations

Schematic direction shown; the algorithm includes generalized Cauchy and subspace steps. Bounds do not make a nonconvex problem globally solvable.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Bounded calibrationClassical molecular dynamics (MD) ↗

    Integrates atomic motion under specified interaction forces.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationAb initio molecular dynamics ↗

    Computes interatomic forces from electronic-structure calculations during motion.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationLennard–Jones potential ↗

    Combines short-range repulsion with an inverse-sixth-power attraction.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationMorse potential ↗

    Represents an anharmonic bond with a finite dissociation energy.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationEmbedded-atom method (EAM) ↗

    Combines pair interactions with an embedding energy dependent on local electron density.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationModified embedded-atom method (MEAM) ↗

    Extends embedding models with angular information.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationTersoff bond-order potential ↗

    Makes bond strength depend on the local bonding environment.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationStillinger–Weber potential ↗

    Uses two-body and three-body terms to favor local tetrahedral structure.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationReaxFF reactive force field ↗

    Uses variable bond orders and charge equilibration to represent chemical reactions.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationAMBER force-field family ↗

    Uses parameterized bonded and nonbonded interactions for biomolecules.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationCHARMM force-field family ↗

    Models biomolecular interactions with chemistry-specific parameter sets.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationOPLS force-field family ↗

    Uses parameterized molecular interactions developed for condensed phases.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationSPC/E water model ↗

    Approximates water using a rigid three-site classical model.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationTIP4P water-model family ↗

    Uses a four-site geometry with an off-oxygen charge site.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationDrude polarizable model ↗

    Uses auxiliary charged particles to represent induced polarization.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationMachine-learned interatomic potential ↗

    Fits atomic energies and forces from reference data using statistical learning.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationShockley diode model ↗

    Approximates diode current with an exponential voltage relation.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationEbers–Moll transistor model ↗

    Models coupled junction currents in a bipolar transistor.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationMOSFET square-law model ↗

    Approximates long-channel transistor current from terminal voltages.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationBSIM compact-model family ↗

    Uses detailed parameterized MOS transistor relations.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationNernst equilibrium potential ↗

    Relates electrochemical equilibrium potential to species activities.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationButler–Volmer kinetics ↗

    Relates interfacial current to electrochemical overpotential.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationTafel approximation ↗

    Approximates high-overpotential behavior of Butler–Volmer kinetics.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationGaussian-process surrogate ↗

    Predicts responses with a probabilistic function model fitted to samples.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationPhysics-informed neural network (PINN) ↗

    Trains a neural approximation using data and governing-equation residuals.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

  • Bounded calibrationNeural operator ↗

    Learns a map between function-valued inputs and outputs.

    For smooth parameter fitting with physical bounds; local minima and identifiability still require assessment.

Relationships to other techniques

Optimize → Optimization & inverse problems

Specific connections

  • Limited-memory bound-constrained relative of BFGS quasi-Newton

    Stores curvature pairs and incorporates bound-aware steps.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Optimization & inverse problems041

Gauss-Newton least squares

Linearizes residuals to solve a nonlinear least-squares problem.

OptimizeNumerical technique
Formulation & short derivation

Representative numerical formulation

min⁡x12∥r(x)∥2\min_x\frac12\lVert r(x)\rVert^2JTJp=−JTrJ^{\mathsf T}Jp=-J^{\mathsf T}r

Derivation / construction sketch

  1. Expand each residual to first order.
  2. Minimize the squared norm of the linearized residual.
  3. Update and repeat, preferably solving the least-squares system by stable factorization.

Symbols & assumptions

Normal equations show the concept but square the condition number; QR/SVD can be preferable. Best near a suitable small-residual solution.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Optimization & inverse problems

Practical use

Calibration of sensors and constitutive parameters.

Engineering application examples

Calibration of sensors and constitutive parameters.

Software product or implementation route

SciPy optimize.least_squares ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Calibration of sensors and constitutive parameters.

Method limitations

Normal equations show the concept but square the condition number; QR/SVD can be preferable. Best near a suitable small-residual solution.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • CalibrationMachine-learned interatomic potential ↗

    Fits atomic energies and forces from reference data using statistical learning.

    For differentiable residual-based parameter fitting, especially near a suitable small-residual solution.

Relationships to other techniques

Optimize → Optimization & inverse problems

Specific connections

  • Can use stable linear solve by QR factorization

    QR avoids explicitly forming the normal equations.

  • Can be stabilized by Levenberg-Marquardt

    Damping limits the correction when the local least-squares model is unreliable.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Optimization & inverse problems042

Levenberg-Marquardt

Regularizes a Gauss-Newton step to balance stability and progress.

OptimizeNumerical technique
Formulation & short derivation

Representative numerical formulation

(JTJ+λI)p=−JTr(J^{\mathsf T}J+\lambda I)p=-J^{\mathsf T}r

Derivation / construction sketch

  1. Construct the linearized residual objective.
  2. Add a penalty on the correction size.
  3. Adjust damping according to agreement between predicted and actual reduction.

Symbols & assumptions

Identity-scaled schematic form; practical implementations use scaling and trust-region logic.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Optimization & inverse problems

Practical use

Nonlinear curve fitting and experimental model calibration.

Engineering application examples

Nonlinear curve fitting and experimental model calibration.

Software product or implementation route

SciPy optimize.least_squares ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Nonlinear curve fitting and experimental model calibration.

Method limitations

Identity-scaled schematic form; practical implementations use scaling and trust-region logic.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • CalibrationNRTL activity model ↗

    Uses local-composition parameters to describe nonideal liquid mixtures.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationUNIQUAC activity model ↗

    Combines molecular size, shape and interaction contributions.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationDebye–Hückel model ↗

    Approximates ionic activity using screened electrostatic interactions.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationArrhenius rate model ↗

    Relates a rate coefficient to temperature through an activation energy.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationTransition-state theory ↗

    Estimates reaction rates from a free-energy barrier.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationMichaelis–Menten kinetics ↗

    Approximates enzyme reaction rates with substrate saturation.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationLangmuir adsorption isotherm ↗

    Models adsorption on equivalent sites with finite occupancy.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationLangmuir–Hinshelwood kinetics ↗

    Models surface reactions involving adsorbed reactants.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationMaxwell viscoelastic model ↗

    Combines an elastic spring and viscous dashpot in series.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationKelvin–Voigt model ↗

    Combines an elastic spring and viscous dashpot in parallel.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationStandard linear solid ↗

    Combines elastic and viscoelastic branches.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationNorton creep law ↗

    Relates creep rate to a power of stress.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationLinear elastic fracture mechanics (LEFM) ↗

    Uses crack-tip intensity parameters in an elastic body.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationCohesive-zone model ↗

    Uses traction-separation relations across a fracture process zone.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationParis fatigue crack-growth law ↗

    Relates cyclic crack-growth rate to stress-intensity-factor range.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationMiner cumulative damage rule ↗

    Adds fractions of fatigue life consumed by load cycles.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationArchard wear model ↗

    Relates wear volume to load, sliding distance and hardness.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationMagnetic-circuit model ↗

    Uses reluctance and magnetomotive force in lumped magnetic paths.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationJiles–Atherton hysteresis model ↗

    Represents path-dependent magnetization with phenomenological parameters.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationGaussian beam model ↗

    Represents a paraxial beam with a Gaussian transverse profile.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationDrude–Lorentz optical model ↗

    Represents free-carrier and bound-charge contributions to permittivity.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationEquivalent-circuit battery model ↗

    Uses fitted electrical elements to approximate terminal behavior.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationForchheimer model ↗

    Adds inertial resistance to porous flow.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • Calibrationvan Genuchten retention model ↗

    Relates water saturation to pressure head with fitted parameters.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationRainfall–runoff model ↗

    Converts precipitation and catchment storage into streamflow.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationEnergy-balance climate model ↗

    Balances incoming and outgoing energy in a simplified climate system.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationSix-degree-of-freedom flight model ↗

    Evolves vehicle translation and rotation using aerodynamic and propulsion forces.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationBicycle vehicle model ↗

    Combines left and right wheels into a planar steering model.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationQuarter-car suspension model ↗

    Represents one wheel assembly and a fraction of vehicle body mass.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationPacejka tire model ↗

    Uses empirical nonlinear formulas for tire forces.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationSwing-equation generator model ↗

    Represents rotor-angle dynamics from mechanical-electrical power imbalance.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationBuilding thermal-zone model ↗

    Balances heat gains, losses and storage within building zones.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationHodgkin–Huxley membrane model ↗

    Uses voltage-dependent ion-channel conductances.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationFitzHugh–Nagumo model ↗

    Simplifies excitation and recovery into two dynamical variables.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationHill muscle model ↗

    Represents muscle mechanics with active and passive elements.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationWindkessel circulation model ↗

    Represents vascular resistance and compliance with lumped elements.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationMonod growth model ↗

    Relates microbial growth to a limiting substrate.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationPhysiologically based compartment model ↗

    Represents exchange between anatomically motivated compartments.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationDigital twin framework ↗

    Links an evolving model of a specific asset with observations.

    For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationStokes-Einstein diffusion relation ↗

    Connects Brownian translational diffusion to temperature, solvent viscosity, and hydrodynamic particle radius.

    Fit a positive hydrodynamic radius or viscosity to diffusion measurements; parameters may be unidentifiable if fitted together.

  • CalibrationEinstein crystal heat-capacity model ↗

    Treats crystal vibrations as independent quantum oscillators at a single frequency.

    Use fitting for an Einstein temperature, interpolation for tabulated responses, or refinement to check derived quantities; the displayed formula itself is explicit. For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationDebye phonon model ↗

    Approximates acoustic phonons by a continuum spectrum with a mode-count cutoff.

    Use quadrature for the full Debye integral and fitting for a Debye temperature; the cubic law is only a low-temperature asymptote. For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

  • CalibrationSommerfeld free-electron model ↗

    Describes conduction electrons as a degenerate, noninteracting Fermi gas.

    Integrate the free-electron density of states with Fermi occupations at finite temperature, or fit a low-temperature heat-capacity coefficient; preserve electron number. For nonlinear least-squares fitting with damping; standard LM does not handle arbitrary constraints.

Relationships to other techniques

Optimize → Optimization & inverse problems

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Optimization & inverse problems043

Sequential quadratic programming

Solves a sequence of locally quadratic constrained subproblems.

OptimizeNumerical technique
Formulation & short derivation

Representative numerical formulation

min⁡p12pTBkp+∇fkTp\min_p\frac12p^{\mathsf T}B_kp+\nabla f_k^{\mathsf T}pck+Jkp=0c_k+J_kp=0

Derivation / construction sketch

  1. Approximate the Lagrangian curvature.
  2. Linearize constraints and solve a quadratic subproblem.
  3. Globalize and update primal and multiplier estimates.

Symbols & assumptions

Equality form shown; inequality constraints require active-set or related treatment and constraint qualifications.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Optimization & inverse problems

Practical use

Constrained geometry and operating-point optimization.

Engineering application examples

Constrained geometry and operating-point optimization.

Software product or implementation route

SciPy SLSQP ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Constrained geometry and operating-point optimization.

Method limitations

Equality form shown; inequality constraints require active-set or related treatment and constraint qualifications.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Constrained designGibbs-energy minimization ↗

    Finds equilibrium by minimizing free energy under conservation constraints.

    For smooth constrained parameter/design optimization with derivatives and constraint regularity.

  • Constrained designCALPHAD model ↗

    Combines assessed phase free energies to predict equilibria.

    For smooth constrained parameter/design optimization with derivatives and constraint regularity.

  • Constrained designModel predictive control ↗

    Optimizes future actions using a predictive model and constraints.

    For smooth constrained parameter/design optimization with derivatives and constraint regularity.

Relationships to other techniques

Optimize → Optimization & inverse problems

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Optimization & inverse problems044

Interior-point optimization

Approaches inequality-constrained solutions through barrier subproblems.

OptimizeNumerical technique
Formulation & short derivation

Representative numerical formulation

min⁡xf(x)−μ∑iln⁡si(x)\min_x f(x)-\mu\sum_i\ln s_i(x)si(x)>0s_i(x)>0

Derivation / construction sketch

  1. Represent inequalities by positive slack functions.
  2. Penalize approach to the feasible boundary with a logarithmic barrier.
  3. Solve a sequence of primal-dual or barrier systems as the parameter decreases.

Symbols & assumptions

Feasibility, scaling, and conditioning near active constraints need care; nonconvex problems lack a general global guarantee.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Optimization & inverse problems

Practical use

Large constrained engineering design problems.

Engineering application examples

Large constrained engineering design problems.

Software product or implementation route

SciPy trust-constr ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Large constrained engineering design problems.

Method limitations

Feasibility, scaling, and conditioning near active constraints need care; nonconvex problems lack a general global guarantee.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Constrained designGibbs-energy minimization ↗

    Finds equilibrium by minimizing free energy under conservation constraints.

    For appropriately formulated inequality-constrained design or control problems; scale constraints and verify feasibility.

  • Constrained designCALPHAD model ↗

    Combines assessed phase free energies to predict equilibria.

    For appropriately formulated inequality-constrained design or control problems; scale constraints and verify feasibility.

  • Constrained designModel predictive control ↗

    Optimizes future actions using a predictive model and constraints.

    For appropriately formulated inequality-constrained design or control problems; scale constraints and verify feasibility.

Relationships to other techniques

Optimize → Optimization & inverse problems

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Approximation & reduction045

Lagrange interpolation

Passes a polynomial through prescribed distinct data points.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

p(x)=∑iyiℓi(x)p(x)=\sum_i y_i\ell_i(x)ℓi(x)=∏j≠ix−xjxi−xj\ell_i(x)=\prod_{j\ne i}\frac{x-x_j}{x_i-x_j}

Derivation / construction sketch

  1. Construct a basis that is one at its own node and zero at all others.
  2. Weight each basis function by its data value.
  3. Sum to reproduce the nodal data.

Symbols & assumptions

High-degree equispaced interpolation can oscillate badly; noisy data may need fitting instead.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Approximation & reduction

Practical use

Tabulated material-property interpolation.

Engineering application examples

Tabulated material-property interpolation.

Software product or implementation route

SciPy interpolate.lagrange ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Tabulated material-property interpolation.

Method limitations

High-degree equispaced interpolation can oscillate badly; noisy data may need fitting instead.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Tabulated dataArrhenius rate model ↗

    Relates a rate coefficient to temperature through an activation energy.

    For small, well-chosen interpolation grids; avoid high-degree equispaced interpolation and extrapolation.

Relationships to other techniques

Approximate & sample → Approximation & reduction

Specific connections

  • Can be evaluated using Barycentric interpolation

    The same interpolation polynomial is evaluated with precomputed barycentric weights.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Approximation & reduction046

Barycentric interpolation

Evaluates the interpolation polynomial using precomputed weights.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

p(x)=∑iwiyi/(x−xi)∑iwi/(x−xi)p(x)=\frac{\sum_i w_i y_i/(x-x_i)}{\sum_i w_i/(x-x_i)}wi=1∏j≠i(xi−xj)w_i=\frac1{\prod_{j\ne i}(x_i-x_j)}

Derivation / construction sketch

  1. Factor common products from the Lagrange basis.
  2. Cancel the common factor between numerator and denominator.
  3. At a node, return the stored nodal value directly.

Symbols & assumptions

A stable evaluation form cannot eliminate poor node choice or intrinsic interpolation conditioning.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Approximation & reduction

Practical use

Repeated evaluations of polynomial surrogate tables.

Engineering application examples

Repeated evaluations of polynomial surrogate tables.

Software product or implementation route

SciPy interpolate.BarycentricInterpolator ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Repeated evaluations of polynomial surrogate tables.

Method limitations

A stable evaluation form cannot eliminate poor node choice or intrinsic interpolation conditioning.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Tabulated dataForchheimer model ↗

    Adds inertial resistance to porous flow.

    For repeated evaluation of a polynomial interpolant with suitable nodes; this is not itself a physical solver.

  • Tabulated datavan Genuchten retention model ↗

    Relates water saturation to pressure head with fitted parameters.

    For repeated evaluation of a polynomial interpolant with suitable nodes; this is not itself a physical solver.

Relationships to other techniques

Approximate & sample → Approximation & reduction

Specific connections

  • Evaluation form of Lagrange interpolation

    The same interpolation polynomial is evaluated with precomputed barycentric weights.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Approximation & reduction047

Cubic spline interpolation

Joins piecewise cubic polynomials with continuity constraints.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

Si(x)=ai+biξ+ciξ2+diξ3S_i(x)=a_i+b_i\xi+c_i\xi^2+d_i\xi^3S,S′,S′′ are continuousS,S',S''\text{ are continuous}

Derivation / construction sketch

  1. Fit a cubic on each interval.
  2. Match values and first two derivatives at interior knots.
  3. Add endpoint conditions to close the system.

Symbols & assumptions

Endpoint choices affect the result; ordinary cubic splines need not preserve positivity or monotonicity.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Approximation & reduction

Practical use

Smooth material curves and measured trajectory reconstruction.

Engineering application examples

Smooth material curves and measured trajectory reconstruction.

Software product or implementation route

SciPy interpolate.CubicSpline ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Smooth material curves and measured trajectory reconstruction.

Method limitations

Endpoint choices affect the result; ordinary cubic splines need not preserve positivity or monotonicity.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Tabulated dataPacejka tire model ↗

    Uses empirical nonlinear formulas for tire forces.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataGeometrically scaled physical model ↗

    Reproduces a system's shape at another size.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataWind-tunnel model ↗

    Uses a controlled air stream around a physical specimen.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataHydraulic flume model ↗

    Uses physical water flow with selected similarity conditions.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataShake-table structural model ↗

    Excites a physical structure with controlled base motion.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataPhotoelastic model ↗

    Uses stress-induced optical birefringence to visualize stress patterns.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataElectrical analog model ↗

    Maps another physical system onto an electrical network.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataHardware-in-the-loop model ↗

    Couples actual hardware to simulated parts of a system.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataDimensional-analysis similarity model ↗

    Uses dimensionless groups to relate tests across scales.

    For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

  • Tabulated dataStokes-Einstein diffusion relation ↗

    Connects Brownian translational diffusion to temperature, solvent viscosity, and hydrodynamic particle radius.

    Interpolate tabulated solvent viscosity versus temperature before evaluating diffusivity; avoid extrapolation and preserve positive viscosity.

  • Tabulated dataEinstein crystal heat-capacity model ↗

    Treats crystal vibrations as independent quantum oscillators at a single frequency.

    Use fitting for an Einstein temperature, interpolation for tabulated responses, or refinement to check derived quantities; the displayed formula itself is explicit. For smooth interpolation of coefficients or responses; ordinary splines do not guarantee positivity or monotonicity.

Relationships to other techniques

Approximate & sample → Approximation & reduction

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Approximation & reduction048

Polynomial least squares

Fits basis coefficients by minimizing data residuals.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

c∗=arg⁡min⁡c∥Ac−y∥22c^*=\arg\min_c\lVert Ac-y\rVert_2^2Aij=ϕj(xi)A_{ij}=\phi_j(x_i)

Derivation / construction sketch

  1. Choose a polynomial basis and assemble a design matrix.
  2. Project the data onto its column space.
  3. Use QR or SVD to obtain a stable least-squares solution.

Symbols & assumptions

Scale coordinates and avoid unnecessarily high degree; outliers and correlated errors require additional treatment.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Approximation & reduction

Practical use

Calibration curves and response-surface approximations.

Engineering application examples

Calibration curves and response-surface approximations.

Software product or implementation route

SciPy linalg.lstsq ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Calibration curves and response-surface approximations.

Method limitations

Scale coordinates and avoid unnecessarily high degree; outliers and correlated errors require additional treatment.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Data fittingHagen–Poiseuille model ↗

    Predicts fully developed laminar flow in a circular pipe.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingDarcy–Weisbach model ↗

    Relates pipe pressure loss to friction factor and flow speed.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingNon-Newtonian power-law fluid ↗

    Relates shear stress to a power of shear rate.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingBingham plastic model ↗

    Represents a material with a yield stress and post-yield viscosity.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingHerschel–Bulkley model ↗

    Combines yield stress with nonlinear post-yield flow.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingPolynomial chaos expansion ↗

    Represents uncertain responses with polynomial functions of random inputs.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingGeometrically scaled physical model ↗

    Reproduces a system's shape at another size.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingWind-tunnel model ↗

    Uses a controlled air stream around a physical specimen.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingHydraulic flume model ↗

    Uses physical water flow with selected similarity conditions.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingShake-table structural model ↗

    Excites a physical structure with controlled base motion.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingPhotoelastic model ↗

    Uses stress-induced optical birefringence to visualize stress patterns.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingElectrical analog model ↗

    Maps another physical system onto an electrical network.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingHardware-in-the-loop model ↗

    Couples actual hardware to simulated parts of a system.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingDimensional-analysis similarity model ↗

    Uses dimensionless groups to relate tests across scales.

    For fitting a low-dimensional response or constitutive curve; scale variables and validate independently.

  • Data fittingCarnahan-Starling hard-sphere equation of state ↗

    Approximates the compressibility factor of a monodisperse hard-sphere fluid from its packing fraction.

    Fit a limited-range surrogate to EOS evaluations only when repeated calls require it; verify the fit against the explicit formula.

Relationships to other techniques

Approximate & sample → Approximation & reduction

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Approximation & reduction049

Chebyshev approximation

Uses Chebyshev bases and clustered nodes to approximate smooth functions.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

pN(x)=∑k=0NakTk(x)p_N(x)=\sum_{k=0}^Na_kT_k(x)Tk(cos⁡θ)=cos⁡(kθ)T_k(\cos\theta)=\cos(k\theta)

Derivation / construction sketch

  1. Map the domain to the standard interval.
  2. Sample or project using Chebyshev structure.
  3. Truncate the expansion after checking coefficient decay or error.

Symbols & assumptions

Smoothness controls convergence; discontinuities cause Gibbs-type behavior.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Approximation & reduction

Practical use

Accurate one-dimensional property and response surrogates.

Engineering application examples

Accurate one-dimensional property and response surrogates.

Software product or implementation route

Chebfun ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Accurate one-dimensional property and response surrogates.

Method limitations

Smoothness controls convergence; discontinuities cause Gibbs-type behavior.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Smooth surrogateArrhenius rate model ↗

    Relates a rate coefficient to temperature through an activation energy.

    For smooth responses on a bounded interval; check coefficient decay and interpolation error.

Relationships to other techniques

Approximate & sample → Approximation & reduction

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Approximation & reduction050

Proper orthogonal decomposition

Extracts dominant modes from a snapshot matrix.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

X=UΣVTX=U\Sigma V^{\mathsf T}Xr=UrΣrVrTX_r=U_r\Sigma_rV_r^{\mathsf T}

Derivation / construction sketch

  1. Collect representative, consistently scaled snapshots.
  2. Compute their singular value decomposition.
  3. Retain dominant modes and project the dynamics or data onto that subspace.

Symbols & assumptions

Optimal low-rank snapshot error does not guarantee predictive accuracy or stable reduced dynamics.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Approximation & reduction

Practical use

Reduced fluid, structural, and thermal models.

Engineering application examples

Reduced fluid, structural, and thermal models.

Software product or implementation route

pyMOR ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Reduced fluid, structural, and thermal models.

Method limitations

Optimal low-rank snapshot error does not guarantee predictive accuracy or stable reduced dynamics.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Reduced-order modelingNavier–Stokes model ↗

    Conserves mass and momentum for a viscous continuum fluid.

    Build a reduced basis from representative snapshots; check predictive accuracy, conservation, and stability after projection.

Relationships to other techniques

Approximate & sample → Approximation & reduction

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Approximation & reduction051

Dynamic mode decomposition

Fits a linear evolution map between successive snapshots.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

A≈YX+A\approx YX^+A~=UrTYVrΣr−1\widetilde A=U_r^{\mathsf T}YV_r\Sigma_r^{-1}

Derivation / construction sketch

  1. Arrange consecutive snapshots in paired matrices.
  2. Fit a least-squares map from present to future snapshots.
  3. Project the map and analyze its eigenvalues and modes.

Symbols & assumptions

Sampling, noise, rank truncation, and nonlinearity affect interpretation and extrapolation.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Approximation & reduction

Practical use

Flow-pattern analysis and data-driven dynamics prediction.

Engineering application examples

Flow-pattern analysis and data-driven dynamics prediction.

Software product or implementation route

PyDMD ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Flow-pattern analysis and data-driven dynamics prediction.

Method limitations

Sampling, noise, rank truncation, and nonlinearity affect interpretation and extrapolation.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Dynamics analysisDigital twin framework ↗

    Links an evolving model of a specific asset with observations.

    For time-resolved snapshots; modes describe a fitted evolution map and need not capture all nonlinear behavior.

Relationships to other techniques

Approximate & sample → Approximation & reduction

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Approximation & reduction052

Gaussian-process regression

Predicts a response using a covariance model and observed data.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

m∗=k∗T(K+σn2I)−1ym_*=k_*^{\mathsf T}(K+\sigma_n^2I)^{-1}yv∗=k∗∗−k∗T(K+σn2I)−1k∗v_*=k_{**}-k_*^{\mathsf T}(K+\sigma_n^2I)^{-1}k_*

Derivation / construction sketch

  1. Specify a prior mean and covariance kernel.
  2. Condition the joint Gaussian distribution on observations.
  3. Compute the posterior mean and variance at a query point.

Symbols & assumptions

Zero-mean formula shown; uncertainty is conditional on the kernel/noise assumptions, not a guarantee of physical accuracy.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Approximation & reduction

Practical use

Expensive-simulation surrogates and uncertainty-aware design exploration.

Engineering application examples

Expensive-simulation surrogates and uncertainty-aware design exploration.

Software product or implementation route

scikit-learn GaussianProcessRegressor ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Expensive-simulation surrogates and uncertainty-aware design exploration.

Method limitations

Zero-mean formula shown; uncertainty is conditional on the kernel/noise assumptions, not a guarantee of physical accuracy.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Surrogate / uncertaintyDigital twin framework ↗

    Links an evolving model of a specific asset with observations.

    For an expensive-simulation or data surrogate; predictive uncertainty is conditional on kernel and noise assumptions.

  • Surrogate / uncertaintyWind-tunnel model ↗

    Uses a controlled air stream around a physical specimen.

    For an expensive-simulation or data surrogate; predictive uncertainty is conditional on kernel and noise assumptions.

Relationships to other techniques

Approximate & sample → Approximation & reduction

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Quadrature & sampling053

Composite trapezoidal rule

Integrates sampled data by joining neighboring values with straight lines.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

Ih=h[12f0+∑i=1n−1fi+12fn]I_h=h[\tfrac12 f_0+\sum_{i=1}^{n-1}f_i+\tfrac12 f_n]

Derivation / construction sketch

  1. Interpolate each interval linearly.
  2. Integrate that line exactly.
  3. Sum the interval contributions.

Symbols & assumptions

Uniform-spacing formula shown; second-order error for sufficiently smooth functions.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Quadrature & sampling

Practical use

Integrating measured heat-flow, force, or current histories.

Engineering application examples

Integrating measured heat-flow, force, or current histories.

Software product or implementation route

SciPy integrate.trapezoid ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Integrating measured heat-flow, force, or current histories.

Method limitations

Uniform-spacing formula shown; second-order error for sufficiently smooth functions.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Integral / post-processingFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    For sampled loads, fluxes, or response histories with enough resolution; account for nonsmooth events.

Relationships to other techniques

Approximate & sample → Quadrature & sampling

Specific connections

  • Underlies time-stepping formula Crank-Nicolson

    Trapezoidal integration of the time derivative produces the endpoint-average method.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Quadrature & sampling054

Composite Simpson rule

Integrates pairs of intervals using quadratic interpolation.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

Ih=h3[f0+4∑i oddfi+2∑i even, 0<i<nfi+fn]I_h=\frac h3[f_0+4\sum_{i\text{ odd}}f_i+2\sum_{i\text{ even},\,0<i<n}f_i+f_n]

Derivation / construction sketch

  1. Fit a quadratic across each consecutive pair of intervals.
  2. Integrate the polynomial exactly.
  3. Add the repeated weights across the domain.

Symbols & assumptions

Uniform grid with even n shown; fourth-order convergence requires smoothness.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Quadrature & sampling

Practical use

Integrating smooth sampled loads and response curves.

Engineering application examples

Integrating smooth sampled loads and response curves.

Software product or implementation route

SciPy integrate.simpson ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Integrating smooth sampled loads and response curves.

Method limitations

Uniform grid with even n shown; fourth-order convergence requires smoothness.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Integral / post-processingFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    For smooth sampled responses with compatible spacing; do not apply its uniform-grid error order blindly.

  • Integral / post-processingHydraulic flume model ↗

    Uses physical water flow with selected similarity conditions.

    For smooth sampled responses with compatible spacing; do not apply its uniform-grid error order blindly.

Relationships to other techniques

Approximate & sample → Quadrature & sampling

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Quadrature & sampling055

Gaussian quadrature

Chooses nodes and weights to integrate high-degree polynomials efficiently.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

∫−11f(x) dx≈∑i=1nwif(xi)\int_{-1}^1 f(x)\,dx\approx\sum_{i=1}^nw_if(x_i)

Derivation / construction sketch

  1. Choose nodes as roots of the relevant orthogonal polynomial.
  2. Determine weights by exactness conditions.
  3. Map the rule to the integration interval or element.

Symbols & assumptions

Gauss-Legendre is exact through degree 2n-1; singular or nonsmooth integrands need special treatment.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Quadrature & sampling

Practical use

Element stiffness and load integrals in finite-element analysis.

Engineering application examples

Element stiffness and load integrals in finite-element analysis.

Software product or implementation route

SciPy integrate.fixed_quad ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Element stiffness and load integrals in finite-element analysis.

Method limitations

Gauss-Legendre is exact through degree 2n-1; singular or nonsmooth integrands need special treatment.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Integral assemblySchrödinger model ↗

    Evolves a nonrelativistic quantum state using a Hamiltonian.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyDirac model ↗

    Describes relativistic spin-half particles with a spinor wave equation.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyBorn–Oppenheimer approximation ↗

    Separates electronic motion from slower nuclear motion.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyHartree–Fock model ↗

    Approximates a many-electron wavefunction by one self-consistent Slater determinant.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyDensity functional theory (DFT) ↗

    Uses electron density to determine ground-state properties with an approximate exchange-correlation functional.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyTight-binding model ↗

    Represents electronic states with localized orbitals and hopping parameters.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyHubbard model ↗

    Models competition between particle hopping and local electron interactions.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyQuantum harmonic oscillator ↗

    Describes a quantum degree of freedom in a quadratic potential.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyParticle-in-a-box model ↗

    Confines a quantum particle within idealized boundaries.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyRadiative transfer equation ↗

    Tracks radiation intensity through emission, absorption and scattering.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblySurface-to-surface radiosity model ↗

    Balances diffuse radiation exchange between surfaces.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyScalar diffraction model ↗

    Uses a scalar wave approximation for light diffraction.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyLifting-line model ↗

    Approximates finite-wing lift using a spanwise circulation distribution.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyBlade-element momentum model ↗

    Combines blade-section loads with momentum balances.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyNeutron transport model ↗

    Tracks neutron angular and energy-dependent transport with interactions.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyHomogenization ↗

    Derives effective properties or equations from smaller-scale structure.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyRepresentative volume element (RVE) ↗

    Uses a finite microstructural sample to estimate bulk response.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyFE² computational homogenization ↗

    Solves microscale problems within a macroscale finite-element calculation.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyPolynomial chaos expansion ↗

    Represents uncertain responses with polynomial functions of random inputs.

    For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyDebye phonon model ↗

    Approximates acoustic phonons by a continuum spectrum with a mode-count cutoff.

    Use quadrature for the full Debye integral and fitting for a Debye temperature; the cubic law is only a low-temperature asymptote. For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblySommerfeld free-electron model ↗

    Describes conduction electrons as a degenerate, noninteracting Fermi gas.

    Integrate the free-electron density of states with Fermi occupations at finite temperature, or fit a low-temperature heat-capacity coefficient; preserve electron number. For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyTight-binding electronic model ↗

    Builds crystal electronic bands from localized orbitals and intersite hopping.

    Use SVD on H(k)-E I to check null states, condition analysis to assess sensitivity, and quadrature for Brillouin-zone averages. General bands require a Hermitian eigensolver, not SVD singular values interpreted as signed energies. For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

  • Integral assemblyNearly-free-electron model ↗

    Predicts band gaps by perturbing free electrons with a weak periodic potential.

    Use SVD for null-state diagnostics, sensitivity checks for small gaps, and quadrature for band averages. Diagonalize the Hermitian plane-wave Hamiltonian for actual energies. For smooth element, energy, or moment integrals; singular or discontinuous integrands need special rules.

Relationships to other techniques

Approximate & sample → Quadrature & sampling

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Quadrature & sampling056

Adaptive quadrature

Subdivides intervals according to local integration-error estimates.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

I≈∑KQKI\approx\sum_KQ_K∑KeK≤max⁡(εabs,εrel∣I∣)\sum_K e_K\le\max(\varepsilon_{\mathrm{abs}},\varepsilon_{\mathrm{rel}}|I|)

Derivation / construction sketch

  1. Compare paired integration rules to estimate local error.
  2. Refine intervals with the largest estimated contribution.
  3. Stop when the global estimated tolerance is met.

Symbols & assumptions

Estimates can miss singularities or narrow features; supply known breakpoints and examine diagnostics.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Quadrature & sampling

Practical use

Accurate one-dimensional response and probability integrals.

Engineering application examples

Accurate one-dimensional response and probability integrals.

Software product or implementation route

SciPy integrate.quad ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Accurate one-dimensional response and probability integrals.

Method limitations

Estimates can miss singularities or narrow features; supply known breakpoints and examine diagnostics.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Integral evaluationPopulation balance model ↗

    Tracks the distribution of particle sizes or other internal properties.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationIdeal gas equation of state ↗

    Relates pressure, volume and temperature for a dilute noninteracting gas.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationVan der Waals equation of state ↗

    Adds molecular attraction and excluded volume to an ideal gas model.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationPeng–Robinson equation of state ↗

    Uses a cubic equation of state for real-fluid behavior.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationSoave–Redlich–Kwong equation of state ↗

    Uses a temperature-dependent attraction correction in a cubic fluid model.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationVirial equation of state ↗

    Represents nonideal behavior as a density or pressure expansion.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationStefan–Boltzmann surface model ↗

    Relates idealized surface radiant emission to the fourth power of temperature.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationNorton creep law ↗

    Relates creep rate to a power of stress.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationLinear elastic fracture mechanics (LEFM) ↗

    Uses crack-tip intensity parameters in an elastic body.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationCohesive-zone model ↗

    Uses traction-separation relations across a fracture process zone.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationParis fatigue crack-growth law ↗

    Relates cyclic crack-growth rate to stress-intensity-factor range.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationMiner cumulative damage rule ↗

    Adds fractions of fatigue life consumed by load cycles.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationArchard wear model ↗

    Relates wear volume to load, sliding distance and hardness.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationMagnetic-circuit model ↗

    Uses reluctance and magnetomotive force in lumped magnetic paths.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationJiles–Atherton hysteresis model ↗

    Represents path-dependent magnetization with phenomenological parameters.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationGeometrical optics ↗

    Approximates light propagation as rays.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationScalar diffraction model ↗

    Uses a scalar wave approximation for light diffraction.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationGaussian beam model ↗

    Represents a paraxial beam with a Gaussian transverse profile.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationDrude–Lorentz optical model ↗

    Represents free-carrier and bound-charge contributions to permittivity.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationStellar structure model ↗

    Couples hydrostatic balance, energy transport and energy generation.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationFLRW cosmological model ↗

    Assumes a homogeneous and isotropic expanding spacetime.

    For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationOrnstein-Zernike equation ↗

    Relates total and direct pair correlations in a homogeneous liquid, linking microscopic structure to scattering.

    Evaluate radial correlation integrals or Fourier-Bessel transforms with controlled truncation and oscillatory-integration error.

  • Integral evaluationPercus-Yevick closure ↗

    Closes the liquid integral equation using an approximate relation between pair correlations and interactions.

    Evaluate radial correlation integrals or Fourier-Bessel transforms with controlled truncation and oscillatory-integration error.

  • Integral evaluationHypernetted-chain (HNC) closure ↗

    Approximates liquid pair structure by neglecting bridge diagrams in the exact closure.

    Evaluate radial correlation integrals or Fourier-Bessel transforms with controlled truncation and oscillatory-integration error.

  • Integral evaluationGreen-Kubo viscosity relation ↗

    Obtains equilibrium shear viscosity from the time integral of microscopic shear-stress fluctuations.

    Integrate a smooth fitted stress-autocorrelation function; account separately for finite sampling and the unobserved long-time tail.

  • Integral evaluationDebye phonon model ↗

    Approximates acoustic phonons by a continuum spectrum with a mode-count cutoff.

    Use quadrature for the full Debye integral and fitting for a Debye temperature; the cubic law is only a low-temperature asymptote. For deterministic low-dimensional integrals; identify singularities and verify error estimates.

  • Integral evaluationSommerfeld free-electron model ↗

    Describes conduction electrons as a degenerate, noninteracting Fermi gas.

    Integrate the free-electron density of states with Fermi occupations at finite temperature, or fit a low-temperature heat-capacity coefficient; preserve electron number. For deterministic low-dimensional integrals; identify singularities and verify error estimates.

Relationships to other techniques

Approximate & sample → Quadrature & sampling

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Quadrature & sampling057

Monte Carlo integration

Estimates an integral by averaging independent random samples.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

I^N=1N∑i=1Nf(Xi)\widehat I_N=\frac1N\sum_{i=1}^Nf(X_i)SE⁡(I^N)≈sN\operatorname{SE}(\widehat I_N)\approx\frac s{\sqrt N}

Derivation / construction sketch

  1. Express the target as an expectation under a chosen sampling distribution.
  2. Generate independent samples and average the integrand.
  3. Estimate uncertainty from sample variability.

Symbols & assumptions

Finite-variance independent-sample form; dimension-independent rate can still have a large variance constant.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Quadrature & sampling

Practical use

Uncertain-load propagation and probabilistic engineering estimates.

Engineering application examples

Uncertain-load propagation and probabilistic engineering estimates.

Software product or implementation route

NumPy Generator ↗

Sampling building block; construct the integrand estimator and its uncertainty analysis in application code.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Uncertain-load propagation and probabilistic engineering estimates.

Method limitations

Finite-variance independent-sample form; dimension-independent rate can still have a large variance constant.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Statistical estimationHeisenberg spin model ↗

    Represents interacting localized magnetic moments.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationIsing model ↗

    Represents discrete spins with interaction energies.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationClassical molecular dynamics (MD) ↗

    Integrates atomic motion under specified interaction forces.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationAb initio molecular dynamics ↗

    Computes interatomic forces from electronic-structure calculations during motion.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationLennard–Jones potential ↗

    Combines short-range repulsion with an inverse-sixth-power attraction.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationMorse potential ↗

    Represents an anharmonic bond with a finite dissociation energy.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationEmbedded-atom method (EAM) ↗

    Combines pair interactions with an embedding energy dependent on local electron density.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationModified embedded-atom method (MEAM) ↗

    Extends embedding models with angular information.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationTersoff bond-order potential ↗

    Makes bond strength depend on the local bonding environment.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationStillinger–Weber potential ↗

    Uses two-body and three-body terms to favor local tetrahedral structure.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationReaxFF reactive force field ↗

    Uses variable bond orders and charge equilibration to represent chemical reactions.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationAMBER force-field family ↗

    Uses parameterized bonded and nonbonded interactions for biomolecules.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationCHARMM force-field family ↗

    Models biomolecular interactions with chemistry-specific parameter sets.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationOPLS force-field family ↗

    Uses parameterized molecular interactions developed for condensed phases.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationSPC/E water model ↗

    Approximates water using a rigid three-site classical model.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationTIP4P water-model family ↗

    Uses a four-site geometry with an off-oxygen charge site.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationDrude polarizable model ↗

    Uses auxiliary charged particles to represent induced polarization.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationMachine-learned interatomic potential ↗

    Fits atomic energies and forces from reference data using statistical learning.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationCoarse-grained molecular model ↗

    Groups atoms into effective interaction sites.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationMartini coarse-grained model ↗

    Uses mapped molecular beads and parameterized interactions.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationDissipative particle dynamics (DPD) ↗

    Combines conservative, dissipative and random pair forces.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationBrownian dynamics ↗

    Uses overdamped stochastic motion for particles in a surrounding medium.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationLangevin dynamics ↗

    Adds friction and random forces to a dynamical model.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationKinetic Monte Carlo ↗

    Samples transitions between states using event rates.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationPotts grain-growth model ↗

    Represents grain orientations as discrete lattice states.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationRadiative transfer equation ↗

    Tracks radiation intensity through emission, absorption and scattering.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationDiscrete-event simulation ↗

    Advances a system through scheduled events.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationAgent-based physical-system model ↗

    Represents interacting entities following local rules.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationMarkov state model ↗

    Represents probabilistic transitions between a finite set of states.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationNeutron transport model ↗

    Tracks neutron angular and energy-dependent transport with interactions.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationSmoothed particle hydrodynamics (SPH) ↗

    Approximates continuum fields through moving particles and kernels.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationDiscrete element method (DEM) ↗

    Evolves contacting discrete bodies with contact laws.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationLattice Boltzmann method (LBM) ↗

    Evolves discrete velocity populations to recover suitable macroscopic flow equations.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationDirect simulation Monte Carlo (DSMC) ↗

    Samples particle motion and collisions in a rarefied gas.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationParticle-in-cell (PIC) ↗

    Couples moving computational particles to fields on a mesh.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationMaterial point method (MPM) ↗

    Transfers particle-carried material state to a computational grid.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationMonte Carlo transport ↗

    Samples particle histories and interactions statistically.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationPhysics-informed neural network (PINN) ↗

    Trains a neural approximation using data and governing-equation residuals.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationNeural operator ↗

    Learns a map between function-valued inputs and outputs.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationGeometrically scaled physical model ↗

    Reproduces a system's shape at another size.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationWind-tunnel model ↗

    Uses a controlled air stream around a physical specimen.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationHydraulic flume model ↗

    Uses physical water flow with selected similarity conditions.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationShake-table structural model ↗

    Excites a physical structure with controlled base motion.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationPhotoelastic model ↗

    Uses stress-induced optical birefringence to visualize stress patterns.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationElectrical analog model ↗

    Maps another physical system onto an electrical network.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationHardware-in-the-loop model ↗

    Couples actual hardware to simulated parts of a system.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationDimensional-analysis similarity model ↗

    Uses dimensionless groups to relate tests across scales.

    For expectation or uncertainty calculations with a stated sampling law; this does not replace a specialized stochastic time integrator.

  • Statistical estimationStokes-Einstein diffusion relation ↗

    Connects Brownian translational diffusion to temperature, solvent viscosity, and hydrodynamic particle radius.

    Propagate uncertainty in temperature, viscosity, and radius through D=kBT/(6πηR); no numerical integration is needed for the nominal value.

Relationships to other techniques

Approximate & sample → Quadrature & sampling

Specific connections

  • Can use variance reduction from Importance sampling

    Change the proposal and correct with density-ratio weights.

  • Has low-discrepancy alternative Quasi-Monte Carlo

    Both average integrand values, but their sampling designs and uncertainty analyses differ.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Quadrature & sampling058

Quasi-Monte Carlo

Uses low-discrepancy points to cover an integration domain evenly.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

I≈1N∑i=1Nf(ui)I\approx\frac1N\sum_{i=1}^Nf(u_i)ui∈[0,1]du_i\in[0,1]^d

Derivation / construction sketch

  1. Map the integral to a unit cube.
  2. Generate a low-discrepancy design such as Sobol points.
  3. Average values; use independent randomized scramblings for error assessment.

Symbols & assumptions

Performance depends on smoothness and effective dimension; deterministic points do not give an IID standard error.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Quadrature & sampling

Practical use

High-dimensional uncertainty propagation and design integration.

Engineering application examples

High-dimensional uncertainty propagation and design integration.

Software product or implementation route

SciPy integrate.qmc_quad ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

High-dimensional uncertainty propagation and design integration.

Method limitations

Performance depends on smoothness and effective dimension; deterministic points do not give an IID standard error.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Uncertainty integrationDiscrete-event simulation ↗

    Advances a system through scheduled events.

    For well-behaved parameter integrals where low-discrepancy coverage helps; use randomized replicates for uncertainty assessment.

  • Uncertainty integrationAgent-based physical-system model ↗

    Represents interacting entities following local rules.

    For well-behaved parameter integrals where low-discrepancy coverage helps; use randomized replicates for uncertainty assessment.

  • Uncertainty integrationMonte Carlo transport ↗

    Samples particle histories and interactions statistically.

    For well-behaved parameter integrals where low-discrepancy coverage helps; use randomized replicates for uncertainty assessment.

  • Uncertainty integrationBayesian model calibration ↗

    Updates uncertain parameters using observations and a statistical likelihood.

    For well-behaved parameter integrals where low-discrepancy coverage helps; use randomized replicates for uncertainty assessment.

  • Uncertainty integrationPolynomial chaos expansion ↗

    Represents uncertain responses with polynomial functions of random inputs.

    For well-behaved parameter integrals where low-discrepancy coverage helps; use randomized replicates for uncertainty assessment.

Relationships to other techniques

Approximate & sample → Quadrature & sampling

Specific connections

  • Low-discrepancy alternative to Monte Carlo integration

    Both average integrand values, but their sampling designs and uncertainty analyses differ.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Quadrature & sampling059

Importance sampling

Changes the sampling distribution to focus on influential regions.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

I=∫f(x)p(x) dx=Eq[f(X)p(X)/q(X)]I=\int f(x)p(x)\,dx=\mathbb E_q[f(X)p(X)/q(X)]

Derivation / construction sketch

  1. Choose a proposal density covering the target contribution.
  2. Sample from that proposal.
  3. Weight each value by the target-to-proposal density ratio.

Symbols & assumptions

Support coverage and finite weight variance are essential; extreme weights can destroy efficiency.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Quadrature & sampling

Practical use

Rare-event estimation and reliability calculations.

Engineering application examples

Rare-event estimation and reliability calculations.

Software product or implementation route

Custom importance-sampling estimator ↗

Implementation route: use a sampling library and compute density-ratio weights; verify support and weight diagnostics. The review explains the construction and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Rare-event estimation and reliability calculations.

Method limitations

Support coverage and finite weight variance are essential; extreme weights can destroy efficiency.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Rare-event / expectation estimationCoarse-grained molecular model ↗

    Groups atoms into effective interaction sites.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationMartini coarse-grained model ↗

    Uses mapped molecular beads and parameterized interactions.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationDissipative particle dynamics (DPD) ↗

    Combines conservative, dissipative and random pair forces.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationBrownian dynamics ↗

    Uses overdamped stochastic motion for particles in a surrounding medium.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationLangevin dynamics ↗

    Adds friction and random forces to a dynamical model.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationKinetic Monte Carlo ↗

    Samples transitions between states using event rates.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationPotts grain-growth model ↗

    Represents grain orientations as discrete lattice states.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationDiscrete-event simulation ↗

    Advances a system through scheduled events.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationAgent-based physical-system model ↗

    Represents interacting entities following local rules.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationNeutron transport model ↗

    Tracks neutron angular and energy-dependent transport with interactions.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationMonte Carlo transport ↗

    Samples particle histories and interactions statistically.

    For a known target and proposal with correct support and controlled weight variance.

  • Rare-event / expectation estimationBayesian model calibration ↗

    Updates uncertain parameters using observations and a statistical likelihood.

    For a known target and proposal with correct support and controlled weight variance.

Relationships to other techniques

Approximate & sample → Quadrature & sampling

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Quadrature & sampling060

Metropolis-Hastings sampling

Builds a Markov chain with a desired stationary density.

Approximate & sampleNumerical technique
Formulation & short derivation

Representative numerical formulation

α(x,y)=min⁡(1,π(y)q(x∣y)π(x)q(y∣x))\alpha(x,y)=\min\left(1,\frac{\pi(y)q(x\mid y)}{\pi(x)q(y\mid x)}\right)

Derivation / construction sketch

  1. Propose a move using a chosen transition distribution.
  2. Accept according to the target/proposal ratio.
  3. Repeat and diagnose mixing and convergence.

Symbols & assumptions

Samples are correlated; initialization, ergodicity, effective sample size, and multimodality matter.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Quadrature & sampling

Practical use

Bayesian inverse problems and posterior uncertainty exploration.

Engineering application examples

Bayesian inverse problems and posterior uncertainty exploration.

Software product or implementation route

PyMC Metropolis ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Bayesian inverse problems and posterior uncertainty exploration.

Method limitations

Samples are correlated; initialization, ergodicity, effective sample size, and multimodality matter.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Equilibrium / posterior samplingHeisenberg spin model ↗

    Represents interacting localized magnetic moments.

    For a specified target distribution; diagnose mixing and correlation. Samples do not generally represent physical time.

  • Equilibrium / posterior samplingIsing model ↗

    Represents discrete spins with interaction energies.

    For a specified target distribution; diagnose mixing and correlation. Samples do not generally represent physical time.

  • Equilibrium / posterior samplingBayesian model calibration ↗

    Updates uncertain parameters using observations and a statistical likelihood.

    For a specified target distribution; diagnose mixing and correlation. Samples do not generally represent physical time.

Relationships to other techniques

Approximate & sample → Quadrature & sampling

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Error analysis & verification061

Richardson extrapolation

Cancels a leading discretization-error term using two resolutions.

Verify & assessVerification / analysis
Formulation & short derivation

Representative numerical formulation

u∗≈uh/r+uh/r−uhrp−1u_*\approx u_{h/r}+\frac{u_{h/r}-u_h}{r^p-1}

Derivation / construction sketch

  1. Assume a leading error term proportional to h to the power p.
  2. Write that expansion at two consistently refined resolutions.
  3. Eliminate the leading coefficient.

Symbols & assumptions

Requires an asymptotic convergence regime and a valid order estimate; shocks and inconsistent grids can invalidate the assumption.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Error analysis & verification

Practical use

Checking mesh and time-step convergence.

Engineering application examples

Checking mesh and time-step convergence.

Software product or implementation route

Post-processing refinement study ↗

Implementation route: compare systematically refined solver runs and evaluate the error model; no particular solver is certified by the estimate.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Checking mesh and time-step convergence.

Method limitations

Requires an asymptotic convergence regime and a valid order estimate; shocks and inconsistent grids can invalidate the assumption.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • VerificationTime-dependent DFT (TDDFT) ↗

    Evolves electron density to approximate excited-state response.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationDissipative particle dynamics (DPD) ↗

    Combines conservative, dissipative and random pair forces.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationBrownian dynamics ↗

    Uses overdamped stochastic motion for particles in a surrounding medium.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationLangevin dynamics ↗

    Adds friction and random forces to a dynamical model.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationKinetic Monte Carlo ↗

    Samples transitions between states using event rates.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationPotts grain-growth model ↗

    Represents grain orientations as discrete lattice states.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationDiscrete dislocation dynamics ↗

    Tracks line defects and their interactions.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationGeometrical optics ↗

    Approximates light propagation as rays.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationPoint reactor kinetics ↗

    Approximates time-dependent neutron population with delayed-neutron groups.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationBateman decay-chain model ↗

    Evolves coupled radioactive parent and daughter populations.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationNewtonian gravitational N-body model ↗

    Evolves masses under mutual inverse-square attraction.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationFinite element method (FEM / FEA) ↗

    Approximates fields with basis functions over elements.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationFinite volume method (FVM) ↗

    Discretizes conservation laws using fluxes across control-volume boundaries.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationFinite difference method (FDM) ↗

    Approximates derivatives with differences on a grid.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationBoundary element method (BEM) ↗

    Recasts suitable field problems as boundary integral equations.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationSpectral method ↗

    Represents fields with global or element-wise high-order basis expansions.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationSmoothed particle hydrodynamics (SPH) ↗

    Approximates continuum fields through moving particles and kernels.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationDiscrete element method (DEM) ↗

    Evolves contacting discrete bodies with contact laws.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationLattice Boltzmann method (LBM) ↗

    Evolves discrete velocity populations to recover suitable macroscopic flow equations.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationDirect numerical simulation (DNS) ↗

    Resolves turbulence without a turbulence closure for the selected flow equations.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationDirect simulation Monte Carlo (DSMC) ↗

    Samples particle motion and collisions in a rarefied gas.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationParticle-in-cell (PIC) ↗

    Couples moving computational particles to fields on a mesh.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationFinite-difference time-domain (FDTD) ↗

    Advances discretized electromagnetic fields in time.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationMaterial point method (MPM) ↗

    Transfers particle-carried material state to a computational grid.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationQM/MM coupling ↗

    Combines quantum mechanics in a selected region with molecular mechanics around it.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationAtomistic–continuum coupling ↗

    Connects particle-level and continuum descriptions.

    For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationGreen-Kubo viscosity relation ↗

    Obtains equilibrium shear viscosity from the time integral of microscopic shear-stress fluctuations.

    Check time-step or correlation-integration refinement where a regular error expansion holds; it does not remove statistical trajectory noise.

  • VerificationEinstein crystal heat-capacity model ↗

    Treats crystal vibrations as independent quantum oscillators at a single frequency.

    Use fitting for an Einstein temperature, interpolation for tabulated responses, or refinement to check derived quantities; the displayed formula itself is explicit. For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

  • VerificationHarmonic lattice dynamics ↗

    Computes phonon modes from a quadratic expansion of crystal potential energy.

    Fit force constants with SVD, cross-check real-time harmonic motion with Verlet, and test displacement/time-step refinement. Obtain phonon frequencies with a Hermitian dynamical-matrix eigensolver. For systematically refined computations in an established asymptotic error regime; use consistent geometry and boundary data.

Relationships to other techniques

Verify & assess → Error analysis & verification

Specific connections

  • Provides error model for Grid convergence index

    The Richardson-style estimate is multiplied by a documented safety factor.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Error analysis & verification062

Grid convergence index

Reports a safety-factored estimate of discretization uncertainty.

Verify & assessVerification / analysis
Formulation & short derivation

Representative numerical formulation

GCIfine=Fs∣(uf−uc)/uf∣rp−1\mathrm{GCI}_{\mathrm{fine}}=F_s\frac{|(u_f-u_c)/u_f|}{r^p-1}

Derivation / construction sketch

  1. Estimate the observed convergence order using a systematic refinement study.
  2. Use a Richardson-style fine-grid error estimate.
  3. Apply the chosen safety factor and report assumptions.

Symbols & assumptions

Relative form fails near zero output; nonmonotonic convergence and coupled grid/time errors need special analysis.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Error analysis & verification

Practical use

Documenting numerical uncertainty in CFD outputs.

Engineering application examples

Documenting numerical uncertainty in CFD outputs.

Software product or implementation route

Grid-convergence post-processing ↗

Implementation route: compute the index from a documented systematic refinement study; the report explains numerical verification practices.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Documenting numerical uncertainty in CFD outputs.

Method limitations

Relative form fails near zero output; nonmonotonic convergence and coupled grid/time errors need special analysis.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Discretization uncertaintyNavier–Stokes model ↗

    Conserves mass and momentum for a viscous continuum fluid.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyEuler flow model ↗

    Neglects viscous stresses in compressible or incompressible flow.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyStokes creeping-flow model ↗

    Neglects inertial terms relative to viscosity.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyBoundary-layer model ↗

    Resolves thin near-wall regions with scale-based simplifications.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyLubrication approximation ↗

    Simplifies viscous flow in thin gaps.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyOldroyd-B model ↗

    Combines solvent viscosity with an elastic polymer stress.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintySaint-Venant shallow-water model ↗

    Depth-averages mass and momentum in free-surface flow.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyKinematic-wave routing ↗

    Simplifies flow routing by approximating dominant slope and friction balance.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyGroundwater flow model ↗

    Combines water conservation with porous-flow relations.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyAdvection–dispersion groundwater model ↗

    Represents contaminant transport and spreading through an aquifer.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyNumerical weather prediction ↗

    Evolves atmospheric dynamics and thermodynamics from an analyzed initial state.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyGeneral circulation model (GCM) ↗

    Represents large-scale atmospheric or oceanic circulation.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyEarth system model (ESM) ↗

    Couples atmosphere, ocean, land, ice and biogeochemical processes.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyOcean circulation model ↗

    Evolves ocean momentum, temperature and salinity.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintySea-ice thermodynamic-dynamic model ↗

    Couples freezing, melting and ice motion.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

  • Discretization uncertaintyElastic seismic-wave model ↗

    Propagates elastic disturbances through Earth materials.

    For a systematic grid-refinement study with a justified observed order and safety factor; it is not physical validation.

Relationships to other techniques

Verify & assess → Error analysis & verification

Specific connections

  • Builds uncertainty estimate from Richardson extrapolation

    The Richardson-style estimate is multiplied by a documented safety factor.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Error analysis & verification063

Method of manufactured solutions

Tests a PDE implementation using a constructed exact solution.

Verify & assessVerification / analysis
Formulation & short derivation

Representative numerical formulation

Lu=fLu=fu=um ⇒ fm=Lumu=u_m\ \Rightarrow\ f_m=Lu_m

Derivation / construction sketch

  1. Choose a smooth analytic manufactured solution.
  2. Substitute it into the governing operator to compute forcing and boundary data.
  3. Solve and verify that error decreases at the expected rate.

Symbols & assumptions

Tests implementation accuracy, not whether the governing physical model describes reality.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Error analysis & verification

Practical use

Code verification for PDE and multiphysics solvers.

Engineering application examples

Code verification for PDE and multiphysics solvers.

Software product or implementation route

MOOSE MMS tools ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Code verification for PDE and multiphysics solvers.

Method limitations

Tests implementation accuracy, not whether the governing physical model describes reality.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Code verificationCahn–Hilliard model ↗

    Evolves a conserved composition field through chemical-potential gradients.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationAllen–Cahn model ↗

    Evolves a nonconserved order parameter toward lower free energy.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationPhase-field crystal model ↗

    Uses a periodic density-like field to represent crystalline ordering.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationFickian diffusion ↗

    Relates diffusive flux to concentration gradients.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationMaxwell–Stefan diffusion ↗

    Represents multicomponent diffusion through interspecies friction.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationAdvection–diffusion–reaction model ↗

    Combines bulk transport, diffusion and reaction sources.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationPennes bioheat model ↗

    Adds perfusion and metabolic heat to tissue heat transfer.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationReaction–diffusion morphogenesis model ↗

    Couples reacting substances with diffusion.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationFinite element method (FEM / FEA) ↗

    Approximates fields with basis functions over elements.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationFinite volume method (FVM) ↗

    Discretizes conservation laws using fluxes across control-volume boundaries.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationFinite difference method (FDM) ↗

    Approximates derivatives with differences on a grid.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationBoundary element method (BEM) ↗

    Recasts suitable field problems as boundary integral equations.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationSpectral method ↗

    Represents fields with global or element-wise high-order basis expansions.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationDirect numerical simulation (DNS) ↗

    Resolves turbulence without a turbulence closure for the selected flow equations.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

  • Code verificationFinite-difference time-domain (FDTD) ↗

    Advances discretized electromagnetic fields in time.

    For an accessible differential operator, construct compatible forcing and boundaries; this tests implementation rather than physical realism.

Relationships to other techniques

Verify & assess → Error analysis & verification

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Error analysis & verification064

Residual-based error estimation

Uses equation and interface residuals to guide error assessment.

Verify & assessVerification / analysis
Formulation & short derivation

Representative numerical formulation

ηK2=hK2∥f+∇⋅(k∇uh)∥K2+∑e⊂∂Khe∥Je∥e2\eta_K^2=h_K^2\lVert f+\nabla\cdot(k\nabla u_h)\rVert_K^2+\sum_{e\subset\partial K}h_e\lVert J_e\rVert_e^2

Derivation / construction sketch

  1. Measure how strongly the approximate solution violates the PDE inside cells.
  2. Add flux-jump contributions on interfaces.
  3. Combine with problem-dependent weights and refine where indicators are large.

Symbols & assumptions

Schematic elliptic estimator; reliability constants and boundary terms depend on assumptions and discretization.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Error analysis & verification

Practical use

Adaptive refinement in finite-element simulation.

Engineering application examples

Adaptive refinement in finite-element simulation.

Software product or implementation route

MFEM error-estimator/adaptivity framework ↗

MFEM supports estimator-driven adaptivity; the displayed residual form is schematic and is not a claim about the exact estimator selected by an example.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Adaptive refinement in finite-element simulation.

Method limitations

Schematic elliptic estimator; reliability constants and boundary terms depend on assumptions and discretization.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Spatial verificationFourier heat conduction ↗

    Relates conductive heat flux to temperature gradient.

    For a suitable PDE discretization with an estimator derived for its operator; a small solver residual alone is not a full error bound.

Relationships to other techniques

Verify & assess → Error analysis & verification

Specific connections

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Error analysis & verification065

Von Neumann stability analysis

Tests Fourier-mode amplification for linear grid schemes.

Verify & assessVerification / analysis
Formulation & short derivation

Representative numerical formulation

ujn=Gneijθu_j^n=G^n e^{ij\theta}∣G(θ)∣≤1|G(\theta)|\le1

Derivation / construction sketch

  1. Insert a Fourier mode into a linear constant-coefficient difference scheme.
  2. Solve for its amplification factor.
  3. Require bounded amplification for every resolvable wave number.

Symbols & assumptions

Typically periodic or infinite uniform grids; boundary effects and nonlinear stability need separate analysis.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Error analysis & verification

Practical use

Selecting stable explicit diffusion and wave time steps.

Engineering application examples

Selecting stable explicit diffusion and wave time steps.

Software product or implementation route

Symbolic or scripted stability analysis ↗

Implementation route: derive the amplification factor and evaluate its magnitude over wave numbers; this is an analysis procedure rather than a solver product.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Selecting stable explicit diffusion and wave time steps.

Method limitations

Typically periodic or infinite uniform grids; boundary effects and nonlinear stability need separate analysis.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Scheme analysisTransient heat equation ↗

    Balances thermal storage, conduction and heat sources.

    For a linearized constant-coefficient uniform-grid subproblem; boundaries and nonlinear effects need separate checks.

Relationships to other techniques

Verify & assess → Error analysis & verification

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Error analysis & verification066

Condition-number analysis

Measures how perturbations in inputs can affect a computed solution.

Verify & assessVerification / analysis
Formulation & short derivation

Representative numerical formulation

κ(A)=∥A∥∥A−1∥\kappa(A)=\lVert A\rVert\lVert A^{-1}\rVert∥δx∥∥x∥≤κ(A)∥δb∥∥b∥\frac{\lVert\delta x\rVert}{\lVert x\rVert}\le\kappa(A)\frac{\lVert\delta b\rVert}{\lVert b\rVert}

Derivation / construction sketch

  1. Perturb a nonsingular linear system with fixed A.
  2. Bound the solution change using operator norms.
  3. Compare relative input and output perturbations.

Symbols & assumptions

The bound shown concerns right-hand-side perturbations; conditioning is a problem property, distinct from algorithm stability.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Error analysis & verification

Practical use

Diagnosing sensitive inverse problems and poorly scaled matrix systems.

Engineering application examples

Diagnosing sensitive inverse problems and poorly scaled matrix systems.

Software product or implementation route

NumPy linalg.cond ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Diagnosing sensitive inverse problems and poorly scaled matrix systems.

Method limitations

The bound shown concerns right-hand-side perturbations; conditioning is a problem property, distinct from algorithm stability.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Sensitivity diagnosisTight-binding model ↗

    Represents electronic states with localized orbitals and hopping parameters.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisHubbard model ↗

    Models competition between particle hopping and local electron interactions.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisHeisenberg spin model ↗

    Represents interacting localized magnetic moments.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisIsing model ↗

    Represents discrete spins with interaction energies.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisAC power-flow model ↗

    Balances complex power on an electrical network.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisDC power-flow approximation ↗

    Linearizes active-power flow under restrictive grid assumptions.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisState-space model ↗

    Represents system evolution with internal states, inputs and outputs.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisTransfer-function model ↗

    Relates linear time-invariant input and output in the transform domain.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisBond-graph model ↗

    Represents energy exchange across mechanical, electrical and other domains.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisSystem-dynamics stock-flow model ↗

    Represents accumulated quantities and their rates of change.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisMarkov state model ↗

    Represents probabilistic transitions between a finite set of states.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisKalman state estimator ↗

    Combines a dynamical model with noisy observations using covariance updates.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisFinite element method (FEM / FEA) ↗

    Approximates fields with basis functions over elements.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisFinite volume method (FVM) ↗

    Discretizes conservation laws using fluxes across control-volume boundaries.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisFinite difference method (FDM) ↗

    Approximates derivatives with differences on a grid.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisBoundary element method (BEM) ↗

    Recasts suitable field problems as boundary integral equations.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisSpectral method ↗

    Represents fields with global or element-wise high-order basis expansions.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisSmoothed particle hydrodynamics (SPH) ↗

    Approximates continuum fields through moving particles and kernels.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisDiscrete element method (DEM) ↗

    Evolves contacting discrete bodies with contact laws.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisLattice Boltzmann method (LBM) ↗

    Evolves discrete velocity populations to recover suitable macroscopic flow equations.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisDirect numerical simulation (DNS) ↗

    Resolves turbulence without a turbulence closure for the selected flow equations.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisDirect simulation Monte Carlo (DSMC) ↗

    Samples particle motion and collisions in a rarefied gas.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisParticle-in-cell (PIC) ↗

    Couples moving computational particles to fields on a mesh.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisFinite-difference time-domain (FDTD) ↗

    Advances discretized electromagnetic fields in time.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisMaterial point method (MPM) ↗

    Transfers particle-carried material state to a computational grid.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisGaussian-process surrogate ↗

    Predicts responses with a probabilistic function model fitted to samples.

    For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisTight-binding electronic model ↗

    Builds crystal electronic bands from localized orbitals and intersite hopping.

    Use SVD on H(k)-E I to check null states, condition analysis to assess sensitivity, and quadrature for Brillouin-zone averages. General bands require a Hermitian eigensolver, not SVD singular values interpreted as signed energies. For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

  • Sensitivity diagnosisNearly-free-electron model ↗

    Predicts band gaps by perturbing free electrons with a weak periodic potential.

    Use SVD for null-state diagnostics, sensitivity checks for small gaps, and quadrature for band averages. Diagonalize the Hermitian plane-wave Hamiltonian for actual energies. For assembled linear systems or fitting matrices; separate problem conditioning from algorithm stability.

Relationships to other techniques

Verify & assess → Error analysis & verification

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Error analysis & verification067

Complex-step differentiation

Estimates an analytic derivative without real subtractive cancellation.

Verify & assessVerification / analysis
Formulation & short derivation

Representative numerical formulation

f′(x)≈Im⁡f(x+ih)hf'(x)\approx\frac{\operatorname{Im}f(x+ih)}h

Derivation / construction sketch

  1. Expand an analytic function in a complex Taylor series.
  2. Extract the imaginary part.
  3. Divide by the perturbation to recover the derivative with second-order truncation error.

Symbols & assumptions

Requires holomorphic operations and complex-compatible code; absolute values, branches, or discarded imaginary parts can break it.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Error analysis & verification

Practical use

Checking gradients in smooth engineering optimization code.

Engineering application examples

Checking gradients in smooth engineering optimization code.

Software product or implementation route

SciPy optimize.least_squares ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Checking gradients in smooth engineering optimization code.

Method limitations

Requires holomorphic operations and complex-compatible code; absolute values, branches, or discarded imaginary parts can break it.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Derivative verificationPacejka tire model ↗

    Uses empirical nonlinear formulas for tire forces.

    Only for smooth analytic, complex-compatible code paths; clipping, absolute values, and phase switches can invalidate it.

Relationships to other techniques

Verify & assess → Error analysis & verification

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider
Error analysis & verification068

Adjoint sensitivity analysis

Computes gradients of scalar outputs with respect to many parameters.

Verify & assessVerification / analysis
Formulation & short derivation

Representative numerical formulation

R(u,p)=0R(u,p)=0RuTλ=JuTR_u^{\mathsf T}\lambda=J_u^{\mathsf T}dJdp=Jp−λTRp\frac{dJ}{dp}=J_p-\lambda^{\mathsf T}R_p

Derivation / construction sketch

  1. Differentiate the state residual and objective.
  2. Introduce an adjoint variable to eliminate the expensive state sensitivity.
  3. Solve one adjoint system per scalar objective and assemble parameter gradients.

Symbols & assumptions

Discrete formulation shown; consistent boundary conditions, differentiation, and solver tolerances are essential.

The sketch summarizes algorithm construction, not a complete convergence proof or implementation. Consult the references for stopping criteria, stability requirements, and variants.

Practical applications & products

Application area

Error analysis & verification

Practical use

Aerodynamic shape optimization and inverse design.

Engineering application examples

Aerodynamic shape optimization and inverse design.

Software product or implementation route

dolfin-adjoint ↗

Named software documentation describes this method or a directly relevant implementation component; verify the selected routine and its assumptions.

Application examples illustrate where a technique is useful. They do not assert an undocumented manufacturer workflow or certify a software implementation.

Example, limitations & references

In practice

Aerodynamic shape optimization and inverse design.

Method limitations

Discrete formulation shown; consistent boundary conditions, differentiation, and solver tolerances are essential.

References & further reading

Suitable physical-model applications

Conditional starting points, not automatic solver selections. Check the actual equations, constraints, stiffness, and matrix structure. Some linked atlas entries are methods or frameworks rather than physical laws. See the linked entries’ assumptions and references.

  • Many-parameter gradientsNavier–Stokes model ↗

    Conserves mass and momentum for a viscous continuum fluid.

    For a differentiable discretized state problem and scalar objectives; use consistent derivatives, boundary conditions, and solver tolerances.

  • Many-parameter gradientsModel predictive control ↗

    Optimizes future actions using a predictive model and constraints.

    For a differentiable discretized state problem and scalar objectives; use consistent derivatives, boundary conditions, and solver tolerances.

Relationships to other techniques

Verify & assess → Error analysis & verification

Specific connections

  • Can supply gradients to Newton optimization

    An adjoint supplies objective sensitivities; Hessian information requires additional work.

Other entries in this discipline

Editorial connections and subject grouping, not a complete dependency graph. Assumptions matter; consult the linked entries’ formulations and references.

↑ Find this technique in the tree
Search Google ↑ Go back to the slider

FROM A QUESTION TO A MODEL

Choose the method that fits
your equations and error budget.

01

Define the outcome

Choose the quantities, length and time scales, operating conditions, and accuracy you need. Distinguish discretization error, algebraic error, roundoff, and uncertainty in the physical inputs.

02

Check the assumptions

Check matrix structure, stiffness, smoothness, constraints, conservation, and available derivative information before selecting a numerical method.

03

Verify, then validate

Check units, conservation, numerical convergence and limiting cases. Compare predictions with independent measurements and report uncertainty and model limitations.

ABOUT THIS ATLAS

A starting point for discovery.

Numerical Modeling brings together discretizations, solvers, approximations, sampling techniques, and verification tools. A numerical method approximates or solves a mathematical problem; it does not establish whether the underlying physical model describes reality.

This release includes 68 entries across 8 subject areas. Methods may have many variants, and software implementations can differ. A converged iteration is not by itself proof of discretization accuracy or physical validity.

References at the point of use

Open any entry’s “Example, limitations & references” section for linked papers, author-written textbooks, or official technical documentation. A technical manual can support several related entries; use the model name to find the relevant section. Some publisher-hosted papers require access.

Descriptions and examples are concise editorial summaries. This educational catalog is not a simulation service, engineering certification, or substitute for validating a design. Follow the licensing and citation requirements of each original source.

Catalog release: September 23, 2026 · Version 1.0

THE NUMERICAL MODELING COMMUNITY

Cumulative visitor map

84Site visitors
Distinct daily visits, cumulative

Counting since 2026-09-23 (UTC). As on To the Infinity, one IP address counts once per UTC day; this is not a count of registered accounts or unique people across all time. Shared networks and automated traffic can affect totals.

Santa Clara, US: 12 cumulative daily visits12Los Angeles, US: 8 cumulative daily visits8New York, US: 7 cumulative daily visits7Tracy, US: 7 cumulative daily visits7Hong Kong, HK: 6 cumulative daily visits6Ashburn, US: 6 cumulative daily visits6Mountain View, US: 4 cumulative daily visits4São Paulo, BR: 3 cumulative daily visits3Frankfurt am Main, DE: 3 cumulative daily visits3Council Bluffs, US: 3 cumulative daily visits3Livermore, US: 3 cumulative daily visits3Your current approximate location (shown only to you)
84 geolocated cumulative daily visits · 11 locations meet the public threshold

Cumulative location countYour approximate current location

Your approximate current location: Columbus, US. This blue marker is shown only in your response.

Locations are approximate IP-derived city/region estimates, never GPS. Public markers require at least three cumulative daily visits. Numerical Modeling stores aggregate locations and short-lived keyed daily identifiers outside the public site, not raw IP addresses. Map boundaries are for orientation, not a statement of jurisdiction.

Visitor locations & map sources
  • Santa Clara, US — 12 cumulative daily visits
  • Los Angeles, US — 8 cumulative daily visits
  • New York, US — 7 cumulative daily visits
  • Tracy, US — 7 cumulative daily visits
  • Hong Kong, HK — 6 cumulative daily visits
  • Ashburn, US — 6 cumulative daily visits
  • Mountain View, US — 4 cumulative daily visits
  • São Paulo, BR — 3 cumulative daily visits
  • Frankfurt am Main, DE — 3 cumulative daily visits
  • Council Bluffs, US — 3 cumulative daily visits
  • Livermore, US — 3 cumulative daily visits

Basemap: World Atlas, derived from Natural Earth public-domain map data. IP geolocation uses the same local database as To the Infinity.