SEARCH · Engineering Papers

Results for “simulation analytics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Analytical Assessment of Simultaneous Parallel Approach Feasibility from Total System Error

In a simultaneous paired approach to closely-spaced parallel runways, a pair of aircraft flies in close proximity on parallel approach paths. The aircraft pair must maintain a longitudinal separation within a range that avoids wake encounters and, if one of the aircraft blunders, avoids collision. Wake avoidance defines the rear gate of the longitudinal separation. The lead aircraft generates a wake vortex that, with the aid of crosswinds, can travel laterally onto the path of the trail aircraft. As runway separation decreases, the wake has less distance to traverse to reach the path of the trail aircraft. The total system error of each aircraft further reduces this distance. The total system error is often modeled as a probability distribution function. Therefore, Monte-Carlo simulations are a favored tool for assessing a "safe" rear-gate. However, safety for paired approaches typically requires that a catastrophic wake encounter be a rare one-in-a-billion event during normal operation. Using a Monte-Carlo simulation to assert this event rarity with confidence requires a massive number of runs. Such large runs do not lend themselves to rapid turn-around during the early stages of investigation when the goal is to eliminate the infeasible regions of the solution space and to perform trades among the independent variables in the operational concept. One can employ statistical analysis using simplified models more efficiently to narrow the solution space and identify promising trades for more in-depth investigation using Monte-Carlo simulations. These simple, analytical models not only have to address the uncertainty of the total system error but also the uncertainty in navigation sources used to alert an abort of the procedure. This paper presents a method for integrating total system error, procedure abort rates, avionics failures, and surveillance errors into a statistical analysis that identifies the likely feasible runway separations for simultaneous paired approaches.

Madden, Michael M.↗

NASA NTRS1974

MIA computer simulation test results report

Results of the first noise susceptibility computer simulation tests of the complete MIA receiver analytical model are presented. Computer simulation tests were conducted with both Gaussian and pulse noise inputs. The results of the Gaussian noise tests were compared to results predicted previously and were found to be in substantial agreement. The results of the pulse noise tests will be compared to the results of planned analogous tests in the Data Bus Evaluation Laboratory at a later time. The MIA computer model is considered to be fully operational at this time.

Unger, G. E.↗

NASA NTRS2004

Crash Simulation of a Vertical Drop Test of a Commuter-Class Aircraft

A finite element model of an ATR42-300 commuter-class aircraft was developed and a crash simulation was executed. Analytical predictions were correlated with data obtained from a 30-ft/s (9.14-m/s) vertical drop test of the aircraft. The purpose of the test was to evaluate the structural response of the aircraft when subjected to a severe, but survivable, impact. The aircraft was configured with seats, dummies, luggage, and other ballast. The wings were filled with 8,700 lb. (3,946 kg) of water to represent the fuel. The finite element model, which consisted of 57,643 nodes and 62,979 elements, was developed from direct measurements of the airframe geometry. The seats, dummies, luggage, fuel, and other ballast were represented using concentrated masses. The model was executed in LS-DYNA, a commercial code for performing explicit transient dynamic simulations. Predictions of structural deformation and selected time-history responses were generated. The simulation was successfully validated through extensive test-analysis correlation.

Jackson, Karen E.↗

DOE OSTI2026

Consequential improvement acquisition function for efficient multi-fidelity Bayesian optimization

Abstract Surrogate-based Bayesian optimization has been widely applied in design optimization to increase sampling efficiency. However, the cost for each evaluation of the objective function can still be very high when physical experiments or large-scale simulations are involved. Multi-fidelity Bayesian optimization is the new approach to further improve the sampling efficiency by reducing the number of expensive samples at the highest fidelity level and supplementing them with less expensive ones at low-fidelity levels. In this paper, a new consequential improvement (CI) acquisition function is proposed to allow for the simultaneous selection of the solution and the fidelity level in problems with a known hierarchy of fidelity levels. The new CI acquisition function incorporates the consequential effectiveness of objective improvement with the considerations of cost, accuracy, and validity differences between high- and low-fidelity samples in engineering practice. The new method of multi-fidelity Bayesian optimization based on the CI is demonstrated with several analytical and simulation-based design examples. In the simulation-based design optimization example, the results show that the CI acquisition function has a decisive advantage in the sampling efficiency over the other methods of multi-fidelity Bayesian optimization with simultaneous selection. The results indicate that the proposed method is particularly advantageous in solving high-dimensional problems and when large cost ratios between high- and low-fidelity evaluations exist and high-fidelity validation is mandatory. Furthermore, the method robustly avoids the prevalent issue of over sampling at low-fidelity levels.

