Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Virtual Reality Simulation of the International Space Welding Experiment

Virtual Reality (VR) is a set of breakthrough technologies that allow a human being to enter and fully experience a 3-dimensional, computer simulated environment. A true virtual reality experience meets three criteria: (1) It involves 3-dimensional computer graphics; (2) It includes real-time feedback and response to user actions; and (3) It must provide a sense of immersion. Good examples of a virtual reality simulator are the flight simulators used by all branches of the military to train pilots for combat in high performance jet fighters. The fidelity of such simulators is extremely high -- but so is the price tag, typically millions of dollars. Virtual reality teaching and training methods are manifestly effective, and we have therefore implemented a VR trainer for the International Space Welding Experiment. My role in the development of the ISWE trainer consisted of the following: (1) created texture-mapped models of the ISWE's rotating sample drum, technology block, tool stowage assembly, sliding foot restraint, and control panel; (2) developed C code for control panel button selection and rotation of the sample drum; (3) In collaboration with Tim Clark (Antares Virtual Reality Systems), developed a serial interface box for the PC and the SGI Indigo so that external control devices, similar to ones actually used on the ISWE, could be used to control virtual objects in the ISWE simulation; (4) In collaboration with Peter Wang (SFFP) and Mark Blasingame (Boeing), established the interference characteristics of the VIM 1000 head-mounted-display and tested software filters to correct the problem; (5) In collaboration with Peter Wang and Mark Blasingame, established software and procedures for interfacing the VPL DataGlove and the Polhemus 6DOF position sensors to the SGI Indigo serial ports. The majority of the ISWE modeling effort was conducted on a PC-based VR Workstation, described below.

Phillips, James A.↗
NASA NTRS1993

Probabilistic assessment of composite structures

A general purpose methodology for integrated probabilistic assessment of composite structures is demonstrated using aircraft fuselage type structures with rectangular cutouts. The computational simulation took into account buckling loads, vibration frequencies, global displacements, and local stresses. The procedure was embedded into the computer code IPACS (Integrated Probabilistic Assessment of Composite Structures) that calculates the simulated scatter, the sensitivity of the composite structural behavior to all the primitive variables that affect the structural behavior. The probabilistic assessment for the composite structure with rectangular cutouts indicates that the uncertainty in the longitudinal ply stress is mainly due to uncertainty in the laminate thickness.

Shiao, Michael C.↗
NASA NTRS2019

SI IV Resonance Line Emission during Solar Flares: Non-LTE, Non-equilibrium, Radiation Transfer Simulations

The Interface Region Imaging Spectrograph routinely observes the Si IV resonance lines. When analyzing quiescent observations of these lines, it has typically been assumed that they form under optically thin conditions. This is likely valid for the quiescent Sun, but this assumption has also been applied to the more extreme flaring scenario. We used 36 electron-beam-driven radiation hydrodynamic solar flare simulations, computed using the RADYN code, to probe the validity of this assumption. Using these simulated atmospheres, we solved the radiation transfer equations to obtain the non-LTE, nonequilibrium populations, line profiles, and opacities for a model silicon atom, including charge exchange processes. This was achieved using the "minority species" version of RADYN. The inclusion of charge exchange resulted in a substantial fraction of Si IV at cooler temperatures than those predicted by ionization equilibrium. All simulations with an injected energy flux F> 5 x 10(exp 10) erg cm(exp -2) s(exp -1) resulted in optical depth effects on the Si IV emission, with differences in both intensity and line shape compared to the optically thin calculation. Weaker flares (down to F ~ 10(exp 9) erg cm(exp -2) s(exp -1)) also resulted in Si IV emission forming under optically thick conditions, depending on the other beam parameters. When opacity was significant, the atmospheres generally had column masses in excess of 5 x 10-6 g cm-2 over the temperature range 40-100 kK, and the Si IV formation temperatures were between 30 and 60 kK. We urge caution when analyzing Si IV flare observations, or when computing synthetic emission without performing a full radiation transfer calculation.

Kerr, Graham S.↗
DOE OSTI2022

