Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer simulations”

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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Advances in magnetospheric physics by numerical simulations - A critical analysis by an outsider

This is a selective, somewhat editorialized assessment, based on a literature survey, of computer simulation as it exists today in magnetospheric physics. Both large (MHD) and small scale dynamical simulations are described and considered from the perspective of what they are trying to do and with what success. Several specific problem areas where simulations are being carried out are called for commentary: global magnetospheric structure, magnetotail and magnetopause reconnection, Kelvin-Helmholtz instability, hydrodynamic expansion of ionospheric ions, electric double layers, heavy ion heating, and auroral kilometric radiation. It is concluded that simulation is a necessary tool for understanding magnetospheric physics and that significant progress has been made in simulation development. However, results should be evaluated knowing that many factors, some real and physical, others structural, may contribute to such output. A plea is made for greater intercalibration among different simulators working in parallel areas, so that facts can be distinguished from artifacts.

Birmingham, Thomas J.↗
NASA NTRS1995

Transferring ecosystem simulation codes to supercomputers

Many ecosystem simulation computer codes have been developed in the last twenty-five years. This development took place initially on main-frame computers, then mini-computers, and more recently, on micro-computers and workstations. Supercomputing platforms (both parallel and distributed systems) have been largely unused, however, because of the perceived difficulty in accessing and using the machines. Also, significant differences in the system architectures of sequential, scalar computers and parallel and/or vector supercomputers must be considered. We have transferred a grassland simulation model (developed on a VAX) to a Cray Y-MP/C90. We describe porting the model to the Cray and the changes we made to exploit the parallelism in the application and improve code execution. The Cray executed the model 30 times faster than the VAX and 10 times faster than a Unix workstation. We achieved an additional speedup of 30 percent by using the compiler's vectoring and 'in-line' capabilities. The code runs at only about 5 percent of the Cray's peak speed because it ineffectively uses the vector and parallel processing capabilities of the Cray. We expect that by restructuring the code, it could execute an additional six to ten times faster.

Skiles, J. W.↗
NASA NTRS1998

Probabilistic Modeling of Ceramic Matrix Composite Strength

