Computational Engineering & Modelling

Numerical Methods & Scientific Computing

Selection, implementation review and verification of numerical schemes for differential and algebraic problems, including stability, accuracy and performance considerations.

Last updated:

  • DivisionComputational Engineering & Modelling
  • Scope items5 defined
  • Typical inputs3 to share
  • Agreed deliverables3 outputs

Service overview

Mathematical modelling, numerical methods, multiphysics, uncertainty analysis and digital twins. This service is part of that division and is scoped individually around what you need to learn, decide or publish.

Numerical methods are the engine under every simulation. Choosing a scheme and a solver, and knowing its stability and accuracy properties, determines whether a calculation is trustworthy, efficient or merely plausible. Research codes often implement a scheme without a clear picture of its order of accuracy, its stability limits or the way errors accumulate, which becomes a problem when results are published.

This service reviews, selects and verifies numerical schemes for ordinary and partial differential equations and algebraic systems. It suits researchers who write their own codes, those comparing algorithms in a methods paper and those who need to justify the numerics in a thesis. The aim is evidence, in the form of tests, that the code solves its equations as accurately as claimed.

If you plan to publish the outcome, see manuscript & publication support, for example computational & simulation manuscripts or original research article writing.

What this service covers

Problems this service can address

  • Numerical instability or non-physical oscillations
  • Unknown order of accuracy of a home-grown code
  • Slow solution times that limit parameter studies

Methods and analytical approaches

  • Finite-difference, finite-volume and finite-element scheme review
  • Stability and convergence analysis
  • Linear and nonlinear solver selection
  • Code verification using manufactured solutions

Applications

  • Research code verification
  • Algorithm comparison studies
  • Methods papers in computational science

Scope, inputs and deliverables

Scope of work

  • Definition of the physics, geometry and quantities of interest
  • Model set-up with verification and, where data allow, validation
  • Parametric or sensitivity studies agreed at scoping
  • Results and documentation suitable for reports or publication
  • Feasibility check on: equations and discretisation used

Typical inputs required from you

  • Equations and discretisation used
  • Existing code or pseudocode where relevant
  • Benchmarks or known solutions

Deliverables

  • Verification test suite description and results
  • Convergence-order evidence
  • Recommendations for scheme or solver changes

Highly specialised requests are first checked for feasibility from this information. Please do not send confidential data through the public form; describe it first and agree secure transfer afterwards.

Project workflow and method

  1. 01Enquiry and requirement assessment
  2. 02Technical scoping and proposal
  3. 03Execution
  4. 04Quality review and delivery, with an agreed round of revisions

Read the full project workflow.

In depth: how the work is done

  1. Choosing a discretisation

    Finite-difference methods are simple on structured grids; finite-volume methods enforce conservation locally and suit transport problems; finite-element methods handle complex geometries and variational formulations; spectral and other high-order methods offer rapid convergence for smooth solutions. The choice depends on the physics, geometry, smoothness of the solution and conservation requirements. We discuss the trade-offs of accuracy, robustness and implementation effort and recommend a scheme in light of your objectives.

  2. Stability, consistency and convergence

    A scheme should be consistent with the continuous equation, stable so that errors do not grow without bound, and therefore convergent. Stability analysis, such as von Neumann analysis for linear problems, and conditions such as Courant-type limits for explicit schemes, determine permissible time steps. Stiff problems require implicit or specially designed integrators. We check whether stability limits are respected and whether artificial dissipation or limiting used for robustness degrades the accuracy that matters.

  3. Solvers for algebraic systems

    Discretisation often leads to large sparse linear or nonlinear systems whose solution dominates computing cost. Direct solvers are robust but memory-hungry; iterative Krylov methods with preconditioning scale better but depend on problem structure. Nonlinear systems use Newton-type methods, requiring good initial guesses and careful Jacobian treatment. Convergence criteria and residual monitoring are specified and reported, since loose tolerances can masquerade as physical behaviour.

  4. Code verification with manufactured solutions

    The method of manufactured solutions chooses an analytical function, substitutes it into the equations to derive a source term and then checks that the code reproduces the function with the expected order of convergence as the mesh is refined. This isolates implementation errors from model errors and is far more discriminating than eyeballing results. Additional tests include conservation checks, symmetry tests, comparison with known analytical solutions and regression tests to protect against later changes.

  5. Reporting algorithms in methods papers

    Methods papers carry a higher burden of proof. They should state the algorithm clearly, give pseudo-code or equations in sufficient detail to implement, specify test problems and metrics, and compare against appropriate existing methods under fair conditions. We help structure the numerical experiments, choose test cases that reveal strengths and weaknesses, and present convergence and cost results honestly, including cases where the proposed method does not perform well.

