Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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 OSTI2019

Continuum Modeling of Inductor Magnetic Hysteresis and Eddy Currents in Resonant Circuits

This paper presents a high-fidelity finite-element modeling technique for magnetic hysteresis and eddy current losses in toroid inductors. The method is based on the separation of ferromagnetic loss characteristics into two components: a quasi-static hysteresis component and a dynamic eddy current component. The Preisach model is used to describe the quasi-static magnetic hysteresis behavior of the core, providing strong guarantees on the reproducibility of the experimentally measured characteristics. This model is used to represent the magnetic field constitutive relationships within a finite-element framework combining the effects of hysteresis and eddy currents in a unified dynamic simulation. The finite-element model of the toroid is used as a high-order inductor model coupled to a resonant circuit simulation. The modeling technique is validated through experimental measurements on two different series RLC circuits. The first circuit is based on an M19 electrical steel toroid having resonant frequency near 200 Hz. The second circuit is based on a T38 ferrite toroid having a resonant frequency near 10 kHz. The models agree closely with the measured voltages, currents, and losses. The models also successfully predict discontinuities in the measured frequency responses due to the existence of bistable operating regimes.

42 ENGINEERING↗
NASA NTRS2011

Pattern Recognition for a Flight Dynamics Monte Carlo Simulation

The design, analysis, and verification and validation of a spacecraft relies heavily on Monte Carlo simulations. Modern computational techniques are able to generate large amounts of Monte Carlo data but flight dynamics engineers lack the time and resources to analyze it all. The growing amounts of data combined with the diminished available time of engineers motivates the need to automate the analysis process. Pattern recognition algorithms are an innovative way of analyzing flight dynamics data efficiently. They can search large data sets for specific patterns and highlight critical variables so analysts can focus their analysis efforts. This work combines a few tractable pattern recognition algorithms with basic flight dynamics concepts to build a practical analysis tool for Monte Carlo simulations. Current results show that this tool can quickly and automatically identify individual design parameters, and most importantly, specific combinations of parameters that should be avoided in order to prevent specific system failures. The current version uses a kernel density estimation algorithm and a sequential feature selection algorithm combined with a k-nearest neighbor classifier to find and rank important design parameters. This provides an increased level of confidence in the analysis and saves a significant amount of time.

Restrepo, Carolina↗
NASA NTRS1981

Differential absorption lidar measurements of atmospheric temperature and pressure profiles

The theory and methodology of using differential absorption lidar techniques for the remote measurement of atmospheric pressure profiles, surface pressure, and temperature profiles from ground, air, and space-based platforms are presented. Pressure measurements are effected by means of high resolution measurement of absorption at the edges of the oxygen A band lines where absorption is pressure dependent due to collisional line broadening. Temperature is assessed using measurements of the absorption at the center of the oxygen A band line originating from a quantum state with high ground state energy. The population of the state is temperature dependent, allowing determination of the temperature through the Boltzmann term. The results of simulations of the techniques using Voigt profile and variational analysis are reported for ground-based, airborne, and Shuttle-based systems. Accuracies in the 0.5-1.0 K and 0.1-0.3% range are projected.

Korb, C. L.↗
DOE OSTI2018

Stellarator microinstabilities and turbulence at low magnetic shear

Gyrokinetic simulations of drift waves in low-magnetic-shear stellarators reveal that simulation domains comprised of multiple turns can be required to properly resolve critical mode structures important in saturation dynamics. Marginally stable eigenmodes important in saturation of ion temperature gradient modes and trapped electron modes in the Helically Symmetric Experiment (HSX) stellarator are observed to have two scales, with the envelope scale determined by the properties of the local magnetic shear and an inner scale determined by the interplay between the local shear and magnetic field-line curvature. Properly resolving these modes removes spurious growth rates that arise for extended modes in zero-magnetic-shear approximations, enabling use of a zero-magnetic-shear technique with smaller simulation domains and attendant cost savings. Analysis of subdominant modes in trapped electron mode (TEM)-driven turbulence reveals that the extended marginally stable modes play an important role in the nonlinear dynamics, and suggests that the properties induced by low magnetic shear may be exploited to provide another route for turbulence saturation.

Faber, B. J. (ORCID:000000034934400X)↗
DOE OSTI2021

