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 829 records · Page 46

Dynamics of Permanent-Magnet Biased Active Magnetic Bearings

Active magnetic radial bearings are constructed with a combination of permanent magnets to provide bias forces and electromagnets to generate control forces for the reduction of cost and the operating energy consumption. Ring-shaped permanent magnets with axial magnetization are attached to a shaft and share their magnet stators with the electromagnets. The magnet cores are made of solid iron for simplicity. A simplified magnetic circuit of the combined magnet system is analyzed with linear circuit theory by approximating the characteristics of permanent magnets with a linear relation. A linearized dynamical model of the control force is presented with the first-order approximation of the effects of eddy currents. Frequency responses of the rotor motion to disturbance inputs and the motion for impulsive forces are tested in the non-rotating state. The frequency responses are compared with numerical results. The decay of rotor speed due to magnetic braking is examined. The experimental results and the presented linearized model are similar to those of the all-electromagnetic design.

Fukata, Satoru↗

Development, Analysis and Testing of the High Speed Research Flexible Semispan Model

This report presents the work performed by Lockheed Martin Engineering and Sciences (LMES) in support of the High Speed Research (HSR) Flexible Semispan Model (FSM) wind-tunnel test. The test was conducted in order to assess the aerodynamic and aeroelastic character of a flexible high speed civil transport wing. Data was acquired for the purpose of code validation and trend evaluation for this type of wing. The report describes a number of activities in preparing for and conducting the wind-tunnel test. These included coordination of the design and fabrication, development of analytical models, analysis/hardware correlation, performance of laboratory tests, monitoring of model safety issues, and wind-tunnel data acquisition and reduction. Descriptions and relevant evaluations associated with the pretest data are given in sections 1 through 6, followed by pre- and post-test flutter analysis in section 7, and the results of the aerodynamics/loads test in section 8. Finally, section 9 provides some recommendations based on lessons learned throughout the FSM program.

Schuster, David M.↗

Electric Field Measurements to Investigate Sheath Formation in a Nanosecond Pulsed Discharge

The objective of the research was to study the temporally and spatially resolved dynamics of sheath formation in a nanosecond pulsed plasma generated in a parallel plate electrode geometry at near atmospheric pressures. We have performed LIF dip measurements at the Sandia Low-Temperature Plasma Research Facility showing the capability to measure electric fields in near atmospheric pressure plasmas although with insufficient spatial resolution near the electrodes. Stark polarization spectroscopy implemented at the University of Minnesota was shown to be able to probe the formation of the sheaths in nanosecond pulsed plasmas. Within 30 ns, a reduction of the sheath thickness to about 250 μm is observed, coinciding with a gradual increase of the discharge current and proportional increase in electric field at the cathode. The electric field evolution as obtained by the fluid model is in excellent agreement with the measurements and shows that an enhanced ionization near the cathode is causing the space charge formation responsible for the increase in electric field. High resolution Stark polarization spectroscopy showed sheaths with dimensions of the order of 50 μm underlining the need of micrometer spatial resolution for sheath characterization in near atmospheric pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A platform to measure isentropes from proton-heated warm dense matter on short pulse laser facilities

We describe the development of an experimental platform that measures the release isentrope of materials heated isochorically to temperatures of a few electron volts, using short-pulse laser-produced protons to heat the sample and long-pulse laser-produced x rays to perform streaked x-ray radiography. The density profiles derived from the radiography data are integrated to generate pressure–density isentropes, independent of prior knowledge of the equation of state of the sample material. In order to understand the sensitivities of isentrope extraction from radiography data, we analyze synthetic radiographs generated by a radiation hydrodynamics code. Noise reduction and high spatial resolution are critical for isentrope reconstruction, as demonstrated by the analysis of a proof-of-principle shot day on the OMEGA-EP facility. In conclusion, the data demonstrate the feasibility of the platform for characterizing isentropes, and we discuss the necessary improvements to enhance precision in differentiating between equation-of-state models.

Equations of state↗

Role of the edge stochastic layer in density pump-out by resonant magnetic perturbations

