Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Low-dimensional Space”

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.

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

At least 73 records · Page 4

Advancing the State-of-the-Practice for Liquid Rocket Engine Injector Design

Current shortcomings in both the overall injector design process and its underlying combustion stability assessment methodology are rooted in the use of empirically based or low fidelity representations of complex physical phenomena and geometry details that have first order effects on performance, thermal environments and combustion stability. The result is a design and analysis capability that is often inadequate to reliably arrive at a suitable injector design in an efficient manner. Specifically, combustion instability has been particularly difficult to predict and mitigate. Large hydrocarbon-fueled booster engines have been especially problematic in this regard. Where combustion instability has been a problem, costly and time-consuming redesign efforts have often been an unfortunate consequence. This paper presents an overview of a recently completed effort at NASA Marshall Space Flight Center to advance the state-of-the-practice for liquid rocket engine injector design. Multiple perturbations of a gas-centered swirl coaxial (GCSC) element that burned gaseous oxygen and RP-1 were designed, assessed for combustion stability, and tested. Three designs, one stable, one marginally unstable and one unstable, were used to demonstrate both an enhanced overall injector design process and an improved combustion stability assessment process. High-fidelity results from state-of-the-art computational fluid dynamics CFD simulations were used to substantially augment and improve the injector design methodology. The CFD results were used to inform and guide the overall injector design process. They were also used to upgrade selected empirical or low-dimensional quantities in the ROCket Combustor Interactive Design (ROCCID) stability assessment tool. Hot fire single element injector testing was used to verify both the overall injector designs and the stability assessments. Testing was conducted at the Air Force Research Laboratory and at Purdue University. Companion papers provide details of the overall injector design process, full- and sub-scale testing, ROCCID-based stability assessments and the CFD simulations.

Tucker, P. K.↗

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗

Laboratory and theoretical models of planetary-scale instabilities and waves

Research work is proceeding in theoretical, numerical, and experimental geophysical fluid dynamics leading up to a reflight of the GFFC (Geophysical Fluid Flow Cell Experiment) on USML-2. The work is intended not only to generate ideas for future space experiments, but to provide fundamental results concerned with nonlinear and chaotic properties of thermal convection and baroclinic waves in terrestrial and planetary atmospheres. The major efforts are focussed on thermal convection in a rapidly rotating annulus relevant to Jovian atmospheric dynamics, and on the chaotic behavior of baroclinic waves relevant to the Earth's atmosphere. The approach, in preparation for USML-2, is primarily theoretical and numerical. Mechanistic process models are solved numerically in order to identify physical mechanisms that may be observed in the GFFC, and which are important in real geophysical applications. The results from numerical simulations of geophysical fluid flow (subject to rotation and stratification) are compared with previous GFFC experiments on Spacelab-3 and with existing and proposed terrestrial laboratory experiments of various types. Pattern recognition algorithms have been employed to generate low-dimensional descriptions of the highly nonlinear and turbulent numerical simulations. Such empirically truncated descriptions provide for simplified but robust physical interpretations of the dynamics, as well as yielding highly efficient computations of these chaotic flows.

Hart, John E.↗

Pressure‐Induced Structural and Magnetic Evolution in Layered Antiferromagnet YbMn 2 Sb 2

Electronic states under pressure exhibit unconventional spin and charge dynamics that provide a powerful route to uncover exotic phases in quantum materials. Here, we present the structural, magnetic, and electronic evolution of YbMn 2 Sb 2 under pressure. Single-crystal X-ray diffraction reveals a pressure-induced structural transition from the space group trigonal P$\overline{3}$m1 to the monoclinic P2 1 /m phase near 3.5 GPa, which remains stable up to 10 GPa. Magnetization measurements display an anomalously weak net magnetic moment and the absence of Curie–Weiss behavior up to 400 K, suggesting the formation of short-range Mn moment pairs that cancel macroscopically and subsequently evolve into long-range order upon cooling. Temperature-dependent resistivity shows semiconducting behavior with a transition at ∼119 K at ambient pressure, while pressure induces a dramatic suppression of resistance and the emergence of metallic-like temperature dependence, stabilized beyond 5 GPa. This pressure-driven semiconductor-metal transition is consistent with our density functional theory calculations, confirming the closing of the band gap under compression. Neutron diffraction under pressure identifies an incommensurate magnetic structure with antiparallel correlations between paired spins. Together, these results demonstrate how pressure-driven structural tuning and competing exchange interactions stabilize unconventional magnetic states in this low-dimensional magnetic semiconductor.

