Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model order reduction”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Application of CFD to aerothermal heating problems

Numerical solutions of the compressible Navier-Stokes equations by an alternating direction implicit scheme, applied to two experimental investigations are presented. The first is cooling by injection of a gas jet through the nose of an ogive-cone, and the second is the aerothermal environment in the gap formed by the wing and elevon section of a test model of the space shuttle. The simulations demonstrate that accurate pressure calculations are easily obtained on a coarse grid, while convergence is obtained after the residual reduces by four orders of magnitude. Accurate heating rates, however, require a fine grid solution, with convergence requiring at least a reduction of six orders of magnitude in the residual. The effect of artificial dissipation on numerical results is also assessed.

Macaraeg, M. G.↗
NASA NTRS1986

Application of CFD to aerothermal heating problems

Numerical solutions of the compressible Navier-Stokes equations by an Alternating Direction Implicit scheme, applied to two experimental investigations at NASA Langley Research Center are presented. The first is cooling by injection of a gas jet through the nose of an ogive-cone, and the second is the aerothermal environment in the gap formed by the wing and elevon section of a test model of the Space Shuttle. The simulations demonstrate that accurate pressure calculations are easily obtained on a coarse grid, witih convergence being obtained after the residual reduces by four orders of magnitude. However, accurate heating rates require a fine grid solution, with convergence requiring at least a reduction of six orders of magnitude in the residual. The effect of artificial dissipation on numerical results is also assessed.

Macaraeg, M. G.↗
NASA NTRS1986

Hot, diamagnetic cavities upstream from the earth's bow shock

On eight occasions the ISEE 1 and 2 spacecraft registered peculiar plasma structures upstream of the earth's bow shock. The events exhibit a temporary, strong reduction in the magnitude of the magnetic field and strong enhancements of the field strength bordering the reduction zone. The low field strength regions featured temperatures from 1-10 million k and pressure an order of magnitude greater than the solar wind. The pressure gradients exceeded the magnetic tension around the structures, although the field of the cavities may be a closed structure. A model is proposed of hot, expanding diamagnetic plasma cavities with scales on the order of a few earth radii. Speculations on the interaction and origin or impetus for the cavities within the bow shock, foreshock, the magnetosphere and the solar wind are discussed. Similarities between the phenomena detected and signatures obtained with the AMPTE releases of chemicals in the solar wind are noted.

Thomsen, M. F.↗
NASA NTRS1988

Finite-dimensional approximation for optimal fixed-order compensation of distributed parameter systems

In controlling distributed parameter systems it is often desirable to obtain low-order, finite-dimensional controllers in order to minimize real-time computational requirements. Standard approaches to this problem employ model/controller reduction techniques in conjunction with LQG theory. In this paper we consider the finite-dimensional approximation of the infinite-dimensional Bernstein/Hyland optimal projection theory. This approach yields fixed-finite-order controllers which are optimal with respect to high-order, approximating, finite-dimensional plant models. The technique is illustrated by computing a sequence of first-order controllers for one-dimensional, single-input/single-output, parabolic (heat/diffusion) and hereditary systems using spline-based, Ritz-Galerkin, finite element approximation. Numerical studies indicate convergence of the feedback gains with less than 2 percent performance degradation over full-order LQG controllers for the parabolic system and 10 percent degradation for the hereditary system.

Bernstein, Dennis S.↗
NASA NTRS2002

Measurements of Turbulence Attenuation by a Dilute Dispersion of Solid Particles in Homogeneous Isotropic Turbulence

This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Eaton, John↗
DOE OSTI2026

C–C Bond Formation during Electrochemical CO 2 Reduction on Pristine Cu(100) Unlikely to Involve Adsorbed CO at Any Potential