A full-ring variable-aperture cadmium zinc telluride system for whole-body single photon emission computed tomography: realistic simulations with phantoms

Single photon emission computed tomography (SPECT) is an imaging modality that has demonstrated its utility in a number of clinical indications. Despite this progress, a high sensitivity, high spatial resolution, multi-tracer SPECT with a large field of view suitable for whole-body imaging of a broad range of radiotracers for theranostics is not available. Purpose We have designed a cadmium zinc telluride (CZT) variable-aperture full-ring SPECT scanner instrumented with a broad-energy tungsten collimator intended to fill this technological gap. The final purpose is to provide a multi-tracer solution for brain and whole-body imaging. Our static SPECT scanner breaks the paradigm of the standard dual- and triple-head rotational SPECT systems, utilizing a larger detector area in each scan increasing the sensitivity. We provide a demonstration of the performance of our design using a realistic model of our detector with simulated body-sized 99mTc phantoms. Methods We developed a realistic model of our detector by using a combination of a Geant4 Monte Carlo simulation and a CZT detector response model based on a finite element model. Our approach models the characteristic low-energy tail effect in CZT that noticeably affects the sensitivity and the quality of the scatter correction in CZT detectors. We implement a modified dual energy window scatter correction adapted to include the CZT low-energy tail effect. A dedicated correction is also developed to eliminate the undesirable truncation observed in images given the presence of detector edges and gaps between detectors, due to the non-rotational nature of our device. Corrections for the attenuation, detector response and the presence of collimators are also included. The images are reconstructed using the maximum-likelihood expectation-maximization algorithm implemented in the reconstruction open software STIR. Detector and reconstruction performance are characterized with a Derenzo phantom and a body-sized National Electrical Manufacturers Association (NEMA) Image Quality (IQ) phantom containing 99mTc. Results Our SPECT design can resolve 6.4mm rods in a Derenzo phantom and obtain a good image contrast with the IQ phantom. Explicit testing of the gap and edge correction is provided, showing an excellent performance in eliminating the image truncation artifacts. Our modified scatter correction shows no overestimation of the contrast-recovery ratio for our realistic CZT detector model, as opposed to the cases without correction and with a standard dual-energy window scatter correction. Conclusions In this paper, we further demonstrate the performance of our design for whole-body imaging purposes. This adds to our previous demonstration of improved qualitative and quantitative 99mTc imaging for brain perfusion and 123I imaging for dopamine transport with respect to state-of-the-art NaI dual-head cameras. We show that our design performs similarly to the VERITON SPECT from Spectrum Dynamics, a commercial full-ring CZT SPECT camera, with the potential advantage of the broader energy range of application given by our custom-design tungsten collimators. Here, our device combines high sensitivity and image resolution with a broad-energy imaging application for the purpose of clinical imaging and theranostics of emerging radionuclides.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2007

Progressive Fracture of Laminated Composite Stiffened Plate

Laminated fiber-reinforced composite stiffened plate with [0/90/plus or minus 45]s plies made of S-Glass/epoxy are evaluated via computational simulation to study damage and fracture progression. The loads are pressure and temperature which varies from 21 to 65.5 C (case I) and from 143.3 to 21 C (case II). An integrated computer code is used for the simulation of the damage progression. Results show that damage initiation begins at low load level, with matrix cracking at the 0 deg. (bottom and top) plies, fiber fracture at the bottom (0 deg.) ply and interply delamination at the top (0 deg. ) ply. Increasing the applied pressure, the damage growth is expended resulting in fracture through the thickness of the structure. At this stage, 90 percent of the plies damage at applied pressure 15.306 MPa for the case I and 15.036 MPa for the case II. After this stage the cracks propagate rapidly and the structure collapses.

Gotsis, Pascalis K.↗
NASA NTRS1999

A One Dimensional, Time Dependent Inlet/Engine Numerical Simulation for Aircraft Propulsion Systems

