Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduced Order Models”

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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Estimation of Timing Performance with Pulsed Laser Links for Small Satellite PNT

Satellite constellations are growing in size and in utility for applications as diverse as providing reliable, low-latency internet service to rural areas orEarth science missions.This expansion of satellite constellations brings attention to the need for responsive and reliable satellite communication.Current research in the PrecisionSpace Systems Laboratory at the University of Florida focuses on using pulsed laser communications as the method of relaying time signals instead of the radio frequency (RF)signals.The upcoming CLICK (CubeSat Laser Infrared CrosslinK) missions (Massachusetts Institute of Technology, University of Florida, NASA Ames Research Center) will demonstrate a spacecraft-to-spacecraft laser communication link and time-transfer.For pre-flight technology assessment of missions such as CLICK, a simulation of a spacecraft constellation in Earth orbit has been created to estimate clock synchronization and precision orbit determination based on measured instrumentation performance.We have developed a novel MATLAB-based numerical simulator to model spacecraft-to-spacecraft laser time-transfer and estimate the offset between the spacecraft clocks over time.This simulation includes timing errors associated with laser pulse generation and detection, as well as non-Gaussian clock drift models.The two on-board clocks modeled are a cesium-based Chip-Scale Atomic Clock (CSAC) and a rubidium-based Miniature Atomic Clock (MAC), both produced by Microchip.An example case of two spacecraft in a circular, low-Earth orbit receiving GPS position, GPS timing, and laser pulse time-of-flight measurements is simulated.The positions and velocities of the two spacecraft at a reference epoch and the constant clock model coefficients are estimated.Polynomial models of different orders were used as clock models. The effect of clock model order on the root-mean-square (RMS) of the clock error is apparent in the case using GPS and lasing measurements, showing that the clock model improves with increasing clock model order. Results compare the estimated clock model of a mission operation that only uses GPS measurements and a mission operation that uses both GPS and laser pulse time-of-flight measurements between spacecraft referenced to their on-board CSACs or on-board MACs. Including lasing measurements reduces the RMS clock model error by approximately 80% of the RMS of the cases with only GPS measurements.This simulation tool can be used to optimize the lasing operations schedule based on mission timing performance objectives. During this presentation, I will discuss the details of the simulation and clock model estimation, as well as discuss the results.I will compare the results of incorporating lasing timing measurements to using only GPS measurements and the outcomes of using different orders of polynomials.

Timing Performance↗
DOE OSTI2024

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A multifidelity method for a nonlocal diffusion model

Nonlocal models feature a finite length scale, referred to as the horizon, such that points separated by a distance smaller than the horizon interact with each other. Such models have proven to be useful in a variety of settings. However, due to the reduced sparsity of discretizations, they are also generally computationally more expensive compared to their local differential equation counterparts. In this work, we introduce a multifidelity Monte Carlo method that combines the high-fidelity nonlocal model of interest with surrogate models that use coarser grids and/or smaller horizons and thus have lower fidelities and lower costs. Using the multifidelity method, the overall computational cost of uncertainty quantification is reduced without compromising accuracy. It is shown for a one-dimensional nonlocal diffusion example that speedups of up to two orders of magnitude can be achieved using the multifidelity method to estimate the expectation of an output of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗
DOE OSTI2020

Accelerated version of NUBEAM capabilities in DIII-D using neural networks

A neural network model of the effects of neutral beam injection on DIII-D has been developed. The training and testing data used by the model have been generated by the NUBEAM module of TRANSP for experimental discharges from the 2018 DIII-D campaign. Using a principle component analysis to reduce the dimensionality of profile data, the model has been shown to reproduce the results of the Monte Carlo code NUBEAM with a high level of accuracy and an execution time orders of magnitude faster than the execution time of NUBEAM. Furthermore, this makes the neural network model uniquely suited to applications in model-based scenario planning (off-line) and active control (on-line), where a large number of simulation runs are required by the associated optimization tasks that need to be performed before and during the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Enhanced analysis of experimental x-ray spectra through deep learning

X-ray spectroscopic data from high-energy-density laser-produced plasmas has long required thorough, time-consuming analysis to extract meaningful source conditions. There are often confounding factors due to rapidly evolving states and finite spatial gradients (e.g., the existence of multi-temperature, multi-density, multi-ionization states, etc.) that make spectral measurements and analysis difficult. Here, in this paper, we demonstrate how deep learning can be applied to enhance x-ray spectral data analysis in both speed and intricacy. Neural networks (NNs) are trained on ensemble atomic physics simulations so that they can subsequently construct a model capable of extracting plasma parameters directly from experimental spectra. Through deep learning, the models can extract temperature distributions as opposed to single or dual temperature/density fits from standard trial-and-error atomic modeling at a significantly reduced computational cost compared to traditional trial-and-error methods. These NNs are envisioned to be deployed with high repetition rate x-ray spectrometers in order to provide detailed real-time analysis of experimental spectra.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Finite-size effects in the static structure factor S ( k ) and S ( 0 ) for a two-dimensional Yukawa liquid