Formation of hydrocarbons containing two or more carbon atoms (C 2+ ) during heterogeneous electrochemical CO and CO 2 reduction (ECOR and ECO 2 R) only occurs, among pure metals, on Cu electrodes. Moreover, the activity and selectivity is facet dependent, with Cu(100) generally preferentially forming ethylene over methane. Previously, we found via quantum-mechanics-based modeling that, unlike standard density functional theory, more accurate correlated wavefunction methods predict that non-electroactive coupling pathways involving two adsorbed COs (*CO) or a *CO and a *COH to form C–C bonds on Cu(100) are kinetically inhibited, with the former also thermodynamically unfavorable. Here, we extend that embedded complete active space second order perturbation theory (ECASPT2) study, further showing that electrochemical coupling of two *COs to form an anionic dimer [OC*–*CO] (1+δ)– , followed by protonation to form [OC*–*COH] δ− , is not kinetically competitive with the reduction of *CO to *COH at relevant ECO/CO 2 R potentials. Our simulations therefore suggest that the ability of Cu(100) to electrochemically synthesize C 2+ molecules from CO and CO 2 is unlikely to be via *CO, at least on pristine Cu(100). Instead, hydrogenated CO species (*COH, *CH x OH, or *CH x ) are most likely to be the key intermediates in C–C bond formation.