antiferromagnetic semiconductor↗

Fault-tolerant resource comparison of qudit and qubit encodings for diagonal quadratic operators

Finite local Hilbert-space truncations arise naturally in quantum simulations of lattice field theories and motivate qudit encodings, but their fault-tolerant advantage over qubit encodings remains unclear. We compare the non-Clifford cost of implementing quadratic diagonal evolutions, exemplified by 𝑈 = 𝑒$^{−𝑖⁢𝑡⁢𝜙^2_𝑥}$ in a uniform field-amplitude discretization of a real scalar field, using either one logical 𝑑-level qudit or 𝑛 𝑏 = ⌈log 2⁡ 𝑑⌉ logical qubits. We analyze two standard settings: product-formula simulation and linear combination of unitaries (LCU) per block encoding, taking the resource metric to be the number of non-Clifford gates after synthesis into a discrete logical gate set. Because tight synthesis bounds for general single-qudit rotations are not known, we express the qudit constructions in terms of embedded two-level SU⁡(2) rotations and derive explicit finite-𝑑 break-even conditions for their synthesis cost; these serve as compiler targets for when qudit encodings can outperform the qubit baseline. Within the constructive models studied here, product-formula implementations would require an exponentially stronger per-primitive synthesis advantage for qudits to win asymptotically, while in the LCU setting the qubit encoding is asymptotically cheaper in 𝑑. Nevertheless, the finite-𝑑 threshold analysis identifies low-dimensional regions in which qudits can yield meaningful constant-factor savings, particularly for LCU-based implementations. As a secondary analysis of the LCU construction, we use an idealized negligible-overhead qubit-qudit code-switching model to give an absolute 𝑇-count comparison and reinterpret the savings as an allowable per-switch overhead budget.

