Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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,279 records · Page 71

Dynamic Modeling of Starting Aerodynamics and Stage Matching in an Axi-Centrifugal Compressor

A DYNamic Turbine Engine Compressor Code (DYNTECC) has been modified to model speed transients from 0-100% of compressor design speed. The impetus for this enhancement was to investigate stage matching and stalling behavior during a start sequence as compared to rotating stall events above ground idle. The model can simulate speed and throttle excursions simultaneously as well as time varying bleed flow schedules. Results of a start simulation are presented and compared to experimental data obtained from an axi-centrifugal turboshaft engine and companion compressor rig. Stage by stage comparisons reveal the front stages to be operating in or near rotating stall through most of the start sequence. The model matches the starting operating line quite well in the forward stages with deviations appearing in the rearward stages near the start bleed. Overall, the performance of the model is very promising and adds significantly to the dynamic simulation capabilities of DYNTECC.

Wilkes, Kevin↗

Calculation of a circular jet in crossflow with a multiple-time-scale turbulence model

Numerical calculation of a three dimensional turbulent flow of a jet in a crossflow using a multiple time scale turbulence model is presented. The turbulence in the forward region of the jet is in a stronger inequilibrium state than that in the wake region of the jet, while the turbulence level in the wake region is higher than that in the front region. The calculated flow and the concentration fields are in very good agreement with the measured data, and it indicated that the turbulent transport of mass, concentration, and momentum is strongly governed by the inequilibrium turbulence. The capability of the multiple time scale turbulence model to resolve the inequilibrium turbulence field is also discussed.

Kim, S.-W.↗

Calculation of a circular jet in crossflow with a multiple-time-scale turbulence model

Numerical calculation of a three dimensional turbulent flow of a jet in a crossflow using a multiple time scale turbulence model is presented. The turbulence in the forward region of the jet is in a stronger inequilibrium state than that in the wake region of the jet, while the turbulence level in the wake region is higher than that in the front region. The calculated flow and the concentration fields are in very good agreement with the measured data, and it indicated that the turbulent transport of mass, concentration, and momentum is strongly governed by the inequilibrium turbulence. The capability of the multiple time scale turbulence model to resolve the inequilibrium turbulence field is also discussed.

Kim, S.-W.↗

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↗

Numerical modeling of laser thermal propulsion flows

An review of the problems associated with modeling laser thermal propulsion flows, a synopsis of the status of such models, and the attributes of a successful model are presented. The continuous gaseous hydrogen laser-supported combustion wave (LSCW) thruster, for which a high-energy laser system (preferably space-based) should exist by the time the propulsion technology is developed, is considered in particular. The model proposed by Raizer (1970) is based on the assumptions of one-dimensional flow at constant pressure with heat conduction as the principal heat transfer mechanism. Consideration is given to subsequent models which account for radiative transfer into the ambient gas; provide a two-dimensional generalization of Raizer's analysis for the subsonic propagation of laser sparks in air; include the effect of forward plasma radiation in a one-dimensional model; and attempt a time-dependent (elliptic) solution of the full Navier-Stokes equations for the flow in a simple axisymmetric thruster. Attention is also given to thruster and nozzle flow models and thermodynamic and transport properties.

Mccay, T. D.↗

Design of a simple active controller to suppress helicopter air resonance

A coupled rotor/fuselage helicopter analysis with the important effects of blade torsional flexibility, unsteady aerodynamics, and forward flight is presented. Using this mathematical model, a nominal configuration is selected that experiences an air resonance instability throughout most of its flight envelope. A simple multivariable compensator using conventional swashplate inputs and a single body roll rate measurement is then designed. The controller design is based on a linear estimator in conjunction with optimal feedback gains, and the design is done in the frequency domain using the Loop Transfer Recovery method. The controller is shown to suppress the air resonance instability throughout wide range helicopter loading conditions and forward flight speeds.

Takahashi, M. D.↗

Quadratic synthesis of integrated active controls for an aeroelastic forward-swept-wing-aircraft