Martirez, John Mark P. [Princeton Plasma Physics L↗
DOE OSTI2025

Integration of renewable energy generation and storage systems for emissions reduction in an islanded campus microgrid

Microgrid connected building communities are projected to play an integral part in the clean energy transition. These types of systems, when integrated with distributed energy resources (DERs) such as combined heat and power (CHP), district heating and cooling, renewable generation, and energy storage, can provide clean, reliable power to critical facilities and vulnerable communities. The intermittent nature of renewable generation is a challenge when integrating renewables into any grid system, but particularly in islanded microgrids. The University of Texas at Austin (UT) operates an islanded microgrid powered by a CHP plant, while also utilizing district heating and cooling systems and thermal energy storage (TES). High fidelity operating data was used to develop a validated reduced order model of UT’s integrated campus energy systems to serve as a testbed for use in a case study. Hypothetical renewable energy installations on land owned by UT in west Texas were modeled and integrated into the validated campus models along with battery energy storage (BES). Simulations showed that a combination of renewable energy from wind, and optimally controlled 24-hour thermal and battery storage systems could reduce carbon dioxide emissions on campus by 45.4%. The additional retrofit of burner systems to utilize hydrogen natural gas blends resulted in an overall annual emissions reduction of 54.7%. Carbon capture and storage eliminated the majority of the remaining emissions with increased plant energy expenditure. The presented simulations display the practical limitations of a CHP system complemented by renewable generation and short-term storage in eliminating emissions. Furthermore, results highlight the need for further research and development in long duration storage technologies and hydrogen fueled turbines to increase penetration of renewable energy and reduce emissions.

CHP↗
DOE OSTI2020

A review of black carbon in snow and ice and its impact on the cryosphere

Black carbon (BC) has emerged as an important short-lived climate forcer. Due to its light absorption properties, BC can darken the snow/ice surface, affect the energy balance, and further lead to acceleration of the cryospheric melting. By reviewing the recent published literatures, we present an overview of the historical changes, spatial distribution of BC in snow/ice, and how these changes are related to the melting of the cryosphere (e.g., glaciers, snow cover, and sea ice). Ice core records show an abrupt increase of BC concentrations that began in the 1850s and continued throughout the 20th century, which is consistent with an increase in BC emissions owing to industrialization. A decrease in the amount of BC since 1970s in Arctic and European ice cores has been partially attributed to the Clean Air Act. However, in the Himalayas, BC records show a continuous increase during this period. BC concentrations in snow and ice in the mid-latitude regions are one to two orders of magnitude higher than those in the polar regions. In particular, BC concentrations in aged snow and granular ice in ablation areas of mountain glaciers are one to three orders of magnitude higher than those in fresh snow or snowpits in the glacier accumulation areas due to BC accumulation during melting season. BC in the surface snow/ice is responsible for about 20% of the albedo reduction in the Tibetan Plateau during glacier melt season. Globally, observations and modeling results indicate that radiative forcing (RF) induced by BC in snow and ice is highest in the mid-latitudes, ranging from several W m–2 in fresh snow to hundreds of W m–2 in aged snow and granular ice in glacier ablation areas. The large BC-in-snow RF and associated snow albedo feedback lead to an acceleration in the total glacier melt (approximately 20%) and/or a reduction in the duration of the snow cover by several days, resulting in an increase of glacier discharge. Given our limited understanding of the role of BC in cryospheric melting, it is important to synthesize the existing research on the multi-scale processes related to BC in snow and ice to identify the gaps in our understanding of these processes and to propose a path forward to improve the quality of our observations of the aforementioned phenomena to fill these gaps.

Kang, Shichang↗
DOE OSTI2021

Atomic Origin of the Autocatalytic Reduction of Monoclinic CuO in a Hydrogen Atmosphere

Reducibility is key for the use of bulk metal oxides in chemical transformation involving redox reactions but probing microscopic processes of oxide reduction is challenging. This is because the insulating nature of bulk oxides restricts ion and electron spectroscopic measurements of oxide surfaces. Herein, using a combination of environmental transmission electron microscopy and atomistic modeling, we report direct in-situ atomic-scale observations of the surface and subsurface dynamics and show that the hydrogen induced CuO reduction occurs via the surface decay of Cu-O/Cu bilayer atomic steps, formation of partially reduced CuO superstructure by the self-ordering of O vacancies in the subsurface, and collapse of Cu-O layers in the bulk. In this work, all these substeps can be traced back to the progressively increased concentration and activity of O vacancies in the surface and subsurface of the oxide, thereby leading to the self-accelerated oxide reduction. These results demonstrate the microscopic details that may have broader applicability in modulating various redox processes.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Interactive software for spectral assignment

A new interactive computer software package for eigenvalue/eigenvector assignment using constant state feedback is described. The package consists of ten subprograms, each associated with a specific design objective, accessible from a main control program. Using this package, primary design objectives of assigning eigenvalues and approximating eigenvectors are first achieved. Then secondary design objectives, including modification of specific eigenvector components, reduction in specified elements of the feedback gain matrix, and reduction in eigensystem sensitivity to changes in plant parameters, are addressed. These secondary objectives are achieved by a systematic modification of the assigned eigenvectors in a small region about the initial assignment. In addition, the program implements the use of spectral assignment procedures with reduced-order system models. Program modes are described and illustrated by numerical examples.

Mielke, R. R.↗
NASA NTRS1989

Reliability-based failure analysis of brittle materials

The reliability of brittle materials under a generalized state of stress is analyzed using the Batdorf model. The model is modified to include the reduction in shear due to the effect of the compressive stress on the microscopic crack faces. The combined effect of both surface and volume flaws is included. Due to the nature of fracture of brittle materials under compressive loading, the component is modeled as a series system in order to establish bounds on the probability of failure. A computer program was written to determine the probability of failure employing data from a finite element analysis. The analysis showed that for tensile loading a single crack will be the cause of total failure but under compressive loading a series of microscopic cracks must join together to form a dominant crack.

Powers, Lynn M.↗
NASA NTRS2020

Computational Study of NASA's Quadrotor Urban Air Taxi Concept

High-fidelity computational fluid dynamics simulations have been carried out in order to analyze NASA's quadrotor urban air taxi concept for urban air mobility, also know as on-demand mobility applications. High-order accurate schemes, dual-time stepping, and the delayed detached-eddy simulation model have been employed. The ow solver has been loosely coupled with a rotorcraft comprehensive analysis code. The vehicle simulated is a six-passenger quadrotor for air taxi operations. A study of power reduction as a function of the rear-rotor to front-rotors vertical separation has been performed, for a quad-rotor without the airframe, in cruise flight conditions. Then, the quadrotor without the airframe has been simulated in hover. The airloads and wake geometries are analyzed. To finish the study the complete quadrotor vehicle is presented. NASA's quadrotor air taxi concept is one of the many concepts being developed by NASA in support of aircraft development for vertical take-o and landing air taxi operations.

Ventura Diaz, Patricia↗
NASA NTRS2020

Nonbinary Systems: Looking Towards the Future of Gender Equity in Planetary Science

Gender equity remains a major issue facing the field of planetary science, and there is broad interest in addressing gender disparities within space science and related disciplines. Many studies of these topics have been performed by professional planetary scientists who are relatively unfamiliar with research in fields such as gender studies and sociology. As a result, they adopt a normative view of gender as a binary choice of ‘male’ or ‘female,’ leaving planetary scientists whose genders do not fit within that model out of such research entirely. Reductive frameworks of gender and an overemphasis on quantification as an indicator of gendered phenomena are harmful to people of marginalized genders, especially those who live at the intersections of multiple axes of marginalization such as race, disability, and socioeconomic status. In order for the planetary science community to best serve its marginalized members as we move into the next decade, a new paradigm must be established. This paper aims to address the future of gender equity in planetary science by recommending better survey practices and institutional policies based on a more profound approach to gender.

Beck E Strauss↗
DOE OSTI2022

Calibration of cloud and aerosol related parameters for solar irradiance forecasts in WRF-solar

Model parameters are a major source of uncertainty in numerical weather prediction. Recently, the Weather Research and Forecasting model with Solar extensions (WRF-Solar) has been upgraded by enhancing the treatment of sub-grid scale cloud and aerosols with augmentations of a sub-grid scale cloud scheme (CLD3) and an upgraded aerosol-aware Thompson-Eidhammer scheme (TE14). However, the value of model parameters associated with these parameterizations are assigned based on limited measurements or theoretical calculations. Calibrating the most sensitive parameters has the potential to improve solar irradiance predictions. Here, we adopted a multiobjective surrogate-based optimization (SBO) framework to calibrate nine parameters used in CLD3 and TE14 that lead to the largest sensitivity in simulated irradiance. The normalized mean-absolute-error (NMAE) of global horizontal irradiance (GHI) and direct normal irradiance (DNI) are minimized by calibrating WRF-Solar over two regions including the Southern Great Plains (SGP) and Central California, in order to focus on parameter calibration under cloudy conditions with different aerosol loading. The results show that generalized linear model (GLM)-based surrogate models approximate physical models well, particularly when the third order and three-way interaction terms are considered. The SBO framework efficiently searches the parameter space for optimal solutions with less computational costs than directly calibrating the physical model. We first calibrate CLD3 parameters over the less-polluted SGP region. Optimized CLD3 parameters alone result in NMAE reduction by 14% for the site-mean and up to 33% for individual cases over the SGP region. With further calibration of TE14 parameters over the Central California during active fire periods, the optimized parameters lead to over 20% reductions of NMAE. Our investigation reveals, however, that optimizing TE14 has a limited impact on irradiance simulations under less-polluted conditions in the SGP.

14 SOLAR ENERGY↗
DOE OSTI2026

Field-tailoring quantum materials via magneto-synthesis: metastable metallic and magnetically suppressed phases in a trimer iridate

We demonstrate that applying modest magnetic fields (< 0.1 T) during high-temperature crystal growth can profoundly alter the structure and ground state of a spin-orbit-coupled, antiferromagnetic trimer lattice. Using BaIrO₃ as a model system, whose ground state is intricately dictated by the trimer lattice, we show that magneto-synthesis , a field-assisted synthesis approach, stabilizes a structurally compressed, metastable metallic and magnetically suppressed phases inaccessible via conventional methods. These effects include a 0.85% reduction in unit cell, 4-order-of-magnitude decrease in resistivity, a 10-fold enhancement of the Sommerfeld coefficient, and the collapse of long-range magnetic order -- all intrinsic and bulk in origin. First-principles calculations confirm that the field-stabilized structure lies substantially above the ground state in energy, highlighting its metastable character. These large, coherent and correlated changes across multiple bulk properties, unlike those caused by dilute impurities, defects or off-stoichiometry, point to an intrinsic field-induced mechanism. The findings establish magneto-synthesis as a powerful new pathway for accessing non-equilibrium quantum phases in strongly correlated materials.

magneto-synthesis↗
DOE OSTI2024

Sparsified time-dependent Fourier neural operators for fusion simulations

This paper presents a sparsified Fourier neural operator for coupled time-dependent partial differential equations (ST-FNO) as an efficient machine learning surrogate for fluid and particle-based fusion codes such as NIMROD (Non-Ideal Magnetohydrodynamics with Rotation - Open Discussion) and GTC (Gyrokinetic Toroidal Code). ST-FNO leverages the structures in the governing equations and utilizes neural operators to represent Green's function-like numerical operators in the corresponding numerical solvers. Once trained, ST-FNO can rapidly and accurately predict dynamics in fusion devices compared with first-principle numerical algorithms. In general, ST-FNO represents an efficient and accurate machine learning surrogate for numerical simulators for multi-variable nonlinear time-dependent partial differential equations, with the proposed architectures and loss functions. The efficacy of ST-FNO has been demonstrated using quiescent H-mode simulation data from NIMROD and kink-mode simulation data from GTC. The ST-FNO H-mode results show orders of magnitude reduction in memory and central processing unit usage in comparison with the numerical solvers in NIMROD when computing fields over a selected poloidal plane. The ST-FNO kink-mode results achieve a factor of 2 reduction in the number of parameters compared to baseline FNO models without accuracy loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Training custom light curve models of SN Ia subpopulations selected according to host galaxy properties

ABSTRACT Type Ia supernova (SN Ia) cosmology analyses include a luminosity step function in their distance standardization process to account for an observed yet unexplained difference in the post-standardization luminosities of SNe Ia originating from different host galaxy populations [e.g. high-mass ($M \gtrsim 10^{10} \, {\rm M}_{\odot }$) versus low-mass galaxies]. We present a novel method for including host-mass correlations in the SALT3 (Spectral Adaptive Light curve Template 3) light curve model used for standardizing SN Ia distances. We split the SALT3 training sample according to host-mass, training independent models for the low- and high-host-mass samples. Our models indicate that there are different average Si ii spectral feature strengths between the two populations, and that the average spectral energy distribution of SNe from low-mass galaxies is bluer than the high-mass counterpart. We then use our trained models to perform an SN cosmology analysis on the 3-yr spectroscopically confirmed Dark Energy Survey SN sample, treating SNe from low- and high-mass host galaxies as separate populations throughout. We find that our mass-split models reduce the Hubble residual scatter in the sample, albeit at a low statistical significance. We do find a reduction in the mass-correlated luminosity step but conclude that this arises from the model-dependent re-definition of the fiducial SN absolute magnitude rather than the models themselves. Our results stress the importance of adopting a standard definition of the SN parameters (x0, x1, c) in order to extract the most value out of the light curve modelling tools that are currently available and to correctly interpret results that are fit with different models.

Taylor, G. (ORCID:0000000157563259)↗
NASA NTRS1980

Geophysical and astronomical models applied in the analysis of very long baseline interferometry

Very long baseline interferometry presents an opportunity to measure at the centimeter level such geodetic parameters as baseline length and instantaneous pole position. In order to achieve such precision, the geophysical and astronomical models used in data analysis must be as accurate as possible. The Mark-3 interactive data analysis system includes a number of refinements beyond conventional practice in modeling precession, nutation, diurnal polar motion, UT1, solid Earth tides, relativistic light deflection, and reduction to solar system barycentric coordinates. The algorithms and their effects on the recovered geodetic, geophysical, and astrometric parameters are discussed.

Ma, C.↗