SEARCH · Search NASA

Results for “low order model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗

DOE OSTI2020

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Magic-angle semimetals with chiral symmetry

Here, we construct and solve a two-dimensional, chirally symmetric model of Dirac cones subjected to a quasiperiodic modulation. In real space, this is realized with a quasiperiodic hopping term. This hopping model, as we show, at the Dirac node energy has a rich phase diagram with a semimetal-to-metal phase transition at intermediate amplitude of the quasiperiodic modulation, and a transition to a phase with a diverging density of states (DOS) and subdiffusive transport when the quasiperiodic hopping is strongest. We further demonstrate that the semimetal-to-metal phase transition can be characterized by the multifractal structure of eigenstates in momentum space and can be considered as a unique “unfreezing” transition. This unfreezing transition in momentum space generates flat bands with a dramatically renormalized bandwidth in the metallic phase similar to the phenomena of the band structure of twisted bilayer graphene at the magic angle. We characterize the nature of this transition numerically as well as analytically in terms of the formation of a band of topological zero modes. For pure quasiperiodic hopping, we provide strong numerical evidence that the low-energy DOS develops a divergence and the eigenstates exhibit Chalker (quantum-critical) scaling despite the model not being random. At particular commensurate limits the model realizes higher-order topological insulating phases. We discuss how these systems can be realized in experiments on ultracold atoms and metamaterials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Cost Efficient and Highly Weather-Resistant Solar Panel Backsheet Produced through Continuous Co-Extrusion Processing

The efficient generation and delivery of ‘clean energy” through photovoltaic (PV) technology relies heavily on the cost, performance and reliability of solar panels, along with the individual components used in assembly of those panels. In order to continue towards the Department of Energy’s (DOE) goals of driving costs down, the lifetime of high-output power generation equipment (i.e. solar modules) must be extended to reduce costs associated with power output degradation and replacement of failed or degraded panels. One of the major sources of module failures in the past 30 years has been the failure of the solar panel backsheet. Failure of this component causes severe output losses due to oxidation and yellowing of the module while creating safety concerns as the major electrical insulation of the module is compromised. Tomark-Worthen LLC was formed in 2012 with the goal of designing and manufacturing novel PV backsheets and encapsulants that would be produced in the United States of America and exceed expectations of domestic and international module producers. These backsheets would improve upon known backsheet design failures (i.e. Isovoltaic’s AAA backsheet) and create domestic manufacturing jobs in the PV sector. This goal became a reality with a first-generation product in 2018 and now, with the support of the DOE, Tomark-Worthen has launched a second-generation product known as PhotoMark® Reflections™ 205-3 and 360-3. This 2+ year effort began in late 2017 with a screening of potential polymeric materials and backsheet structural designs. Performance evaluations of each material provided the necessary data needed to select the proper material combinations and move forward into the design of the manufacturing process. These decisions were made by teaming with a local university, well-respected government laboratories, and private research organizations who provided the tools necessary to predict long-term performance of these materials in much shorter periods of time. Once the go/no-go criteria were met at the mid-point of the project, the effort transitioned to predicting 30+ year performance, achieving certifications, and improving manufacturing efficiencies designed to lower costs in order to be competitive in a cost-driven market. Major project accomplishments include recognition of the new backsheet by Underwriters Laboratories (UL) and international certification for 1000V and 1500V modules by TUV Sud. Both 1000V and 1500V products have been proven on several module manufacturing lines and initial customer orders have been received by both domestic and international module producers. Accelerated exposure testing has shown that the performance of this backsheet exceeds the prior polyamide-based backsheet known to fail in the field. Cost-models and initial customer orders have shown the ability to produce this backsheet at a cost acceptable to many domestic manufacturers, while International costs for standard backsheet remains extremely low. Tomark-Worthen remains one of the only U.S. companies manufacturing backsheet domestically and this has drawn interest from many domestic module manufacturers looking to improve their material supply logistics by avoiding long lead times and custom’s headaches while avoiding current tariffs on Chinese-made components. Along with the launch of a new PV backsheet, this project has provided a critical deliverable to the PV R&D community. This deliverable is the knowledge that not all polyamide-containing backsheets are doomed to failure as the earlier versions did. There has been a fear amongst the PV community to consider these materials due to the failure of one design. This effort has brought a new level of interest to the scientific community as it has been shown that polyamides can be used successfully in a well-designed backsheet. As we move forward from this project, we are excited to see how backsheet technology and manufacturing in the U.S. will continue to progress and grow.

