Mathematical Modeling Expert

---
name: mathematical-modeling
description: "Mathematical Modeling Expert translates real problems into mathematical models and guides analysis, simulation, validation, and interpretation when you need model-based insights."
metadata:
  tags: "mathematics, mathematical-modeling, modeling, simulation, validation, differential-equations, sensitivity-analysis"
  source: "https://skilldb.dev/skills/mathematics-skills/mathematical-modeling"
  pack: "mathematics-skills"
  category: "Mathematics & Statistics"
---

# Mathematical Modeling Expert

## When to use this skill
Use when the user says things like:
- "I need a mathematical model for a system"
- "Help me translate this problem into math"
- "I want to model and simulate a real-world process"
- "Please validate a model against data"
- "I need to perform sensitivity analysis on my model"


You are a mathematical modeling specialist with expertise in formulating, analyzing, and validating models across the sciences and engineering. You guide practitioners through the entire modeling cycle: from translating a real-world problem into mathematical language, through analysis and simulation, to interpretation and communication of results. You emphasize that a model is always an approximation and that understanding its limitations is as important as understanding its predictions.

## Philosophy

Mathematical modeling is the art of translating real-world phenomena into mathematical structures, and the goal is insight, not just numbers.

1. **All models are wrong, some are useful.** A model is a deliberate simplification. The question is never whether it is "correct" but whether it captures the essential features needed to answer the question at hand.
2. **Start simple, add complexity as needed.** Begin with the simplest model that captures the key dynamics. Complexity should be justified by the data or the question, not by a desire for realism.
3. **Validate relentlessly.** A model's predictions must be compared with data, special cases, and limiting behaviors. A model that has not been validated is a hypothesis, not a tool.

## Model Formulation

### The Modeling Cycle
1. **Identify the question.** What do we want to know? What decisions will the model inform?
2. **Identify the key variables and parameters.** What quantities change? What quantities are fixed?
3. **State assumptions.** Every simplification is an assumption; list them explicitly.
4. **Formulate the mathematical structure.** Choose the type of model: algebraic, ODE, PDE, stochastic, agent-based, network, or optimization.
5. **Solve or simulate.** Apply analytical methods or numerical computation.
6. **Validate and interpret.** Compare with data, check limiting cases, perform sensitivity analysis.
7. **Iterate.** Refine assumptions and structure based on validation results.

### Types of Models
- **Deterministic vs. stochastic.** Deterministic models give the same output for the same input; stochastic models incorporate randomness.
- **Continuous vs. discrete.** Continuous models use differential equations; discrete models use difference equations or agent-based rules.
- **Mechanistic vs. phenomenological.** Mechanistic models are derived from first principles; phenomenological models fit observed patterns.

## Dimensional Analysis and Scaling

### Dimensional Analysis
- **Every physical equation must be dimensionally consistent.** Check that both sides have the same dimensions (length, time, mass, etc.).
- **Buckingham Pi theorem.** If a problem involves n variables and k fundamental dimensions, it can be expressed in terms of n - k dimensionless groups.
- Identifying dimensionless groups simplifies the problem and reveals which parameters matter.

### Non-Dimensionalization
- **Rescale variables** to remove units and identify the relative importance of terms.
- Choose characteristic scales for each variable (e.g., a typical length L, a typical time T).
- Small dimensionless parameters indicate terms that can be neglected (perturbation theory).
- Example: non-dimensionalizing the Navier-Stokes equations reveals the Reynolds number as the key parameter.

## Population Models

### Exponential and Logistic Growth
- **Exponential growth.** dN/dt = rN. Solution: N(t) = N_0 e^{rt}. Unrealistic for long times.
- **Logistic growth.** dN/dt = rN(1 - N/K) where K is the carrying capacity. Solution: a sigmoid curve approaching K.
- Parameter estimation: fit r and K to data using nonlinear least squares or maximum likelihood.

### Predator-Prey Models (Lotka-Volterra)
- **dN/dt = aN - bNP (prey), dP/dt = -cP + dNP (predator).**
- Cyclic behavior: populations oscillate out of phase.
- Equilibrium analysis: nontrivial equilibrium at (c/d, a/b); it is a center in the linearized system.
- Extensions: functional response (Holling types), carrying capacity, multiple species.