State variable representations of flexible aircraft are obtained using a modeling approach which first defines a mean-reference axis coordinate system. The active system design approach for forward-swept-wing aircraft also expresses structural deformations of the wing in terms of free normal vibration modes relative to this mean reference axis. Calculations which assume quasi-steady incompressible aerodynamics are performed to obtain generalized force expressions for flight conditions. The use of this quadratic synthesis technique shows that increased performance and redundancy over decentralized approaches can be achieved using integrated active longitudinal stability augmentation and aeroelastic stabilization.

Gilbert, M. G.↗

Centrality determination with a forward detector in the RHIC Beam Energy Scan

Recently, Chatterjee et al [1] used a hadronic transport model to estimate the resolution with which various experimental quantities select the impact parameter of relativistic heavy ion collisions at collision energies relevant to the Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider (RHIC). Measures based on particle multiplicity at forward rapidity were found to be significantly worse than those based on midrapidity multiplicity. Using the same model, we show that a slightly more sophisticated measure greatly improves the resolution based on forward rapidity particles; this improvement persists even when the model is filtered through a realistic simulation of a recent upgrade detector to the STAR experiment. Furthermore, these results highlight the importance of optimizing centrality measures based on particles detected at forward rapidity, especially for experimental studies that search for a critical point in the QCD phase diagram. Such measurements usually focus on proton multiplicity fluctuations at midrapidity, hence selecting events based on multiplicity at midrapidity raises the possibility of nontrivial autocorrelations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Probabilistic Groundwater Modeling of the RDX Plume at Los Alamos National Laboratory to Support Risk Assessment - 20359

A plume of the contaminant hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) with concentrations greater than the New Mexico tap water drinking standard (9.66 ppb) is present in the regional aquifer near the southwestern boundary of Los Alamos National Laboratory (LANL). A risk assessment is performed for exposure to regional aquifer groundwater, with long-term predictions of RDX concentrations provided by a calibrated, probabilistic, numerical fate and transport model. The structure of the model is hierarchical, with the RDX Regional Aquifer (RA) groundwater model acting as the primary tool for analysis of downgradient RDX concentrations. The RA model is deeply informed by the conceptual site model (CSM) and is calibrated using site RDX concentration data and hydraulic head measurements, along with other analyses. The model is calibrated using data through December 2019. Where data are scarce other lines of evidence are used to inform inputs, including the multiphase RDX Vadose Zone (VZ) model and Pipe and Disk analytical screening tool. Model inputs are described with informative prior distributions using a robust approach to development that incorporates all available lines of evidence for every parameter used as an input to the RA model. Calibration is performed using a classical nonlinear optimization routine, which is then used to initialize a Bayesian calibration. The Bayesian calibration constrains the uncertainty in the classical calibration, ultimately providing posterior distributions for all model parameters. The challenges of the calibration include high-dimensional parameter space, including spatially heterogeneous hydraulic conductivities, comparatively sparse data, and low RDX concentrations. Posterior parameter distributions developed in the calibration process are then used for stochastic predictive model runs into the future. The result of the forward runs is spatially and temporally explicit estimates of head and concentration with uncertainty at all points in space and time. The probabilistic modeling approach presented here includes innovative computational and statistical methods that leverage high-performance computing (HPC) resources. It makes use of prior modeling work performed at the chromium plume site in the central LANL area, with extensive updates. The risk assessment will ultimately be used to support decision-making at the site, using multiple appropriately weighted sources of information, as well as uncertainty, sensitivity, and value-of-information analyses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effects of Nose Radius and Extreme Cooling on Boundary-Layer Transition for Three Smooth 15°-Cone-Cylinders in Free Flight at Mach Numbers to 8.50

