Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source and propagation 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 109 records · Page 6

Dynamics of Intra-Cell Thermal Front Propagation in Lithium-Ion Battery Safety Issues

Thermal runaway (TR), a critical failure mode in lithium-ion batteries (LIBs), poses significant safety risks and hinders wider application of LIBs. TR typically begins at a localized heat source and spreads across the cell. Understanding thermal front propagation (TFP) characteristics, such as front and velocity, is crucial for assessing energy release and temperature distribution for battery hazardous estimation. Recent studies assume that TR within cells propagates at a near-constant velocity, based on the reaction kinetics and thermal properties. Here, in this study, an intra-battery TR model is further proposed and it indicates that TFP velocity stabilizes when the front is distanced from the heat source. Theoretical estimates for propagation velocity and front are developed and validated through numerical simulations and experimental tests from the NREL Battery Failure Databank. The energy release rate during TFP and the impact of preheating based on a point heat source are explored. This work clarifies the long-standing clouds of the thermal font propagation behaviors within the single cell, highlights the power and beauty of mathematics modeling to describe the complicated thermal behaviors, and provides important guidelines for thermal hazardous understanding for next-generation batteries.

25 ENERGY STORAGE↗

Phased Microphone Array on Aircraft Fuselage

This paper studies some challenging features for phased microphone arrays on aircraft fuselage, including the effects of sound scattering such as refraction by nonuniform boundary layer flows, reflection by curved surfaces, and diffraction by sharp edges and smooth geometry. The objective is to show if and how source locations and source amplitudes can be accurately captured by the array data analysis, in the presence of these features. The study starts with canonical problems to separately examine the individual features, making use of their respective analytical solutions, which provide exact results for all quantities needed for the array data processing. To simulate realistic flight environments, a generic aircraft geometry is then utilized, involving multiple scattering mechanisms and multiple noise sources to model aircraft engine noise. For this application, numerical solutions are computed to provide noise propagation and scattering from the sources to the microphones, as well as the quantities needed for array data analysis. The study reveals the effects of various features on fuselage-mounted arrays and provides guidelines to account for these effects in array data analysis.

Yueping Guo↗

Using Deep Learning to Develop a High Resolution Planetary Boundary Layer Model for Infrasound Propagation

Infrasound, with frequencies less than 20 Hz, is generated by both natural and anthropogenic sources. When one of these sources exerts a force on the atmosphere, infrasonic waves are generated. The propagation of these waves largely depends on temperature, wind speed, and wind direction. Previous work has used deep learning to accurately predict atmospheric specifications to altitudes of ~40 km. However, this model breaks down for local distances because it is too low resolution. Here we use a high-resolution meteorological dataset collected in Las Vegas, Nevada, USA to develop a deep learning model that can predict temperature, wind speed, and wind direction. Predictions are compared to ground truth observations to show that the model performs well at predicting temperature and wind direction but struggles with prediction wind speed. Model limitations and improvements are also discussed.

54 ENVIRONMENTAL SCIENCES↗

Nuclear data uncertainty propagation and modeling uncertainty impact evaluation in neutronics core simulation

Uncertainty analysis is a critical requirement in reactor simulation as it is used to quantify the reliability of best-estimate calculation. A comprehensive uncertainty analysis should characterize all sources of uncertainties in a computationally-feasible and scientifically-defendable manner. Here we employ a well-established reduced order modeling (ROM) based uncertainty quantification methodology to propagate uncertainties throughout neutronic calculations. ROM relies on recent advances in randomized data mining techniques applied to large data streams. In our proposed implementation, the nuclear data uncertainties are first propagated from multi-group level through lattice physics calculation to generate few-group parameter uncertainties, described using a vector of mean values and a covariance matrix. Employing an ROM-based compression of the covariance matrix, the few-group uncertainties are then propagated through downstream core simulation in a computationally efficient manner. This straightforward approach, albeit efficient as compared to brute force forward and/or adjoint-based methods, often employs a number of assumptions that have been unquestioned in the literature of neutronic uncertainty analysis. This manuscript argues that these assumptions could introduce another source of uncertainty referred to as modeling uncertainties, whose magnitude needs to be quantified in tandem with nuclear data uncertainties. Thus, our primary goal is to explore the interactions between these two uncertainty sources in order to assess whether modeling uncertainties have an impact on parameter uncertainties. To explore this endeavor, the impact of a number of modeling assumptions on core attributes uncertainties is quantified. The study employs a CANDU reactor model, with Serpent and NEWT as lattice physics solvers and NESTLE-C as core simulator. The modeling assumptions investigated include those related with the uncertainty propagation method employed, e.g., deterministic vs. stochastic, the few-group energy structure employed to represent the cross-sections, the resonance treatment in lattice physics calculation, the reference values for the cross-section, and the number of samples employed to render ROM compression. Results indicate that some of the modeling assumptions could have a non-negligible impact on the core responses propagated uncertainties, highlighting the need for a more comprehensive approach to combine parameter and modeling uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Detecting Large Explosions With Machine Learning Models Trained on Synthetic Infrasound Data

