Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation Use Case”

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 433 records · Page 24

Molecular Modes Elucidate the Nuclear Magnetic Resonance Relaxation of Viscous Fluids

The Bloembergen, Purcell, and Pound (BPP) theory of nuclear magnetic resonance (NMR) relaxation in fluids dating back to 1948 continues to be the linchpin in interpreting NMR relaxation data in applications ranging from characterizing fluids in porous media to medical imaging (MRI). The BPP theory is founded on assuming molecules are hard spheres with 1 H– 1 H dipole pairs reorienting randomly; assumptions that are severe in light of modern understanding of liquids. Nevertheless, it is intriguing to this day that the BPP theory was consistent with the original experimental data for glycerol, a hydrogen-bonding molecular fluid for which the hard-sphere-rigid-dipole assumption is inapplicable. To better understand this incongruity, atomistic molecular simulations are used to compute 1 H NMR T 1 relaxation dispersion (i.e., frequency dependence) in two contrasting cases: glycerol, and a (non hydrogen-bonding) viscosity standard. At high viscosities, simulations predict distinct functional forms of T 1 for glycerol compared to the viscosity standard, in agreement with modern measurements, yet both in contrast to BPP theory. The cause of these departures from BPP theory is elucidated, without assuming any relaxation models and without any free parameters, by decomposing the simulated T 1 response into dynamic molecular modes for both intramolecular and intermolecular interactions. The decomposition into dynamic molecular modes provides an alternative framework to understand the physics of NMR relaxation for viscous fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Perspective on Quantum Computing Applications in Quantum Chemistry Using 25-100 Logical Qubits

The intersection of quantum computing and quantum chemistry represents a promising frontier for achieving quantum utility in domains of both scientific and societal relevance. Owing to the exponential growth of classical resource requirements for simulating quantum systems, quantum chemistry has long been recognized as a natural candidate for quantum computation. This perspective focuses on identifying scientifically meaningful use cases where early fault-tolerant quantum computers, which are considered to be equipped with approximately 25-100 logical qubits, could deliver tangible impact. While recent advances in classical computing have pushed the boundaries of tractable simulations to unprecedented scales, this logical-qubit regime represents the first window where quantum devices can pursue qualitatively distinct strategies, such as polynomial-scaling phase estimation, direct simulation of quantum dynamics, and active-space embedding, that remain challenging for classical solvers, such as multireference charge-transfer and conical-intersection states central to photochemistry and materials design. We highlight near-term opportunities in algorithm and software design, discuss representative chemical problems suited for quantum acceleration, and propose strategic roadmaps and collaborative pathways for advancing practical quantum utility in quantum chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating Matrix Crack and Delamination Interaction in a Clamped Tapered Beam

Blind predictions were conducted to validate a discrete crack methodology based on the Floating Node Method to simulate matrix-crack/delamination interaction. The main novel aspects of the approach are: (1) the implementation of the floating node method via an 'extended interface element' to represent delaminations, matrix-cracks and their interaction, (2) application of directional cohesive elements to infer overall delamination direction, and (3) use of delamination direction and stress state at the delamination front to determine migration onset. Overall, good agreement was obtained between simulations and experiments. However, the validation exercise revealed the strong dependence of the simulation of matrix-crack/delamination interaction on the strength data (in this case transverse interlaminar strength, YT) used within the cohesive zone approach applied in this work. This strength value, YT, is itself dependent on the test geometry from which the strength measurement is taken. Thus, choosing an appropriate strength value becomes an ad-hoc step. As a consequence, further work is needed to adequately characterize and assess the accuracy and adequacy of cohesive zone approaches to model small crack growth and crack onset. Additionally, often when simulating damage progression with cohesive zone elements, the strength is lowered while keeping the fracture toughness constant to enable the use of coarser meshes. Results from the present study suggest that this approach is not recommended for any problem involving crack initiation, small crack growth or multiple crack interaction.

De Carvalho, N. V.↗

Laser wavelength dependence of particle acceleration mechanisms in high intensity laser–solid density plasma interactions