Uncertainties associated with the primitive random variables such as manufacturing process (processing temperature, fiber volume ratio, void volume ratio), constituent properties (fiber, matrix and interface), and geometric parameters (ply thickness, interphase thickness) have been simulated to quantify the scatter in the first matrix cracking strength (FMCS) and the ultimate tensile strength of SCS-6/RBSN (SiC fiber (SCS-6) reinforced reaction-bonded silicon nitride composite) ceramic matrix composite laminate at room temperature. Cumulative probability distribution function for the FMCS and ultimate tensile strength at room temperature (RT) of (0)(sub 8), (0(sub 2)/90(sub 2), and (+/-45(sub 2))(sub S) laminates have been simulated and the sensitivity of primitive variables to the respective strengths have been quantified. Computationally predicted scatter of the strengths for a uniaxial laminate have been compared with those from limited experimental data. Also the experimental procedure used in the tests has been described briefly. Results show a very good agreement between the computational simulation and the experimental data. Dominating failure modes in (0)(sub 8), (0/90)(sub s) and (+/-45)(sub S) laminates have been identified. Results indicate that the first matrix cracking strength for the (0)(sub S), and (0/90)(sub S) laminates is sensitive to the thermal properties, modulus and strengths of both the fiber and matrix whereas the ultimate tensile strength is sensitive to the fiber strength and the fiber volume ratio. In the case of a (+/-45)(sub S), laminate, both the FMCS and the ultimate tensile strengths have a small scatter range and are sensitive to the fiber tensile strength as well as the fiber volume ratio.

Shan, Ashwin R.↗
NASA NTRS2005

Transient Three-Dimensional Analysis of Nozzle Side Load in Regeneratively Cooled Engines

Three-dimensional numerical investigations on the start-up side load physics for a regeneratively cooled, high-aspect-ratio nozzle were performed. The objectives of this study are to identify the three-dimensional side load physics and to compute the associated aerodynamic side load using an anchored computational methodology. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics formulation, and a transient inlet condition based on an engine system simulation. Computations were performed for both the adiabatic and cooled walls in order to understand the effect of boundary conditions. Finite-rate chemistry was used throughout the study so that combustion effect is always included. The results show that three types of shock evolution are responsible for side loads: generation of combustion wave; transitions among free-shock separation, restricted-shock separation, and simultaneous free-shock and restricted shock separations; along with oscillation of shocks across the lip. Wall boundary conditions drastically affect the computed side load physics: the adiabatic nozzle prefers free-shock separation while the cooled nozzle favors restricted-shock separation, resulting in higher peak side load for the cooled nozzle than that of the adiabatic nozzle. By comparing the computed physics with those of test observations, it is concluded that cooled wall is a more realistic boundary condition, and the oscillation of the restricted-shock separation flow pattern across the lip along with its associated tangential shock motion are the dominant side load physics for a regeneratively cooled, high aspect-ratio rocket engine.

Wang, Ten-See↗
NASA NTRS2004

Space Shuttle UHF Communications Performance Evaluation

An extension boom is to be installed on the starboard side of the Space Shuttle Orbiter (SSO) payload bay for thermal tile inspection and repairing. As a result, the Space Shuttle payload bay Ultra High Frequency (UHF) antenna will be under the boom. This study is to evaluate the Space Shuttle UHF communication performance for antenna at a suitable new location. To insure the RF coverage performance at proposed new locations, the link margin between the UHF payload bay antenna and Extravehicular Activity (EVA) Astronauts at a range distance of 160 meters from the payload bay antenna was analyzed. The communication performance between Space Shuttle Orbiter and International Space Station (SSO-ISS) during rendezvous was also investigated. The multipath effects from payload bay structures surrounding the payload bay antenna were analyzed. The computer simulation tool based on the Geometrical Theory of Diffraction method (GTD) was used to compute the signal strengths. The total field strength was obtained by summing the direct fields from the antennas and the reflected and diffracted fields from the surrounding structures. The computed signal strengths were compared to the signal strength corresponding to the 0 dB link margin. Based on the results obtained in this study, RF coverage for SSO-EVA and SSO- ISS communication links was determined for the proposed payload bay antenna UHF locations. The RF radiation to the Orbiter Docking System (ODS) pyros, the payload bay avionics, and the Shuttle Remote Manipulator System (SRMS) from the new proposed UHF antenna location was also investigated to ensure the EMC/EMI compliances.

Hwu, Shian U.↗
NASA NTRS2009

Presentation Extensions of the SOAP

A set of extensions of the Satellite Orbit Analysis Program (SOAP) enables simultaneous and/or sequential presentation of information from multiple sources. SOAP is used in the aerospace community as a means of collaborative visualization and analysis of data on planned spacecraft missions. The following definitions of terms also describe the display modalities of SOAP as now extended: In SOAP terminology, View signifies an animated three-dimensional (3D) scene, two-dimensional still image, plot of numerical data, or any other visible display derived from a computational simulation or other data source; a) "Viewport" signifies a rectangular portion of a computer-display window containing a view; b) "Palette" signifies a collection of one or more viewports configured for simultaneous (split-screen) display in the same window; c) "Slide" signifies a palette with a beginning and ending time and an animation time step; and d) "Presentation" signifies a prescribed sequence of slides. For example, multiple 3D views from different locations can be crafted for simultaneous display and combined with numerical plots and other representations of data for both qualitative and quantitative analysis. The resulting sets of views can be temporally sequenced to convey visual impressions of a sequence of events for a planned mission.

Carnright, Robert↗
DOE OSTI2021

Potential design problems for ITER fusion device

