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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Influence of Source Propagation Direction and Shear Flow Profile in Impedance Eduction of Acoustic Liners

The acoustic impedance of liners is a key parameter for their design, and depends on the flow conditions, i.e., the sound pressure level and the presence of a grazing flow. The surface impedance of a locally reacting liner is defined as a local intrinsic property relating the acoustic pressure to the normal acoustic particle velocity at the liner surface. Impedance eduction techniques are now widely used to retrieve the impedance of liners in aeroacoustic facilities in the presence of a shear grazing flow. While surface impedance is intrinsic by definition, the educed impedance has recently been shown to depend on the direction of the incident waves relative to the mean flow. Different studies have investigated this issue by considering different acoustic propagation models used in the education process in the hope of matching the educed values. The purpose of the present work is to continue the previous investigations by evaluating the influence of the shear flow profile on the educed impedance, while considering a Bayesian inference process in order to evaluate the uncertainty on the educed values. The identified uncertainties were not able to totally account for the observed discrepancies between educed impedances.

Roncen, R.↗
NASA NTRS2018

Hydrodynamic Modeling of Solar Acoustic Wave Propagation

Acoustic waves in the sun are separated into two categories: global and local. Global oscillations are excited through convection, and resonate in the solar interior. As with all resonant media, propagation of these acoustic waves is limited to specific frequencies for a given wavenumber. These modes are easily distinguished in 1-v diagrams, such as in figure 1. Global modes dominate the spectrum in the 2-5 mHz range. Local oscillations, so-called "sunquakes", are usually excited by particularly strong flares and can be observed as ring-shaped waves propagating away from the point of impact. The propagation is best displayed in a time-distance diagram, where a slit passing from the source to some distance away. Distance is displayed on the x-axis, and the time-series is plotted along the y-axis; an example is shown in figure 1, showing the sunquake induced by the Sept 2017 X-class flare.

Hydrodynamic↗
DOE OSTI2024

Evolving Outer Heliosphere: Tracking Solar Wind Transients from 1 au to the VLISM with IBEX and Voyager 1

Abstract Interstellar Boundary Explorer (IBEX) observations of energetic neutral atom (ENA) fluxes from the heliosphere have greatly enriched our understanding of the interaction of the solar wind (SW) with the local interstellar medium (LISM). However, there has been recent controversy surrounding the inability of most ENA models to produce as high an intensity of ∼0.5–6 keV ENAs as IBEX observes at 1 au, especially as a function of time. In our previous study (E. J. Zirnstein et al.), we introduced a new model that utilizes a data-driven magnetohydrodynamic simulation of the SW–LISM interaction to propagate pickup ions through the heliosheath (HS) after they are nonadiabatically heated at the heliospheric termination shock. E. J. Zirnstein et al. only simulated and analyzed IBEX observations from the direction of Voyager 2. In this study, we expand our model to include fluxes from the direction of Voyager 1, as well as in the low-latitude part (middle) of the ribbon (10° below the ecliptic plane). We show that the model results at Voyager 1 are consistent with E. J. Zirnstein et al.’s results at Voyager 2 in terms of a secondary ENA source contribution of ≲20% from both directions. Our results in the middle of the ribbon also reproduce the data, when including a time-dependent secondary ENA source. Finally, we demonstrate with our simulation that three large pressure waves likely merged in the VLISM and were observed by Voyager 1 as “pf2,” while at least one of the wave’s effects in the HS was observed by IBEX as a brief enhancement in ENA flux in early 2016.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Numerical simulation of hyperbolic heat conduction with convection boundary conditions and pulse heating effects

The paper describes the numerical simulation of hyperbolic heat conduction with convection boundary conditions. The effects of a step heat loading, a sudden pulse heat loading, and an internal heat source are considered in conjunction with convection boundary conditions. Two methods of solution are presened for predicting the transient behavior of the propagating thermal disturbances. In the first method, MacCormack's predictor-corrector method is employed for integrating the hyperbolic system of equations. Next, the transfinite element method, which employs specially tailored elements, is used for accurately representing the transient response of the propagating thermal wave fronts. The agreement between the results of various numerical test cases validate the representative behavior of the thermal wave fronts. Both methods represent hyperbolic heat conduction behavior by effectively modeling the sharp discontinuities of the propagating thermal disturbances.