The role of the edge stochastic layer on particle transport is addressed in DIII-D plasmas with applied resonant magnetic perturbations (RMPs) causing density pump-out (a term used to denote the density reduction at the top of the pedestal due to the applied RMPs). Using an analytical model that adds the stochastic parallel transport of electrons (Harvey et al 1981 Phys. Rev. Lett. 47 102) to the fluid equations, the ambipolar radial electric field and particle flux are calculated simultaneously. In this model the nonambipolar electron particle flux, driven by the stochastic magnetic field, is predominantly balanced by the nonambipolar perpendicular ion flux, driven by toroidal viscosity, across a narrow stochastic layer of order two percent of the plasma radius (0.98 < ψn < 1). The model reproduces the level of density pump-out and its dependence on RMP amplitude near the separatrix for three plasma discharges for which the pedestal foot is at medium to high collisionality $(2\leqslant {\nu }_{\text{sep}}^{\ast }\leqslant 30)$. In this work, the experimentally observed and numerically tested inverse density dependence of density pump-out (Hu et al 2020 Nucl. Fusion 60 076001) is accurately captured by the model: an increase in collisionality with density in the pedestal foot results in a decrease in stochastic diffusivity, and hence a decrease in the level of pump-out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduction of intensity variations on a photovoltaic array with compound parabolic concentrators

The reduction of nonuniformities in the intensity distribution of light focused on a photovoltaic array by a compound parabolic concentrator is investigated. The introduction of small distortions into the surfaces of the reflector in order to diffuse the incident collimated light to fill the angular acceptance of the concentrator is calculated by means of ray tracing to decrease the irradiance nonuniformity at the cost of a lowered effective concentration of the concentrator. Measurements of the intensity distribution on a scale test model in terrestrial sunlight with corrugated aluminized mylar reflectors are shown to be in good agreement with the ray tracing results. A two-stage concentrator consisting of a focusing primary and a nonimaging secondary is also shown to result in a fairly uniform intensity distribution except in the case of a 4-deg incidence angle, which may be corrected by the introduction of distortions into one or both concentration stages.

Greenman, P.↗

Exposure time considerations in high temperature low cycle fatigue

The Conventional Strainrange Partitioning (CSRP) method for High-Temperature, Low-Cycle Fatigue (HTLCF) life prediction has its origins in the modeling of first-order, creep-fatigue waveform effects while treating as second-order effects, the influence of metallurgical or environmental time dependencies. Procedures are proposed to include the latter explicitly in the inelastic strainrange--life relations. For brevity, only the CP life relation will be presented in detail. The exposure-time effect within the CP inelastic strainrange (tensile creep reversed by compressive plasticity) was determined by tensile stresshold-time experiments for 316 SS at 816 C. Reductions in CP cyclic life of a factor of about two were observed with an increase in exposure time or a corresponding decrease in creep rate by a factor of about 100. The CP life relation has been modified to be expressed in terms of either Steady State Creep Rate (SSCR) or Exposure Time (ET). The applicability and accuracy of the time-dependent CP life relations is demonstrated by conducting verification experiments involving complex hysteresis loops. Metallographic examination revealed time-dependent degradation attributable to oxide formation and precipitation of carbides along grain boundaries.

Kalluri, S.↗

A reduced version of the NMC DERF 2 data set

The National Meterological Center (NMC) Dynamical Extended Range Forecast (DERF 2) data represents a major computational effort to better ascertain the potential for extended range forecasts and to develop a strategy for performing operational extended range forecasts using dynamical models. A major stumbling block for using this data has been the sheer volume of data that must be processed to perform even simple calculations. The product of the data reduction described is a manageable data set that fits comfortably on five magnetic tapes or on one compact disc. The document outlines the data reduction process of the second phase of DERF data. It contains the description of the fields and the resolution of both the original and final fields. In order to assist the users of this data set, maps of selected fields, using both the original truncation at rhomboidal 30 and the truncation of the final data at triangular 20, are displayed.

Schubert, Siegfried↗

Integrated active and passive control design methodology for the LaRC CSI evolutionary model