Three highly polished 15°- included- angle cone- cylinders with hemispherical tips of several diameters ( 2, 3, and 4 in.) have been flown in order to obtain boundary- layer transition data at very low wall to local stream temperature ratios, and heat- transfer data. All surfaces had a 2-microinch average roughness height. Laminar flow existed over the entire hemispherical nose of the 2- and 3-inch-tip- diameter models throughout the complete flight history. Extreme cooling to wall to local stream temperature ratios at the sonic point as low as 0.20 did not cause transition on the nose for diameters as large as 3 inches. However, extreme cooling did cause early transition on the 4-inch model where it appears probable that transition occurred forward of the 45° station at a wall to local stream temperature ratio of about 0.26. Variations in tip diameter influenced transition downstream of the nose under conditions of extreme cooling. The 2-inch- tip model was laminar at all cone- cylinder stations at temperature ratios as low as 0.32 whereas the 3- and 4-inch-tip models were turbulent at the same local flow conditions but at higher wall to local temperature ratios. Transition on the cone and cylinder of the 3- and 4-inch- tip bodies appeared to be sensitive to local Mach number, and occurred at higher local temperature ratios when values of local Mach number were higher. Increasing the nose diameter from 2 to 3 inches significantly changed the local flow conditions for which laminar flow existed on the cone- cylinder afterbody. However, a further increase in tip size t o a 4-inch diameter had no discernable effect on the local flow conditions at transition. The transition results of the 3- and 4-inch-nose-diameter smooth bodies are similar to those observed on a 7/8-inch-nose-diameter body with roughened surfaces. Turbulent boundary layers resulted in both cases at very low wall to local stream temperature ratios. Both laminar and turbulent heat-transfer data were in good agreement with theoretical Stanton numbers when heat-transfer reduction due to tip blunting was considered.

CYLINDRICAL AFTERBODY↗

Concerning the Perturbations Generated by Flux Transfer Events

The interaction of the solar wind with the Earth's magnetosphere is often highly unsteady. Bursts of magnetic reconnection at multiple locations on the dayside equatorial magnetopause generate flux transfer events or FTEs: twisted ropes of interconnected magnetosheath and magnetospheric magnetic field lines. Once formed, the events move antisunward, displacing and perturbing the ambient media. This talk explores the perturbations predicted by both global numerical simulations and analytical models. Results from global hybrid code simulation confirm the predictions of analytical models indicating that the events generate standing forward slow mode waves as their speeds relative to the magnetosheath flow approach the Alfven velocity. Geometric considerations lead to the conclusion that events generated by component reconnection on the dayside magnetopause move poleward and exhibit strong signatures during intervals of southward IMF orientation, but move towards the flanks and exhibit weak signatures during intervals of northward IMF orientation. Changing event orientations and magnetosheath/magnetospheric magnetic field orientation can enhance the amplitudes of events reaching the flanks. Although the orientations of events on the flanks inferred from multispacecraft timing techniques are consistent with the predictions of the component reconnection model, occurrence patterns versus latitude and IMF orientation require an explanation in terms of both the component and antiparallel reconnection models.

Sibeck, David G.↗

Titan Haze

The Titan haze exerts a dominating influence on surface visibility and atmospheric radiative heating at optical and near-infrared wavelengths and our desire to understand surface composition and atmospheric dynamics provides a strong motivation to study the properties of the haze. Prior to the Cassini/Huygens missions the haze was known to be global in extent, with a hemispheric contrast asymmetry, with a complicated structure in the polar vortex region poleward of about 55 deg latitude, and with a distinct layer near 370 km altitude outside of the polar vortex at the time of the Voyager 2 flyby. The haze particles measured by the Pioneer and Voyager spacecraft were both highly polarizing and strongly forward scattering, a combination that seems to require an aggregation of small (several tens of nm radius) primary particles. These same properties were seen in the Cassini orbiter and Huygens Probe data. The most extensive set of optical measurements were made inside the atmosphere by the Descent Imager/Spectral Radiometer (DISR) instrument on the Huygens Probe. At the probe location as determined by the DISR measurements the average haze particle contained about 3000 primary particles whose radius is about 40 nm. Three distinct vertical regions were seen in the DISR data with differing particle properties. Refractive indices of the particles in the main haze layer resemble those reported by Khare et al. between O.3S and about 0.7 micron but are more absorbing than the Khare et al. results between 0.7 micron and the long-wavelength limit of the DISR spectra at 1.6 micron. These and other results are described by Tomasko et al., and a broader summary of results was given by Tomasko and West,. New data continue to stream in from the Cassini spacecraft. New data analyses and new laboratory and model results continue to move the field forward. Titan's 'detached' haze layer suffered a dramatic drop in altitude near equinox in 2009 with implications for the circulation and seasonal change in the stratosphere. The book chapter associated with this talk will also present new material on thermal-infrared data analysis and on new developments in laboratory work and haze microphysical modeling.