Glass, David E.↗
NASA NTRS2016

Analysis and Simulations of Space Radiation Induced Single Event Effects and Transients

Spacecraft electronic are affected by the space radiation environment. Among the different types of radiation effects that can affect spacecraft electronics is the single event transient (SET). The space environment is responsible for many of the single event transients which can upset the performance of the spacecraft avionics hardware. In this chapter we first explore the origins of single event transients, then explore the modeling of a single event transient in digital and analog circuit. The chapter also addresses the concept of crosstalk that could develop among digital circuits in the present of a SET event. The chapter also provides a discussion of SET hardening. We then provide a discussion concerning propagation of a single event transient event at the local, subsystem, and system level in a spacecraft using two different models, one of the models developed by the author, known as the state transition model. The final goal of the chapter is to provide a qualitatively methodology for assessing single event transients and its effects so that spacecraft avionics engineers can develop either hardware or software countermeasures in their designs. SET is not a form of electromagnetic interference (EMI) in its origin, but semantically SET is very similar to EMI because they are both caused a current source not previously accounted for. SET has the same effects as EMI and it can cause interference problems in electronic circuits via multiple coupling mechanisms similar to EMI, and therefore makes such circuits incompatible.

Perez, Reinaldo J.↗
DOE OSTI2022

Proper orthogonal decomposition based reduced-order modeling of flux-Limited gray thermal radiation

Here, in this work, a proper orthogonal decomposition (POD) based reduced-order model (ROM) is developed to solve gray, flux-limited thermal radiation diffusion. We focus on the variable opacity radiation penetration benchmark posed by Olson, Auer, and Hall. The T -3 relationship for opacity in conjunction with high-temperature radiation penetrating an initially cold material produces a strong thermal radiation shock. This class of problems is particularly challenging for standard POD-based reduced-order modeling due to the nonlinearities presented by 1) the T 4 source term, and 2) flux-limited diffusion operator. To address these challenges and develop a cost competitive ROM, we employ a “hyper-reduction” technique through discrete empirical interpolation (DEIM) and allow for adaptive reduced-order projections through principal interval decomposition (PID). Performance of the proposed methodology is quantified by comparing the cost savings and accuracy relative to a full-order computation. Reference solutions and snapshot data are obtained through a full-order calculation performed by the University of Chicago maintained astrophysics code, FLASH. For consistency and potential extensibility, the developed ROM is also implemented in FLASH. We find that in the initialization regime, where the thermal radiation wave is initially created by the warming of the material, this class of problems is highly reducible and suitable for POD-based ROMs. However, the strong convective nature of the wave propagation regime is less reducible and more challenging to create an efficient ROM.

42 ENGINEERING↗
NASA NTRS1984

Fan noise reduction achieved by removing tip flow irregularities behind the rotor - forward arc test configurations

The noise source caused by the interaction of the rotor tip flow irregularities (vortices and velocity defects) with the downstream stator vanes was studied. Fan flow was removed behind a 0.508 meter (20 in.) diameter model turbofan through an outer wall slot between the rotor and stator. Noise measurements were made with far-field microphones positioned in an arc about the fan inlet and with a pressure transducer in the duct behind the stator. Little tone noise reduction was observed in the forward arc during flow removal; possibly because the rotor-stator interaction noise did not propagate upstream through the rotor. Noise reductions were maded in the duct behind the stator and the largest decrease occurred with the first increment of flow removal. This result indicates that the rotor tip flow irregularity-stator interaction is as important a noise producing mechanism as the normally considered rotor wake-stator interaction.

Dittmar, J. H.↗
NASA NTRS1992

Automatic differentiation of advanced CFD codes for multidisciplinary design

