Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Lattice simulations of inflation

The scalar field theory of cosmological inflation constitutes nowadays one of the preferred scenarios for the physics of the early universe. In this paper we aim at studying the inflationary universe making use of a numerical lattice simulation. Various lattice codes have been written in the last decades and have been extensively used for understating the reheating phase of the universe, but they have never been used to study the inflationary phase itself far from the end of inflation (i.e. about 50 e-folds before the end of inflation). In this paper we use a lattice simulation to reproduce the well-known results of some simple models of single-field inflation, particularly for the scalar field perturbation. The main model that we consider is the standard slow-roll inflation with an harmonic potential for the inflaton field. We explore the technical aspects that need to be accounted for in order to reproduce with precision the nearly scale invariant power spectrum of inflaton perturbations. We also consider the case of a step potential, and show that the simulation is able to correctly reproduce the oscillatory features in the power spectrum of this model. Even if a lattice simulation is not needed in these cases, that are well within the regime of validity of linear perturbation theory, this sets the basis to future work on using lattice simulations to study more complicated models of inflation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1977

A real-time subsynoptic, meso, and microscale severe storm forecast system

A system of differential equations is integrated numerically in space and time over several different matrices in an effort to simulate the atmospheric wave structures which organize severe local storms. Preliminary results with case studies indicate that dynamical fields produced by the numerical simulations can be translated into very fine scale space and time zones where severe storm forecast indices can be developed. These fine scale indices are now available in real-time when run on Langley's STAR 100 computer system.

Kaplan, M. L.↗
DOE OSTI2021

Grid-Forming Frequency Shaping Control for Low-Inertia Power Systems

As power systems transit to a state of high renewable penetration, little or no presence of synchronous generators makes the prerequisite of well-regulated frequency for grid-following inverters unrealistic. Thus, there is a trend to resort to grid-forming inverters which set frequency directly. We propose a novel grid-forming frequency shaping control that is able to shape the aggregate system frequency dynamics into a first-order one with the desired steady-state frequency deviation and Rate of Change of Frequency (RoCoF) after a sudden power imbalance. The no overshoot property resulting from the first-order dynamics allows the system frequency to monotonically move towards its new steady-state without experiencing frequency Nadir, which largely improves frequency security. We prove that our grid-forming frequency-shaping control renders the system internally stable under mild assumptions. The performance of the proposed control is verified via numerical simulations on a modified Icelandic Power Network test case.

control systems↗
NASA NTRS1991

Model studies on the role of moist convection as a mechanism for interaction between the mesoscales

A three year research effort is described which had as its goal the development of techniques to improve the numerical prediction of cumulus convection on the meso-beta and meso-gamma scales. Two MESO models are used, the MASS (mesoscale) and TASS (cloud scale) models. The primary meteorological situation studied is the 28-29 Jun. 1986 Cooperative Huntsville Meteorological Experiment (COHMEX) study area on a day with relatively weak large scale forcing. The problem of determining where and when convection should be initiated is considered to be a major problem of current approaches. Assimilation of moisture data from satellite, radar, and surface data is shown to significantly improve mesoscale simulations. The TASS model is shown to reproduce some observed mesoscale features when initialized with 3-D observational data. Convection evolution studies center on comparison of the Kuo and Fritsch-Chappell cumulus parameterization schemes to each other, and to cloud model results. The Fritsch-Chappell scheme is found to be superior at about 30 km resolution, while the Kuo scheme does surprisingly well in simulating convection down to 10 km in cases where convergence features are well-resolved by the model grid. Results from MASS-TASS interaction experiments are presented and discussed. A discussion of the future of convective simulation is given, with the conclusion that significant progress is possible on several fronts in the next few years.

Waight, Kenneth T., III↗
NASA NTRS1976

Monte Carlo calculations of diatomic molecule gas flows including rotational mode excitation

The direct simulation Monte Carlo method was used to solve the Boltzmann equation for flows of an internally excited nonequilibrium gas, namely, of rotationally excited homonuclear diatomic nitrogen. The semi-classical transition probability model of Itikawa was investigated for its ability to simulate flow fields far from equilibrium. The behavior of diatomic nitrogen was examined for several different nonequilibrium initial states that are subjected to uniform mean flow without boundary interactions. A sample of 1000 model molecules was observed as the gas relaxed to a steady state starting from three specified initial states. The initial states considered are: (1) complete equilibrium, (2) nonequilibrium, equipartition (all rotational energy states are assigned the mean energy level obtained at equilibrium with a Boltzmann distribution at the translational temperature), and (3) nonequipartition (the mean rotational energy is different from the equilibrium mean value with respect to the translational energy states). In all cases investigated the present model satisfactorily simulated the principal features of the relaxation effects in nonequilibrium flow of diatomic molecules.