Anderson, Carrie M.↗

Interlaced Characterization and Calibration (ICC) for Improved Computational Simulation Credibility

Accurate material characterization and model calibration are pivotal for simulations used for high-consequence engineering decisions. Current characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data is collected for a specific model of interest, (3) provide only mean parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work developed a new paradigm—coined Interlaced Characterization and Calibration (ICC)—which drives forward the state-of-the-art in model calibration by bringing together recent advancements into one improved workflow. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) provides uncertainty metrics on the calibrated model parameters, and (4) incorporates these advances into a quasi real-time feedback loop. The ICC framework was validated synthetically with both low-fidelity and high-fidelity simulations paired with several different elastoplastic material models, and was also demonstrated experimentally with an aluminum 6061 cruciform exemplar specimen. Results showed that the ICC framework—in which Bayesian optimal experimental design actively guided the experiment— resulted in calibrations with similar or better accuracy than predetermined experiments based on subject matter expertise. Moreover, the ICC framework produced a complete model calibration— with quantified uncertainties on model parameters—in 1 week, a 5 - 10× increase in efficiency over traditional approaches. Thus, the ICC paradigm improves both the calibration process and quality, by (1) improving efficiency, which increases agility of solid mechanics modeling and enables utilization of computational simulation (CompSim) at earlier stages of the design cycle and (2) providing quantified, and in some cases reduced, parameter uncertainties, which increases confidence in model predictions and supports credible decision making.

97 MATHEMATICS AND COMPUTING↗

Optimum waiting time for acquisition of return link PN signals

The subject of this paper is a model that takes into account the effect that acquisition of a PN signal on the forward link has on the acquisition of a PN signal on the return link. The model is used to determine how long the start of the PN search on the return link should be delayed in order to minimize the combined acquisition time (delay + acquisition time) as a function of desired acquisition probability. It is assumed that the return link a-priori epoch uncertainty information is available. Software has been developed that models the return link PN receiver and incorporates in it the effect of forward link PN acquisition statistics in assessing performance.

Kolata, W.↗

Integrated Instrument Simulator Suites for Earth Science

The NASA Earth Observing System Simulators Suite (NEOS3) is a modular framework of forward simulations tools for remote sensing of Earth's Atmosphere from space. It was initiated as the Instrument Simulator Suite for Atmospheric Remote Sensing (ISSARS) under the NASA Advanced Information Systems Technology (AIST) program of the Earth Science Technology Office (ESTO) to enable science users to perform simulations based on advanced atmospheric and simple land surface models, and to rapidly integrate in a broad framework any experimental or innovative tools that they may have developed in this context. The name was changed to NEOS3 when the project was expanded to include more advanced modeling tools for the surface contributions, accounting for scattering and emission properties of layered surface (e.g., soil moisture, vegetation, snow and ice, subsurface layers). NEOS3 relies on a web-based graphic user interface, and a three-stage processing strategy to generate simulated measurements. The user has full control over a wide range of customizations both in terms of a priori assumptions and in terms of specific solvers or models used to calculate the measured signals.This presentation will demonstrate the general architecture, the configuration procedures and illustrate some sample products and the fundamental interface requirements for modules candidate for integration.

emission↗

Integrated Instrument Simulator Suites for Earth Science