Automated multidisciplinary design of aircraft and other flight vehicles requires the optimization of complex performance objectives with respect to a number of design parameters and constraints. The effect of these independent design variables on the system performance criteria can be quantified in terms of sensitivity derivatives which must be calculated and propagated by the individual discipline simulation codes. Typical advanced CFD analysis codes do not provide such derivatives as part of a flow solution; these derivatives are very expensive to obtain by divided (finite) differences from perturbed solutions. It is shown that sensitivity derivatives can be obtained accurately and efficiently using the ADIFOR source translator for automatic differentiation. In particular, it is demonstrated that the 3-D, thin-layer Navier-Stokes, multigrid flow solver called TLNS3D is amenable to automatic differentiation in the forward mode even with its implicit iterative solution algorithm and complex turbulence modeling. It is significant that by using computational differentiation, consistent discrete nongeometric sensitivity derivatives have been obtained from an aerodynamic 3-D CFD code in a relatively short time, e.g., O(man-week) not O(man-year).

Bischof, C.↗
NASA NTRS2000

GENOA-PFA: Progressive Fracture in Composites Simulated Computationally

GENOA-PFA is a commercial version of the Composite Durability Structural Analysis (CODSTRAN) computer program that simulates the progression of damage ultimately leading to fracture in polymer-matrix-composite (PMC) material structures under various loading and environmental conditions. GENOA-PFA offers several capabilities not available in other programs developed for this purpose, making it preferable for use in analyzing the durability and damage tolerance of complex PMC structures in which the fiber reinforcements occur in two- and three-dimensional weaves and braids. GENOA-PFA implements a progressive-fracture methodology based on the idea that a structure fails when flaws that may initially be small (even microscopic) grow and/or coalesce to a critical dimension where the structure no longer has an adequate safety margin to avoid catastrophic global fracture. Damage is considered to progress through five stages: (1) initiation, (2) growth, (3) accumulation (coalescence of propagating flaws), (4) stable propagation (up to the critical dimension), and (5) unstable or very rapid propagation (beyond the critical dimension) to catastrophic failure. The computational simulation of progressive failure involves formal procedures for identifying the five different stages of damage and for relating the amount of damage at each stage to the overall behavior of the deteriorating structure. In GENOA-PFA, mathematical modeling of the composite physical behavior involves an integration of simulations at multiple, hierarchical scales ranging from the macroscopic (lamina, laminate, and structure) to the microscopic (fiber, matrix, and fiber/matrix interface), as shown in the figure. The code includes algorithms to simulate the progression of damage from various source defects, including (1) through-the-thickness cracks and (2) voids with edge, pocket, internal, or mixed-mode delaminations.

Murthy, Pappu L. N.↗
NASA NTRS1989

Are supernova remnants quasi-parallel or quasi-perpendicular accelerators

Observations of shock waves in the solar system which show a pronounced difference in the plasma wave and particle environment depending on whether the shock is propagating along or perpendicular to the interplanetary magnetic field are discussed. Theories for particle acceleration developed for quasi-parallel and quasi-perpendicular shocks, when extended to the interstellar medium suggest that the relativistic electrons in radio supernova remnants are accelerated by either the Q parallel or Q perpendicular mechanisms. A model for the galactic magnetic field and published maps of supernova remnants were used to search for a dependence of structure on the angle Phi. Results show no tendency for the remnants as a whole to favor the relationship expected for either mechanism, although individual sources resemble model remnants of one or the other acceleration process.

Spangler, S. R.↗
NASA NTRS2012

Gravity Wave Variances and Propagation Derived from AIRS Radiances