Parametrized Hamiltonian simulation using quantum optimal control

Analog quantum simulation offers a hardware-specific approach to studying quantum dynamics, but mapping a model Hamiltonian onto the available device parameters requires matching the hardware dynamics. In this work, we introduce a paradigm for quantum Hamiltonian simulation that leverages digital decomposition techniques and optimal control to perform analog simulation. We validate this approach by constructing the optimal analog controls for a superconducting transmon device to emulate the dynamics of an extended Bose-Hubbard model. We demonstrate the role of control time, digital error, and pulse complexity, and we explore the accuracy and robustness of these controls. We conclude by discussing the opportunity for implementing this paradigm in near-term quantum devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

An evaluation of thematic mapper simulator data for mapping forest cover

Computer-aided analysis techniques applied to Thematic Mapper Simulator (TMS) data were evaluated for the purpose of mapping forest cover types. Classification results obtained using a supervised set of training statistics and various combinations of three and four channel subsets of the seven available TMS channels are compared for the L2 (Minimum Euclidean Distance), GML (Gaussian Maximum Likelihood), and SECHO (Supervised Extraction and Classification of Homogeneous Objects) classification algorithms. SECHO performed significantly better than either of the two per-point classifiers for the untransformed data. Overall classification results of the Karhunen-Loeve transformation increased for the L2 algorithm, but decreased for both the GML and SECHO algorithms.

Dean, M. E.↗
NASA NTRS1989

Man-in-the-control-loop simulation of manipulators

A method to achieve man-in-the-control-loop simulation is presented. Emerging real-time dynamics simulation suggests a potential for creating an interactive design workstation with a human operator in the control loop. The recursive formulation for multibody dynamics simulation is studied to determine requirements for man-in-the-control-loop simulation. High speed computer graphics techniques provides realistic visual cues for the simulator. Backhoe and robot arm simulations are implemented to demonstrate the capability of man-in-the-control-loop simulation.

Chang, J. L.↗
NASA NTRS2003

Simulation of Combustion Systems with Realistic g-jitter

