uncertainty quantification
/ you-kew /
Every simulation eats inputs and produces an output: feed in material properties, loads, and parameters, and out comes a predicted stress, temperature, or price. But the inputs are never known exactly — a measured stiffness has error bars, a parameter is only estimated, the weather forecast feeding a model is itself uncertain. Uncertainty quantification (UQ) is the discipline of carrying that input uncertainty through the simulation and reporting it on the output, so that instead of a bare number you produce an honest statement: 'the answer is 42, plus or minus 5, with this distribution'. A computed prediction without an uncertainty attached is, in this field's view, only half a result.
The core task is forward propagation: given a distribution over the uncertain inputs, what is the resulting distribution over the output? The simplest robust tool is Monte Carlo: draw many random samples of the inputs from their distributions, run the full simulation on each, and look at the spread of outputs — its mean estimates the expected answer, its standard deviation estimates the uncertainty. Monte Carlo is wonderfully general (it treats the simulation as a black box) but its error falls only like 1/sqrt(N), so cutting the error in half costs four times the samples — though, crucially, that rate is independent of the number of uncertain inputs, which is why Monte Carlo wins in high dimensions. When the simulation is smooth and the uncertain inputs are few, polynomial chaos and other stochastic-expansion methods can converge far faster by building a surrogate that captures how the output depends on the inputs. Alongside propagation sits SENSITIVITY ANALYSIS — which inputs actually drive the output's uncertainty — telling you where to spend effort reducing it.
UQ is what turns a number into evidence, and it connects to the broader landscape of forward and inverse problems. A FORWARD problem pushes known inputs to an output; UQ's forward propagation puts error bars on it. An INVERSE problem runs the other way — inferring uncertain parameters from observed data — and its outputs are naturally uncertain too, often described by a posterior distribution in the Bayesian framing. The honest cautions: UQ propagates the uncertainties you MODELLED, not the ones you forgot — an unmodelled bias or a wrong assumption about the input distribution will not show up in the error bars, however carefully computed. Garbage assumptions in, confident-looking garbage error bars out. And UQ quantifies uncertainty; it does not, by itself, certify that the simulation is verified or the model validated — those must be established first for the error bars to mean anything.
A beam's predicted deflection depends on its elastic modulus E, known only as 200 GPa plus or minus 10 GPa. Monte Carlo UQ draws 10000 random E values from that distribution, runs the deflection model on each, and reports the output spread: deflection 2.0 mm with standard deviation 0.1 mm. Want half that output uncertainty bar from sampling noise alone? You need roughly four times the samples — the 1/sqrt(N) tax.
Monte Carlo pushes a distribution of inputs through the model to get a distribution of outputs.
UQ propagates only the uncertainties you actually modelled. An unmodelled bias or a wrong input distribution gives confident-looking error bars that are themselves wrong — and UQ presumes the code is already verified and the model validated.