As the first gravity wave (GW) climatology study using nadir-viewing infrared sounders, 50 Atmospheric Infrared Sounder (AIRS) radiance channels are selected to estimate GW variances at pressure levels between 2-100 hPa. The GW variance for each scan in the cross-track direction is derived from radiance perturbations in the scan, independently of adjacent scans along the orbit. Since the scanning swaths are perpendicular to the satellite orbits, which are inclined meridionally at most latitudes, the zonal component of GW propagation can be inferred by differencing the variances derived between the westmost and the eastmost viewing angles. Consistent with previous GW studies using various satellite instruments, monthly mean AIRS variance shows large enhancements over meridionally oriented mountain ranges as well as some islands at winter hemisphere high latitudes. Enhanced wave activities are also found above tropical deep convective regions. GWs prefer to propagate westward above mountain ranges, and eastward above deep convection. AIRS 90 field-of-views (FOVs), ranging from +48 deg. to −48 deg. off nadir, can detect large-amplitude GWs with a phase velocity propagating preferentially at steep angles (e.g., those from orographic and convective sources). The annual cycle dominates the GW variances and the preferred propagation directions for all latitudes. Indication of a weak two-year variation in the tropics is found, which is presumably related to the Quasi-biennial oscillation (QBO). AIRS geometry makes its out-tracks capable of detecting GWs with vertical wavelengths substantially shorter than the thickness of instrument weighting functions. The novel discovery of AIRS capability of observing shallow inertia GWs will expand the potential of satellite GW remote sensing and provide further constraints on the GW drag parameterization schemes in the general circulation models (GCMs).

AIRS↗
NASA NTRSDate not supplied

Anomaly Detection in Flight Operational Data Using Deep Learning

In this session, we demonstrate two recently developed deep learning models for anomaly detection in flight operational data by the Data Sciences Group at NASA Ames Research Center. The first model is Convolutional Variational Auto-Encoder (CVAE) [1], which is an unsupervised deep encoder-decoder model, designed specifically for finding anomalies in heterogeneous multivariate time series data. We will demonstrate its application to finding anomalies in streaming data from NASA’s Digital Information Platform’s Fuser source. CVAE identifies data instances that are not representative of expected nominal behavior as anomalous. Since it is an unsupervised approach, the flagged anomalies will need to be reviewed by the subject matter experts (SMEs) for validation and labeling and is designed to assist with vulnerability discovery within Safety Monitoring System programs. The second model is Robust and Explainable Semi-supervised Anomaly Detection (RESAD) model [2], which builds on CVAE to allow learning from both minimally labeled data (previously reviewed by the SMEs) as well as majority unlabeled data. RESAD takes advantage of graph theoretic techniques to propagate the labels from the labeled data to the unlabeled data based on a pre-defined similarity metric and structures the learned feature space from flight time-series so that data of the same class would cluster tightly together. This model characteristic is enabled by training with an augmented loss function and allows learning of a more informative feature space for down-stream tasks such as search and active learning. We demonstrate RESAD using data from the NASA DASHlink project [3].

anomaly detection↗
NASA NTRS1995

Reflection and transmission of GMIR shock at the heliopause and their relation to the 2- and 3-kHz radio emissions

We use Voyager 2 plasma and magnetic field data together with a one-fluid MHD model to study the interactions of the 1991 Global Merged Interaction Region (GMIR) shock with the heliopause. The 1991 GMIR is an extraordinarily large global solar wind structure in radial, longitudinal and latitudinal extents. It has a strong shock at the leading edge. After its penetration through the termination shock, the GMIR shock first propagates through the subsonic solar wind, then interacts with the heliopause. The interaction produces a transmitted shock propagating outward in the interstellar medium, and a reflected shock propagating backward in the subsonic solar wind. We identify the reflected shock and the transmitted shock as the possible source of the radio noise detected at Voyagers. The plasma frequency behind the reflected and the transmitted shock can be, respectively, responsible for the 2- and 3-kHz radio emissions. The two bands of radio noise are emitted from sources on both sides of the heliopause starting at about the same time. If the emission is generated by f(sub p)-radiation then the heliopause is located at R approximately 130 AU. If the emission is generated by 2f(sub p)-radiation the n R approximately 150 AU. Because the relative speed of the interstellar plasma with respect to the sun appears to be sub-Alfvenic, it is very unlikely there is a fast-bow shock of the heliosphere.

Whang, Y. C.↗
NASA NTRS1996

Ultraviolet Opacity and Fluorescence in Supernova Envelopes