The NASA Lewis Research Center (LeRC) and the Arnold Engineering Development Center (AEDC) have developed a closely coupled computer simulation system that provides a one dimensional, high frequency inlet/engine numerical simulation for aircraft propulsion systems. The simulation system, operating under the LeRC-developed Application Portable Parallel Library (APPL), closely coupled a supersonic inlet with a gas turbine engine. The supersonic inlet was modeled using the Large Perturbation Inlet (LAPIN) computer code, and the gas turbine engine was modeled using the Aerodynamic Turbine Engine Code (ATEC). Both LAPIN and ATEC provide a one dimensional, compressible, time dependent flow solution by solving the one dimensional Euler equations for the conservation of mass, momentum, and energy. Source terms are used to model features such as bleed flows, turbomachinery component characteristics, and inlet subsonic spillage while unstarted. High frequency events, such as compressor surge and inlet unstart, can be simulated with a high degree of fidelity. The simulation system was exercised using a supersonic inlet with sixty percent of the supersonic area contraction occurring internally, and a GE J85-13 turbojet engine.

Doug Garrard↗
NASA NTRS1989

Development and simulation study of a new inverse-pinch high Coulomb transfer switch

The inverse-pinch plasma switch was studied using a computer simulation code. The code was based on a 2-D, 2-temperature magnetohydrodynamic (MHD) model. The application of this code was limited to the disk-type inverse-pinch plasma switch. The results of the computer analysis appear to be in agreement with the experimental results when the same parameters are used. An inverse-pinch plasma switch for closing has been designed and tested for high-power switching requirements. An azimuthally uniform initiation of breakdown is a key factor in achieving an inverse-pinch current path in the switch. Thus, various types of triggers, such as trigger pins, wire-brush, ring trigger, and hypocycloidal-pinch (HCP) devices have been tested for uniform breakdown. Recently, triggering was achieved by injection of a plasma-ring (plasma puff) that is produced separately with hypocycloidal-pinch electrodes placed under the cathode of the main gap. The current paths at switch closing, initiated by the injection of a plasma-ring from the HCP trigger are azimuthally uniform, and the local current density is significantly reduced, so that damage to the electrodes and the insulator surfaces is minimized. The test results indicate that electron bombardment on the electrodes and the insulator surfaces is minimized. The test results indicate that electron bombardment on the electrodes is four orders of magnitude less than that of a spark-gap switch for the same switching power. Indeed, a few thousand shots with peak current exceeding a mega-ampere and with hold-off voltage up to 20 kV have been conducted without showing measurable damage to the electrodes and insulators.

Choi, Sang H.↗
NASA NTRS2016

An Efficient Ray-Tracing Method for Determining Terrain Intercepts in EDL Simulations

The calculation of a ray's intercept from an arbitrary point in space to a prescribed surface is a common task in computer simulations. The arbitrary point often represents an object that is moving according to the simulation, while the prescribed surface is fixed in a defined frame. For detailed simulations, this surface becomes complex, taking the form of real-world objects such as mountains, craters or valleys which require more advanced methods to accurately calculate a ray's intercept location. Incorporation of these complex surfaces has commonly been implemented in graphics systems that utilize highly optimized graphics processing units to analyze such features. This paper proposes a simplified method that does not require computationally intensive graphics solutions, but rather an optimized ray-tracing method for an assumed terrain dataset. This approach was developed for the Mars Science Laboratory mission which landed on the complex terrain of Gale Crater. First, this paper begins with a discussion of the simulation used to implement the model and the applicability of finding surface intercepts with respect to atmosphere modeling, altitude determination, radar modeling, and contact forces influencing vehicle dynamics. Next, the derivation and assumptions of the intercept finding method are presented. Key assumptions are noted making the routines specific to only certain types of surface data sets that are equidistantly spaced in longitude and latitude. The derivation of the method relies on ray-tracing, requiring discussion on the formulation of the ray with respect to the terrain datasets. Further discussion includes techniques for ray initialization in order to optimize the intercept search. Then, the model implementation for various new applications in the simulation are demonstrated. Finally, a validation of the accuracy is presented along with the corresponding data sets used in the validation. A performance summary of the method will be shown using the analysis from the Mars Science Laboratory's terminal descent sensing model. Alternate uses will also be shown for determining horizon maps and orbiter set times.

Shidner, Jeremy D.↗
DOE OSTI2023