Yoshikawa, K. K.↗
DOE OSTI2020

Bipolar planetary nebulae from outflow collimation by common envelope evolution

ABSTRACT The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and the dense toroidal-shaped ejecta left over from common envelope (CE) evolution. Here we use the 3D hydrodynamic adaptive mesh refinement (AMR) code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the smoothed particle hydrodynamics (SPH) code phantom serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta pre-PN outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post-CE winds may be aspherical, our models show that collimation via ‘inertial confinement’ will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock-focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC-driven bipolar lobes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Balancing Impedance and Controllability in Response Reconstruction with $\mathrm{TS-IMMAT}$

In a prior work the authors proposed a variant on the Impedance Matched Multi-Axis Test (IMMAT) in which a fixture is defined, called the Transmission Simulator (TS), and the desired environment is matched at a set of sensors on the TS. If the motion of the TS is matched then the response of the rest of the component will also match, provided that the attached component has the same dynamics as it did when the environment was measured. Hence, one would like the TS to be flexible so that it reproduces the boundary conditions that the component of interest experiences during flight, but the more flexible the TS, the more shakers might be needed to control its response. This work presents a derivation that gives expressions for these two potential error sources in TS-IMMAT. Then, various case studies are presented, both on simulated and real hardware, to understand the importance of each error term in practical testing. The theory explains the phenomena that were observed when using measurements from a component that flew on a sounding rocket. The environmental response was measured and then various fixtures were attached, each comprising more of the next assembly, or the hardware to which the component was attached in flight. MIMO testing was repeated with each fixture and the results were compared to seek to understand the role of the impedance match in this type of testing. Here the results show that the number of modes that are active in the transmission simulator is also very important, and so the best solution balances these two considerations. An improved method of simulating the MIMO test is then proposed, so simulations can be used to predict what fixture, or transmission simulator, will give the best results in a TS-IMMAT test.

42 ENGINEERING↗
DOE OSTI2023

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗
NASA NTRS1983

Precise thermal control test demonstration on simulated space telescope main ring

A simulated section of the main ring of the Space Telescope, (ST) the part that holds the primary mirror, was manufactured and tested for thermal stability characteristics in a complex structural assembly. The other telescope subassemblies are also mounted on the main ring, which is constructed of a titanium alloy. The entire ring is wrapped in heater elements around the outside diameter to reduce the external temperature gradient effects by two orders of magnitude. A thermal vacuum test was performed using a stainless steel alloy tailored to represent the titanium ring's thermal characteristics. The tests examined the thermal performance and controllability of the test ring, the effects of a 24-hr slew manuever that produced worst-case environment changes, hot and cold simulation of predicted ST orbital thermal environments, and the effects of heater failure. Pointing and thermal control features were verified, but a heater failure on the 'on' side was found to quickly lead to exceeding the thermal design limits.

Bettini, R. G.↗
NASA NTRSDate not supplied

Development of an Airspace Simulation and Modeling Tool for Enhanced Spectrum Management

The emergence of new aerial vehicles into the National Airspace System creates an increased demand for aeronautical communications to support aviation operations. However, the issue of spectrum scarcity remains an ever-present concern, and the growing demand cannot be supported using existing spectrum allocation strategies. As a result, a new spectrum management approach is required, and the National Aeronautics and Space Administration (NASA) is investigating advanced concepts to modernize the management and use of aviation spectrum by leveraging the latest advancements in wireless communications, big data and machine learning. This research proposes an autonomous spectrum allocation concept, which allocates communications resources, such as spectrum and power, based on the predicted communications and air traffic demands throughout the airspace, as opposed to the use of fixed allocations as is done today. This approach will result in improved spectrum utilization efficiency and enhanced airspace capacity. The autonomous spectrum allocation concept decomposes into three research areas: demand prediction, resource allocation, and use case evaluation. As part of the use case evaluation effort, a modeling and simulation capability is currently under development. This simulation capability includes the implementation of various features, including visualization of both live or virtually-generated airspace traffic, simulation scenario development, simulation management with data collection, and flight plan creation with corresponding trajectory generation. This modeling and simulation capability will continue to evolve as new and advanced airspace applications are introduced into existing and emerging operational environments.