A general design methodology to integrate active control with passive damping was demonstrated on the NASA LaRC CSI Evolutionary Model (CEM), a ground testbed for future large, flexible spacecraft. Vibration suppression controllers designed for Line-of Sight (LOS) minimization were successfully implemented on the CEM. A frequency-shaped H2 methodology was developed, allowing the designer to specify the roll-off of the MIMO compensator. A closed loop bandwidth of 4 Hz, including the six rigid body modes and the first three dominant elastic modes of the CEM was achieved. Good agreement was demonstrated between experimental data and analytical predictions for the closed loop frequency response and random tests. Using the Modal Strain Energy (MSE) method, a passive damping treatment consisting of 60 viscoelastically damped struts was designed, fabricated and implemented on the CEM. Damping levels for the targeted modes were more than an order of magnitude larger than for the undamped structure. Using measured loss and stiffness data for the individual damped struts, analytical predictions of the damping levels were very close to the experimental values in the (1-10) Hz frequency range where the open loop model matched the experimental data. An integrated active/passive controller was successfully implemented on the CEM and was evaluated against an active-only controller. A two-fold increase in the effective control bandwidth and further reductions of 30 percent to 50 percent in the LOS RMS outputs were achieved compared to an active-only controller. Superior performance was also obtained compared to a High-Authority/Low-Authority (HAC/LAC) controller.

Voth, Christopher T.↗

Lightning Strike Distance Distribution Beyond a Preexisting Lightning Area

The 45th Weather Squadron (45 WS) asked the Applied Meteorology Unit (AMU) to review the 30-year-old, lightning stand-off distances of 5 nautical miles (nmi) for applicability to today's operations. This was based on the realization that previous lightning strike distance studies did not match how 45 WS issues lightning warnings (Roeder, 2008). The previous lightning distance studies were from the point of origin of the lightning or the average starting location that would tend to be in the core of the thunderstorm. However, the 45 WS issues lightning warnings based on the edge of a preexisting lightning area. Before beginning the AMU project, it took several years to develop a method to calculate a distance distribution beyond a preexisting area (Roeder, 2015). The AMU pulled Lightning Detection and Ranging (LDAR) sensor data from 1/1/2013 to 12/31/2013. This dataset consisted of 37 million individual source data points from the LDAR sensors. Only sources within 50 km north, south, east or west of the LDAR grid center were included in the dataset. This limited the use of LDAR data to that with the greatest accuracy of source detection and increased data processing. Points were grouped into flashes based on spatial and temporal criteria. Based on the sensitivity analysis the AMU performed on the flash clustering algorithm, a time value of 0.3 seconds was found to model flashes adequately. Distance parameters were tested from 1,500 to 7,500 meters (m) in 500 m increments. Distance parameters of both 3,000 m and 4,000 m produced results in the plotting tool that were most representative of the physical behavior of lightning. Thus statistics were gathered for the most representative of these spatial and temporal criteria on the flash size and the polygon expansion distance in order to find the correct data distributions. The best fit curves for the LDAR polygon expansion frequency vs. distances for both the 3 kilometer (km) and 4 km distance threshold values were exponential decay functions and had R2 values of > 0.998, indicating good model fits. The equations of the best fit curves were then used to calculate a desired safety radius of 4 nmi for either 3 km or 4 km distance threshold criteria. The AMU analysis concludes the safe reduction of the 5 nmi lightning warning circles to 4 nmi should improve the operational impact by 36% if based on distance from the center of the property area being protected. If based on the edge of the property being protected, then the reduction is 4.5 nmi to 4 nmi and the operational impact is decreased by 21%. For the 6 nmi lightning warning circles, the recommended 4 nmi stand-off distance will result in a safe reduction of operational impact of 31% if based on the center of the area being protected, or 16% if based on the edge of the property.

Flash clustering algorithm↗

Implementation of Stochastic Gradient Descent in an Automated Glow Peak Identification Software for Multiple Thermoluminescent Dosimeter Types