The international thermonuclear experimental reactor (ITER) is a worldwide project currently being built in France for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. ITER is of the tokamak based design using strong magnetic fields to confine the very hot plasma needed to induce the fusion reaction. Tokamak devices are currently the front leading designs. Building a successful magnetic fusion device for energy production is of great challenge. A key obstacle to such design is the performance during abnormal events including plasma disruptions and so-called edge-localized modes (ELMs). In these events, a massive and sudden release of energy occurs quickly, due to loss of full or partial plasma confinement, leading to very high transient power loads on the reactor surface boundaries. A successful reactor design should tolerate several of these transient events without serious damages such as melting and vaporization of the structure. This paper highlights, through comprehensive state-of-the-art computer simulation of the entire ITER interior design during such transient events, e.g., ELMs occurring at normal operation and disruptions during abnormal operation, potential serious problems with current plasma facing components (PFCs) design. The HEIGHTS computer package is used in these simulations. The ITER reactor design was simulated in full and exact 3D geometry including all known relevant physical processes involved during these transient events. The current ITER divertor design may not work properly and may requires significant modifications or new innovative design to prevent serious damage and to ensure successful operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Efficiency in nonequilibrium molecular dynamics Monte Carlo simulations

Hybrid algorithms combining nonequilibrium molecular dynamics and Monte Carlo (neMD/MC) offer a powerful avenue for improving the sampling efficiency of computer simulations of complex systems. These neMD/MC algorithms are also increasingly finding use in applications where conventional approaches are impractical, such as constant-pH simulations with explicit solvent. However, selecting an optimal nonequilibrium protocol for maximum efficiency often represents a non-trivial challenge. This work evaluates the efficiency of a broad class of neMD/MC algorithms and protocols within the theoretical framework of linear response theory. The approximations are validated against constant pH-MD simulations and shown to provide accurate predictions of neMD/MC performance. An assessment of a large set of protocols confirms (both theoretically and empirically) that a linear work protocol gives the best neMD/MC performance. Lastly, a well-defined criterion for optimizing the time parameters of the protocol is proposed and demonstrated with an adaptive algorithm that improves the performance on-the-fly with minimal cost.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

User's manual for a computer program for the emulation/simulation of a space station Environmental Control and Life Support System (ESCM)

This manual describes how to use the Emulation Simulation Computer Model (ESCM). Based on G189A, ESCM computes the transient performance of a Space Station atmospheric revitalization subsystem (ARS) with CO2 removal provided by a solid amine water desorbed subsystem called SAWD. Many performance parameters are computed some of which are cabin CO2 partial pressure, relative humidity, temperature, O2 partial pressure, and dew point. The program allows the user to simulate various possible combinations of man loading, metabolic profiles, cabin volumes and certain hypothesized failures that could occur.

Yanosy, James L.↗
DOE OSTI2018

An Event Detection Framework for Virtual Observation System: Anomaly Identification for an ACME Land Simulation

Based on previous work on in-situ data transfer infrastructure and compiler-based software analysis, we have designed a virtual observation system for real time computer simulations. This paper presents an event detection framework for a virtual observation system. By using signal processing and detection approaches to the memory-based data streams, this framework can be reconfigured to capture high-frequency events and low-frequency events. These approaches used in the framework can dramatically reduce the data transfer needed for in-situ data analysis (between distributed computing nodes or between the CPU/GPU nodes). In the paper, we also use a terrestrial ecosystem system simulation within the Earth System Model to demonstrate the practical values of this effort.

Yao, Cindy↗
NASA NTRS1989

Numerical simulation of controlled large space structures