Eric J. Knoblock↗
NASA NTRS1993

Accuracy requirements

Satellite and surface measurements, if they are to serve as a climate monitoring system, must be accurate enough to permit detection of changes of climate parameters on decadal time scales. The accuracy requirements are difficult to define a priori since they depend on unknown future changes of climate forcings and feedbacks. As a framework for evaluation of candidate Climsat instruments and orbits, we estimate the accuracies that would be needed to measure changes expected over two decades based on theoretical considerations including GCM simulations and on observational evidence in cases where data are available for rates of change. One major climate forcing known with reasonable accuracy is that caused by the anthropogenic homogeneously mixed greenhouse gases (CO2, CFC's, CH4 and N2O). Their net forcing since the industrial revolution began is about 2 W/sq m and it is presently increasing at a rate of about 1 W/sq m per 20 years. Thus for a competing forcing or feedback to be important, it needs to be of the order of 0.25 W/sq m or larger on this time scale. The significance of most climate feedbacks depends on their sensitivity to temperature change. Therefore we begin with an estimate of decadal temperature change. Presented are the transient temperature trends simulated by the GISS GCM when subjected to various scenarios of trace gas concentration increases. Scenario B, which represents the most plausible near-term emission rates and includes intermittent forcing by volcanic aerosols, yields a global mean surface air temperature increase Delta Ts = 0.7 degrees C over the time period 1995-2015. This is consistent with the IPCC projection of about 0.3 degrees C/decade global warming (IPCC, 1990). Several of our estimates below are based on this assumed rate of warming.

Delgenio, Anthony↗
DOE OSTI2020

Data-driven selection of stiff chemistry ODE solver in operator-splitting schemes

Most computational fluid dynamics simulations of practical combustion applications employ operator-splitting schemes, where chemistry and transport are separated and integrated with distinct numerical methods. The changes in composition due to chemistry are evaluated by solving ordinary differential equations (ODE) in each cell of the computational domain, which typically dominates the computational cost when detailed chemistry is considered. In this work, a data-driven approach for the selection of chemistry ODE solvers in operator-splitting schemes is presented. Neural networks are used to predict the ODE solvers CPU times and errors for a given thermochemical state. This allows the selection of an optimal ODE solver on a cell-by-cell, timestep-by-timestep basis. The models are trained using a wide set of thermochemical states generated through partially-stirred reactors and flames simulations. The methodology is validated by quantifying the prediction errors, the classification accuracy, and the computational speedup. The model predicts the optimal ODE solver for 70 to 95% of the validation cases and decreases the computional cost by a factor of 3 or more. The generalizability of the methodology to different chemical mechanisms and different fuels is assessed and it is shown that the model’s performance is only slightly degraded and its applicability is significantly enhanced if the inputs to the neural networks are restricted to a small set of thermochemical state variables present in most chemical mechanisms. In conclusion, the models are used in an homogeneous reactor case and a multi-dimensional CFD simulation of a diesel spray at high pressure where a speedup of more than 3 is achieved.

42 ENGINEERING↗
NASA NTRS2011

Crush Test Abuse Stand

The purpose of this system is to simulate an internal short on battery cells by causing deformation (a crushing force) in a cell without penetration. This is performed by activating a hydraulic cylinder on one side of a blast wall with a hydraulic pump located on the other. The operator can control the rate of the crush by monitoring a local pressure gauge connected to the hydraulic cylinder or a load cell digital display located at the hydraulic pump control area. The internal short simulated would be considered a worst-case scenario of a manufacturer fs defect. This is a catastrophic failure of a cell and could be a very destructive event. Fully charged cells are to have an internal short simulated at the center of the length of the cell (away from terminals). The crush can be performed with a .- to 1-in. (.0.6- to 2.5-cm) rod placed crossways to the cell axis, causing deformation of the cell without penetration. The OCV (open-circuit voltage) and temperature of the cells, as well as the pressure and crushing force, are recorded during the operation. Occurrence of an internal short accompanied by any visible physical changes such as venting, fires, or explosions is reported. Typical analytical data examined after the test would be plots of voltage, temperature, and pressure or force versus time. The rate of crushing force can be increased or decreased based on how fast the operator pumps the hydraulic pump. The size of cylinder used to compress the battery cell can be easily changed by adding larger or smaller fittings onto the end of the hydraulic cylinder based on the battery/cell size being tested. The cell is crushed remotely and videotaped, allowing the operator to closely monitor the situation from a safe distance.

Collins, Jacob↗
NASA NTRS1991

Numerical study of axial turbulent flow over long cylinders

The effects of transverse curvature are investigated by means of direct numerical simulations of turbulent axial flow over cylinders. Two cases of Reynolds number of about 3400 and layer-thickness-to-cylinder-radius ratios of 5 and 11 were simulated. All essential turbulence scales were resolved in both calculations, and a large number of turbulence statistics were computed. The results are compared with the plane channel results of Kim et al. (1987) and with experiments. With transverse curvature the skin friction coefficient increases and the turbulence statistics, when scaled with wall units, are lower than in the plane channel. The momentum equation provides a scaling that collapses the cylinder statistics, and allows the results to be interpreted in light of the plane channel flow. The azimuthal and radial length scales of the structures in the flow are of the order of the cylinder diameter. Boomerang-shaped structures with large spanwise length scales were observed in the flow.

Neves, J. C.↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
NASA NTRS2022

The Effect of Solar Wind Turbulence on Parallel and Oblique Firehose Instabilities

We consider the firehose instability coexisting with the omnipresent ambient solar wind turbulence. The characteristic temporal and spatial scales of the turbulence are comparable to those of the instability. Therefore, turbulence may violate the common assumption of a uniform and stationary background used to describe instabilities and make the properties of the instabilities different. To investigate this effect, we perform threedimensional hybrid simulations with particle-in-cell ions and a quasi-neutralizing electron fluid. We find that the turbulence significantly reduces the growth rates and saturation levels of both instabilities. Comparing the cases with and without turbulence, the former results in a higher temperature anisotropy in the asymptotic marginally stable state at large times. In the former case, the distribution function averaged over the simulation box is also closer to the initial one.

S A Markovskii↗
DOE OSTI2020

Dark Lock-in Thermography Identifies Solder Bond Failure as the Root Cause of Series Resistance Increase in Fielded Solar Modules

Correlating the electrical performance of photovoltaic modules to the spatially resolved photoluminescence, electroluminescence, and dark lock-in thermography (DLIT) images is an important long-term goal for developing solar cell technology. These images offer highly sophisticated and detailed information about the spatial distribution and (if imaged at successive time intervals) temporal degradation of local series and shunt resistances. There have been extensive studies to correlate these imaging techniques to local characteristics at the cell level. However, it has been difficult to extract and quantify module-level information from the techniques. In this study, we interpret module current-voltage (I-V) measurements along with corresponding DLIT images within a module-scale simulation framework, demonstrating that series resistance degradation of the module I-V characteristics, in this case, can be attributed to solder bond failures. Our simulations highlight how current crowding associated with a failed solder joint (or a section of the solder pad) translates to the characteristic point-like (asymmetric doublet) heating pattern in neighboring solder joints (or the neighboring regions of the solder pad). Overall, the correlation between series resistance and solder joint degradation could inform expedited diagnosis of field degradation of solar modules.

42 ENGINEERING↗
NASA NTRS2004

Theory and Simulation of Real and Ideal Magnetohydrodynamic Turbulence

Incompressible, homogeneous magnetohydrodynamic (MHD) turbulence consists of fluctuating vorticity and magnetic fields, which are represented in terms of their Fourier coefficients. Here, a set of five Fourier spectral transform method numerical simulations of two-dimensional (2-D) MHD turbulence on a 512(sup 2) grid is described. Each simulation is a numerically realized dynamical system consisting of Fourier modes associated with wave vectors k, with integer components, such that k = |k| less than or equal to k(sub max). The simulation set consists of one ideal (non-dissipative) case and four real (dissipative) cases. All five runs had equivalent initial conditions. The dimensions of the dynamical systems associated with these cases are the numbers of independent real and imaginary parts of the Fourier modes. The ideal simulation has a dimension of 366104, while each real simulation has a dimension of 411712. The real runs vary in magnetic Prandtl number P(sub M), with P(sub M) is a member of {0.1, 0.25, 1, 4}. In the results presented here, all runs have been taken to a simulation time of t = 25. Although ideal and real Fourier spectra are quite different at high k, they are similar at low values of k. Their low k behavior indicates the existence of broken symmetry and coherent structure in real MHD turbulence, similar to what exists in ideal MHD turbulence. The value of PM strongly affects the ratio of kinetic to magnetic energy and energy dissipation (which is mostly ohmic). The relevance of these results to 3-D Navier-Stokes and MHD turbulence is discussed.

Shebalin, John V.↗