By the time the expanding envelope of a Type 2 supernova becomes transparent in the optical continuum, most of the gamma-ray luminosity produced by radioactive Fe/Co/Ni clumps propagates into the hydrogen/helium envelope and is deposited there, if at all. The resulting fast electrons excite He 1 and H 1, the two- photon continua of which are the dominant internal sources of ultraviolet radiation. The UV radiation is blocked by scattering in thousands of resonance lines of metals and converted by fluorescence into optical and infrared emission lines that escape freely. We describe results of Monte Carlo calculations that simulate non-LTE scattering and fluorescence in more than five million allowed lines of Ca, Sc, Ti, V, Cr, Mn, Fe, Co, and Ni. For a model approximating conditions in the envelope of SN 1987A, the calculated emergent spectrum resembles the observed one. For the first 2 yr after explosion, the ultraviolet radiation (lambda less than or approximately equals 3000) is largely blocked and converted into a quasi continuum of many thousands of weak optical and infrared emission lines and some prominent emission features, such as the Ca 2 lambdalambda8600 triplet. Later, as the envelope cools and expands, it becomes more transparent, and an increasing fraction of the luminosity emerges in the UV band.

Li, Hongwei↗
NASA NTRS2019

Assessment of Mixer-Ejector Nozzle with Thermal Acoustic Shield for Jet Noise Reduction

A tendency for excessive exhaust jet mixing noise from low bypass ratio turbofan engines is recognized as a key challenge in the design of commercial supersonic aircraft. In this work we investigate a unique combination of two noise mitigation methods as a novel strategy to reduce jet mixing noise. First, a thermal acoustic shield (TAS) is used to reflect high frequency acoustic waves at small angles to the jet axis; second, a mixer-ejector (ME) nozzle is used to mechanically shield noise propagating at large angles to the axis. The ME shroud also provides a convenient location for a TAS nozzle and improves TAS effectiveness by limiting the downstream extent of high frequency noise generation. In an additional benefit for a velocity-matched TAS stream, the ME allows a reduction in strength of the TAS outer shear layer which could serve as a secondary noise source. The present work provides a quantitative assessment of the ME-TAS concept, using a combination of RANS CFD simulations, acoustic analogy calculations for the farfield Green's function, and surrogate-based modeling and parameter space exploration. We first evaluate a subscale configuration, then use scaling arguments to apply subscale results to the systems-level analysis of a flight configuration; the latter configuration includes a generic low bypass ratio turbofan engine with an engine-driven electric generator for supplementary heating of the TAS stream. Additional RANS CFD calculations are performed for a notional ME-TAS geometry based on the full scale configuration, and various modeling assumptions and operational characteristics are evaluated. The ME-TAS concept is shown to provide effective shielding for high frequency jet noise, and should enable comparable noise suppression to a stand-alone ME of considerably greater length, weight and drag. In addition to investigating the integrated ME-TAS system, the present work differs from previous research into TAS and related fluidic shield concepts through the inclusion of modern numerical analysis tools and the systematic numerical examination of various design parameters.

thermal acoustic shield↗
DOE OSTI2025