Finite-size effects in the static structure factor S(k) are analyzed for an amorphous substance. As the number of particles is reduced, S(0) increases greatly, up to an order of magnitude. Meanwhile, there is a decrease in the height of the first peak S peak. These finite-size effects are modeled accurately by the Binder formula for S(0) and our empirical formula for S peak. Procedures are suggested to correct for finite-size effects in S(k) data and in the hyperuniformity index H ≡ S(0)/S peak. These principles generally apply to S(k) obtained from particle positions in noncrystalline substances. As a result, the amorphous substance we simulate is a two-dimensional liquid, with a soft Yukawa interaction modeling a dusty plasma experiment.

2-dimensional systems↗
NASA NTRS1973

Jupiter's radiation belts: Can Pioneer 10 survive?

Model calculations of Jupiter's electron and proton radiation belts indicate that the Galilean satellites can reduce particle fluxes in certain regions of the inner magnetosphere by as much as six orders of magnitude. Average fluxes should be reduced by a factor of 100 or more along the Pioneer 10 trajectory through the heart of Jupiter's radiation belts in early December. This may be enough to prevent serious radiation damage to the spacecraft.

Hess, W. N.↗
NASA NTRS1984

A Study of the Adequacy of Quasi-geostrophic Dynamics for Modeling the Effect of Frontal Cyclones on the Larger Scale Flow

The major objectives of this study are to test the validity of quasi-geostrophic (QG) dynamics, compared to primitive equation (PE) dynamics for modeling the effect of cyclone waves on the larger scale flow, and to study the formation of frontal cyclones and the dynamics of occluded frontogenesis. In order to extend the realism of the initial conditions over those used in earlier work of this nature and in order to include strong surface frontal zones in the initial conditions, a horizontal resolution on the order of 100 km between gridpoints is needed. In order to allow for integrations of up to two weeks simulated time and to keep the logistics and cost of the runs feasible, the horizontal domain of the models must be reduced. Based on work concerning polar lows, the tropospheric static stability was reduced and was also allowed to vary in different tropospheric layers. Previous work used a 'standard atmosphere' temperature lapse rate. Stratospheric stability remains unchanged. The use of one of these 'reduced stability' cases, with no surface frontal zone present resulted in the discovery of normal mode perturbations.

Mudrick, S.↗
NASA NTRS1987

The effect of decay of the amplitude of oscillation on random process models for QPO X-ray stars

Random process models provide a useful mathematical framework for analysis of the quasi-periodic oscillations (QPO) recently discovered in the X-ray emission from certain luminous Galactic X-ray stars. In this paper, consideration is given to the effects of the power spectrum of the decay of the amplitude of oscillation that is expected in some physical models for the QPO phenomena. The resulting changes in the power spectrum depend on the ratio of the decay time of the amplitude of oscillation to the lifetime of the shot envelope. For values of this ratio less than about 10, it is found that incoherent QPO peaks in the power spectrum are significantly reduced in height and total power and are broadened in width, making them more difficult to detect. Coherent terms in the power spectrum, if present, are reduced in height and total power, but their width is unaffected. The red noise component remains unchanged. The present results are applied to the beat-frequency modulated accretion model in order to set constraints on the physical properties of the boundary layer between the inner edge of the accretion disk and the neutron star magnetosphere.

Shibazaki, N.↗
NASA NTRS1990

Analysis and design of algorithm-based fault-tolerant systems

An important consideration in the design of high performance multiprocessor systems is to ensure the correctness of the results computed in the presence of transient and intermittent failures. Concurrent error detection and correction have been applied to such systems in order to achieve reliability. Algorithm Based Fault Tolerance (ABFT) was suggested as a cost-effective concurrent error detection scheme. The research was motivated by the complexity involved in the analysis and design of ABFT systems. To that end, a matrix-based model was developed and, based on that, algorithms for both the design and analysis of ABFT systems are formulated. These algorithms are less complex than the existing ones. In order to reduce the complexity further, a hierarchical approach is developed for the analysis of large systems.

Nair, V. S. Sukumaran↗
NASA NTRS1995

Analysis of uncertainty in force balance calibration

In order to reduce errors encountered in the measurement and identification of loads using force balances, a third-order forward polynomial relation between loads and output voltages is proposed. This full third-order model represents an alternative to the second-order model currently in use at NASA Langley Research Center (LaRC) and many other installations worldwide. The new model requires the identification of 84 coefficients (including the 28 used presently) for each of the six outputs. The existing LaRC calibration loading sequence is insufficient for identification of many of these 504 coefficients because critical three-load combinations are absent. Accordingly, a new loading sequence that permits the identification of all 504 coefficients has been developed and is described fully. It is apparent from numerical tests that the new third-order model is clearly superior to second-order models in the presence of small amounts of random measurement noise, assuming that there are indeed higher-order interactions between loads. As the amount of noise increases, however, a third-order model becomes less attractive due to its ability to match the noise itself more faithfully than a second-order model. Numerical results suggest that the transition occurs when the magnitude of random noise becomes of the same order as that of the physical higher-order interactions.