14 SOLAR ENERGY↗

DOE OSTI2021

Variational quantum solver employing the PDS energy functional

In our previous work (J. Chem. Phys. 2020, 153, 201102),we reported a new class of quantum algorithms that are based on the quantum computation of the connected moment expansion to find the ground and excited state energies. In particular, the Peeters-Devreese-Soldatov (PDS) formulation is found variational and bearing the potential for further combining with the existing variational quantum infrastructure. Following this direction, here we propose a variational quantum solver employing the PDS energy gradient. In comparison with the usual variational quantum eigensolver (VQE) and the original static PDS approach, the proposed variational quantum solver offers an effective approach to achieve high ac-curacy at finding the ground state and its energy through the rotation of the trial wave function of modest quality guided by the low order PDS energy gradients, thus improves the ac-curacy and efficiency of the quantum simulation. We demonstrate the performance of the proposed variational quantum solver for toy models, H2molecule, and strongly correlated planar H4system in some challenging situations. In all the case studies, the proposed variational quantum approach out-performs the usual VQE and static PDS calculations even at the lowest order.

Peng, Bo↗

DOE OSTI2021

Neat, Simple, and Wrong: Debunking Electrostatic Fallacies Regarding Noncovalent Interactions

Multipole moments such as charge, dipole, and quadrupole are often invoked to rationalize intermolecular phenomena, but a low-order multipole expansion is rarely a valid description of electrostatics at the length scales that characterize nonbonded interactions. This is illustrated by examining several common misunderstandings rooted in erroneous electrostatic arguments. First, the notion that steric repulsion originates in Coulomb interactions is easily disproved by dissecting the interaction potential for Ar2. Second, the Hunter-Sanders model of π–π interactions, which is based on quadrupolar electrostatics, is shown to have no basis in accurate calculations. Third, curved “buckybowls” exhibit unusually large dipole moments, but these are ancillary to the forces that control their intermolecular interactions, as illustrated by two examples involving corannulene. Finally, the assumption that interactions between water and small anions are dictated by the dipole moment of H 2 O is shown to be false in the case of binary halide–water complexes. Furthermore, these examples present a compelling case that electrostatic explanations based on low-order multipole moments are very often counterfactual for nonbonded interactions at close range and should not be taken seriously in the absence of additional justification.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Multifractally-enhanced superconductivity in thin films

The multifractal superconducting state originates from the interplay of Anderson localization and interaction effects. In this article we overview the recent theory of the superconductivity enhancement by multifractality and extend it to describe the spectral properties of superconductors on the scales of the order of the superconducting gap. Specifically, using the approach based on renormalization group within the nonlinear sigma model, we develop the theory of a multifractal superconducting state in thin films. We derive a modified Usadel equation that incorporates the interplay of disorder and interactions at energy scales larger than the spectral gap and study the effect of such an interplay on the low-energy physics. We determine the spectral gap at zero temperature which occurs to be proportional to the multifractally enhanced superconducting transition temperature. The modified Usadel equation results in the disorder-averaged density of states that, near the spectral gap, resembles the one obtained in the model of a spatially random superconducting order parameter. We reveal strong mesoscopic fluctuations of the local density of states in the superconducting state. Such strong mesoscopic fluctuations imply that the interval of energies in which the superconducting gap establishes is parametrically large in systems with multifractally-enhanced superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

DOE OSTI2023

Advances in 3D Geologic Modeling of Alluvial Basins with a Focus on Facies and Property Modeling