A glow-curve analysis code was previously developed in C++ to analyze thermoluminescent dosimeter glow curves using automated peak detection while applying a first-order kinetics model. A newer version of this code was implemented to improve the automated peak detection and curve fitting models. The Stochastic Gradient Descent Algorithm was introduced to replace the prior approach of taking first and second-order derivatives for peak detection. Additionally, early stopping mechanisms were invoked to improve the previously used Levenberg-Marquardt Algorithm employed for curve fitting. The two software versions were compared through glow curve analysis of different thermoluminescent dosimeter materials and calculation of the corresponding figures of merit. Altogether improvements were shown, namely an increase in the number of peaks detected and a reduction of the mean figure of merit by approximately 46%.

137Cs↗

Application of Modified Meshgraphnets for Subsurface Prediction during CO2 Sequestration

In the face of the increasingly dire consequences of anthropogenic climate change, capturing and storing carbon dioxide is paramount. However, several impediments exist to the safe and effective subsurface storage of CO2, such as cost of transport, identification of suitable sites for subsurface storage, and assessment of long-term risk from storage in subsurface aquifers. Accurate subsurface modeling is necessary to ensure that CO2 storage is both safe and effective. Still, such modeling has traditionally required either substantial time and computational power (numerical simulation) or a substantial amount of pre-existing data for training (machine learning models). Additionally, these models lack flexibility in dealing with both changes in discretization of the input data and generalizability beyond the data on which they are trained. In order to address these issues, this research applies graph neural networks (GNNs) to predict subsurface saturation and pressure during CO₂ injection in a model of the Illinois Basin-Decatur Project (IBDP). GNNs provide a flexible, intuitive method for representing and manipulating complex unstructured data, which is often found in many practical domain problems such as fluid flow and subsurface characterization. These unstructured grids are easily represented in GNNs by representing spatially-localized features such as permeability, porosity, saturation, and pressure as nodes in a graph and relationships between these properties as edges connecting these nodes. This research applies a specific GNN model called MeshGraphNets (MGN) to model the change in CO2 saturation and pressure over a 50-month time period (36 months of injection, 14 months post-injection). The MGN model leverages a message passing process that allows the network to learn both the spatial and temporal dynamics of this system simultaneously. Additionally, training on a limited dataset (64 realizations, 20 time points each) resulted in a high degree of accuracy in saturation prediction both within the same timeframe as the training (20 months, 0.039 average RMSE) and when projecting out to the end of injection (36 months, 0.053 average RMSE). Temporal predictions such as those generated by MGNs and other similar models are prone to accumulated error over time; in order to address this, a multi-step rollout (MSR) training process was applied to calculate training loss. This method mimics the forward prediction during inference by “rolling out” multiple time points in a single training step using the previous prediction as input to the MGN model. By calculating the loss several time steps forward from the current prediction, the model is forced to find a more stable state over time. Application of MSR to the MGN model resulted in an average 15% reduction in inference error over time during forward prediction. This study showcases the immense potential of GNNs as a game-changing methodology for predicting pressure and saturation evolution in CCS projects, ultimately paving the way for more sustainable and effective carbon storage solutions. Presentation prepared for the 2024 AiChE Annual Meeting, October 27 to November 1 2024, San Diego, CA.

Holcomb, Paul↗

Radar backscattering from a sea having an anisotropic large-scale surface, part 2

A two scale scattering model was derived that combines specular reflections from sea waves and Bragg scattering in a manner consistent with energy conservation. The effect of the tilting of the small scale roughness by the large scale roughness was included, which accounted for the reduction of reflected power. The special case of backscattering for which the transmitted polarization equaled the received polarization was considered. An anisotropic large scale surface was used to specify the probability density function of the large scale surface normal. In order to isolate the azimuthal variation of the normalized radar cross section produced by the anisotropic probability density function, an isotropical small scale spectrum was assumed.

Wentz, F. J.↗

Exposure time considerations in high temperature low cycle fatigue