Bursal, Faruk H.↗
NASA NTRS2004

Inducer Hydrodynamic Forces in a Cavitating Environment

Marshall Space Flight Center has developed and demonstrated a measurement device for sensing and resolving the hydrodynamic loads on fluid machinery. The device - a derivative of the six-component wind tunnel balance - senses the forces and moments on the rotating device through a weakened shaft section instrumented with a series of strain gauges. This rotating balance was designed to directly measure the steady and unsteady hydrodynamic loads on an inducer, thereby defining the amplitude and frequency content associated with operating in various cavitation modes. The rotating balance was calibrated statically using a dead-weight load system in order to generate the 6 x 12 calibration matrix later used to convert measured voltages to engineering units. Structural modeling suggested that the rotating assembly first bending mode would be significantly reduced with the balance s inclusion. This reduction in structural stiffness was later confirmed experimentally with a hammer-impact test. This effect, coupled with the relatively large damping associated with the rotating balance waterproofing material, limited the device s bandwidth to approximately 50 Hertz Other pre-test validations included sensing the test article rotating assembly built-in imbalance for two configurations and directly measuring the assembly mass and buoyancy while submerged under water. Both tests matched predictions and confirmed the device s sensitivity while stationary and rotating. The rotating balance was then demonstrated in a water test of a full-scale Space Shuttle Main Engine high-pressure liquid oxygen pump inducer. Experimental data was collected a scaled operating conditions at three flow coefficients across a range of cavitation numbers for the single inducer geometry and radial clearance. Two distinct cavitation modes were observed symmetric tip vortex cavitation and alternate-blade cavitation. Although previous experimental tests on the same inducer demonstrated two additional cavitation modes at lower inlet pressures, these conditions proved unreachable with the rotating balance installed due to the intense dynamic environment. The sensed radial load was less influenced by flow coefficient than by cavitation number or cavitation mode although the flow coefficient range was relatively narrow. Transition from symmetric tip vortex to alternate-blade cavitation corresponded to changes in both radial load magnitude and radial load orientation relative to the inducer. Sensed moments indicated that the effective load center moved downstream during this change in cavitation mode. An occurrence of "higher+rdex cavitation" was also detected in both the stationary pressures and the rotating balance data although the frequency of the phenomena was well above the reliable bandwidth of the rotating balance. In summary the experimental tests proved both the concept and device s capability despite the limitations and confirmed that hydrodynamically-induced forces and moments develop in response to the unbalanced pressure field, which is, in turn, a product of the cavitation environment.

Skelley, Stephen E.↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
NASA NTRS2016

Biological Potential in Serpentinizing Systems

Generation of the microbial substrate hydrogen during serpentinization, the aqueous alteration of ultramafic rocks, has focused interest on the potential of serpentinizing systems to support biological communities or even the origin of life. However the process also generates considerable alkalinity, a challenge to life, and both pH and hydrogen concentrations vary widely across natural systems as a result of different host rock and fluid composition and differing physical and hydrogeologic conditions. Biological potential is expected to vary in concert. We examined the impact of such variability on the bioenergetics of an example metabolism, methanogenesis, using a cell-scale reactive transport model to compare rates of metabolic energy generation as a function of physicochemical environment. Potential rates vary over more than 5 orders of magnitude, including bioenergetically non-viable conditions, across the range of naturally occurring conditions. In parallel, we assayed rates of hydrogen metabolism in wells associated with the actively serpentinizing Coast Range Ophiolite, which includes conditions more alkaline and considerably less reducing than is typical of serpentinizing systems. Hydrogen metabolism is observed at pH approaching 12 but, consistent with the model predictions, biological methanogenesis is not observed.

the potential of serpentinizing systems to support↗
NASA NTRS1992

Integrating aerodynamics and structures in the minimum weight design of a supersonic transport wing

An approach is presented for determining the minimum weight design of aircraft wing models which takes into consideration aerodynamics-structure coupling when calculating both zeroth order information needed for analysis and first order information needed for optimization. When performing sensitivity analysis, coupling is accounted for by using a generalized sensitivity formulation. The results presented show that the aeroelastic effects are calculated properly and noticeably reduce constraint approximation errors. However, for the particular example selected, the error introduced by ignoring aeroelastic effects are not sufficient to significantly affect the convergence of the optimization process. Trade studies are reported that consider different structural materials, internal spar layouts, and panel buckling lengths. For the formulation, model and materials used in this study, an advanced aluminum material produced the lightest design while satisfying the problem constraints. Also, shorter panel buckling lengths resulted in lower weights by permitting smaller panel thicknesses and generally, by unloading the wing skins and loading the spar caps. Finally, straight spars required slightly lower wing weights than angled spars.

Barthelemy, Jean-Francois M.↗