Aydogdu, Ibrahim [Georgia Institute of Technology,↗

NASA NTRS1973

Continuous and discrete describing function analysis of the LST system

A describing function of the control moment gyros (CMG) frictional nonlinearity is derived using the analytic torque equation. Computer simulation of the simplified Large Space Telescope (LST) system with the analytic torque expression is discussed along with the transfer functions of the sampled-data LST system, and the discrete describing function of the GMC frictionality.

Kuo, B. C.↗

NASA NTRS1991

Hypersonic waverider analysis - A comparison of numerical and experimental results

In this study two waverider configurations are analyzed numerically and experimentally. The geometries are designed using an inverse design and optimization code with the optimization procedure driven by maximized L/D and limited by geometric constraints. Euler and Navier-Stokes numerical simulations are performed with a partially flux-split finite-difference code developed at NASA Ames, F3D. Experimental data is obtained at NASA Langley in the Unitary Plan Wind Tunnel (UPWT) and the 20-Inch Mach 6 wind tunnel for flow visualization, force/moment, and surface pressure measurements at both on- and off-design conditions. Euler simulations are compared with analytic results obtained from conical flow theory, and Navier-Stokes simulations are compared with experimental results both qualitatively and quantitatively for on- and off-design conditions.

Jones, K. D.↗

NASA NTRS2013

Environmentally Friendly Rotorcraft Concepts: 2011-2014: US/France MoA Helicopter Aeromechanics - Task 5

Motivation: Greenhouse gas concentrations have increased well beyond their pre-industrial levels due to human activities; Aircraft emissions are becoming regulated in industrialized nations, so future rotorcraft will need to be designed for minimal environmental impact. Objectives: Determine suitable metrics to measure environmental impact of rotorcraft, focusing first on air pollution; Integrate performance metrics with rotorcraft design codes; Generate rotorcraft designs with minimal environmental impact. Outcomes: Obtain an understanding of how targeting reduced emissions affects the rotorcraft design process; Advance the capability of existing rotorcraft design tools to include environmental performance metrics. Tools: Evaluation and predesign - NDARC (NASA Design and Analysis of Rotorcraft) and CREATION (Concepts of Rotorcraft Enhanced Assessment Through Integrated Optimization Network); Flight simulation - CAMRAD II (Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics) and HOST (Helicopter Overall Simulation Tool).

Environmentally Friendly↗

DOE OSTI2021

Bayesian force fields from active learning for simulation of inter-dimensional transformation of stanene

Abstract We present a way to dramatically accelerate Gaussian process models for interatomic force fields based on many-body kernels by mapping both forces and uncertainties onto functions of low-dimensional features. This allows for automated active learning of models combining near-quantum accuracy, built-in uncertainty, and constant cost of evaluation that is comparable to classical analytical models, capable of simulating millions of atoms. Using this approach, we perform large-scale molecular dynamics simulations of the stability of the stanene monolayer. We discover an unusual phase transformation mechanism of 2D stanene, where ripples lead to nucleation of bilayer defects, densification into a disordered multilayer structure, followed by formation of bulk liquid at high temperature or nucleation and growth of the 3D bcc crystal at low temperature. The presented method opens possibilities for rapid development of fast accurate uncertainty-aware models for simulating long-time large-scale dynamics of complex materials.

Chemistry↗

DOE OSTI2023

Mechanical behavior analysis using small-size graphite disc compression testing with digital image correlation methods

The ASTM International Standard Test Method for Tensile Stress-Strain of Carbon and Graphite (ASTM C749) requires dog-bone shaped specimens. Unfortunately, specimen sizes required in ASTM C749 are much too large for irradiation within a test reactor and are incompatible with irradiation capsule volumes, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed as an alternative that provides a means for testing tensile properties using smaller specimens when geometric constraints are present. However, the important stress–strain relation throughout the test loading history that could be recorded using ASTM C749 is lost by using ASTM D8289 because of the difficulties involved with using traditional strain gauges or extensometers on such a small specimen. In this work, a digital image correlation (DIC) method, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 splitting tensile test on eight small samples of graphite grade Mersen 2114 to track their deformation during compressive loading. Initial results showed that the DIC technique can measure the surface displacement/strain on these small (Ø6 × 3 mm) graphite specimens with high accuracy and good repeatability, indicating that the DIC technique, combined with the small-disc splitting tensile test in ASTM D8289, can capture the strain evolution in the test. In conclusion, further comparison of measured strains with analytical and numerical simulation results provides additional understanding of the mechanical behaviour of this nuclear graphite material in complex stress states.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Upscaling Reactive Transport and Clogging in Shale Microcracks by Deep Learning

Fracture networks in shales exhibit multiscale features. A rock system may contain a few main fractures and thousands of microcracks, whose length and aperture are orders of magnitude smaller than the former. It is computationally prohibitive to resolve all the fractures explicitly for such multiscale fracture networks. One traditional approach is to model the small-scale features (e.g., microcracks in shales) as an effective medium. Although this fracture-matrix conceptualization significantly reduces the problem complexity, there are classes of physical processes that cannot be accurately upscaled by effective medium approximations, e.g., microcrack clogging during mineral reactions. Here, we employ deep learning in place of effective medium theory to upscale physical processes in small-scale features. Specifically, we consider reactive transport in a fracture-microcrack network where microcracks can be clogged by precipitation. A deep learning multiscale algorithm is developed, in which the microcracks are upscaled as a wall boundary condition of the main fractures. The wall boundary condition is constructed by recurrent neural networks, which take concentration histories as input and predict the solute transport from main fractures to microcracks. The deep learning multiscale algorithm is firstly employed in specific scenarios, then a general model is developed which can work under various conditions. The new approach is validated against fully resolved simulations and an analytical solution, providing a reliable and efficient solution for problems that cannot be upscaled by effective medium models.

58 GEOSCIENCES↗

DOE OSTI2023

Measuring the giant radio galaxy length distribution with the LoTSS

Many massive galaxies launch jets from the accretion disk of their central black hole, but only ~10 3 instances are known in which the associated outflows form giant radio galaxies (GRGs, or giants): luminous structures of megaparsec extent that consist of atomic nuclei, relativistic electrons, and magnetic fields. Large samples are imperative to understanding the enigmatic growth of giants, and recent systematic searches in homogeneous surveys constitute a promising development. For the first time, it is possible to perform meaningful precision statistics with GRG lengths, but a framework to do so is missing. We measured the intrinsic GRG length distribution by combining a novel statistical framework with a LOFAR Two-metre Sky Survey (LoTSS) sample of freshly discovered giants. In turn, this allowed us to answer an array of questions on giants. For example, we can now assess how rare a 5 Mpc giant is compared with one of 1 Mpc, and how much larger – given a projected length – the corresponding intrinsic length is expected to be. Notably, we can now also infer the GRG number density in the Local Universe. We assumed the intrinsic GRG length distribution to be Paretian (i.e. of power-law form) with tail index ξ, and predicted the observed distribution by modelling projection and selection effects. To infer ξ, we also systematically searched the LoTSS for hitherto unknown giants and compiled the largest catalogue of giants to date. We show that if intrinsic GRG lengths are Pareto distributed with index ξ, then projected GRG lengths are also Pareto distributed with index ξ. Selection effects induce curvature in the observed projected GRG length distribution: angular length selection flattens it towards the lower end, while surface brightness selection steepens it towards the higher end. We explicitly derived a GRG’s posterior over intrinsic lengths given its projected length, laying bare the ξ dependence. We also discovered 2060 giants within LoTSS DR2 pipeline products; our sample more than doubles the known population. Spectacular discoveries include the largest, second-largest, and fourth-largest GRG known (l p = 5.1 Mpc, l p = 5.0 Mpc, and l p = 4.8 Mpc), the largest GRG known hosted by a spiral galaxy (l p = 2.5 Mpc), and the largest secure GRG known beyond redshift 1 (l p = 3.9 Mpc). We increase the number of known giants whose angular length exceeds that of the Moon from 10 to 23; among the discoveries is the angularly largest known radio galaxy in the Northern Sky, which is also the angularly largest known GRG ($\phi$ = 2°). Combining theory and data, we determined that intrinsic GRG lengths are well described by a Pareto distribution, and measured the index ξ = –3.5 ± 0.5. This implies that, given its projected length, a GRG’s intrinsic length is expected to be just 15% larger. Finally, we determined the comoving number density of giants in the Local Universe to be n GRG = 5 ± 2(100 Mpc) –3 . We developed a practical mathematical framework that elucidates the statistics of giant radio galaxy lengths. Through a LoTSS search, we also discovered 2060 new giants. By combining both advances, we determined that intrinsic GRG lengths are well described by a Pareto distribution with index ξ = –3.5 ± 0.5, and that giants are truly rare in a cosmological sense: most clusters and filaments of the Cosmic Web are not currently home to a giant. Thus, our work yields new observational constraints for analytical models and simulations featuring radio galaxy growth.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

BEYONDPLANCK X. Planck Low Frequency Instrument frequency maps with sample-based error propagation

We present Planck Low Frequency Instrument (LFI) frequency sky maps derived within the BEYONDPLANCK framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical, and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples, each of which corresponds to one possible realization of the various modeled instrumental systematic corrections, including correlated noise, time-variable gain, as well as far sidelobe and bandpass corrections. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncertainties into higher-level science products, including astrophysical component maps, angular power spectra, and cosmological parameters. We show that the two dominant sources of LFI instrumental systematic uncertainties are correlated noise and gain fluctuations, and the products presented here support – for the first time – full Bayesian error propagation for these effects at full angular resolution. We compared our posterior mean maps with traditional frequency maps delivered by the Planck Collaboration, and find generally good agreement. The most important quality improvement is due to significantly lower calibration uncertainties in the new processing, as we find a fractional absolute calibration uncertainty at 70 GHz of Δg 0 /g 0 = 5 x 10 -5 , which is nominally 40 times smaller than that reported by Planck 2018. However, we also note that the original Planck 2018 estimate has a nontrivial statistical interpretation, and this further illustrates the advantage of the new framework in terms of producing self-consistent and well-defined error estimates of all involved quantities without the need of ad hoc uncertainty contributions. We describe how low-resolution data products, including dense pixel-pixel covariance matrices, may be produced from the posterior samples directly, without the need for computationally expensive analytic calculations or simulations. We conclude that posterior-based frequency map sampling provides unique capabilities in terms of low-level systematics modeling and error propagation, and may play an important role for future Cosmic Microwave Background (CMB) B-mode experiments aiming at nanokelvin precision.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Deep learning surrogate model for kinetic Landau-fluid closure with collision

In this work, the kinetic Landau-fluid (LF) closure with collision and periodic boundary condition is used in the development of the deep learning (DL) surrogate model. A classical neural network, namely, feedforward neural network or sometimes termed multilayer perceptron, is constructed and trained to learn the kinetic LF closure in the static limit and arbitrary mean free path in configuration space. The preliminary relation between best hyperparameters and critical parameters for data generation is found. Compared with the numerical approach (non-Fourier method) of the LF closure, the deep learning surrogate model shows an order of magnitude of improvement in terms of accuracy. Perhaps most importantly, the surrogate model closure has been integrated for the first time with fluid simulations. Our DL-enabled fluid simulations, for the first time, give the correct Landau damping rate for a wide range of wave vectors, while the Hammett–Perkins closure cannot produce the correct damping rate. We correctly connect the collisionless Hammett–Perkins closure and collisional Braginskii closure to reproduce the intrinsic nonlocal feature of the heat flux with DL techniques. We address the most concerning error accumulation problem and find that simulations with the deep learning surrogate model are as good as, if not better than, simulations with the analytic closure in terms of long-term numerical stability in the linear Landau damping test.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Thickness dependence of piezo-bimorph adaptive mirror bending

A new generation of adaptive x-ray optics (AXO) is being installed on high-coherent-flux x-ray beamlines worldwide to correct and control the optical wavefront with sub-nm precision. These ultra-smooth mirrors achieve high reflectivities at glancing angles of incidence and can be hundreds of mm long. One type of adaptive x-ray mirror relies on piezoelectric ceramic strips which are segmented into channels and actuated to induce local, longitudinal bending, generating one-dimensional shape changes in the mirror substrate. A recently described mirror model uses a three-layer geometry with parallel actuators on the front and back surfaces of a thicker mirror substrate. By analogy to a solved problem in the thermal actuation of a tri-metal strip, we show that the achievable bending radius varies approximately as the square of the substrate thickness. We provide an analytic solution and simulate bending using a finite-element model.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

What next: Further implosion space exploration on the path to NIF extended yield capability

We present quantitative motivations and assessments of various proposed and ongoing directions to further improving yields and target gain of igniting indirect-drive implosions at the National Ignition Facility (NIF). These include increasing compression and confinement time, improving hohlraum and ablator efficiency, and further increasing peak power and laser energy. 1D hydroscaled simulations, augmented by analytic 1D theory, have been used to project yield improvements for each of these implosion optimization tracks, normalized to the best current performing 4 MJ shot. At current NIF capabilities of 2.2 MJ, 450 TW, we project several paths could reach 15 MJ yield levels. We also expect several key implosion physics questions will be addressed in attempting to reach this yield level. These include demonstrating to what extent lower adiabat designs leading to higher compression will increase gain and efficiency, and whether we can reduce residual kinetic energy and ablator-fuel mix that is probably limiting the current burn-up fraction. For an envisaged NIF upgrade to EL = 3 MJ at fixed 450 TW peak power, scaling capsule size and fuel thicknesses faster than pure hydroscaling should allow for yields that could reach up to 60–80 MJ, depending on the efficiency gains realized in increasing deuterium-tritium fuel thickness, reducing hohlraum losses, and switching to lower Z ablators. The laser-plasma instability and beam transmission scaling in these larger hohlraums is shown to be favorable if the spot size is increased with hohlraum scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗

DOE OSTI2020

A Lagrangian perturbation theory in the presence of massive neutrinos

We develop a Lagrangian Perturbation Theory (LPT) framework to study the clustering of cold dark matter (CDM) in cosmologies with massive neutrinos. We follow the trajectories of CDM particles with Lagrangian displacements fields up to third order in perturbation theory. Once the neutrinos become non-relativistic, their density fluctuations are modeled as being proportional to the CDM density fluctuations, with a scale-dependent proportionality factor. This yields a gravitational back-reaction that introduces additional scales to the linear growth function, which is accounted for in the higher order LPT kernels. Through non-linear mappings from Eulerian to Lagrangian frames, we ensure that our theory has a well behaved large scale behavior free of unwanted UV divergences, which are common when neutrino and CDM densities are not treated on an equal footing, and in resummation schemes that manifestly break Galilean invariance. We use our theory to construct correlation functions for both the underlying matter field, as well as for biased tracers using ConvolutionLPT. Redshift-space distortions effects are modeled using the Gaussian Streaming Model. When comparing our analytical results to simulated data from the Quijote 1 simulation suite, we find good accuracy down to r = 20 Mpc h -1 at redshift z = 0.5, for the real space and redshift space monopole particle correlation functions with no free parameters. The same accuracy is reached for the redshift space quadrupole if we additionally consider an effective field theory parameter that shifts the pairwise velocity dispersion. For modeling the correlation functions of tracers we adopt a simple Lagrangian biasing scheme with only density and curvature operators, which we find sufficient to reach down to r = 20 Mpc h -1 when comparing to simulated halos.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Diluted axion star collisions with neutron stars

Diluted axion star, a self-gravitating object with the quantum pressure balancing gravity, has been predicted in many models with a QCD axion or axion-like particle. It can be formed in the early universe and composes a sizable fraction of dark matter. One could detect the transient radio signals when it passes by a magnetar with the axion particle converted into photon in the magnetic field. Using both numerical and semi-analytic approaches, we simulate the axion star's dynamic evolution and estimate the fraction of axion particles that can have a resonance conversion during such a collision event. We have found that both self-gravity and quantum pressure are not important after the diluted axion star enters the Roche radius. A free-fall approximate can capture individual particle trajectories very well. Finally, with some optimistic cosmological and astrophysical assumptions, the QCD axion parameter space can be probed from detecting such a collision event by radio telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