Explosions produce low-frequency acoustic (infrasound) waves capable of propagating globally, but the spatio-temporal variability of the atmosphere makes detecting events difficult. Machine learning (ML) is well-suited to identify the subtle and nonlinear patterns in explosion infrasound signals, but a previous lack of ground-truth data inhibited training of generalized models. We introduce a physics-based method that propagates infrasound sources through realistic atmospheres to create 28,000 synthetic events, which are used to train ML classifiers. A simple artificial neural network and modern temporal convolutional network discriminate synthetic events from background noise with >90% accuracy and, more importantly, successfully identify the majority of real-world explosion signals recorded during the Humming Road Runner experiment. ML models trained entirely on physics-based synthetics advance explosion detection capabilities and make ML more viable to related fields lacking training data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Numerical modelling of turbulent flow in a combustion tunnel

A numerical technique is presented for the analysis of turbulent flow associated with combustion. The technique uses Chorin's random vortex method (rvm), an algorithm capable of tracing the action of elementary turbulent eddies and their cumulative effects without imposing any restriction upon their motion. In the past, the rvm has been used with success to treat nonreacting turbulent flows, revealing in particular the mechanics of large-scale flow patterns, the so-called coherent structures. Introduced here is a flame propagation algorithm, also developed by Chorin, in conjunction with volume sources modelling the mechanical effects of the exothermic process of combustion. As an illustration of its use, the technique is applied to flow in a combustion tunnel where the flame is stabilized by a back-facing step. Solutions for both nonreacting and reacting flow fields are obtained which satisfactorily describe the essential features of turbulent combustion in a lean propane-air mixture that were observed in the laboratory by means of high speed Schlieren photography.

Ghoniem, A. F.↗

A theoretical model of the wave particle interaction of plasma in space

A theoretical model, based on the kinetic theory for the perturbation of plasma in the magnetosphere, is proposed to study the observed disturbances which are caused by both natural and artificial sources that generate wave-like perturbations propagating around the globe. The proposed model covers the wave propagation through a media of transitional (from collisional to collisionless) fully ionized magnetoactive plasma. A systematic formulation of the problem is presented and the method of solution for the transitional model of magnetosphere is discussed. The possible emission of hydromagnetic waves in the magnetosphere during the quiet and disturbed time are also discussed.

Wu, S. T.↗

Gravity wave momentum flux in the lower stratosphere over convection

This work describes a method for estimating vertical fluxes of horizontal momentum carried by short horizontal scale gravity waves (lambda(sub x) = 10-100 km) using aircraft measured winds in the lower stratosphere. We utilize in situ wind vector and pressure altitude measurements provided by the Meteorological Measurement System (MMS) on board the ER-2 aircraft to compute the momentum flux vectors at the flight level above deep convection during the tropical experiment of the Stratosphere Troposphere Exchange Project (STEP-Tropical). Data from Flight 9 are presented here for illustration. The vertical flux of horizontal momentum these observations points in opposite directions on either side of the location of a strong convective updraft in the cloud shield. This property of internal gravity waves propagating from a central source compares favorably with previously described model results.

Alexander, M. Joan↗

Element Abundances in the Galactic Cosmic Rays with Atomic Number (Z) in the Interval 30 is less than or equal to Z is less than or equal to 40

We present new results on the elemental abundances of galactic cosmic rays with atomic number, Z, greater than 30, and comparison of these observations with abundances expected from galactic propagation of various suggested models of the cosmic-ray source. We combine preliminary results from the 2003-04 flight of the Trans-Iron Galactic Element Recorder (TIGER) cosmic-ray detector with previously reported results from the 2001-02 flight. This instrument flew over Antarctica for nearly 32 days at a mean atmospheric depth of 5.2 mb in December 2001 - January 2002. At the time of submission of this abstract, January 8, 2004, TIGER was again in the air over Antarctica having completed 22 days of an expected 30day flight at a mean atmospheric depth of about 4 nb, Data from the first flight demonstrated excellent resolution of individual elements, and we expect similar resolution from the second flight.

Barbier, Louis↗

Density and white light brightness in looplike coronal mass ejections - Temporal evolution