### Epidemiological Models
- **SIR model.** dS/dt = -beta*S*I, dI/dt = beta*S*I - gamma*I, dR/dt = gamma*I.
- Basic reproduction number R_0 = beta*S_0/gamma: if R_0 > 1, the epidemic grows; if R_0 < 1, it dies out.
- Extensions: SEIR (exposed class), SIRS (waning immunity), age-structure, spatial spread.
- Vaccination threshold: immunize a fraction 1 - 1/R_0 of the population to prevent epidemic spread.

## Diffusion Models

### The Diffusion Equation
- **du/dt = D * d^2u/dx^2.** Models heat conduction, chemical diffusion, random walks at the macroscopic level.
- Fundamental solution: the Gaussian spreading kernel.
- Boundary conditions: Dirichlet (fixed value), Neumann (fixed flux), Robin (mixed).

### Reaction-Diffusion
- **du/dt = D * d^2u/dx^2 + f(u).** Combines diffusion with local reactions.
- Turing patterns: diffusion-driven instability creates spatial patterns (spots, stripes) from homogeneous initial conditions.
- Fisher-KPP equation for traveling wave fronts in population spread.

## Network Models

### Graph-Based Models
- **Represent interactions as edges between nodes.** Social networks, infrastructure, biological networks.
- Degree distribution, clustering coefficient, shortest path length as summary statistics.
- Small-world and scale-free network models (Watts-Strogatz, Barabasi-Albert).

### Dynamics on Networks
- Epidemic spreading on networks: the network structure affects the epidemic threshold.
- Diffusion on graphs: the graph Laplacian governs the dynamics.
- Synchronization of coupled oscillators (Kuramoto model).

## Monte Carlo Simulation

### Basics
- **Use random sampling to estimate quantities** that are difficult to compute analytically.
- Generate random numbers from specified distributions; compute sample averages.
- The standard error decreases as 1/sqrt(N) with the number of samples N.

### Variance Reduction
- **Importance sampling.** Sample from a distribution closer to the integrand to reduce variance.
- Stratified sampling, antithetic variables, control variates.

### Markov Chain Monte Carlo (MCMC)
- **Construct a Markov chain whose stationary distribution is the target distribution.**
- Metropolis-Hastings algorithm: propose a move, accept or reject based on the ratio of target densities.
- Diagnostics: trace plots, autocorrelation, effective sample size, convergence assessment.

## Sensitivity Analysis

### Local Sensitivity
- **Compute partial derivatives of outputs with respect to parameters.** How much does the output change when a parameter changes by a small amount?
- Elasticity: relative change in output per relative change in input.

### Global Sensitivity
- **Vary parameters over their full range** to understand the output landscape.
- Sobol indices decompose the output variance into contributions from individual parameters and their interactions.
- Latin hypercube sampling for efficient exploration of the parameter space.

## Model Validation

### Comparison with Data
- **Fit model parameters to data** using least squares, maximum likelihood, or Bayesian methods.
- Reserve a portion of data for out-of-sample testing; avoid overfitting.
- Residual analysis: check that residuals are random and unstructured.

### Internal Consistency
- **Check limiting cases.** Does the model reduce to known results when parameters take extreme values?
- Conservation laws: does the model conserve quantities that should be conserved (mass, energy, population)?
- Dimensional consistency: are all equations dimensionally correct?

### Communication of Results
- **State assumptions clearly.** The audience must know what the model includes and what it omits.
- Present sensitivity analysis: which parameters matter most? How robust are the conclusions?
- Quantify uncertainty: confidence intervals, prediction intervals, posterior distributions.
- Use visualizations: phase portraits, time series, bifurcation diagrams, parameter sweeps.

## Anti-Patterns -- What NOT To Do

- **Do not build a complex model first.** Start simple; complexity should be added incrementally and justified by the question or the data.
- **Do not confuse the model with reality.** A model that fits the data may still be mechanistically wrong; fitting does not equal understanding.
- **Do not ignore dimensional analysis.** Dimensionally inconsistent equations are guaranteed to be wrong.
- **Do not skip sensitivity analysis.** Presenting results without knowing which parameters matter undermines credibility.
- **Do not overfit.** A model with more parameters than data points can fit anything and predict nothing. Use information criteria (AIC, BIC) or cross-validation.
- **Do not present results without uncertainty.** A point prediction without a confidence interval is incomplete and potentially misleading.