The unsaturated zone alluvium reference case is one of several geologic systems under consideration by the U.S. Department of Energy Office of Nuclear Energy for hosting repositories for spent nuclear fuel and associated waste (Sevougian et al., 2019). As noted by Mariner et al. (2018), the generic alluvial basin offers positive attributes that merit its consideration as a reference case by the Spent Fuel and Waste Science and Technology (SFWST) campaign. There are hundreds of alluvial basins and sub-basins scattered across the arid western United States (Figure 1-1). Precipitation and infiltration rates are relatively low with high evapotranspiration, resulting in vertical separation between repository and water table and thus longer transport paths to an aquifer. Accumulations of alluvial sediments within these basins are typically on the order of hundreds of meters, and locally may exceed 1,000 m in thickness, as is the case for the Deming sub-basin in southern New Mexico. A thick geologic host medium, which serves as the natural barrier system (NBS) in the conceptual model framework of a geologic disposal system, isolates the waste packages from receptors in the biosphere. Further, alluvial basin fill is typically comprised of stacked playa and lacustrine deposits along the basin axis (Perry et al., 2018). Characterized by low permeability, these layers protect the biosphere above the repository and the groundwater resources below the repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Local positional and spin symmetry breaking as a source of magnetism and insulation in paramagnetic EuTi O 3

Here we consider theoretically the paramagnetic phases of EuTiO 3 that represent configurations created by two sets of microscopic degrees of freedom (m-DOFs): positional symmetry breaking due to octahedral rotations and magnetic symmetry breaking due to spin disorder. The effect of these sets of m-DOFs on the electronic structure and properties of the para phases is assessed by considering sufficiently large (super) cells with the required nominal global average symmetry, allowing, however, the $\mathcal{local}$ positional and magnetic symmetries to be lowered. We find that tendencies for local symmetry breaking can be monitored by following total energy lowering in mean-fieldlike density-functional theory, without recourse for strong correlation effects. While most nominally cubic ABO 3 perovskites are known for their symmetry breaking due to the B-atom sublattice, the case of f-electron magnetism in EuTiO 3 is associated with A-sublattice symmetry breaking and its coupling to structural distortions. We find that (i) paramagnetic $\mathcal{cubic}$ EuTiO 3 has an intrinsic tendency for both magnetic and positional symmetry breaking, while paramagnetic $\mathcal{tetragonal}$ EuTiO 3 has only magnetic symmetry lowering and no noticeable positional symmetry lowering with respect to low-temperature antiferromagnetic tetragonal phase. (ii) Properly modeled paramagnetic tetragonal and cubic EuTiO 3 have a nonzero local magnetic moment on each Eu ion, consistent with the experimental observations of local magnetism in the para phases of EuTiO 3 significantly above the Néel temperature. Interestingly, (iii) the local positional distortion modes in the short-range ordered para phases are inherited from the long-range ordered low-temperature antiferromagnetic ground-state phase.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Topological order in matrix Ising models

We study a family of models for an N_1 \times N_2 N 1 × N 2 matrix worth of Ising spins S_{aB} S a B . In the large N_i N i limit we show that the spins soften, so that the partition function is described by a bosonic matrix integral with a single ‘spherical’ constraint. In this way we generalize the results of to a wide class of Ising Hamiltonians with O(N_1,\mathbb{Z})\times O(N_2,\mathbb{Z}) O ( N 1 , ℤ ) × O ( N 2 , ℤ ) symmetry. The models can undergo topological large N N phase transitions in which the thermal expectation value of the distribution of singular values of the matrix S_{aB} S a B becomes disconnected. This topological transition competes with low temperature glassy and magnetically ordered phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

A machine-learning inverse model framework for rapid forecasting and history matching in unconventional reservoirs

Model-based optimization for real-time forecasting in unconventional reser-voirs requires novel methods and work?ows since the strategies and work?ows used in conventional reservoirs are either inapplicable, or prohibitively expen-sive and time-consuming. Insu?cient site data and computational expense of high-?delity simulations mean that work?ows with high-?delity simulations are not ideal for usage in comprehensive uncertainty quanti?cation stud-ies that require 1000s of forward model runs. We present an alternative, novel work?ow for unconventional reservoirs, based on the interplay between reduced-order models and machine-learning. Our physics-informed machine-learning (PIML) work?ow addresses the challenges to real-time reservoir management in uncoventionals, namely lack of data (the time-frame for which the wells have been producing), and computational expense of high-?delity modeling. We use the machine-learning paradigm of transfer-learning to bind together fast but less accurate reduced-order models with slow, but accurate high-?delity models and circumvent the di?culties inherent in the current state-of-the-art for unconventionals. Such a PIML work?ow, grounded in physics, is a viable candidate for real-time history matching and production forecasting in a fractured shale gas reservoir. The signi?cance of our approach is that while it is developed for a particu-lar well and site in the Marcelus Shale gas reservoir of the Appalachian basin (MSEEL), it is not wedded to it. We expect the same work?ow can be ap-plied to other shale formations (e.g., Woodford, Barnett, Utica, EagleFord) should site-data become available, using the same set of machine-learning techniques from transfer learning. Some ?ne-tuning (or minimal retraining of the neural networks) will be required to transfer knowledge across shale gas sites/formations but it is a clearly superior alternative to developing a new machine-learning model altogether when considering a di?erent site.