Godwood, Samuel [Univ. of Liverpool (United Kingdo↗

Anisotropic light-tailored RKKY interaction in two-dimensional 𝑑-wave altermagnets

Altermagnets are known in spintronics for their intrinsic spin-splitting and unconventional magnetic responses, particularly to magnetic impurities. However, effectively controlling the magnetic exchange interactions in altermagnets is challenging for practical applications. Here, in this work, we propose using circularly polarized light to tune the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction in two-dimensional 𝑑-wave altermagnets. Using the real-space retarded Green's functions approach, our results show that while the Heisenberg and Ising exchanges dominate, a notable Dzyaloshinskii–Moriya (DM) interaction also plays a key role. Furthermore, the inherent strength of altermagnetism imprints chirp-like signatures into the magnetic responses, which can be dynamically tuned via light. We mainly demonstrate that gate-induced Rashba spin-orbit coupling is essential in response to light—light selectively and anisotropically adjusts the DM interaction without affecting the other exchanges. Our findings further indicate that rotating the altermagnet by 45° relative to the light's polarization direction generates a Dirac-like dispersion and different DM interactions. We finally extract critical thresholds where light reverses DM interactions along one axis or balances both in-plane components. The anisotropic light-driven control of RKKY interactions in altermagnets not only highlights their unique properties but also opens new avenues for engineering tailored magnetic characteristics in spintronic applications.

altermagnetism↗

Generalization error guaranteed auto-encoder-based nonlinear model reduction for operator learning

Many physical processes in science and engineering are naturally represented by operators between infinite-dimensional function spaces. The problem of operator learning, in this context, seeks to extract these physical processes from empirical data, which is challenging due to the infinite or high dimensionality of data. An integral component in addressing this challenge is model reduction, which reduces both the data dimensionality and problem size. In this paper, we utilize low-dimensional nonlinear structures in model reduction by investigating Auto-Encoder-based Neural Network (AENet). AENet first learns the latent variables of the input data and then learns the transformation from these latent variables to corresponding output data. Our numerical experiments validate the ability of AENet to accurately learn the solution operator of nonlinear partial differential equations. Furthermore, we establish a mathematical and statistical estimation theory that analyzes the generalization error of AENet. Finally, our theoretical framework shows that the sample complexity of training AENet is intricately tied to the intrinsic dimension of the modeled process, while also demonstrating the robustness of AENet to noise.

Auto-encoder↗

Polynomial chaos expansions on principal geodesic Grassmannian submanifolds for surrogate modeling and uncertainty quantification

In this work we introduce a manifold learning-based surrogate modeling framework for uncertainty quantification in high-dimensional stochastic systems. Our first goal is to perform data mining on the available simulation data to identify a set of low-dimensional (latent) descriptors that efficiently parameterize the response of the high-dimensional computational model. To this end, we employ Principal Geodesic Analysis on the Grassmann manifold of the response to identify a set of disjoint principal geodesic submanifolds, of possibly different dimension, that captures the variation in the data. Since operations on the Grassmann require the data to be concentrated, we propose an adaptive algorithm based on Riemannian K-means and the minimization of the sample Fréchet variance on the Grassmann manifold to identify “local” principal geodesic submanifolds that represent different system behavior across the parameter space. Polynomial chaos expansion is then used to construct a mapping between the random input parameters and the projection of the response on these local principal geodesic submanifolds. Here, the method is demonstrated on four test cases, a toy-example that involves points on a hypersphere, a Lotka-Volterra dynamical system, a continuous-flow stirred-tank chemical reactor system, and a two-dimensional Rayleigh-Bénard convection problem.

42 ENGINEERING↗

Orbital-Selective Instabilities and Spin Fluctuations at the Verge of Superconductivity in Interlayer-Expanded Iron Selenide

Understanding electron correlation-driven instabilities and their coupling to structural phases is essential for deciphering multiorbital pairing in unconventional superconductors. We investigate Li x (C 5 H 5 N) y Fe 2 Se 2 (x ∼ 0.6; y ∼ 0.7−0.9), a tetragonal β-FeSe intercalate with a superconducting transition temperature (T c = 39 K) closely tied to an expanded Fe-layer spacing (∼11.4 Å). High-resolution synchrotron Xray diffraction and core-level absorption spectroscopy reveal subtle lattice distortions on cooling without a symmetry-breaking transition. Instead, the material exhibits negative thermal expansion (NTE) in the two-dimensional Fe network below T S ∼ 70 K, and stiffening of local Se−Fe−Se bond dynamics near T c . The spatially incoherent rearrangement of FeSe 4 tetrahedra and the site-local fluctuations, signal reduced electron correlations compared to those of parent β-FeSe (T c = 8 K). Complementary X-ray emission spectroscopy, a fast local probe of Fe 3d valence states, detects persistent local Fe spin moments below T S , unlike quenching in related systems. These findings indicate that decoupling of Fe planes leads to an electronically driven lattice instability. The latter emerges as NTE induced from weak, orbital-selective localization of in-plane Fe 3d states rather than conventional transverse vibrations. Governed by Hund’s coupling, this selectivity permits coexistence of local spin fluctuations with itinerant d-electrons critical for enhancing T c . These results suggest that intercalation-driven d-orbital differentiation moderates electron correlations, providing a pathway to optimize the superconductivity in low-dimensional quantum materials.

36 MATERIALS SCIENCE↗

Dark Energy Survey Year 3 results: Simulation-based 𝑤CDM inference from weak lensing and galaxy clustering maps with deep learning: Analysis design

Data-driven approaches using deep learning are emerging as powerful techniques to extract non-Gaussian information from cosmological large-scale structure. Here, this work presents the first simulation-based inference (SBI) pipeline that combines weak lensing and galaxy clustering maps in a realistic Dark Energy Survey Year 3 (DES Y3) configuration and serves as preparation for a forthcoming analysis of the survey data. We develop a scalable forward model based on the CosmoGridV1 suite of N-body simulations to generate over one million self-consistent mock realizations of DES Y3 at the map level. Leveraging this large dataset, we train deep graph convolutional neural networks on the full survey footprint in spherical geometry to learn low-dimensional features that approximately maximize mutual information with target parameters. These learned compressions enable neural density estimation of the implicit likelihood via normalizing flows in a ten-dimensional parameter space spanning cosmological 𝑤CDM, intrinsic alignment, and linear galaxy bias parameters, while marginalizing over baryonic, photometric redshift, and shear bias nuisances. To ensure robustness, we extensively validate our inference pipeline using synthetic observations derived from both systematic contaminations in our forward model and independent Buzzard galaxy catalogs. Our forecasts yield significant improvements in cosmological parameter constraints, achieving 2−3× higher figures of merit in the 𝛺 𝑚 − 𝑆 8 plane relative to our implementation of baseline two-point statistics and effectively breaking parameter degeneracies through probe combination. These results demonstrate the potential of SBI analyses powered by deep learning for upcoming Stage-IV wide-field imaging surveys.

Thomsen, A. [Zurich, ETH] (ORCID:0000000203099021)↗

La 4 Co 4 X ( X = Pb , Bi , Sb ) : A demonstration of antagonistic pairs as a route to quasi-low-dimensional ternary compounds

We outline how pairs of strongly immiscible elements, referred to here as antagonistic pairs, can be used to synthesize ternary compounds with low or quasi-reduced-dimensional motifs intrinsically built into their crystal structures. By identifying third elements that are mutually compatible with a given antagonistic pair, ternary compounds can be formed in which the third element segregates the immiscible atoms into spatially separated substructures. Quasi-low-dimensional structural units, such as sheets, chains, or clusters are a natural consequence of the immiscible atoms seeking to avoid close contact in the solid state. Further, as proof of principle, we present the discovery, crystal growth, and basic physical properties of La 4 ⁢Co 4 ⁢$\mathrm{X}$ (X = Pb, Bi, Sb), a family of intermetallic compounds based on the antagonistic pairs Co-Pb and Co-Bi. La 4 ⁢Co 4 ⁢$\mathrm{X}$ adopts an orthorhombic crystal structure (space group Pbam) containing quasi-two-dimensional Co slabs and La-X polyhedra that stack in an alternating manner along the α axis. Consistent with our proposal, the La atoms separate the Co and X substructures, ensuring there are no direct contacts between the members of the immiscible (antagonistic) pair. Within the Co slabs, the atoms occupy the vertices of corner sharing tetrahedra and triangles, and this bonding motif produces narrow electronic bands near the Fermi level that favor magnetism. The Co is moment bearing in each La 4 ⁢Co 4 $\mathrm{X}$ compound studied, and we show that whereas La 4 ⁢Co 4 ⁢Pb behaves as a three-dimensional antiferromagnet with T N =220K, La 4 ⁢Co 4 ⁢Bi and La 4⁢ Co 4 ⁢Sb have behavior consistent with low-dimensional magnetic coupling and ordering, with T N =153K and 143 K, respectively. In addition to the Pb-, Bi-, and Sb-based La 4 ⁢Co 4 ⁢$\mathrm{X}$ compounds, we also were likely able to produce an analogous La 4 ⁢Co 4 ⁢Sn in polycrystalline form, although we were unable to isolate single crystals. We anticipate that identifying and using mutually compatible third elements together with an antagonistic pair represents a generalizable design principle for discovering new materials and new structure types containing low-dimensional substructures.

36 MATERIALS SCIENCE↗

Data-Driven Compression of Electron-Phonon Interactions

First-principles calculations of electron interactions in materials have seen rapid progress in recent years, with electron-phonon ( e − ph ) interactions being a prime example. However, these techniques use large matrices encoding the interactions on dense momentum grids, which reduces computational efficiency and obscures interpretability. For e − ph interactions, existing interpolation techniques leverage locality in real space, but the high dimensionality of the data remains a bottleneck to balance cost and accuracy. Here we show an efficient way to compress e − ph interactions based on singular value decomposition (SVD), a widely used matrix and image compression technique. Leveraging (un)constrained SVD methods, we accurately predict material properties related to e − ph interactions—including charge mobility, spin relaxation times, band renormalization, and superconducting critical temperature—while using only a small fraction (1%–2%) of the interaction data. These findings unveil the hidden low-dimensional nature of e − ph interactions. Furthermore, they accelerate state-of-the-art first-principles e − ph calculations by about 2 orders of magnitude without sacrificing accuracy. Our Pareto-optimal parametrization of e − ph interactions can be readily generalized to electron-electron and electron-defect interactions, as well as to other couplings, advancing quantitative studies of condensed matter. Published by the American Physical Society 2024

Physics↗

Characterizing skyrmion flow phases with principal component analysis

Principal component analysis (PCA) is a powerful method that can identify patterns in large, complex data sets by constructing low-dimensional order parameters from higher-dimensional feature vectors. There are increasing efforts to use space-and-time-dependent PCA to detect transitions in nonequilibrium systems that are difficult to characterize with equilibrium methods. Here, we demonstrate that feature vectors incorporating the position and velocity information of driven skyrmions moving through random disorder permit PCA to resolve different types of disordered skyrmion motion as a function of driving force and the ratio of the Magnus force to the dissipation. Since the Magnus force creates gyroscopic motion and a finite Hall angle, skyrmions can exhibit a greater range of flow phases than what is observed in overdamped driven systems with quenched disorder. We show that in addition to identifying previously known skyrmion flow phases, PCA detects several additional phases, including different types of channel flow, moving fluids, and partially ordered states. Guided by the PCA analysis, we further characterize the disordered flow phases to elucidate the different microscopic dynamics and show that the changes in the PCA-derived order parameters can be connected to features in bulk transport measures, including the transverse and longitudinal velocity-force curves, differential conductivity, topological defect density, and changes in the skyrmion Hall angle as a function of drive. We discuss how asymmetric feature vectors can be used to improve the resolution of the PCA analysis, and how this technique can be extended to find disordered phases in other nonequilibrium systems with time-dependent dynamics.

36 MATERIALS SCIENCE↗

Data-Driven Closures and Assimilation for Stiff Multiscale Random Dynamics

Here, we introduce a data-driven and physics-informed framework for propagating uncertainty in stiff, multiscale random ordinary differential equations (RODEs) driven by correlated (colored) noise. Unlike systems subjected to Gaussian white noise, a deterministic equation for the joint probability density function (PDF) of RODE state variables does not exist in closed form. Moreover, such an equation would require as many phase-space variables as there are states in the RODE system. To alleviate this curse of dimensionality, we instead derive exact, albeit unclosed, reduced-order PDF (RoPDF) equations for low-dimensional observables/quantities of interest. The unclosed terms take the form of state-dependent conditional expectations, which are directly estimated from data at sparse observation times. However, for systems exhibiting stiff, multiscale dynamics, data sparsity introduces regression discrepancies that compound during RoPDF evolution. This is overcome by introducing a kinetic-like defect term to the RoPDF equation, which is learned by assimilating in sparse, low-fidelity RoPDF estimates. Two assimilation methods are considered, namely nudging and deep neural networks, which are successfully tested against Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗

Scalable Bayesian Physics-Informed Kolmogorov-Arnold Networks

Uncertainty quantification (UQ) plays a pivotal role in scientific machine learning, especially when surrogate models are used to approximate complex systems. Although multilayer perceptions (MLPs) are commonly employed as surrogates, they often suffer from overfitting due to their large number of parameters. Kolmogorov-Arnold networks (KANs) offer an alternative solution with fewer parameters. However, gradient-based inference methods, such as Hamiltonian Monte Carlo (HMC), may result in computational inefficiency when applied to KANs, especially for large-scale datasets, due to the high cost of back-propagation. To address these challenges, we propose a novel approach, combining the dropout Tikhonov ensemble Kalman inversion (DTEKI) with Chebyshev KANs. This gradient-free method effectively mitigates overfitting and enhances numerical stability. In addition, we incorporate the active subspace method to reduce the parameter-space dimensionality, allowing us to improve the accuracy of predictions and obtain more reliable uncertainty estimates. Extensive experiments demonstrate the efficacy of our approach in various test cases, including scenarios with large datasets and high noise levels. Our results show that the new method achieves comparable or better accuracy, much higher efficiency as well as stability compared to HMC, in addition to scalability. Moreover, by leveraging the low-dimensional parameter subspace, our method preserves prediction accuracy while substantially reducing further the computational cost.

97 MATHEMATICS AND COMPUTING↗

Uncertainty Quantification and Sensitivity Analysis of Low-Dimensional Manifold via Co-Kurtosis PCA in Combustion Modeling

For multi-scale multi-physics applications e.g., the turbulent combustion code Pele, robust and accurate dimensionality reduction is crucial to solving problems at exascale and beyond. A recently developed technique, Co-Kurtosis based Principal Component Analysis (CoK-PCA) which leverages principal vectors of co-kurtosis, is a promising alternative to traditional PCA for complex chemical systems. To improve the effectiveness of this approach, we employ Artificial Neural Networks for reconstructing thermo-chemical scalars, species production rates, and overall heat release rates corresponding to the full state space. Our focus is on bolstering confidence in this deep learning based non-linear reconstruction through Uncertainty Quantification (UQ) and Sensitivity Analysis (SA). UQ involves quantifying uncertainties in inputs and outputs, while SA identifies influential inputs. One of the noteworthy challenges is the computational expense inherent in both endeavors. To address this, we employ the Monte Carlo methods to effectively quantify and propagate uncertainties in our reduced spaces while managing computational demands. Our research carries profound implications not only for the realm of combustion modeling but also for a broader audience in UQ. By showcasing the reliability and robustness of CoK-PCA in dimensionality reduction and deep learning predictions, we empower researchers and decision-makers to navigate complex combustion systems with greater confidence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Indications of low dimensionality in magnetospheric dynamics

Using three separate but related approaches, the question of whether the dynamic response of the magnetosphere to the solar wind input may be described by a low-order system of equations is examined. First, it is determined that the dimension of the subset (the attractor) in the high-dimensional magnetospheric phase space associated with the westward auroral electrojet (AL) index for some of the data sets compiled by Bargatze et al. (1985) is 4.0 + or - 0.2, seemingly independent of activity level. Second, direct modeling of the magnetosphere, considering the bulk properties of the tail plasma, leads to a system of equations that is similar to those previously reported as a dripping faucet model; here, the focus is specifically on the prediction of a natural frequency in this model. Finally, a peak is identified with the predicted frequency in power spectra of AL computed for intervals with both low and high activity. Peaks at other frequencies also appear in the spectra, and such resonances would be expected for a chaotic nonlinear oscillator. Combining these approaches it is concluded that at least some aspects of magnetospheric dynamics may be meaningfully modeled by low-dimensional sets of equations.

Roberts, D. A.↗

Two Reduced Resolution Filter Approaches to Data Assimilation

In this paper we evaluate the performance of two reduced resolution filter approaches to data assimilation. The main distinction between these approaches is in the manner they propagate error covariances. Both account for error covariances in a space with dimension m smaller than the model's state vector dimension n. In the first approach the m dimensional error covariance matrix is interpolated to the n-dimensional space and propagated with the n-dimensional dynamics. In the second approach the low-dimensional error covariance matrix is propagated by a dynamical operator generated in the m-dimensional space. Our experiments indicate that the first approach provides a more reliable simplified scheme for error covariance propagation than the second approach.

Todling, Ricardo↗