We investigate the generation of relativistic electrons and the subsequent ion acceleration due to target-normal sheath acceleration when ultra-intense (⁠ I > 10 18 W/cm 2 ⁠) short pulse (⁠ τ L < 10ps⁠) lasers are incident onto solid density targets as laser wavelength is varied. Scaling laws for the hot electron temperature, T hot ⁠, and the maximum ion energy, E max ⁠, are recast as a function of laser wavelength. These predictions are compared to results from particle-in-cell computer simulations in a variety of geometries, including cases where realistic plasma density profiles as determined by a radiation hydrodynamics code are used. It is found that the wavelength dependence observed in simulation is less pronounced than what is predicted from the well-established scaling laws. An assessment of how switching to longer laser wavelengths, specifically 2 μm Tm:YLF technology, would impact current high energy density science applications and diagnostics is made.

Electromagnetism↗

Radar Reflectivity Simulated by a 2-D Spectra Bin Model: Sensitivity of Cloud-aerosol Interaction

The Goddard Cumulus Ensemble (GCE) model with bin spectra microphysics is used to simulate mesoscale convective systems.The model uses explicit bins to represent size spectra of cloud nuclei, water drops, ice crystals, snow and graupel. Each hydrometeorite category is described by 33 mass bins. The simulations provide a unique data set of simulated raindrop size distribution in a realistic dynamic frame. Calculations of radar parameters using simulated drop size distribution serve as an evaluation of numerical model performance. In addition, the GCE bin spectra modes is a very useful tool to study uncertainties related to radar observations; all the environmental parameters are precisely known. In this presentation, we concentrate on the discussion of Z-R (ZDR-R) relation in the simulated systems. Due to computational limitations, the spectra bin model has been run in two dimensions with 31 stretched vertical layers and 1026 horizontal grid points (1 km resolution). Two different cases, one in midlatitude continent, the other in tropical ocean, have been simulated. The continental case is a strong convection which lasted for two hours. The oceanic case is a persistent system with more than 10 hours' life span. It is shown that the simulated Z-R (ZDR-R) relations generally agree with observations using radar and rain gauge data. The spatial and temporal variations of Z-R relation in different locations are also analyzed. Impact of aerosols on cloud formation and raindrop size distribution was studied. Both clean (low CCN) and dirty (high CCN) cases are simulated. The Z-R relation is shown to vary considerable in the initial CCN concentrations.

Li, Kiaowen↗

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few a priori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, M.↗

Design of a Laboratory Hall Thruster with Magnetically Shielded Channel Walls, Phase III: Comparison of Theory with Experiment

A proof-of-principle effort to demonstrate a technique by which erosion of the acceleration channel in Hall thrusters of the magnetic-layer type can be eliminated has been completed. The first principles of the technique, now known as "magnetic shielding," were derived based on the findings of numerical simulations in 2-D axisymmetric geometry. The simulations, in turn, guided the modification of an existing 6-kW laboratory Hall thruster. This magnetically shielded (MS) thruster was then built and tested. Because neither theory nor experiment alone can validate fully the first principles of the technique, the objective of the 2-yr effort was twofold: (1) to demonstrate in the laboratory that the erosion rates can be reduced by >order of magnitude, and (2) to demonstrate that the near-wall plasma properties can be altered according to the theoretical predictions. This paper concludes the demonstration of magnetic shielding by reporting on a wide range of comparisons between results from numerical simulations and laboratory diagnostics. Collectively, we find that the comparisons validate the theory. Near the walls of the MS thruster, theory and experiment agree: (1) the plasma potential has been sustained at values near the discharge voltage, and (2) the electron temperature has been lowered by at least 2.5-3 times compared to the unshielded (US) thruster. Also, based on carbon deposition measurements, the erosion rates at the inner and outer walls of the MS thruster are found to be lower by at least 2300 and 1875 times, respectively. Erosion was so low along these walls that the rates were below the resolution of the profilometer. Using a sputtering yield model with an energy threshold of 25 V, the simulations predict a reduction of ~600 at the MS inner wall. At the outer wall ion energies are computed to be below 25 V, for which case we set the erosion to zero in the simulations. When a 50-V threshold is used the computed ion energies are below the threshold at both sides of the channel. Uncertainties, sensitivities and differences between theory and experiment are also discussed.

plasma simulations↗

Statistical 3D 'Atomistic' Simulation of Decanano MOSFETs