MOOSE ProbML: Parallelizable Probabilistic Machine Learning and Uncertainty Quantification Capabilities

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a widely used open- source finite element software for performing multiphysics multiscale simulations in a massively parallel fashion. Recently, the computational team at Idaho National Laboratory (INL) has implemented Probabilistic Machine Learning (ProbML) capabilities in MOOSE—in a parallelized fashion—and enable active learning with large-scale computational models for tasks such as surrogate model development, scale bridging, forward/inverse uncertainty quantification (UQ), Bayesian optimization, etc. This presentation summarizes these developments in MOOSE along with demonstrations on several real applications relevant to nuclear energy. At the fundamental level, samplers like Monte Carlo/Latin Hypercube, variance reduction, parallelized Markov Chain Monte Carlo (MCMC) support uncertainty propagation in both forward and inverse settings. These samplers can be integrated with the Gaussian processes (GP) suite in MOOSE, which offer several variants like scalar GPs, multi-output GPs, and deep GPs, to enable active learning. These GPs can be tuned using gradient-based optimization methods like Adam and its variants or gradient-free methods like the elliptical slice sampler (a variant of MCMC adept under Gaussian settings) for more complex covariance kernels or likelihoods whose gradient computations can be cumbersome. A variety of batch acquisition functions permit parallelized evaluation of the computational model and support different learning objectives with high efficiency like Bayesian inference, global surrogate development, optimization, etc. Furthermore, libtorch integration supports training, evaluation, and re-training of neural networks and other complex machine learning models in active learning settings. The impacts of these developments are shown on several real applications: (1) nuclear fuel inverse UQ and model inadequacy assessment using the Kennedy O’Hagan framework; (2) uncertainty aware surrogate modeling for additive manufacturing to predict field quantities; (3) nuclear reactor rare events analysis; and (4) complex fluid flow prediction using a global surrogate with quantified prediction uncertainty. Finally, the outlook of MOOSE ProbML is discussed for both outer-loop and inner-loop computations in the broad view to accelerate fuels and materials qualification, address gaps in knowledge and data, and assess new reactor/fuel systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2014

The New NASA Orbital Debris Engineering Model ORDEM 3.0

The NASA Orbital Debris Program Office (ODPO) has released its latest Orbital Debris Engineering Model, ORDEM 3.0. It supersedes ORDEM 2.0. This newer model encompasses the Earth satellite and debris flux environment from altitudes of low Earth orbit (LEO) through geosynchronous orbit (GEO). Debris sizes of 10 microns through 1 m in non-GEO and 10 cm through 1 m in GEO are modeled. The inclusive years are 2010 through 2035. The ORDEM model series has always been data driven. ORDEM 3.0 has the benefit of many more hours from existing data sources and from new sources that weren't available to past versions. Returned surfaces, ground tests, and remote sensors all contribute data. The returned surface and ground test data reveal material characteristics of small particles. Densities of fragmentation debris particles smaller than 10 cm are grouped in ORDEM 3.0 in terms of high-, medium-, and lowdensities, along with RORSAT sodium-potassium droplets. Supporting models have advanced significantly. The LEO-to-GEO ENvironment Debris model (LEGEND) includes an historical and a future projection component with yearly populations that include launched and maneuvered intacts, mission related debris (MRD), and explosion and collision fragments. LEGEND propagates objects with ephemerides and physical characteristics down to 1 mm in size. The full LEGEND yearly population acts as an a priori condition for a Bayesian statistical model. Specific, well defined populations are added like the Radar Ocean Reconnaissance Satellite (RORSAT) sodium-potassium (NaK) droplets, recent major accidental and deliberate collision fragments, and known anomalous debris event fragments. For microdebris of sizes 10 microns to 1 mm the ODPO uses an in-house Degradation/Ejecta model in which a MLE technique is used with returned surface data to estimate populations. This paper elaborates on the upgrades of this model over previous versions highlighting the material density splits and consequences of that to the penetration risk to spacecraft.

Krisko, P. H.↗
DOE OSTI2020

Radius-to-frequency Mapping and FRB Frequency Drifts

We build a model of radius-to-frequency mapping in magnetospheres of neutron stars and apply it to frequency drifts observed in fast radio bursts (FRBs). We assume that an emission patch propagates along the dipolar magnetic field lines, producing coherent emission with frequency, direction, and polarization defined by the local magnetic field. The observed temporal evolution of the frequency depends on the relativistic effects of time contraction and the curvature of the magnetic field lines. The model generically produces linear scaling of the drift rate, $\dot{\omega}$ ∝-${\omega}$, matching both numerically and parametrically the rates observed in FBRs; a more complicated behavior of $\dot{\omega}$ is also possible. Fast rotating magnetospheres produce higher drifts rates for similar viewing parameters than the slowly rotating ones. In the case of repeaters, the same source may show variable drift patterns depending on the observing phase. We expect rotational of polarization position angle through a burst, though by smaller amount than in radio pulsars. All of these findings compare favorably with properties of FBRs, strengthening their possible loci in the magnetospheres of neutron stars.

Magnetars↗