Skip to main content
The new Teacher Workspace is here. Your first 3 assignments are free. Try it →

Finite difference method

The finite difference method is a numerical way to solve heat and mass transfer equations by replacing derivatives with differences on a grid of points. In this course, it is used for conduction, diffusion, and time-dependent transport problems.

Last updated July 2026

What is the finite difference method?

The finite difference method is a grid-based numerical technique for solving heat and mass transfer equations when an exact formula is hard or impossible to get. Instead of treating temperature or concentration as a smooth curve everywhere, you calculate values at specific nodes and use those values to approximate derivatives.

That makes it a natural fit for conduction and diffusion problems. Heat transfer equations and mass transfer equations often involve spatial derivatives, and when the geometry, boundary conditions, or time dependence get complicated, you can rewrite the differential equation as a set of algebraic equations at each node. Then you solve the system step by step.

The core move is discretization. You split space into small intervals, and sometimes time into small steps too. A spatial derivative like dT/dx becomes a difference such as (T at one node minus T at the neighboring node) divided by the spacing between them. A second derivative, which shows up in diffusion, is built from nearby node values as well.

For steady-state problems, the temperatures or concentrations do not change with time, so the finite difference method gives you a grid of equations that describe the final distribution across the domain. For unsteady problems, you also march forward in time, updating each node from one time level to the next. That is why the method shows up in topics like multidimensional steady-state diffusion and unsteady-state conduction.

A simple way to picture it is a furnace wall or a cooling plate. You cannot solve every real shape with a clean hand calculation, but you can place nodes across the material, assign boundary conditions at the edges, and compute the interior values. The smaller the grid spacing, the closer the approximation is to the real physical profile, although the math and computation get heavier.

Why the finite difference method matters in Heat and Mass Transfer

Finite difference method is the bridge between the transport equations you write down and the temperature or concentration profile you actually want. In Heat and Mass Transfer, that matters because many real problems stop being friendly the moment the geometry becomes two-dimensional, the boundary conditions become mixed, or the process changes over time.

It also turns the governing equation into something you can solve with a calculator, spreadsheet, or computer. That is why it shows up in numerical methods for conduction problems and unsteady diffusion. Instead of searching for a rare closed-form solution, you build a grid, apply the boundary conditions, and compute the field at each node.

The method also teaches you how the physics behaves locally. A steep temperature gradient, for example, creates a larger difference between neighboring nodes, which shows up directly in the discrete equations. That makes the method useful for checking whether your results make physical sense, not just whether the algebra works.

It becomes even more useful in inverse heat and mass transfer problems, where you work backward from measured data to estimate an unknown heat flux, surface condition, or source term. In those cases, finite difference is often part of the setup that links unknown inputs to observed outputs.

Keep studying Heat and Mass Transfer Unit 7

Official unit cheatsheet

open one-pager

How the finite difference method connects across the course

Discretization

Finite difference method depends on discretization, which means replacing a continuous region with a finite set of points. If your grid is too coarse, the approximation can miss steep temperature or concentration changes. If it is finer, the result usually improves, but you also solve more equations.

Boundary Conditions

The method only works once the edges of the problem are defined. Boundary conditions tell you what happens at the surfaces, such as a fixed temperature, a specified heat flux, or a convective boundary. Those conditions become the starting and ending equations in the finite difference grid.

Grid System

A grid system is the map of nodes where the solution is computed. In one dimension it may be a simple line of points, but in multidimensional conduction or diffusion it becomes a mesh. The layout of that grid affects both accuracy and how easy the equations are to solve.

Neumann Boundary Condition

Neumann conditions specify a derivative at the boundary, usually heat flux in conduction or mass flux in diffusion. Finite difference handles this by converting the derivative into a difference formula at the edge node. That makes Neumann boundaries very common in heat transfer models with insulated or flux-driven surfaces.

Is the finite difference method on the Heat and Mass Transfer exam?

A problem set question will usually ask you to set up the finite difference equations, not just name the method. You may be given a wall, plate, or diffusion slab with boundary conditions and asked to write the node equations, solve for interior temperatures or concentrations, and interpret the pattern.

For unsteady problems, you might need to update the solution from one time step to the next and compare an explicit or implicit scheme. The real skill is matching the equation to the physical setup, then checking whether your grid spacing and boundary conditions make sense.

In inverse problems, expect the method to show up as the numerical backbone for estimating an unknown surface condition from measured data. A good answer usually shows the discretized equations, the boundary setup, and a short interpretation of the result.

The finite difference method vs Finite Element Method

Finite difference method and finite element method both solve transport equations numerically, but they set up the problem differently. Finite difference builds equations from point-to-point differences on a grid, while finite element divides the domain into elements and uses interpolation functions. In Heat and Mass Transfer, finite difference is often the simpler choice for regular geometries.

Key things to remember about the finite difference method

  • Finite difference method replaces derivatives with difference formulas at grid points.

  • It is especially useful for conduction and diffusion problems that are hard to solve exactly.

  • Steady-state problems give you a spatial temperature or concentration field, while unsteady problems also move through time.

  • Boundary conditions are part of the setup, not an afterthought, because they shape the node equations.

  • A finer grid usually improves accuracy, but it also increases the amount of computation.

Frequently asked questions about the finite difference method

What is finite difference method in Heat and Mass Transfer?

It is a numerical method that approximates heat or mass transfer equations by replacing derivatives with differences between values at nearby grid points. In this course, you use it to solve conduction and diffusion problems when an exact solution is too hard to find.

How does finite difference method work for conduction problems?

You divide the material into nodes, write a difference equation at each interior node, and apply the boundary conditions at the edges. The result is a system of algebraic equations that gives you the temperature at each point in the grid.

What is the difference between finite difference and finite element method?

Finite difference uses point-based difference formulas on a grid, while finite element uses elements and interpolation functions. Both are numerical methods, but finite difference is usually easier to set up for simple, regular heat transfer geometries.

Why does grid size matter in finite difference method?

Grid spacing controls how closely the numerical model matches the real gradient in temperature or concentration. Smaller spacing usually means better accuracy, especially where the profile changes quickly, but it also means more equations to solve.

Finite Difference Method | Heat and Mass Transfer | Fiveable