---
title: Multifidelity Trajectory Simulation Estimators
url: https://www.emergentmind.com/topics/multifidelity-trajectory-simulation-estimators
type: topic
---

# Multifidelity Trajectory Simulation Estimators

Multifidelity trajectory simulation estimators are statistical constructions that combine one trusted high-fidelity simulator with cheaper lower-fidelity models to estimate trajectory-dependent quantities of interest under explicit computational budgets. In the literature, the target may be a scalar expectation such as $\mu=\mathbb{E}[f(Y)]$, the full probabilistic response $p(y_H)$, trajectory statistics such as means, standard deviations, and peak magnitudes, stacked multi-output means, variances, and Sobol main effect variances, or path-integral control and safety probabilities [2103.15342; 1808.03379; 2001.02892; 2310.00125; 2510.07756; 2509.16007]. The common premise is that low-fidelity trajectory information, if exploited through regression, Bayesian conditioning, control variates, multilevel telescoping, or hyperparameter tuning, can reduce estimator variance or mean-squared error for a fixed cost.

## 1. Problem classes and mathematical setting

A recurring setup uses one high-fidelity simulator $Y$ with cost $c_0$ and $n$ lower-fidelity surrogates $X^{(i)}$, $i=1,\dots,n$, with costs $c_i$, all sampled subject to a total budget
\[
\sum_{i=0}^n c_i\,n_i \le B.
\]
The objective is to estimate $\mu:=\mathbb{E}[f(Y)]$ under that cost constraint [2103.15342]. In the trajectory-simulation setting, the simulator output may itself be a time series, so the estimator can target either the entire trajectory or a scalar functional extracted from it.

For parameterized ordinary differential equations,
\[
\frac{du}{dt}(t,k)=f\bigl(t,u(t,k),k\bigr),\qquad u(0,k)=u_0(k),
\]
a natural hierarchy is induced by time-integration levels with timesteps $h_j=h/r^{j-1}$ and $N_j=Nr^{j-1}$, giving discrete trajectories $u_j(\cdot,k)$ at progressively finer resolutions [1808.03379]. In this setting, the fidelity notion is numerical rather than statistical: the lower-fidelity models are trajectories computed with larger timesteps.

In uncertainty propagation, the high-fidelity output $y_H(x)$ and low-fidelity output $y_L(x)$ are both driven by a random input $x\in\Omega\subset\mathbb{R}^d$. The target density satisfies
\[
p(y_H)=\int p(y_H\mid y_L)\,p(y_L)\,dy_L,
\]
so the multifidelity task becomes learning the conditional relationship between high- and low-fidelity responses from a small paired dataset while cheaply sampling the low-fidelity model many times [2001.02892].

For multi-output estimators, each model $\ell$ may map $z\in\mathbb{R}^I$ to $f_\ell(z)\in\mathbb{R}^D$, and the target estimator can stack means, variances, and Sobol main effect variances into a single vector-valued quantity [2310.00125]. For stochastic control and risk quantification, the quantity of interest may be a level-dependent approximation $P_l$ to either a path-integral control law or a safety probability, with trajectory samples generated across a hierarchy of time or state discretizations [2510.07756].

A plausible implication is that “multifidelity” in trajectory simulation is not tied to a single construction. The literature treats at least three distinct sources of fidelity variation: timestep or discretization level, model-form simplification, and surrogate hyperparameters such as $\Delta t$ in a coarse simulator [1808.03379; 2510.07756; 2509.16007].

## 2. Linear-surrogate and bandit estimators under a budget

One line of work formulates multifidelity approximation as a linear regression problem over subsets of available surrogates. For any subset $S\subseteq\{1,\dots,n\}$,
\[
f(Y)=X_S^\top\beta_S+\epsilon_S,
\]
where $X_S=(1,X^{(i)}:i\in S)^\top$, $\beta_S\in\mathbb{R}^{|S|+1}$, and $\epsilon_S$ is zero-mean noise with variance $\sigma_S^2$, independent of $X_S$ [2103.15342]. The estimator is built in two phases: an exploration stage with $m$ joint samples $(X_S,f(Y))$ and an exploitation stage in which a chosen surrogate model is sampled alone.