A 3D statistical 'atomistic' simulation technique has been developed to study the effect of the random dopant induced parameter fluctuations in aggressively scaled MOSFETs. Efficient implementation of the 'atomistic' simulation approach has been used to investigate the threshold voltage standard deviation and lowering in the case of uniformly doped MOSFETs, and in fluctuation-resistant architectures utilising epitaxial-layers and delta-doping. The effect of the random doping in the polysilicon gate on the threshold voltage fluctuations has also been thoroughly investigated. The influence of a single-charge trapping on the channel conductivity in decanano MOSFETs is studied in the 'atomistic' framework as well. Quantum effects are taken into consideration in our 'atomistic' simulations using the density gradient formalism.

Asenov, Asen↗

Effect of Impact Location on the Response of Shuttle Wing Leading Edge Panel 9

The objective of this paper is to compare the results of several simulations performed to determine the worst-case location for a foam impact on the Space Shuttle wing leading edge. The simulations were performed using the commercial non-linear transient dynamic finite element code, LS-DYNA. These simulations represent the first in a series of parametric studies performed to support the selection of the worst-case impact scenario. Panel 9 was selected for this study to enable comparisons with previous simulations performed during the Columbia Accident Investigation. The projectile for this study is a 5.5-in cube of typical external tank foam weighing 0.23 lb. Seven locations spanning the panel surface were impacted with the foam cube. For each of these cases, the foam was traveling at 1000 ft/s directly aft, along the orbiter X-axis. Results compared from the parametric studies included strains, contact forces, and material energies for various simulations. The results show that the worst case impact location was on the top surface, near the apex.

Lyle, Karen H.↗

A Study of Heavy Precipitation Events in Taiwan During 10-13 August, 1994: Mesoscale Model Simulations - Part 2

Several heavy precipitation episodes occurred over Taiwan from August 10 to 13, 1994. Precipitation patterns and characteristics are quite different between the precipitation events that occurred from August 10 and I I and from August 12 and 13. In Part I (Chen et al. 2001), the environmental situation and precipitation characteristics are analyzed using the EC/TOGA data, ground-based radar data, surface rainfall patterns, surface wind data, and upper air soundings. In this study (Part II), the Penn State/NCAR Mesoscale Model (MM5) is used to study the precipitation characteristics of these heavy precipitation events. Various physical processes (schemes) developed at NASA Goddard Space Flight Center (i.e., cloud microphysics scheme, radiative transfer model, and land-soil-vegetation surface model) have recently implemented into the MM5. These physical packages are described in the paper, Two way interactive nested grids are used with horizontal resolutions of 45, 15 and 5 km. The model results indicated that Cloud physics, land surface and radiation processes generally do not change the location (horizontal distribution) of heavy precipitation. The Goddard 3-class ice scheme produced more rainfall than the 2-class scheme. The Goddard multi-broad-band radiative transfer model reduced precipitation compared to a one-broad band (emissivity) radiation model. The Goddard land-soil-vegetation surface model also reduce the rainfall compared to a simple surface model in which the surface temperature is computed from a Surface energy budget following the "force-re store" method. However, model runs including all Goddard physical processes enhanced precipitation significantly for both cases. The results from these runs are in better agreement with observations. Despite improved simulations using different physical schemes, there are still some deficiencies in the model simulations. Some potential problems are discussed. Sensitivity tests (removing either terrain or radiative processes) are performed to identify the physical processes that determine the precipitation patterns and characteristics for heavy rainfall events. These sensitivity tests indicated that terrain can play a major role in determining the exact location for both precipitation events. The terrain can also play a major role in determining the intensity of precipitation for both events. However, it has a large impact on one event but a smaller one on the other. The radiative processes are also important for determining, the precipitation patterns for one case but. not the other. The radiative processes can also effect the total rainfall for both cases to different extents.

Tao, Wei Kuo↗

Designing with Ada for satellite simulation: A case study

A FORTRAN oriented and an Ada oriented design for the same system are compared to learn whether an essentially different design was produced using Ada. The designs were produced by an experiment that involves the parallel development of software for a spacecraft dynamics simulator. Design differences are identified in the use of abstractions, system structure, and simulator operations. Although the designs were vastly different, this result may be influenced by some special characteristics discussed.

Agresti, W. W.↗