Code-verification techniques for the method-of-moments implementation of the combined-field integral equation

Code verification plays an important role in establishing the credibility of computational simulations by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, the numerical solution to integral equations incurs multiple interacting sources of numerical error, as well as other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the combined-field integral equation. Here, we demonstrate the effectiveness of these approaches for cases with and without coding errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Transferability of data-driven, many-body models for CO 2 simulations in the vapor and liquid phases

Here, extending on the previous work by Riera et al. [J. Chem. Theory Comput. 16, 2246–2257 (2020)], we introduce a second generation family of data-driven many-body MB-nrg models for CO 2 and systematically assess how the strength and anisotropy of the CO 2 –CO 2 interactions affect the models’ ability to predict vapor, liquid, and vapor–liquid equilibrium properties. Building upon the many-body expansion formalism, we construct a series of MB-nrg models by fitting one-body and two-body reference energies calculated at the coupled cluster level of theory for large monomer and dimer training sets. Advancing from the first generation models, we employ the charge model 5 scheme to determine the atomic charges and systematically scale the two-body energies to obtain more accurate descriptions of vapor, liquid, and vapor–liquid equilibrium properties. Challenges in model construction arise due to the anisotropic nature and small magnitude of the interaction energies in CO 2 , calling for the necessity of highly accurate descriptions of the multidimensional energy landscape of liquid CO 2 . These findings emphasize the key role played by the training set quality in the development of transferable, data-driven models, which, accurately representing high-dimensional many-body effects, can enable predictive computer simulations of molecular fluids across the entire phase diagram.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1973

Multibeam antenna study, Phase 2

An antenna system concept for point to point communication between any two points in the Continental United States was developed using a synchronous orbit satellite. The design and fabrication of a two dimensional artificial dielectric lens contained in a parallel plate waveguide system is presented. The experimental pattern and voltage standing wave ratio results are listed from tests with four different lenses: a plano-convex natural dielectric (Teflon) lens, a plano-convex artificial dielectric lens, a two-point corrected Teflon lens, and a two-point corrected artificial dielectric lens. A computer simulated design of an operational two lens, 24 beam configuration is also presented. The computer design presents overlays of each beam onto the Continental United States from a fixed synchronous orbit satellite with a definition of these areas where beam-to-beam isolation is not sufficient to permit the reuse of frequency.

Byrnes, P. J.↗
NASA NTRS1978

SAMICS marketing and distribution model

A SAMICS (Solar Array Manufacturing Industry Costing Standards) was formulated as a computer simulation model. Given a proper description of the manufacturing technology as input, this model computes the manufacturing price of solar arrays for a broad range of production levels. This report presents a model for computing these marketing and distribution costs, the end point of the model being the loading dock of the final manufacturer.

Source record↗
NASA NTRS1978

Hydrogen maser frequency standard computer model for automatic cavity tuning servo simulations

A computer model of the JPL hydrogen maser frequency standard was developed. This model allows frequency stability data to be generated, as a function of various maser parameters, many orders of magnitude faster than these data can be obtained by experimental test. In particular, the maser performance as a function of the various automatic tuning servo parameters may be readily determined. Areas of discussion include noise sources, first-order autotuner loop, second-order autotuner loop, and a comparison of the loops.

Potter, P. D.↗
NASA NTRS1979

Simulation of auroral double layers

Some basic properties of plasma double layers are deduced from a particle-in-cell computer simulation and related to parallel electric-field structures above the auroral regions. The simulation results on the processes leading to double-layer formation are examined, particularly in relation to the transient stage and double-layer structure and stability. It is concluded that: (1) a large potential difference applied to a finite-length plasma will be concentrated in a shocklike localized region instead of occurring over the entire length of the system; (2) the initial stage in double-layer formation is dominated by a large-potential pulse propagating in the direction of the induced electrostatic drift; (3) the entire potential is dropped over a specific scale length once the double layer has formed; and (4) this scale length is expected to be of the order of 1 km for a double layer above a discrete auroral arc with a potential of 10 kV and the electric-field vector parallel to the magnetic-field vector.

Hubbard, R. F.↗
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↗