Three ambient coronal models suitable for studies of time-dependent phenomena were used to investigate the propagation of coronal mass ejections initiated in each atmosphere by an identical energy source. These models included those of a static corona with a dipole magnetic field, developed by Dryer et al. (1979); a steady polytropic corona with an equatorial coronal streamer, developed by Steinolfson et al. (1982); and Steinolfson's (1988) model of heated corona with an equatorial coronal streamer. The results indicated that the first model does not adequately represent the general characteristics of observed looplike mass ejections, and the second model simulated only some of the observed features. Only the third model, which included a heating term and a streamer, was found to yield accurate simulation of the mess ejection observations.

Steinolfson, R. S.↗

Imaging the Shallow Structure of the Yucca Flat at the Source Physics Experiment Phase II Site with Horizontal-to-Vertical Spectral Ratio Inversion and a Large- N Seismic Array

The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.

58 GEOSCIENCES↗

Structure and Dynamics of the Quiet Solar Chromosphere

The grant supported research on the structure of the quiet, nonmagnetic chromosphere and on wave excitation and propagation in both the nonmagnetic chromosphere and the magnetic network. The work on the structure of the chromosphere culminated in the recognition that between two competing views of the solar chromosphere, older models by Avrett and collaborators (referred to as VAL) and the newer, dynamical model by Carlsson & Stein (referred to as CS), the clear decision is in favor of the older models, and this in spite of the evident lack of physics, which does not include wave motion and oscillations. The contrast between the static VAL models and the dynamical CS model can be stated most succinctly by comparing the temperature variation implied by the VAL models and the temperature fluctuations of the CS model, which are, respectively, of the order of 10% for the VAL model (at heights where hydrogen is 50% ionized) and a factor of 10 (at the upper boundary of their chromospheric model). The huge fluctuations of the CS model have never been observed, whereas the smaller temperature variations of the VAL models are consistent with ground-based and space-based observations. While it should be obvious which model describes the Sun and which one fails, the case is far from settled in the minds of solar physicists. Thus, much educational work remains to be done and, of course, more research to develop arguments that make the case more convincing. The research on waves and oscillations has been based on a unified theory of excitation of acoustic waves in the field-free atmosphere and of transverse and longitudinal waves in magnetic flux tubes located in the magnetic network by noting, first, that impulsive excitation of all these waves in gravitationally stratified media leads to oscillations at the respective cutoff frequencies and, second, that the observed oscillation frequencies in the nonmagnetic and magnetic parts of the chromosphere match corresponding cutoff frequencies in the upper solar photosphere. The dynamical simulations by Carlsson and Stein have been most instructive and of fundamental importance for understanding wave propagation in a stratified medium by their "flaws", the most important of which is an intensity excess at the H2v emission peak in the H line of Ca II, which is surprising since the observed intensity should have been an upper limit to the simulated intensity. The only plausible explanation for a predicted intensity that is higher than observed is that energy is spread horizontally in upward propagation in the Sun, but not in the plane-wave modeling of CS (as well as by almost everyone else). Investigations of the horizontal size of the region disturbed by the upward-propagating shock in the acoustic-wave propagation implies that the waves in H2v bright-point oscillations emanate from a point source with a diameter corresponding to the width of an intergranular lane, about 100 km, and reach a size of about 4000 km in the upper layers of the chromosphere. Linear, analytic modeling of waves emanating from a point source in a stratified atmosphere shows that the upward-expanding propagation channel does not have sharp boundaries and that the shape of the wave front depends on the order of the wave behind the initial pulse. Otherwise, the behavior of the linear waves resembles that of the nonlinear shock waves observed in the Sun. Research that needs to be done to firm up the conclusions reached above concerns the numerical simulations of nonlinear waves and oscillations in a three-dimensional stratified atmosphere with impulsive excitation, and observations linking directly the horizontal size of the disturbed area in upward propagation to individual waves.

Kalkofen, Wolfgang↗

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗

Isotopic composition of cosmic-ray boron and nitrogen

New measurements of the cosmic-ray boron and nitrogen isotopes at earth and of the elemental abundances of boron, carbon, nitrogen, and oxygen are presented. A region of mutually allowed values for the cosmic-ray nitrogen source ratios is determined, and the cosmic-ray escape mean free path is determined as a function of energy using a leaky box model for cosmic-ray propagation in the Galaxy. Relative to O-16, a N-15 source abundance consistent with solar system composition and a N-14 source abundance which is a factor of about three underabundant relative to the solar value are found.

Krombel, K. E.↗

A generalized machine learning workflow to visualize mechanical discontinuity