Modeling a Miniaturized Scanning Electron Microscope Focusing Column - Lessons Learned in Electron Optics Simulation

This presentation discusses work done to assess the design of a focusing column in a miniaturized Scanning Electron Microscope (SEM) developed at the NASA Marshall Space Flight Center (MSFC) for use in-situ on the Moon-in particular for mineralogical analysis. The MSFC beam column design uses purely electrostatic fields for focusing, because of the severe constraints on mass and electrical power consumption imposed by the goals of lunar exploration and of spaceflight in general. The resolution of an SEM ultimately depends on the size of the focused spot of the scanning beam probe, for which the stated goal here is a diameter of 10 nanometers. Optical aberrations are the main challenge to this performance goal, because they blur the ideal geometrical optical image of the electron source, effectively widening the ideal spot size of the beam probe. In the present work the optical aberrations of the mini SEM focusing column were assessed using direct tracing of non-paraxial rays, as opposed to mathematical estimates of aberrations based on paraxial ray-traces. The geometrical ray-tracing employed here is completely analogous to ray-tracing as conventionally understood in the realm of photon optics, with the major difference being that in electron optics the lens is simply a smoothly varying electric field in vacuum, formed by precisely machined electrodes. Ray-tracing in this context, therefore, relies upon a model of the electrostatic field inside the focusing column to provide the mathematical description of the "lens" being traced. This work relied fundamentally on the boundary element method (BEM) for this electric field model. In carrying out this research the authors discovered that higher accuracy in the field model was essential if aberrations were to be reliably assessed using direct ray-tracing. This led to some work in testing alternative techniques for modeling the electrostatic field. Ultimately, the necessary accuracy was attained using a BEM/Fourier series hybrid approach. The presentation will give background remarks about the MSFC mini Lunar SEM concept and electron optics modeling, followed by a description of the alternate field modeling techniques that were tried, along with their incorporation into a ray-trace simulation. Next, the validation of this simulation against commercially available software will be discussed using an example lens as a test case. Then, the efficacy of aberration assessment using direct ray-tracing will be demonstrated, using this same validation case. The discussion will include practical error checks of the field solution. Finally, the ray-trace assessment of the MSFC mini Lunar SEM concept will be shown and discussed. The authors believe this presentation will be of general interest to practitioners of modeling and simulation, as well as those with a general optics background. Because electron optics and photon optics share many basic concepts (e.g., lenses, images, aberrations, etc.), the appeal of this presentation need not be restricted to just those interested in charged particle optics.

Modeling Miniaturizing Scanning↗

Probabilistic Analysis of a Composite Crew Module

An approach for conducting reliability-based analysis (RBA) of a Composite Crew Module (CCM) is presented. The goal is to identify and quantify the benefits of probabilistic design methods for the CCM and future space vehicles. The coarse finite element model from a previous NASA Engineering and Safety Center (NESC) project is used as the baseline deterministic analysis model to evaluate the performance of the CCM using a strength-based failure index. The first step in the probabilistic analysis process is the determination of the uncertainty distributions for key parameters in the model. Analytical data from water landing simulations are used to develop an uncertainty distribution, but such data were unavailable for other load cases. The uncertainty distributions for the other load scale factors and the strength allowables are generated based on assumed coefficients of variation. Probability of first-ply failure is estimated using three methods: the first order reliability method (FORM), Monte Carlo simulation, and conditional sampling. Results for the three methods were consistent. The reliability is shown to be driven by first ply failure in one region of the CCM at the high altitude abort load set. The final predicted probability of failure is on the order of 10-11 due to the conservative nature of the factors of safety on the deterministic loads.

FROM↗

Fine pointing of Laser Beams by using Laser Arrays for Applications to CubeSats