With exploration estimate $\hat\beta_S$, the linear-regression Monte Carlo estimator is
\[
\widehat\mu_S=\frac{1}{N_S}\sum_{\ell=1}^{N_S} X_{S,\ell}^\top\hat\beta_S.
\]
Its conditional mean-squared error is
\[
R_S\mid_{\hat\beta_S}
=\mathbb{E}\bigl[(\widehat\mu_S-\mu)^2\mid \hat\beta_S\bigr]
=(x_S^\top(\hat\beta_S-\beta_S))^2+\frac1{N_S}\hat\beta_S^\top\Sigma_S\hat\beta_S,
\]
and the average conditional MSE has an asymptotic form
\[
\overline R_S(m)\approx
\frac{c_S\,\beta_S^\top\Sigma_S\beta_S}{B-c_{\exp}m}
+\frac{\sigma_S^2\,\mathrm{Tr}(x_Sx_S^\top\Lambda_S^{-1})}{m}.
\]
This decomposition separates exploitation error, which scales with the remaining budget, from exploration error, which decays with $m$ [2103.15342].

The same work casts each subset $S$ as an arm in a modified stochastic bandit. Policies that explore each arm uniformly and then exploit the empirically best arm are termed uniform explore-then-exploit policies. The proposed Adaptive Explore-Then-Commit (AETC) algorithm computes empirical upper bounds
\[
\hat H_S(t)=\min_{m'\ge t}\left[\frac{k_{1,t}(S)}{B-c_{\exp}m'}+\frac{k_{2,t}(S)+\alpha_t}{m'}\right],
\]
with
\[
k_{1,t}(S)=c_S\,\hat\beta_S^\top\hat\Sigma_S\hat\beta_S,\qquad
k_{2,t}(S)=\hat\sigma_S^2\,\mathrm{Tr}(\hat x_S\hat x_S^\top\hat\Lambda_S^{-1}),
\]
and commits when further exploration is no longer preferable [2103.15342].

The theoretical guarantees are explicit. Theorem 3.1 states that under mild sub-Gaussian design assumptions, for fixed $S$ and large $m$, with probability at least $1-O(m^{-2})$,
\[
R_S(m)\le \beta_S^\top\Sigma_S\,\beta_S/N_S + (|S|+1)\sigma_S^2(\ln m)/m,
\]
hence
\[
R_S(m)=O\!\left(\frac1{N_S}\right)+O\!\left(\frac{\ln m}{m}\right).
\]
Theorem 4.1 states that if $m(B),S(B)$ denote AETC’s exploration length and exploitation model, then almost surely as $B\to\infty$,
\[
m(B)/m_{S^*}\to1,\qquad S(B)\to S^*,
\]
where $m_{S^*},S^*$ are the oracle choices minimizing the leading-order MSE. The paper therefore describes AETC as “trajectory-wise optimal” [2103.15342].

The vector-valued extension replaces scalar noise by $\eta_S\sim(0,\Gamma_S)$ in
\[
f(Y)=\beta_S X_S+\eta_S
\]
and measures error through a quadratic $Q$-risk,
\[
R_Q(\hat Y)=\mathbb{E}\bigl[\|Q(\hat Y-\mathbb{E}[f(Y)])\|_2^2\bigr].
\]
The resulting bounds depend only on scalar traces such as $\mathrm{Tr}(Q\Gamma_SQ^\top)$ and $\mathrm{Tr}(\Sigma_S\beta_S^\top Q^\top Q\beta_S)$, so the approach avoids estimating high-dimensional parameters directly [2103.15342].

A common misconception is that multifidelity estimation requires an ordered fidelity hierarchy or prior covariance knowledge. This framework states the opposite: it does not assume that fidelity $i$ is strictly cheaper or less accurate than $j$, does not require any nested hierarchy, and does not require prior knowledge of $\mathrm{Cov}(X^{(i)},X^{(j)})$ or $\mathrm{Cov}(X_S,f(Y))$; these quantities are learned during exploration [2103.15342].

## 3. Time-integration hierarchies and convergence acceleration

For parametric ODE trajectories, one three-level construction begins by defining low-, medium-, and high-fidelity models through timesteps
\[
h_1=h,\qquad h_2=h/r,\qquad h_3=h/r^2,
\]
with discrete solutions computed by any global-order-$p$ one-step or multistep scheme [1808.03379]. The exact trajectory is viewed in the Hilbert space $H=L^2([0,T];\mathbb{R}^M)$, while each discrete trajectory lies in a level-dependent space $H_j=\ell^2(\{0,\dots,N_j\};\mathbb{R}^M)$ equipped with a weighted inner product based on composite Newton–Cotes weights.

The offline stage constructs a reduced basis of size $n\ll \dim(H_j)$ by a greedy sweep on inexpensive level-1 solutions. Using a large training set $\{k^{(1)},\dots,k^{(Q)}\}\subset\mathcal K$, one computes $u_1(\cdot,k^{(\ell)})$ and selects “important” parameters $\mathcal K_n=\{k_1,\dots,k_n\}$ via pivoted Cholesky, QR, or the greedy rule
\[
k_{q+1}=\arg\max_{k\in\{k^{(1)},\dots,k^{(Q)}\}}
\|u_1(\cdot,k)-\mathcal P_{q,1}u_1(\cdot,k)\|_1.
\]
For each selected $k_q$, the expensive trajectories $u_2(\cdot,k_q)$ and $u_3(\cdot,k_q)$ are then computed [1808.03379].

For a new parameter $k$, coefficients $v(k)\in\mathbb{R}^n$ are obtained from the level-1 least-squares problem
\[
v(k)=\arg\min_{w\in\mathbb{R}^n}
\left\|u_1(\cdot,k)-\sum_{q=1}^n w_q\,u_1(\cdot,k_q)\right\|_1,
\]
or equivalently by the normal equations
\[
G_1\,v(k)=f_1(k).
\]
The multifidelity surrogates are
\[
\widehat u_j(\cdot,k)=\sum_{q=1}^n v_q(k)\,u_j(\cdot,k_q),\qquad j=1,2,3.
\]
Under standard $O(h_j^p)$ global-truncation and quadrature-accuracy assumptions, Theorem 3.1 gives
\[
\|\mathcal P_n u(\cdot,k)-\widehat u_j(\cdot,k)\|_\infty
\lesssim C_1 h^p + C_2(h_j)^p,
\]
so the $j$-dependent term decays like $(h_j)^p$ [1808.03379].

The same paper then adds a three-step acceleration procedure. First, the discrete surrogates $\widehat u_j$ are built at all three levels. Second, each level is lifted to continuous time by a global $p$th-order B-spline interpolant $\widehat w_j(t,k)$ defined on a normalized knot vector and fitted by interpolation at the discrete time nodes. Third, Richardson extrapolation is applied pointwise in level index $j$. For fixed $t,k,m$ one sets $x_j=\widehat w_j^{(m)}(t,k)$, estimates the observed order
\[
p^*=\log_r\!\left(\frac{x_1-x_2}{x_2-x_3}\right),
\]
and computes
\[
x^*=\frac{r^{p^*}x_3-x_2}{r^{p^*}-1}=c\,x_3+(1-c)\,x_2.
\]
The resulting accelerated continuous surrogate $\widehat w^*(t,k)$ is reported to be practically $O(h^{p+1})$-accurate in many cases [1808.03379].

The numerical examples establish the trajectory-specific character of the method. For the damped harmonic oscillator with $k\in[5,25]$, $T=5$, $h_1=0.1$, $h_2=0.05$, $h_3=0.025$, and methods RK2, RK3, RK4 and AB2, AB3, AB4, the estimated orders $p^*$ recover the nominal method orders, the spline-accelerated errors $\|\widehat w^*-u\|_{L^\infty_t}$ attain their minima at the predicted Richardson weights, and log-error versus $h$ shows an approximate $p+1$ slope [1808.03379]. When $k\sim U[5,25]$ and 1,000 Monte Carlo samples are used to estimate statistical moments, direct high-fidelity moments have cost $1000\times X$, whereas trifidelity moments via $\widehat w^*$ have cost approximately $295\,X$ and achieve an order-of-magnitude lower $\ell^2_k$ error in mean and standard deviation. In the predator–prey example, the accelerated surrogate again shows roughly $O(h^{p+1})$ decay in $h$ for AB2, AB3, and RK4 [1808.03379].

## 4. Bayesian multifidelity Monte Carlo for trajectory uncertainty propagation

A different construction begins from the identity
\[
p(y_H)=\int p(y_H\mid y_L)\,p(y_L)\,dy_L
\]
and learns the conditional density from a small paired dataset
\[
D=\{(x_i,y_{L,i},y_{H,i})\}_{i=1}^n,
\qquad n\ll \#\text{LF-samples}.
\]
In the Bayesian multi-fidelity Monte-Carlo framework, the conditional relation is modeled as
\[
y_H=f(y_L)+\epsilon,\qquad \epsilon\sim N(0,\sigma_n^2),
\]
with a Gaussian-process prior
\[
f(\cdot)\sim GP(m(\cdot),k(\cdot,\cdot)).
\]
At a new low-fidelity value $y_L^*$, the posterior predictive distribution is
\[
p(y_H\mid y_L^*,D)=N\bigl(y_H;m_*(y_L^*),v_*(y_L^*)+\sigma_n^2\bigr)
\]
[2001.02892].

For any functional $\phi(y_H)$,
\[
I=\mathbb{E}[\phi(y_H)]
=\iint \phi(y_H)\,p(y_H\mid y_L)\,dy_H\,p(y_L)\,dy_L.
\]
Replacing the unknown conditional by the GP posterior estimate yields the BMFMC estimator
\[
\hat I_{\mathrm{BMFMC}}
=\frac1{N_L}\sum_{j=1}^{N_L}
\mathbb{E}_{\hat p(\cdot\mid y_L^{(j)})}[\phi(y_H)].
\]
For the mean, $\phi(y_H)=y_H$, this reduces to
\[
\hat I=\frac1{N_L}\sum_{j=1}^{N_L} m_*(y_L^{(j)}).
\]
Its bias is
\[
\mathrm{Bias}[\hat I]
=\mathbb{E}_x[m_*(y_L(x))]-\mathbb{E}_x[y_H(x)],
\]
and the variance is approximated by
\[
\mathrm{Var}[\hat I]
\approx \frac1{N_L}\mathrm{Var}_{y_L}[m_*(y_L)]
+\frac1{N_L^2}\sum_{j=1}^{N_L}\mathrm{Var}_{f\mid D}[m_*(y_L^{(j)})],
\]
with mean-squared error $\mathrm{MSE}=\mathrm{Bias}^2+\mathrm{Var}[\hat I]$ [2001.02892].

The practical algorithm separates cheap and expensive work. First, one draws $N_L\gg n$ low-fidelity samples. Second, one adds one or two low-dimensional features $\gamma_i(x)$ at no extra cost and forms
\[
z(x)=[y_L(x),\gamma_1(x),\dots,\gamma_p(x)].
\]
Third, one selects $n$ training indices by a space-filling design in $z$-space, evaluates the high-fidelity model there, and fits the GP $f:z\mapsto y_H$. Fourth, one predicts posterior means and variances at all low-fidelity samples and forms the Monte Carlo average [2001.02892].

The error analysis distinguishes model error, due to a misspecified conditional family, from epistemic uncertainty, due to finite $n$. A key inequality is
\[
\mathbb{E}_{y_L,\gamma}[\mathrm{Var}(y_H\mid y_L,\gamma)]
\le
\mathbb{E}_{y_L}[\mathrm{Var}(y_H\mid y_L)],
\]
showing that informative additional features reduce the average conditional variance. The computational cost of GP fitting is $O(n^3)$ and prediction at $N_L$ points costs $O(N_L n^2)$, with total work dominated by cheap low-fidelity evaluations plus $n$ high-fidelity runs [2001.02892].

For trajectory-type outputs, the framework allows at least three representations: applying BMFMC to scalar trajectory statistics such as maxima or time averages, concatenating observations at a small set of time-points into a multi-output GP, or embedding time through delays via feature vectors of the form $z(x)=[y_L(t),y_L(t+\Delta t),\dots]$ [2001.02892]. The paper also states that credible intervals from the GP posterior predictive variance can be used to certify reliability, and that overly wide intervals motivate additional high-fidelity training points or richer features.

The numerical demonstrations are directly trajectory-oriented. In stochastic flow past a cylinder, the high-fidelity model is high-order DG with $k=6$, the low-fidelity model is $k=3$ and is approximately $8\times$ cheaper, the random inputs include cylinder radius, position, viscosity, and an uncertain-inflow Gaussian field of dimension approximately $200$ reduced to $10$ via Karhunen–Loève, and $N_L=10\,000$ low-fidelity runs with $n=150$ high-fidelity runs are sufficient for BMFMC to recover the full PDF of the maximum lift coefficient to within negligible Kullback–Leibler divergence of order $10^{-3}$, with overall speed-up approximately $7\times$ [2001.02892]. In the bending-wall fluid–structure interaction problem, the quantity of interest is the tip deflection trajectory $Q(t)$ or its peak magnitude, $N_L=7000$, $n=50$, and the cheapest low-fidelity model is $28\times$ cheaper; the reported speed-up is greater than $23\times$ while credible intervals reliably bracket the reference [2001.02892].

## 5. Control variates, multi-output covariance structure, and multilevel–multifidelity fusion

Approximate control variates form another major class of multifidelity trajectory estimators. In a general linear estimator,
\[
\hat Q_{MF}(\beta,\alpha;Z)
=
Q_0(z^0)+\sum_{i=1}^M \alpha_i\,[Q_i(\beta_i,z_i^1)-Q_i(\beta_i,z_i^2)],
\]
$Q_0$ is the high-fidelity quantity of interest, $Q_i(\beta_i,\cdot)$ are parametrized low-fidelity models, and the total cost is
\[
C(\beta,N,N_1,\dots,N_M)=Nw_0+\sum_{i=1}^M N_i\,w_i(\beta_i)
\]
[2509.16007]. For fixed sample counts, the optimal weights are
\[
\alpha^*=-C_{LL}^{-1}C_{qq},
\]
where $C_{qq}=\mathrm{Cov}[Q_0(z^0),\Delta Q_{LF}]$ and $C_{LL}=\mathrm{Cov}[\Delta Q_{LF},\Delta Q_{LF}]$, giving
\[
\mathrm{Var}[\hat Q_{MF}]
=
\mathrm{Var}[Q_0]/N - C_{qq}^\top C_{LL}^{-1}C_{qq}.
\]
An equivalent explicit form expresses the variance through $V_0$, correlation vectors $\rho(\beta)$, correlation matrices $P(\beta)$, and sample-allocation matrices $f$ and $F$ [2509.16007].

For multi-output and multi-statistic estimation, the covariance structure becomes the central object. If each model output is $f_\ell(z)\in\mathbb{R}^D$, then the literature derives closed-form covariance expressions for mean–mean, variance–variance, mean–variance, main-effect Sobol, and variance–main-effect estimators. Denoting the high-fidelity stacked estimator by $Q_0$ and the stacked low-fidelity differences by
\[
\Delta=[Q_1(Z_1^*)-Q_1(Z_1);\dots;Q_K(Z_K^*)-Q_K(Z_K)],
\]
the multi-output ACV estimator is
\[
\widetilde Q(\alpha)=Q_0+\alpha \Delta,
\]
with variance
\[
\mathrm{Var}[\widetilde Q]
=
\mathrm{Var}[Q_0]+\alpha\,\mathrm{Var}[\Delta]\,\alpha^\top
+\mathrm{Cov}[Q_0,\Delta]\alpha^\top
+\alpha\,\mathrm{Cov}[\Delta,Q_0].
\]
Minimizing the determinant or trace yields
\[
\alpha^*=-\mathrm{Cov}[Q_0,\Delta]\mathrm{Var}[\Delta]^{-1},
\]
and therefore
\[
\mathrm{Var}[\widetilde Q]
=
\mathrm{Var}[Q_0]
-
\mathrm{Cov}[Q_0,\Delta]\mathrm{Var}[\Delta]^{-1}\mathrm{Cov}[Q_0,\Delta]^\top
\]
[2310.00125].

The sample-allocation problem is then posed under costs $c_0,\dots,c_K$ as
\[
c_0N_0+\sum_{\ell=1}^K c_\ell(N_\ell+N_\ell^*)\le \text{Budget},
\]
with the objective of minimizing $\det \mathrm{Var}[\widetilde Q]$ over integer sample counts [2310.00125]. In simple scalar cases, the ratio
\[
N_\ell^*/N_0 = (|\rho_{0\ell}|/c_\ell)\sqrt{\mathrm{Var}[f_0]/\mathrm{Var}[f_\ell]}
\]
shows explicitly that cheap and highly correlated low-fidelity models should be sampled more heavily [2310.00125].

A closely related but distinct construction fuses multilevel Monte Carlo and multi-fidelity Monte Carlo. With level-dependent approximations $P_l$ to a quantity of interest $P$, MLMC uses the telescoping identity
\[
\mathbb{E}[P_L]
=
\mathbb{E}[P_1]+\sum_{l=2}^L\mathbb{E}[P_l-P_{l-1}],
\]
whereas MFMC combines model outputs through control variates. The integrated MLMF estimator is
\[
Y_{\mathrm{MLMF}}
=
P_L^{(N_L)}+\sum_{l=1}^{L-1} a_l\bigl[P_l^{(N_l)}-P_l^{(N_{l+1})}\bigr],
\qquad a_L\equiv 1.
\]
It can also be rewritten as
\[
Y_{\mathrm{MLMF}}
=
a_1P_1^{(N_1)}+\sum_{l=2}^L\bigl[a_lP_l^{(N_l)}-a_{l-1}P_{l-1}^{(N_l)}\bigr]
\]
[2510.07756].

The theory for MLMF is explicit. The unbiasedness theorem states that $\mathbb{E}[Y_{\mathrm{MLMF}}]=\mathbb{E}[P_L]$. Under the optimal coefficients
\[
a_l^*=\rho_{l,L}\,\frac{\sigma_L}{\sigma_l},
\]
the variance becomes
\[
\mathrm{Var}[Y_{\mathrm{MLMF}}]=\sum_{l=1}^L \frac{R_l}{N_l},
\qquad
R_l=V_L(\rho_{l,L}^2-\rho_{l-1,L}^2),
\]
with $\rho_{0,L}=0$. Under the stated assumptions of descending $|\rho_{l,L}|$ and increasing costs $C_l$, the estimator is unbiased, consistent, and achieves RMS error $\varepsilon$ at cost $O(\varepsilon^{-2})$ [2510.07756].

These two strands are complementary rather than competing. One emphasizes closed-form covariance blocks for complex stacked trajectory statistics; the other emphasizes unbiased telescoping across levels and models with provable $O(\varepsilon^{-2})$ cost scaling [2310.00125; 2510.07756].

## 6. Automated tuning, empirical domains, and recurrent limitations

The effectiveness of multifidelity trajectory estimators depends on the relative correlations and computational costs of the available models, and one recent direction treats low-fidelity hyperparameters themselves as optimization variables [2509.16007]. In that formulation, the target is to choose both model hyperparameters $\beta$ and sample allocations $(N,N_i)$ to minimize $\mathrm{Var}[\hat Q_{MF}]$ subject to a cost constraint:
\[
\beta^*,N^*=\arg\min_{\beta,N,N_i}\mathrm{Var}[\hat Q_{MF}](\beta,N,N_i)
\quad\text{subject to}\quad
C(\beta,N,N_i)\le C_{\text{budget}}.
\]
The equivalent penalized form uses a Lagrangian
\[
L(\beta,N,N_i,\lambda)
=
\mathrm{Var}[\hat Q_{MF}](\beta,N,N_i)+\lambda[C(\beta,N,N_i)-C_{\text{budget}}],
\]
with KKT conditions
\[
\partial \mathrm{Var}/\partial N+\lambda w_0=0,\qquad
\partial \mathrm{Var}/\partial N_i+\lambda w_i(\beta_i)=0,\qquad
\partial \mathrm{Var}/\partial \beta_i+\lambda N_i\,\partial w_i/\partial \beta_i=0
\]
[2509.16007].

The practical algorithm is bi-level. The inner loop solves the multifidelity sample-allocation problem for fixed $\beta$, while the outer loop solves for $\beta^*$ using Efficient Global Optimization, that is, Bayesian optimization with a Gaussian-process surrogate and Expected Improvement [2509.16007]. The online workflow evaluates pilot samples to estimate $\rho(\beta)$, $P(\beta)$, and $\sigma(\beta)$, then solves an approximate control variates allocation for the remaining budget. Three pilot modes are discussed—offline, projection, and online—and the paper uses the online mode to hedge pilot cost against accuracy [2509.16007].

The reported trajectory application is entry, descent, and landing for the ADEPT SR-1 flight test using NASA’s POST2 code. The uncertain inputs comprise approximately $80$ variables. Four models are used: full POST2 with $\Delta t=0.001\,\mathrm{s}$ and cost $w_0=1\times$; reduced-physics POST2 at approximately $1/5\times$; coarse-timestep POST2 at approximately $1/78\times$; and a machine-learning surrogate at approximately $1/310\,000\times$ [2509.16007]. The tunable models vary over $\Delta t\in[0.001,0.25]\,\mathrm{s}$, and the correlation-versus-cost sweep shows that as $\Delta t$ increases, cost decreases monotonically while correlations can be nonmonotonic, including a sharp drop around $\Delta t\approx 0.1\,\mathrm{s}$ [2509.16007]. This directly supports the paper’s recommendation that hand-tuning $\Delta t$ is unreliable and that models that look “very inaccurate” may still be useful if correlations remain high.

The budget dependence of tuning overhead is also explicit. For single-parameter tuning of the time-of-flight quantity of interest, with budgets $C=\{500,1000,2000\}$, pilot sizes $N_{\text{pilot}}\in\{10,50,100\}$, and EGO iterations $N_{\text{iter}}\in\{5,10,20\}$ across $100$ random trials, tuned ACV lies between hand-tuned and oracle ACV for all budgets and strictly beats hand-tuned ACV at $C\ge 1000$ [2509.16007]. For two-parameter tuning up to $C=10\,000$, the tuning overhead falls below $5\%$ once $C\ge 5000$, and for $C\ge 5000$ tuned estimators consistently outperform hand-tuned ones and approach the oracle solution; by contrast, at low $C$, heavy tuning with large $N_{\text{pilot}}$ or $N_{\text{iter}}$ can degrade performance [2509.16007]. A parametric ACV enhancement based on a search over $16$ admissible recursion trees in the Generalized Multifidelity framework yields up to $50\%$ more variance reduction for the landing-latitude quantity of interest at $C=5000$, $N_{\text{iter}}=10$ [2509.16007].

The diversity of empirical domains is summarized below.

| Setting | Estimator mechanism | Reported outcome |
|---|---|---|
| Parametric ODE moments | Three-level surrogate + spline lifting + Richardson extrapolation [1808.03379] | Cost $1000\times X$ vs $\approx 295\,X$ with order-of-magnitude lower $\ell^2_k$ error |
| Flow past a cylinder | BMFMC with DG $k=6$ HF and $k=3$ LF [2001.02892] | KL divergence of order $10^{-3}$ with $n=150$ HF runs; speed-up $\approx 7\times$ |
| Bending-wall FSI | BMFMC on tip deflection trajectory or peak [2001.02892] | $n=50$ HF runs; speed-up $>23\times$ |
| 5-dimensional Ornstein–Uhlenbeck risk estimation | MLMF with feature-level hierarchy [2510.07756] | $5\times$ lower cost at the same error |
| Path-integral control | MLMF with surrogate dynamics [2510.07756] | Optimal cost recovered at one-third the sampling cost |
| EDL trajectory UQ | MOACV and tuned ACV/GMF [2310.00125; 2509.16007] | $\sim 7.5\times$ variance reduction for touchdown latitude; up to $50\%$ more variance reduction with GMF tuning |

Several limitations recur across the literature. Heavy pilot or tuning cost can be detrimental when the total budget is small [2509.16007]. Multi-output ACV requires $\mathrm{Cov}[\Delta]$ to be invertible, so redundant statistics may need to be removed or regularized [2310.00125]. Higher-dimensional hyperparameter spaces are identified as challenging and may require more advanced derivative-free or cost-aware optimizers [2509.16007]. More broadly, a plausible implication is that practical success depends less on any single fidelity hierarchy than on the joint management of correlation structure, estimator covariance, and budget partitioning across exploration, training, and exploitation.

Source: https://www.emergentmind.com/topics/multifidelity-trajectory-simulation-estimators