The NASA Earth Observing System Simulators Suite (NEOS3) is a modular framework of forward simulations tools for remote sensing of Earth's Atmosphere from space. It was initiated as the Instrument Simulator Suite for Atmospheric Remote Sensing (ISSARS) under the NASA Advanced Information Systems Technology (AIST) program of the Earth Science Technology Office (ESTO) to enable science users to perform simulations based on advanced atmospheric and simple land surface models, and to rapidly integrate in a broad framework any experimental or innovative tools that they may have developed in this context. The name was changed to NEOS3 when the project was expanded to include more advanced modeling tools for the surface contributions, accounting for scattering and emission properties of layered surface (e.g., soil moisture, vegetation, snow and ice, subsurface layers). NEOS3 relies on a web-based graphic user interface, and a three-stage processing strategy to generate simulated measurements. The user has full control over a wide range of customizations both in terms of a priori assumptions and in terms of specific solvers or models used to calculate the measured signals.This presentation will demonstrate the general architecture, the configuration procedures and illustrate some sample products and the fundamental interface requirements for modules candidate for integration.

radiative transfer↗

Report on the LLNL Global Full-waveform Inversion Workflow and Progress

LLNL has developed the SPiRaL global seismic tomography model based upon millions of body wave travel times and surface wave dispersion curves in the period range of 25-200 seconds (Simmons et al., 2021). The motivation to construct the SPiRaL model (which stands for S-, P-, Rayleigh, and Love waves) is based upon the guidance provided by a 3-D model workshop held in Berkeley in 2007 and documented in a paper by Zucca et al. (2009). It was concluded in Zucca et al. (2009) that a “practical and obtainable” model is one that predicts body wave travel times for event location while also predicting surface wave velocities for moment tensor estimation. We have largely followed that guidance which has led to the construction of SPiRaL. Looking forward, we now seek to improve the global model through full waveform inversion (FWI). This necessitates the development/customization of an efficient and flexible modeling workflow to automate the process using LLNL’s supercomputing platforms. This report briefly outlines the work that being performed at LLNL to develop the necessary tools and framework to update the SPiRaL global model with full waveform information.

58 GEOSCIENCES↗

On the Impact of Granularity of Space-Based Urban CO2 Emissions in Urban Atmospheric Inversions: A Case Study for Indianapolis, IN

Quantifying greenhouse gas (GHG) emissions from cities is a key challenge towards effective emissions management. An inversion analysis from the INdianapolis FLUX experiment (INFLUX) project, as the first of its kind, has achieved a top-down emission estimate for a single city using CO2 data collected by the dense tower network deployed across the city. However, city-level emission data, used as a priori emissions, are also a key component in the atmospheric inversion framework. Currently, fine-grained emission inventories (EIs) able to resolve GHG city emissions at high spatial resolution, are only available for few major cities across the globe. Following the INFLUX inversion case with a global 1x1 km ODIAC fossil fuel CO2 emission dataset, we further improved the ODIAC emission field and examined its utility as a prior for the city scale inversion. We disaggregated the 1x1 km ODIAC non-point source emissions using geospatial datasets such as the global road network data and satellite-data driven surface imperviousness data to a 3030 m resolution. We assessed the impact of the improved emission field on the inversion result, relative to priors in previous studies (Hestia and ODIAC). The posterior total emission estimate (5.1 MtC/yr) remains statistically similar to the previous estimate with ODIAC (5.3 MtC/yr). However, the distribution of the flux corrections was very close to those of Hestia inversion and the model-observation mismatches were significantly reduced both in forward and inverse runs, even without hourly temporal changes in emissions. EIs reported by cities often do not have estimates of spatial extents. Thus, emission disaggregation is a required step when verifying those reported emissions using atmospheric models. Our approach offers gridded emission estimates for global cities that could serves as a prior for inversion, even without locally reported EIs in a systematic way to support city-level Measuring, Reporting and Verification (MRV) practice implementation.

EIs↗