The Conventional Strainrange Partitioning (CSRP) method for High-Temperature, Low-Cycle Fatigue (HTLCF) life prediction has its origins in the modeling of first-order, creep-fatigue waveform effects while treating as second-order effects, the influence of metallurgical or environmental time dependencies. Procedures are proposed to include the latter explicitly in the inelastic strainrange-life relations. For brevity, only the CP life relation will be presented in detail. The exposure-time effect within the CP inelastic strainrange (tensile creep reversed by compressive plasticity) was determined by tensile stresshold-time experiments for 316 SS at 816 C. Reductions in CP cyclic life of a factor of about two were observed with an increase in exposure time or a corresponding decrease in creep rate by a factor of about 100. The CP life relation has been modified to be expressed in terms of either Steady State Creep Rate (SSCR) or Exposure Time (ET). The applicability and accuracy of the time-dependent CP life relations is demonstrated by conducting verification experiments involving complex hysteresis loops. Metallographic examination revealed time-dependent degradation attributable to oxide formation and precipitation of carbides along grain boundaries.

Kalluri, S.↗

Internal Mixing Studied for GE/ARL Ejector Nozzle

To achieve jet noise reduction goals for the High Speed Civil Transport aircraft, researchers have been investigating the mixer-ejector nozzle concept. For this concept, a primary nozzle with multiple chutes is surrounded by an ejector. The ejector mixes low-momentum ambient air with the hot engine exhaust to reduce the jet velocity and, hence, the jet noise. It is desirable to mix the two streams as fast as possible in order to minimize the length and weight of the ejector. An earlier model of the mixer-ejector nozzle was tested extensively in the Aerodynamic Research Laboratory (ARL) of GE Aircraft Engines at Cincinnati, Ohio. While testing was continuing with later generations of the nozzle, the earlier model was brought to the NASA Lewis Research Center for relatively fundamental measurements. Goals of the Lewis study were to obtain details of the flow field to aid computational fluid dynamics (CFD) efforts and obtain a better understanding of the flow mechanisms, as well as to experiment with mixing enhancement devices, such as tabs. The measurements were made in an open jet facility for cold (unheated) flow without a surrounding coflowing stream.

Zaman, Khairul↗

Assessment of the driving factors of CO 2 mitigation costs of household biogas systems in China: A LMDI decomposition with cost analysis model

China has been placing a substantial focus on biogas for reducing energy consumption and carbon dioxide (CO 2 ) emissions. The operation mode of biogas systems may make the CO 2 reduction target over-optimistic. There is limited research to investigate the influential factors that may be causing the gap between the actual and theoretical CO 2 reduction costs of biogas systems in China. In this research, by using field survey data of 209 biogas users and 489 non-biogas users from 19 villages in 2015, the gap between actual and theoretical unit CO 2 reduction cost is quantified at approximately 156 USD/t CO 2 . By employing the Logarithmic Mean Divisia Index I (LMDI) model, it is found that both the cost effect (48%) and the reduction effect (52%) contribute to the unit CO 2 reduction cost gap. Furthermore, four influential factors–household labor, accessibility to the energy resource, acceptance of biogas technology, and subsidy–significantly narrow the gap between actual and theoretical CO 2 reduction costs, while the levelized subsidy contributes to widening the gap. On average, biogas systems should be operated for at least four years and the substitution rate should be more than 67% in order to keep the gap between actual and theoretical CO 2 reduction costs under 50%.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

An analytic modeling and system identification study of rotor/fuselage dynamics at hover

A combination of analytic modeling and system identification methods have been used to develop an improved dynamic model describing the response of articulated rotor helicopters to control inputs. A high-order linearized model of coupled rotor/body dynamics including flap and lag degrees of freedom and inflow dynamics with literal coefficients is compared to flight test data from single rotor helicopters in the near hover trim condition. The identification problem was formulated using the maximum likelihood function in the time domain. The dynamic model with literal coefficients was used to generate the model states, and the model was parametrized in terms of physical constants of the aircraft rather than the stability derivatives resulting in a significant reduction in the number of quantities to be identified. The likelihood function was optimized using the genetic algorithm approach. This method proved highly effective in producing an estimated model from flight test data which included coupled fuselage/rotor dynamics. Using this approach it has been shown that blade flexibility is a significant contributing factor to the discrepancies between theory and experiment shown in previous studies. Addition of flexible modes, properly incorporating the constraint due to the lag dampers, results in excellent agreement between flight test and theory, especially in the high frequency range.

Hong, Steven W.↗