Large Space Structures do not have much damping, which necessitates the installation of a controller onto the structure. If the controller is improperly designed, the structure may become unstable and be destroyed. Since Large Space Structures are extremely expensive pieces of hardware, new controllers must not be tested first on the structure. They must first be tested in computer simulations. Until now, the usual procedure for simulating controlled Large Space Structures is to compute a reduced order modal representation of the structure and then apply the controller. However, this procedure entails modal truncation error. A new software package which is free from this error is currently under development within the Center for Space Construction. The more accurate finite element representation of the structure is used in the simulation, instead of the less accurate reduced order modal representation. This software also features an efficient matrix storage scheme, which effectively deals with the asymmetric system matrices which occur when control is added to the structure. Also, an integration algorithm was chosen so that the simulation is a reliable indicator of system stability or instability. The software package is fairly general in nature. Linearity of the finite element model and of the controller is the only assumption made. Actuator dynamics, sensor dynamics, noise, and disturbances can be handled by the package. In addition, output feedback of displacement, velocity, and/or acceleration signals can be simulated. Kalman state estimation was also implemented. This software was tested on a finite element model of a real Large Space Structure: The Mini-Mast Truss. Mini-Mast is a testbed at NASA-Langley which is currently under development. A 714 degree of freedom finite element model was computed, and a 19 state controller was designed for it. Torque wheel dynamics were added to the model, and the entire closed loop system was simulated with the software package.

Quan, Ralph↗
NASA NTRS1990

Numerical simulations of interacting disk galaxies

Galaxy-galaxy interactions have long attracted many extragalactic astronomers in various aspects. A number of computer simulations performed in the 1970s have successfully reproduced the peculiar morphologies observed in interacting disk galaxies and clarified that tidal deformation explains most of the observed global peculiarities. However, most of these simulations have used test particles in modelling the disk component. Tidal response of a self-gravitating disk remains to be further clarified. Another topic which is intensely discussed at present is the relation between galaxy-galaxy interactions and activity. Many observations suggest that interactions trigger strong starbursts and possibly active galactic nuclei (AGN). However, the detailed mechanism of triggering is not yet clear. It is vital here to understand the dynamics of interstellar gas. In order to understand various phenomena related to galaxy-galaxy interactions (mainly for disk galaxies), the author performed a series of numerical simulations on close galaxy encounters which includes both interstellar gas and self-gravitating disk components. In these simulations, the galaxy model to be perturbed (target galaxy) consists of a halo and a disk. The halo was treated as a rigid spherical gravitational field which is assumed to remain fixed during the interaction. The disk is composed of stars and gas. The stellar disk was constructed by 20000 collisionless particles of the same mass. Those particles move in the halo gravitational field, interacting with each other and with the perturber. Therefore, the self-graviy of the disk is properly taken into account. Stellar particles were initially given circular velocities with small random motions required to stabilize the disk against local axisymmetric disturbances. The gravitational field of the stellar disk was calculated by the particle-mesh scheme (e.g. Hockney and Eastwood 1981). The gaseous component was modelled by the cloud-particle scheme (e.g. Roberts and Hausman 1984). Here, the authors represent the gas as an ensemble of small spheres (i.e. clouds) and include the creation of an OB star in a cloud-cloud collision and subsequent velocity push on nearby clouds due to a supernova explosion.

Noguchi, Masafumi↗
NASA NTRS1992

Comparison of simulated and actual wind shear radar data products

Prior to the development of the NASA experimental wind shear radar system, extensive computer simulations were conducted to determine the performance of the radar in combined weather and ground clutter environments. The simulation of the radar used analytical microburst models to determine weather returns and synthetic aperture radar (SAR) maps to determine ground clutter returns. These simulations were used to guide the development of hazard detection algorithms and to predict their performance. The structure of the radar simulation is reviewed. Actual flight data results from the Orlando and Denver tests are compared with simulated results. Areas of agreement and disagreement of actual and simulated results are shown.

Britt, Charles L.↗
NASA NTRS1992

Electric and hybrid electric vehicle study utilizing a time-stepping simulation

The applicability of NASA's advanced power technologies to electric and hybrid vehicles was assessed using a time-stepping computer simulation to model electric and hybrid vehicles operating over the Federal Urban Driving Schedule (FUDS). Both the energy and power demands of the FUDS were taken into account and vehicle economy, range, and performance were addressed simultaneously. Results indicate that a hybrid electric vehicle (HEV) configured with a flywheel buffer energy storage device and a free-piston Stirling convertor fulfills the emissions, fuel economy, range, and performance requirements that would make it acceptable to the consumer. It is noted that an assessment to determine which of the candidate technologies are suited for the HEV application has yet to be made. A proper assessment should take into account the fuel economy and range, along with the driveability and total emissions produced.