A method of fine pointing of laser beams by using laser arrays has been developed. The telescope system combines a lens system and a VCSEL/Photodetector Array. It does not use moving parts. In computer simulations, it is applied to CubeSats that use body pointing. Body pointing was used by the Aerospace Corporation for CubeSats in LEO in NASA’s Optical Communications and Sensors Demonstration (OCSD) program. Computer simulations of this fine pointing capability have been applied previously to CubeSats in the OCSD program. In this paper, computer simulations of laser pointing using this telescope design are applied to CubeSats in LLO, at 100 km. These CubeSats could form part of the LunaNet, the lunar communications and navigation network, part of the NASA ARTEMIS Program. With more accurate pointing, a laser beam with smaller divergence can be used. For the case of the AeroCube-7B vehicle that was used in the OCSD program, computer simulations will show, for example, that the divergence of the output beam can be reduced from approximately 0.06° FWHM to 0.014°. For the proposed electro-optical system, reaction times to pointing changes and vibrations are on a nanosecond time scale, much faster than those for fine pointing systems that use moving parts such as fast steering mirrors, including MEMS, or that use quad-cell photodetectors to improve the body pointing of the CubeSat. Other possible applications are to Optical Multiple Access (OMA)for simultaneously communicating with ground stations at different locations and to Wavelength-Division Multiplexing (WDM)for increasing data rate transmission.

SMD↗

Fine Pointing of Laser Beams by Using Laser Arrays for Applications to CubeSats

A method of fine pointing of laser beams by using laser arrays has been developed. The telescope system combines a lens system and a VCSEL/Photodetector Array. It does not use moving parts. In computer simulations, it is applied to CubeSats that use body pointing. Body pointing was used by the Aerospace Corporation for CubeSats in LEO in NASA’s Optical Communications and Sensors Demonstration (OCSD) program. Computer simulations of this fine pointing capability have been applied previously to CubeSats in the OCSD program. In this paper, computer simulations of laser pointing using this telescope design are applied to CubeSats in LLO, at 100 km. These CubeSats could form part of the LunaNet, the lunar communications and navigation network, part of the NASA ARTEMIS Program. With more accurate pointing, a laser beam with smaller divergence can be used. For the case of the AeroCube-7B vehicle that was used in the OCSD program, computer simulations will show, for example, that the divergence of the output beam can be reduced from approximately 0.06° FWHM to 0.014°. For the proposed electro-optical system, reaction times to pointing changes and vibrations are on a nanosecond time scale, much faster than those for fine pointing systems that use moving parts such as fast steering mirrors, including MEMS, or that use quad-cell photodetectors to improve the body pointing of the CubeSat. Other possible applications are to Optical Multiple Access (OMA)for simultaneously communicating with ground stations at different locations and to Wavelength-Division Multiplexing (WDM)for increasing data rate transmission.

HEC↗

Digital simulation of the serpentuator using MARSYAS

Serpentuator is a serpentine teleoperator device for intravehicular and extravehicular activities in space. The serpentuator is simulated using simulation software system MARSYAS and using the Component-Connection Simulation model and the Direct Simulation model. A comparison of the results for the two cases shows that under identical conditions, simulation execution time in the Component-Connection model case is reduced by a factor of the order of 100. A visual display of the serpentuator positions is obtained using the AMTRAN system on the Datacraft DC 6024 computer.

Singh, S., P.↗

Detailed Characterization of CZT Detector Response for Improved Coded-Aperture Imaging Performance

Gamma-ray imaging is a powerful method for locating and quantifying sources of radiation. The coded-aperture technique demonstrates superior angular resolution in comparison to other methods (e.g., Compton reconstruction). In this method, a mask constructed of highly attenuating material encodes the scene as a shadow pattern on a position-sensitive detector; this pattern can then be used to recreate the origin(s) of incident radiation. This is typically done through convolution of the mask and shadow patterns. Iterative methods which attempt to reconstruct the observed shadow pattern using a weighted combination of simulated patterns may also be employed. In either case, errors in event position reconstruction due to detector imperfections alter the shadow pattern and will therefore degrade system performance and may introduce imaging artifacts. These effects can be mitigated with a detailed understanding of such errors – allowing for the generation of representative simulations that include the errors and/or correction of raw imager data to remove the errors. We present a calibration process for a commercially available cadmium zinc telluride (CZT) gamma imager which provides a comprehensive characterization of the spatial and energy dependence of event reconstruction. By illuminating a mask featuring a regular grid of pinholes with a calibration source, the localized response of the detector can be measured with fine granularity. These local responses are combined to generate a full detector response map which can be used to distort simulations in a manner that is representative of the observed detector data. Details of the calibration procedure and an assessment of the impact of its end products on the performance of iterative imaging methods will be presented.

Ziock, Klaus-Peter↗

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