Srinivasan, Shriram↗

DOE OSTI2024

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Re-evaluation of the sequential 3 α reaction in stellar conditions

The formation of carbon via the triple alpha process is a key step in stellar production reactions of light ions and heavily impacts the life-cycle of red giant stars. Here, we propose a novel theoretical approach for sequential mechanism of the helium burning process, based on the Imaginary Time Method (ITM) and semiclassical models. We compare our results to those of the NACRE collaboration, which has set the standard for the past decades. Our calculations differ by several orders of magnitude in the low temperature region, which may seriously affect the early evolution of helium burning stars.

Depastas, Theodoros↗

DOE OSTI2025

Confronting new NICER mass-radius measurements with phase transition in dense matter and twin compact stars

The (re)analysis of data on the X-ray emitting pulsars PSR J0030+0451 and J0740+6620, as well as new results on PSR J0437-4715 and J1231-1411, are confronted with the predictions of the equation of state (EoS) models allowing for strong first-order phase transition for the mass-radius (M-R) diagram. Here, we use models that are based on a covariant density functional (CDF) EoS for nucleonic matter at low densities and a quark matter EoS, parameterized by the speed of sound, at higher densities. To account for the variations in the ellipses for PSR J0030+0451 obtained from different analyses, we examined three scenarios to assess their consistency with our models, focusing particularly on the potential formation of twin stars. We found that in two scenarios, where the ellipses for PSR J0030+0451 and J0437-4715 with masses close to the canonical mass ∼ 1.4 M ⊙ are significantly separated, our models allow for the presence of twin stars as a natural explanation for potential differences in the radii of these stars.

X-ray binaries↗

DOE OSTI2024

Calculating the two-photon exchange contribution to K L → μ + μ − decay

We present a theoretical framework within which both the real and imaginary parts of the complex, two-photon exchange amplitude contributing to K L → μ + μ − decay can be calculated using lattice quantum chromodynamics. The real part of this two-photon amplitude is of approximately the same size as that coming from a second-order weak strangeness-changing neutral-current process. Thus a test of the standard model prediction for this second-order weak process depends on an accurate result of this two-photon amplitude. A limiting factor of our proposed method comes from low-energy three-particle π π γ states. The contribution from these states will be significantly distorted by the finite volume of our calculation—a distortion for which there is no available correction. However, a simple estimate of the contribution of these three-particle states suggests their contribution to be at most a few percent allowing their neglect in a lattice calculation with a 10% target accuracy. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2025

An evaluation of multi-fidelity methods for quantifying uncertainty in projections of ice-sheet mass change

Abstract. This study investigated the computational benefits of using multi-fidelity statistical estimation (MFSE) algorithms to quantify uncertainty in the mass change of Humboldt Glacier, Greenland, between 2007 and 2100 using a single climate change scenario. The goal of this study was to determine whether MFSE can use multiple models of varying cost and accuracy to reduce the computational cost of estimating the mean and variance of the projected mass change of a glacier. The problem size and complexity were chosen to reflect the challenges posed by future continental-scale studies while still facilitating a computationally feasible investigation of MFSE methods. When quantifying uncertainty introduced by a high-dimensional parameterization of the basal friction field, MFSE was able to reduce the mean-squared error in the estimates of the statistics by well over an order of magnitude when compared to a single-fidelity approach that only used the highest-fidelity model. This significant reduction in computational cost was achieved despite the low-fidelity models used being incapable of capturing the local features of the ice-flow fields predicted by the high-fidelity model. The MFSE algorithms were able to effectively leverage the high correlation between each model's predictions of mass change, which all responded similarly to perturbations in the model inputs. Consequently, our results suggest that MFSE could be highly useful for reducing the cost of computing continental-scale probabilistic projections of sea-level rise due to ice-sheet mass change.

54 ENVIRONMENTAL SCIENCES↗