In this project a transient, fully three-dimensional computer simulation code was developed to simulate the effects of realistic g-jitter on a number of combustion systems. The simulation code is capable of simulating flame spread on a solid and nonpremixed or premixed gaseous combustion in nonturbulent flow with simple combustion models. Simple combustion models were used to preserve computational efficiency since this is meant to be an engineering code. Also, the use of sophisticated turbulence models was not pursued (a simple Smagorinsky type model can be implemented if deemed appropriate) because if flow velocities are large enough for turbulence to develop in a reduced gravity combustion scenario it is unlikely that g-jitter disturbances (in NASA's reduced gravity facilities) will play an important role in the flame dynamics. Acceleration disturbances of realistic orientation, magnitude, and time dependence can be easily included in the simulation. The simulation algorithm was based on techniques used in an existing large eddy simulation code which has successfully simulated fire dynamics in complex domains. A series of simulations with measured and predicted acceleration disturbances on the International Space Station (ISS) are presented. The results of this series of simulations suggested a passive isolation system and appropriate scheduling of crew activity would provide a sufficiently "quiet" acceleration environment for spherical diffusion flames.

Mell, William E.↗
NASA NTRS2017

Order-Reduced Solution of the Nonlinear High-Fidelity Generalized Method of Cells Micromechanics Relations

The High-Fidelity Generalized Method of Cells (HFGMC) is one technique for accurately simulating nonlinear composite material behavior. The HFGMC uses a higher-order approximation for the subcell displacement field that allows for a more accurate determination of the subcell stressstrain fields at the cost of some computational efficiency. In order to reduce computational costs associated with the solution of the ensuing system of simultaneous equations, the HFGMC global system of equations for doubly-periodic repeating unit cells with nonlinear constituents was reduced in size through the use of a Petrov-Galerkin-based Proper Orthogonal Decomposition order-reduction scheme. A number of cases were presented that address the computational feasibility of using order-reduction techniques to solve solid mechanics problems involving complex microstructures.

order-reduction↗
NASA NTRS1990

Climate and atmospheric modeling studies. Climate applications of Earth and planetary observations. Chemistry of Earth and environment

The research conducted during the past year in the climate and atmospheric modeling programs concentrated on the development of appropriate atmospheric and upper ocean models, and preliminary applications of these models. Principal models are a one-dimensional radiative-convective model, a three-dimensional global climate model, and an upper ocean model. Principal applications have been the study of the impact of CO2, aerosols and the solar 'constant' on climate. Progress was made in the 3-D model development towards physically realistic treatment of these processes. In particular, a map of soil classifications on 1 degree x 1 degree resolution has been digitized, and soil properties have been assigned to each soil type. Using this information about soil properties, a method was developed to simulate the hydraulic behavior of soils of the world. This improved treatment of soil hydrology, together with the seasonally varying vegetation cover, will provide a more realistic study of the role of the terrestrial biota in climate change. A new version of the climate model was created which follows the isotopes of water and sources of water (or colored water) throughout the planet. Each isotope or colored water source is a fraction of the climate model's water. It participates in condensation and surface evaporation at different fractionation rates and is transported by the dynamics. A major benefit of this project has been to improve the programming techniques and physical simulation of the water vapor budget of the climate model.

Source record↗
DOE OSTI2022

Insight into Pyroelectricity and Phase Transitions in Ferroelectrics from Nonequilibrium Approach: The Case of PbTiO 3

Computational methodology is proposed to simulate time dependent temperature variations via thermostat approach. Here, the methodology is applied to study pyroelectricity under nonequilibrium conditions in prototypical ferroelectric PbTiO 3 . It is found that at room temperature the nonequilibrium effects in pyroelectricity begin manifesting above 1.0 THz and originate from the overlap with the intrinsic soft mode frequency. Therefore, the proposed computational methodology can be used to predict equilibrium pyroelectric response by choosing a frequency of temperature variation below the intrinsic polar mode frequencies for the material. Probing the dynamics of ferroelectric phase transition with the proposed methodology revealed that it takes about a nanosecond for the computational supercell to undergo phase transition. It is concluded that computational techniques capable of simulating picoseconds can be used to compute pyroelectric response from nonequilibrium dynamical approach while those that can reach nanoseconds are suitable for simulating phase transitions in ferroelectrics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The DESI N -body Simulation Project – II. Suppressing sample variance with fast simulations

ABSTRACT Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches $\sim 20\, h^{-3}\, \mathrm{Gpc}^{3}$. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. AbacusSummit is a suite of high-resolution dark-matter-only simulations designed for this purpose, with $200\, h^{-3}\, \mathrm{Gpc}^{3}$ (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-N-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. $k \lt 0.3\, h\, \mathrm{Mpc}^{-1}$ for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in AbacusSummit suite to extend the effective volume ∼20 times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

cosmology: theory↗
DOE OSTI2019

Doppler broadening in Mg 20 ( β p γ ) Ne 19 decay

The $^{15}$$\textbf{O}$$\Big({α,γ}\Big)$$^{19}$$\textbf{Ne}$ bottleneck reaction in Type I x-ray bursts is the most important thermonuclear reaction rate to constrain experimentally, to improve the accuracy of burst light-curve simulations. A proposed technique to determine the thermonuclear rate of this reaction employs the $^{20}$$\textbf{Mg}$$\Big({βpα}\Big)$$^{15}$$\textbf{O}$ decay sequence. The key $^{15}$$\textbf{O}$$\Big({α,γ}\Big)$$^{19}$$\textbf{Ne}$ resonance at an excitation of 4.03 MeV is now known to be fed in $^{20}$$\textbf{Mg}$$\Big({βpγ}\Big)$$^{19}$$\textbf{Ne}$ however, the energies of the protons feeding the 4.03 MeV state are unknown. Knowledge of the proton energies will facilitate future $^{20}$$\textbf{Mg}$$\Big({βpα}\Big)$$^{15}$$\textbf{O}$ measurements. To determine the energy of the proton transition feeding the 4.03 MeV state in $^{19}$$\textbf{Ne}$. A fast beam of $^{20}$$\textbf{Mg}$ was implanted into a plastic scintillator, which was used to detect $β$ particles. 16 high purity germanium detectors were used to detect $γ$ rays emitted following $βp$ decay. A Monte Carlo method was used to simulate the Doppler broadening of $^{19}$$\textbf{Ne}γ$-ray lines and compare to the experimental data. The center of mass energy between the proton and $^{19}$$\textbf{Ne}$, feeding the 4.03 MeV state, is measured to be $1.21^{+0.25}_{-0.22}$ $\textbf{MeV}$, corresponding to a $^{20}$$\textbf{Na}$ excitation energy of $7.44^{+0.25}_{-0.22}$$\textbf{MeV}$. Absolute feeding intensities and $γ$-decay branching ratios of $^{19}$$\textbf{Ne}$ states were determined including the 1615 keV state, which has not been observed before in this decay. A new $γ$ decay branch from the 1536 keV state in $^{19}$$\textbf{Ne}$ to the ground state is reported. The lifetime of the 1507 keV state in $^{19}$$\textbf{Ne}$ is measured to be $4.3^{+1.3}_{-1.1}$ ps resolving discrepancies in the literature. Conflicting $^{20}$$\textbf{Mg}$$\Big({βp}\Big)$ decay schemes in published literature are clarified.The utility of this Doppler broadening technique to provide information on $β$-delayed nucleon emission and excited-state lifetimes has been further demonstrated. In particular, knowledge of the proton energies feeding the 4.03 MeV $^{19}$$\textbf{Ne}$ state in $^{20}$$\textbf{Mg}$ $β$ decay will facilitate future measurements of the $α$-particle branching ratio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

GENFIRE: A generalized Fourier iterative reconstruction algorithm for high-resolution 3D imaging

Tomography has made a radical impact on diverse fields ranging from the study of 3D atomic arrangements in matter to the study of human health in medicine. Despite its very diverse applications, the core of tomography remains the same, that is, a mathematical method must be implemented to reconstruct the 3D structure of an object from a number of 2D projections. Here, we present the mathematical implementation of a tomographic algorithm, termed GENeralized Fourier Iterative REconstruction (GENFIRE), for high-resolution 3D reconstruction from a limited number of 2D projections. GENFIRE first assembles a 3D Fourier grid with oversampling and then iterates between real and reciprocal space to search for a global solution that is concurrently consistent with the measured data and general physical constraints. The algorithm requires minimal human intervention and also incorporates angular refinement to reduce the tilt angle error. We demonstrate that GENFIRE can produce superior results relative to several other popular tomographic reconstruction techniques through numerical simulations and by experimentally reconstructing the 3D structure of a porous material and a frozen-hydrated marine cyanobacterium. As a result, equipped with a graphical user interface, GENFIRE is freely available from our website and is expected to find broad applications across different disciplines.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Satellite remote sensing of H2SO4 aerosol using the 8- to 12-microns window region: Application to Mount Pinatubo

Monitoring stratospheric aerosols containing H2SO4 using the brightness temperature (BT) difference between 11 and 8.3 microns (BT(sub 8)-BT(sub 11)) spectral channels is demonstrated using theoretical calculations and satellite observations. Assuming an aqueous solution of 50% and 75% sulfuric acid, radiative transfer calculations indicate that over oceans an increase in the optical depth of the stratospheric aerosol results in an increase in BT(sub 8)-BT(sub 11). Theoretical simulations suggest that the technique is sensitive to visible optical depths greater than approximately 0.15. The simulations also demonstrate a lack of sensitivity to the particle size distribution. Changes in pre- and post-Pinatubo observations by the High-resolution Infrared Radiation Sounder 2 (HIRS2) on board the NOAA 10 are consistent with observed optical depth measurements and confirm the sensitivity of these channels to the presence of the aerosol. The technique is also applied to cold tropical convective clouds and desert regions where the signal, though evident, is less conclusive. Time series analysis is applied to the NOAA 10 and NOAA 12 combined BT(sub 8-BT(sub 11) observations to detect the periodicity of the spread of the volcanic aerosol. Over a region of the southern Pacific a 18- to 26-day period is present. Model simulations were conducted to demonstrate a trispectral with observations near 8, 11, and 12 microns. The trispectral approach has high potential in that the spectral signature of cirrus, water vapor, and H2SO4 aerosols are different. Observations from NOAA 10 and NOAA 11 are combined to demonstrate the capabilities of these infrared wavelengths of detecting the aerosol. The signal is clearly evident when a region of the South Atlantic is compared for pre- and post-Pinatubo conditions.

Ackerman, Steven A.↗
DOE OSTI2022

A Note on the Effects of Linear Topology Preservation in Monte Carlo Simulations of Knotted Proteins

Monte Carlo simulations are a powerful technique and are widely used in different fields. When applied to complex molecular systems with long chains, such as those in synthetic polymers and proteins, they have the advantage of providing a fast and computationally efficient way to sample equilibrium ensembles and calculate thermodynamic and structural properties under desired conditions. Conformational Monte Carlo techniques employ a move set to perform the transitions in the simulation Markov chain. While accepted conformations must preserve the sequential bonding of the protein chain model and excluded volume among its units, the moves themselves may take the chain across itself. We call this a break in linear topology preservation. In this manuscript, we show, using simple protein models, that there is no difference in equilibrium properties calculated with a move set that preserves linear topology and one that does not. However, for complex structures, such as those of deeply knotted proteins, the preservation of linear topology provides correct equilibrium results but only after long relaxation. In any case, to analyze folding pathways, knotting mechanisms and folding kinetics, the preservation of linear topology may be an unavoidable requirement.

Especial, João N. C. (ORCID:0000000254765370)↗
NASA NTRS1984

Initial research program for the National Transonic Facility

The construction and checkout of the National Transonic Facility (NTF) have been completed, and detailed calibration is now in progress. The initial NTF research program covers a wide range of study areas falling into three major elements: (1) the assessment of Reynolds number sensitivities for a broad range of configurations and flow phenomena; (2) validation of the ability of NTF to simulate full-scale aerodynamics; and (3) the development of test techniques for improved test simulations in existing wind tunnels. This paper, therefore, is a status report on these various elements of the initial NTF research program.

Gloss, B. B.↗
NASA NTRS2001

EFQPSK Versus CERN: A Comparative Study

This report presents a comparative study on Enhanced Feher's Quadrature Phase Shift Keying (EFQPSK) and Constrained Envelope Root Nyquist (CERN) techniques. These two techniques have been developed in recent times to provide high spectral and power efficiencies under nonlinear amplifier environment. The purpose of this study is to gain insights into these techniques and to help system planners and designers with an appropriate set of guidelines for using these techniques. The comparative study presented in this report relies on effective simulation models and procedures. Therefore, a significant part of this report is devoted to understanding the mathematical and simulation models of the techniques and their set-up procedures. In particular, mathematical models of EFQPSK and CERN, effects of the sampling rate in discrete time signal representation, and modeling of nonlinear amplifiers and predistorters have been considered in detail. The results of this study show that both EFQPSK and CERN signals provide spectrally efficient communications compared to filtered conventional linear modulation techniques when a nonlinear power amplifier is used. However, there are important differences. The spectral efficiency of CERN signals, with a small amount of input backoff, is significantly better than that of EFQPSK signals if the nonlinear amplifier is an ideal clipper. However, to achieve such spectral efficiencies with a practical nonlinear amplifier, CERN processing requires a predistorter which effectively translates the amplifier's characteristics close to those of an ideal clipper. Thus, the spectral performance of CERN signals strongly depends on the predistorter. EFQPSK signals, on the other hand, do not need such predistorters since their spectra are almost unaffected by the nonlinear amplifier, Ibis report discusses several receiver structures for EFQPSK signals. It is observed that optimal receiver structures can be realized for both coded and uncoded EFQPSK signals with not too much increase in computational complexity. When a nonlinear amplifier is used, the bit error rate (BER) performance of the CERN signals with a matched filter receiver is found to be more than one decibel (dB) worse compared to the bit error performance of EFQPSK signals. Although channel coding is found to provide BER performance improvement for both EFQPSK and CERN signals, the performance of EFQPSK signals remains better than that of CERN. Optimal receiver structures for CERN signals with nonlinear equalization is left as a possible future work. Based on the numerical results, it is concluded that, in nonlinear channels, CERN processing leads towards better bandwidth efficiency with a compromise in power efficiency. Hence for bandwidth efficient communications needs, CERN is a good solution provided effective adaptive predistorters can be realized. On the other hand, EFQPSK signals provide a good power efficient solution with a compromise in band width efficiency.

Borah, Deva K.↗