Schreiber, Jeffrey G.↗
NASA NTRS2006

Modeling of Transient Flow Mixing of Streams Injected into a Mixing Chamber

Ignition is recognized as one the critical drivers in the reliability of multiple-start rocket engines. Residual combustion products from previous engine operation can condense on valves and related structures thereby creating difficulties for subsequent starting procedures. Alternative ignition methods that require fewer valves can mitigate the valve reliability problem, but require improved understanding of the spatial and temporal propellant distribution in the pre-ignition chamber. Current design tools based mainly on one-dimensional analysis and empirical models cannot predict local details of the injection and ignition processes. The goal of this work is to evaluate the capability of the modern computational fluid dynamics (CFD) tools in predicting the transient flow mixing in pre-ignition environment by comparing the results with the experimental data. This study is a part of a program to improve analytical methods and methodologies to analyze reliability and durability of combustion devices. In the present paper we describe a series of detailed computational simulations of the unsteady mixing events as the cold propellants are first introduced into the chamber as a first step in providing this necessary environmental description. The present computational modeling represents a complement to parallel experimental simulations' and includes comparisons with experimental results from that effort. A large number of rocket engine ignition studies has been previously reported. Here we limit our discussion to the work discussed in Refs. 2, 3 and 4 which is both similar to and different from the present approach. The similarities arise from the fact that both efforts involve detailed experimental/computational simulations of the ignition problem. The differences arise from the underlying philosophy of the two endeavors. The approach in Refs. 2 to 4 is a classical ignition study in which the focus is on the response of a propellant mixture to an ignition source, with emphasis on the level of energy needed for ignition and the ensuing flame propagation issues. Our focus in the present paper is on identifying the unsteady mixing processes that provide the propellant mixture in which the ignition source is to be placed. In particular, we wish to characterize the spatial and temporal mixture distribution with a view toward identifying preferred spatial and temporal locations for the ignition source. As such, the present work is limited to cold flow (pre-ignition) conditions

Voytovych, Dmytro M.↗
DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗
DOE OSTI2021

RKKY coupled local-moment magnetism in NaFe 1 – x Cu x As

A central question in a large class of strongly correlated electron systems, including heavy fermion compounds and iron pnictides, is the identification of different phases and their origins. It has been shown that the antiferromagnetic (AFM) phase in some heavy fermion compounds is induced by Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between localized moments, and that the competition between this interaction and Kondo effect is responsible for quantum criticality. However, conclusive experimental evidence of the RKKY interaction in pnictides is lacking. Here, using high resolution 23 Na NMR measurements on lightly Cu-doped metallic single crystals of NaFe 1–x Cu x As (x ≈ 0.01) and numerical simulation, we show direct evidence of the RKKY interaction in this pnictide system. Aided by computer simulation, we identify the 23 Na NMR satellite resonances with the RKKY oscillations of spin polarization at Fe sites. Our 23 Na spin-lattice and spin-spin relaxation data exhibit a signature of an itinerant and inhomogeneous AFM phase in this system, accompanied by a simulation of Cu-induced perturbation to the ordered moments on the Fe sites. Our NMR results indicate coexistence of local and itinerant magnetism in lightly Cu-doped NaFe 1–x Cux As.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Dakota Technical Support Demonstrations and Software Improvements

Caterpillar, Inc. (Caterpillar) and Sandia National Laboratories (SNL) continued their partnership to improve SNL's open source Dakota software product and effectively deploy it to Caterpillar engineering design and analysis problems. Several new features that are expected to benefit not only Caterpillar but also the entire Dakota user community were introduced. In addition, SNL and Caterpillar cooperated extensively to conduct a demonstration of Dakota's uncertainty quantification capabilities on a computational simulation of relevance to Caterpillar.

97 MATHEMATICS AND COMPUTING↗