Quality control and good practice

Quality-control approach

  • Verification: mesh, time-step and solver sensitivity are examined before results are reported.
  • Validation: results are compared with measurements or reference solutions where these exist, and the absence of such data is stated.
  • Boundary conditions and material or fluid properties are recorded with their sources.
  • Sensitivity to the least certain inputs is shown or explained.

Common pitfalls this service helps avoid

  • Reporting the nominal order of accuracy without testing. Implementation errors often reduce the observed order; only measurement shows what the code achieves.
  • Excessive numerical diffusion. Stabilising techniques can smear features that carry the physics of interest.
  • Loose solver tolerances. Unconverged iterative solutions can introduce errors larger than discretisation error.
  • Comparing algorithms unfairly. Methods differ in cost per step, so comparisons should use equal accuracy or equal computing time.

What good reporting looks like

A numerical-methods report describes the equations, the scheme and its theoretical properties, the implementation, test problems with analytical or manufactured solutions, observed convergence rates, stability limits, solver settings and tolerances, computational cost and limitations. Where algorithms are compared, the criteria and test cases are stated and the comparison is designed to be fair. The report supports the methods section of a manuscript and gives readers confidence in the numerics.

An illustrative scenario

Hypothetical example written for explanation. It is not a client case, a result or a testimonial.

Imagine a doctoral researcher who has written a finite-volume code for a convection-diffusion problem. The code produces smooth profiles, but the observed convergence order from a manufactured-solution test is lower than the scheme's nominal order. Investigating, the researcher finds an error in the treatment of a boundary flux and an unintended first-order limiter in smooth regions. After correction, the observed order matches the theory, and the thesis can cite the verification test as evidence.

Questions to settle before you begin

  • Which properties of the solution, such as shocks or boundary layers, must the scheme preserve?
  • What order of accuracy do I expect, and how will I measure it?
  • What limits the time step?
  • How will I know solvers have converged?
  • How will I compare my method fairly with existing ones?

Key terms

Consistency
The discrete equations approach the continuous ones as the mesh is refined.
Stability
Errors do not grow unboundedly during computation.
Courant number
A dimensionless ratio relating wave speed, time step and cell size.
Stiffness
A property of equations containing widely separated time scales.
Manufactured solution
A prescribed function used to test a code by constructing a matching source term.

Working together on Computational Engineering & Modelling projects

Computational engagements begin by agreeing the question, the quantities of interest and the accuracy that matters, so that modelling effort goes where it changes the answer. We ask for geometry or system descriptions, operating conditions, material or property data and any measurements for validation. Early on, a small pilot case often checks feasibility and cost. Models are verified and, where data allow, validated before parametric results are produced. Software, versions and settings are recorded so results can be reproduced. Where inputs are missing, assumptions are listed and agreed rather than quietly filled in.

Frequently asked questions

Can you write or fix my research code?

Scope is agreed in advance. Reviewing and verifying numerical approaches is a typical task; ownership and licensing of your code remain with you.

Can you review code I wrote myself?

Yes. The review focuses on numerical properties and verification evidence, and on any issues with accuracy, stability and efficiency visible from the code and tests.

Do you write production code?

The service concentrates on method selection, verification and review. Implementation support can be discussed if it fits the scope.

Which languages are supported?

Review is concerned with algorithms and can be applied to code in widely used scientific languages. Specifics are agreed at scoping.

What is a manufactured solution?

An analytical function chosen to test a code by deriving the source term that makes it an exact solution, allowing precise error measurement.

How is this service scoped and priced?

After your enquiry we review what you need, ask for missing information and send a written scope and quotation. Price depends on complexity, deliverables and timeline. See the project workflow.

Discuss your project: Numerical Methods & Scientific Computing

Share your objectives, data situation and target outcome. We will assess the scope and suggest sensible next steps.