Accurate detection and mapping of mechanical discontinuity in materials has widespread industrial and research applications. Herein, we developed a generalized machine-learning framework for visualizing single mechanical discontinuity embedded in material of any composition, velocity, density, porosity, and size with limited data. The proposed visualization of discontinuity requires accurate estimations of the length, location, and orientation of the embedded discontinuity by processing multipoint wave-transmission measurements. k-Wave simulator is used to create a large dataset of elastic waveforms recorded during multi-point wave-transmission measurements through materials containing single mechanical discontinuity. k-Wave simulator considers the wave attenuation, dispersion, and mode conversion in wave motion. Discrete wavelet transform (DWT) and statistical feature extraction are essential for data preprocessing prior to the data-driven model development. DWT also minimizes the effect of noise. Using hyper-parameter tuning and cross validation, gradient boosting regression can visualize the mechanical discontinuity with an accuracy of 0.85, in terms of coefficient of determination. A double-layered neural network-based regression has better performance with an accuracy of 0.95. Use of convolutional neural network converts the predictive task from a waveform processing to an image processing problem. Convolutional neural network achieved a generalization performance of 0.91. The proposed generalized workflow requires robust simulation of wave propagation, signal processing, feature engineering, and model evaluation. Sensors closest to the source and those located opposite the source are the most significant for the desired visualization. Notably, the sensors closest to the source capture the non-linear associations, whereas the sensor on the border opposite to the source capture the linear associations between the measured waveforms and the properties of the mechanical discontinuity.

42 ENGINEERING↗

Propagation of electromagnetic waves parallel to the magnetic field in the nightside Venus ionosphere

The propagation of electromagnetic waves parallel to the magnetic field in the nightside Venus ionosphere is presented in a theoretical and numerical analysis. The model assumes a source of electromagnetic radiation in the Venus atmosphere, such as that produced by lightning. Specifically addressed is wave propagation in the altitude range z = 130-160 km at the four frequencies detectable by the Pioneer Venus Orbiter Electric Field Detector: 100 Hz, 730 Hz, 5.4 kHz, and 30 kHz. Parameterizations of the wave intensities, peak electron density, and Poynting flux as a function of magnetic field are presented. The waves are found to propagate most easily in conditions of low electron density and high magnetic field. The results of the model are consistent with observational data.

Huba, J. D.↗

Time-resolved measurements of HO 2 radical in a heated plasma flow reactor

Time-resolved, absolute HO 2 number density in diluted H 2 –O 2 –Ar, CH 4 –O 2 -Ar, and C 2 H 4 –O 2 –Ar mixtures excited by a repetitive ns pulse discharge in a heated plasma flow reactor is measured by Cavity Ringdown Spectroscopy (CRDS). The experimental results are obtained at $\textit{T}$ = 300-600 K and $\textit{P}$ = 130 Torr, both during the discharge pulse burst and in the afterglow. In this work, the HO 2 number density is inferred from the CRDS data using a spectral model exhibiting good agreement with previous measurements of absolute HO 2 absorption cross sections. In the room-temperature H 2 –O 2 mixture, as well as in CH 4 –O 2 and C 2 H 4 –O 2 mixtures over the entire temperature range studied, HO2 is generated only during the discharge burst and decays in the afterglow. However, in the H 2 –O 2 mixture at elevated temperatures, $\textit{T}$ = 400-600 K, HO 2 persists in the afterglow up to 10 ms after the discharge burst, comparable with the flow residence time in the reactor. Comparison with kinetic modeling shows that the sustained reactivity after the source of radicals is turned off is due to a chain propagation / hydrogen oxidation process, which dominates the radical recombination reactions. The kinetic modeling predictions are in good agreement with the relative HO 2 number density measured in all three mixtures, although the model underpredicts the absolute number densities in H 2 –O 2 at $\textit{T}$ = 400-600 K by up to a factor of two. Detection of the sustained low-temperature reactivity in H 2 –O 2 , initiated by the radical generation in the plasma, suggests that the plasma excitation may also affect kinetics of oxidation and reforming of fuels exhibiting low-temperature chemistry below hot ignition point.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The isotopic composition of cosmic rays with 5 is less than or equal to z which is less than or equal to 26

Results obtained from a high altitude balloon flight from Thompson, Canada in August, 1973 are reported. The instrument consisted of a spark chamber, a Lucite Gerenkov counter and thirteen layers of scintillators. For heavy particles the Cerenkov-range method of analysis was used to determine the mass of particles energetic enough to produce a Cerenkov signal and then stop in the layered scintillators. The data appear to be consistent with current cosmic-ray propagation models. Using a simple exponential path length propagation model this data is extrapolated to the cosmic-ray source and some implications of the data are discussed as to the nature of the source.

Fisher, A. J.↗