Computational engineering / Reproducible scientific computing

From physical laws to reproducible computation.

Connecting physics, mathematical models, numerical methods, and executable code through carefully verified engineering problems.

An equation, made executable.Illustrative translation / Newtonian cooling

01 Energy balance

mcpT˙=−hA(T−T∞)

Newton’s law of cooling in a lumped thermal model.

Spatially uniform temperature. Constant ambient temperature and coefficients. Convection only; no heat generation.

m
mass kg
cp
specific heat J kg−1 K−1
h
heat-transfer coefficient W m−2 K−1
A
surface area m2

02 Temperature ODE

T˙=−β(T−T∞)β=hAmcp

T — body temperature, K
T∞ — ambient temperature, K
β — cooling rate constant, s−1

Dots denote time derivatives.
The initial temperature completes the model.

03 Executable step Python

# T_inf: ambient temperature
beta = h * A / (m * cp)
dT_dt = -beta * (T - T_inf)

# Explicit Euler: use the old T
T_next = T + dt * dT_dt

Explicit Euler · monotone relaxation when 0 < β Δt ≤ 1.

Balance → ODE → DiscretizationConservation of energy with Newtonian convection; both sides of the balance have units of power (W). Assume m, cp, A > 0 and h ≥ 0. In code, cp is cp, T_inf is T∞, and dt = Δt > 0 (s). An illustrative step, not experimental validation.

The method

A traceable path from theory to implementation.

Each investigation keeps its assumptions, equations, algorithms, and evidence connected.

A running program is a beginning.
It is not evidence of accuracy.

  1. Physical Principles

    Identify the system and relevant balance laws. Make units, dimensions, and physical interactions explicit.

  2. Mathematical Modeling

    State assumptions, variables, parameters, and initial or boundary conditions. Define the model’s domain of validity.

  3. Numerical Methods

    Choose a discretization appropriate to the governing problem. Examine consistency, stability, and convergence.

  4. Implementation

    Translate the method into readable Python, with reproducible inputs and time grids consistent with the update rule.

  5. Verification & Interpretation

    Compare with analytical or independent reference solutions. Quantify error and distinguish numerical verification from physical validation.

Verification checks the calculation. Validation evaluates the model against independent physical evidence.

The investigations

Published projects

One problem. An explicit model.
Computation open to inspection.

Computational investigation
  • Heat Transfer
  • ODE
  • Python
  • Numerical Verification

How Long Does a Hot Sphere Take to Cool?

A lumped-capacitance heat-transfer investigation: deriving the cooling model, implementing explicit Euler integration, and checking analytical agreement, step-size convergence, limiting behavior, and discrete energy balance. This is an idealized teaching model, not an experimental measurement.

Areas of study

Connected by computation.

Fields of interest guiding future investigations.

About the project

An independent,
long-term investigation.

Equation to Code is an independent project devoted to translating governing equations into reproducible computational experiments, with particular attention to physical interpretation, numerical accuracy, and engineering relevance.

The detailed work lives in version-controlled notebooks and repositories: assumptions, derivations, implementations, experiments, and their limitations.

About the author