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 793 records · Page 44

Modeling and Simulation Credibility Assessments of Musculoskeletal Computational Models for Simulating Astronaut Injuries Due to a Poor Spacesuit Fit

The musculoskeletal (MS) system of astronauts is subject to physiological changes, potentially leading to injuries due to the exposure to different gravitational environments experienced during spaceflight. These injuries can occur while an astronaut is performing an Extravehicular Activity (EVA) in space, on lunar or planetary surfaces or while wearing a spacesuit during terrestrial training for an EVA. The OpenSim MS modeling software can assess EVA induced MS injury mechanisms such as muscle strains, ligament injuries and joint injuries. One area of concern, since there are only a few different spacesuit sizes with limited adjustability, is the possibility of a poorly fitting spacesuit. This can cause unnatural joint motions and torques resulting in various MS injuries. A credibility assessment of the OpenSim modeling and simulation procedures is performed per NASA-STD-7009A to provide information on the credibility of the model’s use in simulating EVA related injury mechanisms. The credibility assessment evaluated various OpenSim models against the following eight credibility factors: data pedigree, input pedigree, code verification, solution verification, conceptual validation, referent validation, results uncertainty and results robustness (sensitivity). The models evaluated for EVA injuries will require additional credibility factor analysis and upgrades to the model features, such as adding ligaments to a whole-body model, to reliably predict and analyze the EVA injuries expected to occur due to a poor spacesuit fit. The degree of elevation strategy required to increase the credibility assessment scores will depend on the model complexity and the injury mechanism.

Christopher A Gallo↗

Generation of spiral bevel gears with zero kinematical errors and computer aided simulation of their meshing and contact

There is proposed a method for generation of Gleason's spiral bevel gears which provides the following properties of meshing and contact: (1) the contact normal keeps its original direction within the neighborhood of the main contact point; (2) the contact ellipse moves along the gear tooth surface; and (3) the kinematical errors caused by Gleason's method of cutting are almost zero. Computer programs for the simulation of meshing and bearing contact are developed.

Litvin, F. L.↗

The Virtual Blast Furnace - An Integrated High Performance Computing Modeling, Simulation, and Visualization Capability for Steel Manufacturing (Final Report)

Many manufacturing industries require substantial capital and utilize energy intensive processes that involve complex phenomena. One example of such an industry is the steel industry, which is the fourth largest energy consuming industry in the U.S. By harnessing the power of High-Performance Computing (HPC) to enhance current simulation and visualization methods in the steel industry, it should be possible to increase resolution and/or decrease time of these methods by a factor of 1000. In this way, information can be obtained in a time frame that is useful for making business and engineering decisions, optimizing manufacturing processes and, ultimately, improving the completeness of U.S. industries. For example, if coke usage in blast furnaces were optimized such that the average coke rate was reduced from 797 lb/net tonne of hot metal (NTHM) to 604 lb/NTHM, costs could be reduced by $894 million/year. Additionally, members of the steel industry need the flexibility to efficiently operate blast furnaces at a range of production rates in order to meet fluctuating market demands. Large scale parameter studies can be utilized to discover workable operating parameters at a range of production rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Object-Oriented Graphical User Interface for a Reusable Rocket Engine Intelligent Control System

An intelligent control system for reusable rocket engines under development at NASA Lewis Research Center requires a graphical user interface to allow observation of the closed-loop system in operation. The simulation testbed consists of a real-time engine simulation computer, a controls computer, and several auxiliary computers for diagnostics and coordination. The system is set up so that the simulation computer could be replaced by the real engine and the change would be transparent to the control system. Because of the hard real-time requirement of the control computer, putting a graphical user interface on it was not an option. Thus, a separate computer used strictly for the graphical user interface was warranted. An object-oriented LISP-based graphical user interface has been developed on a Texas Instruments Explorer 2+ to indicate the condition of the engine to the observer through plots, animation, interactive graphics, and text.

Litt, Jonathan S.↗

Modulating factors of hydrologic exchanges in a large-scale river reach: Insights from three-dimensional computational fluid dynamics simulations

Hyporheic exchange is a critical mechanism shaping hydrological and biogeochemical processes along a river corridor. Recent studies on quantifying the hyporheic exchange were mostly limited to local scales due to field inaccessibility, computational demand, and complexity of geomorphology and subsurface geology. Surface flow conditions and subsurface physical properties are well known factors on modulating the hyporheic exchange, but quantitative understanding of their impacts on the strength and direction of hyporheic exchanges at reach scales is absent. In this study, a high resolution computational fluid dynamics (CFD) model that couples surface and subsurface flow and transport is employed to simulate hyporheic exchanges in a 7-km long reach along the main-stem of the Columbia River. Assuming that the hyporheic exchange does not affect surface water flow conditions due to its negligible magnitude compared to the volume and velocity of river water, we developed a one-way coupled surface and subsurface water flow model using the commercial CFD software STAR-CCM+. The model integrates the Reynolds-averaged Navier-Stokes (RANS) equation solver with a realizable k-e two-layer turbulence model, a two-layer all y^+ wall treatment, and the volume of fluid (VOF) method, and is used to simulate hyporheic exchanges by tracking the free water-air interface as well as flow in the river and the subsurface porous media. The model is validated against measurements from acoustic Doppler current profiler (ADCP) in the stream water and hyporheic fluxes derived from a set of temperature profilers installed across the riverbed. The validated model is then employed to systematically investigate how hyporheic exchanges are influenced by surface water fluid dynamics strongly regulated by upstream dam operations, as well as subsurface structures (e.g. thickness of riverbed and subsurface formation layers) and hydrogeological properties (e.g. permeability). The results suggest that the thickness of riverbed alluvium layer is the dominant factor for reach-scale hyporheic exchanges, followed by the alluvium permeability, the depth of the underlying impermeable layer, and the assumption of hydrostatic pressure.

Bao, Jie↗

Real-time simulation of F100-PW-100 turbofan engine using the hybrid computer

A real-time hybrid computer simulation of the F100-PW-100 augmented turbofan is presented. The digital portion of the hybrid computer is used to perform the bivariate function generation associated with modeling the performance of the engine's rotating components. The remaining calculations are performed on the analog computer. Steady state simulation data along with sea level, static, transient data are presented to show that the real-time simulation matches baseline digital simulation results over a wide range of power settings and flight conditions. Steady state simulation data are compared with sea level, experimental data to show that the real-time hybrid and baseline digital simulations do adequately predict the performance of the actual engine. FORTRAN listings and analog patching diagrams are provided.

Szuch, J. R.↗

Digital quantum simulation of cavity quantum electrodynamics: insights from superconducting and trapped ion quantum testbeds

We explore the potential for hybrid development of quantum hardware where currently available quantum computers simulate open cavity quantum electrodynamical (CQED) systems for applications in optical quantum communication, simulation and computing. Our simulations make use of a recent quantum algorithm that maps the dynamics of a singly excited open Tavis–Cummings model containing N atoms coupled to a lossy cavity. We report the results of executing this algorithm on two noisy intermediate-scale quantum computers: a superconducting processor and a trapped ion processor, to simulate the population dynamics of an open CQED system featuring N = 3 atoms. By applying technology-specific transpilation and error mitigation techniques, we minimize the impact of gate errors, noise, and decoherence in each hardware platform, obtaining results which agree closely with the exact solution of the system. These results can be used as a recipe for efficient and platform-specific quantum simulation of cavity–emitter systems on contemporary and future quantum computers.

cavity QED↗

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sampling Size Optimization for Bioburden Density Estimation in Planetary Protection

Planetary protection (PP) is a discipline that focuses on minimizing the biological contamination of spacecraft to ensure compliance with international policy. Precise estimation of bioburden - the total number of microbes in or on spacecraft hardware – and the bioburden density are of utmost importance for PP. Such estimation is the way concordance with requirements is demonstrated, and it is critical for quantifying the potential risk of inadvertently contaminating other planetary bodies. Although a suite of molecular techniques have been used to thoroughly characterize and profile the microbiome of various cleanroom environments and spacecraft, the gold standard remains the physical enumeration of microbes via culturing of samples directly taken from spacecraft and associated surfaces. However, due to technical, budgetary, and programmatic constraints, only a manageable portion (around 10%) of the entire spacecraft surface is directly sampled with cotton swabs or wipes. To generate the bioburden current best estimate (CBE) for components not directly verifiable, the accepted approach is to apply a NASA-defined bioburden estimate based on the components’ manufacturing or assembly environment. This approach utilizes a prespecified bioburden density estimation that applies a maximum value across the total surface area of the specified component. For hardware components that underwent similar assembly processes, an implied bioburden is adopted for all components, based on a direct verification of a representative component within the same lot. Once all components have a CBE, the bioburden estimates are generated. In previous publication [ 1], we have shown that statistical risks quantifying the accuracy of the estimates for sampled, prespecified, and implied components can be derived and ranked. For mean squared error (MSE) function, the risks are available analytically and hence a cost function can be obtained to optimize the risks with respect to the sampling area and sampling cost. Since the sampling area and sampling cost are two complimentary variables, their sum will have a well-defined minimum. This paper presents the multivariate optimization of the integrated risk of an empirical Bayes estimator to determine the optimal sampling schedule for a given number of components. It is assumed that given a number of components, N, the bioburden density for each component can either be sampled, implied, or prespecified. The multivariate optimization searches through different options to sample, imply or prespecify the bioburden density for a component, and account for the component’s surface area and cost of sampling. The idea of the optimization is based on the observation that the statistical risk of using an estimator is a monotonically decreasing function of the sampled area. The larger the sampled area, the lower the risk of using the estimator as the estimator becomes more and more accurate as the sampling area increases. On the other hand, the cost of sampling is monotonically increasing as the sampled surface grows. This makes the risk and total cost of sampling complimentary variables which can be counterbalanced to achieve an optimal overall value with respect to the sampled surface. In this paper, the integrated risk has been used to quantify the accuracy of the estimator. This risk has been selected because it depends on neither the true value of the parameter nor on the collected data. The cost of each sample was also available to obtain the total cost of sampling of N components. The paper will present the results based on computer-simulated data as well as the data collected during the InSight mission. The computer-simulated data have N components with randomly generated total areas and each component assigned to one of the three categories according to the method of estimating of bioburden density: sampled, implied, or prespecified. The cost of sampling is also available. The cost of sampling is estimated based on a cost model provided by the planetary protection group at JPL. For this paper, the overall cost was assumed to be a linear function of exposure. The optimization process finds the allocation of the components to the three categories that minimizes the tradeoff between integrated risk and total cost. For the InSight data, a set of components is selected representing all three categories, and optimization is performed to determine if the performed allocation was optimal or if a better allocation could have been obtained. To the best of our knowledge, this work is the first attempt not only perform an accurate estimation of bioburden density but also do it in an optimal way.

97 - MATHEMATICS AND COMPUTING↗

Real-time simulation of an automotive gas turbine using the hybrid computer

A hybrid computer simulation of an Advanced Automotive Gas Turbine Powertrain System is reported. The system consists of a gas turbine engine, an automotive drivetrain with four speed automatic transmission, and a control system. Generally, dynamic performance is simulated on the analog portion of the hybrid computer while most of the steady state performance characteristics are calculated to run faster than real time and makes this simulation a useful tool for a variety of analytical studies.

Costakis, W.↗

Noncoherent DTTLs for Symbol Synchronization

Noncoherent data-transition tracking loops (DTTLs) have been proposed for use as symbol synchronizers in digital communication receivers. [Communication- receiver subsystems that can perform their assigned functions in the absence of synchronization with the phases of their carrier signals ( carrier synchronization ) are denoted by the term noncoherent, while receiver subsystems that cannot function without carrier synchronization are said to be coherent. ] The proposal applies, more specifically, to receivers of binary phase-shift-keying (BPSK) signals generated by directly phase-modulating binary non-return-to-zero (NRZ) data streams onto carrier signals having known frequencies but unknown phases. The proposed noncoherent DTTLs would be modified versions of traditional DTTLs, which are coherent. The symbol-synchronization problem is essentially the problem of recovering symbol timing from a received signal. In the traditional, coherent approach to symbol synchronization, it is necessary to establish carrier synchronization in order to recover symbol timing. A traditional DTTL effects an iterative process in which it first generates an estimate of the carrier phase in the absence of symbol-synchronization information, then uses the carrier-phase estimate to obtain an estimate of the symbol-synchronization information, then feeds the symbol-synchronization estimate back to the carrier-phase-estimation subprocess. In a noncoherent symbol-synchronization process, there is no need for carrier synchronization and, hence, no need for iteration between carrier-synchronization and symbol- synchronization subprocesses. The proposed noncoherent symbolsynchronization process is justified theoretically by a mathematical derivation that starts from a maximum a posteriori (MAP) method of estimation of symbol timing utilized in traditional, coherent DTTLs. In that MAP method, one chooses the value of a variable of interest (in this case, the offset in the estimated symbol timing) that causes a likelihood function of symbol estimates over some number of symbol periods to assume a maximum value. In terms that are necessarily oversimplified to fit within the space available for this article, it can be said that the mathematical derivation involves a modified interpretation of the likelihood function that lends itself to noncoherent DTTLs. The proposal encompasses both linear and nonlinear noncoherent DTTLs. The performances of both have been computationally simulated; for comparison, the performances of linear and nonlinear coherent DTTLs have also been computationally simulated. The results of these simulations show that, among other things, the expected mean-square timing errors of coherent and noncoherent DTTLs are relatively insensitive to window width. The results also show that at high signal-to-noise ratios (SNRs), the performances of the noncoherent DTTLs approach those of their coherent counterparts at, while at low SNRs, the noncoherent DTTLs incur penalties of the order of 1.5 to 2 dB.

Simon, Marvin↗

Universal Solvent Viscosity Reduction via Hydrogen Bonding Disruptors

Liquid Ion Solutions LLC (DBA RoCo Global) in partnership with Carnegie Mellon University and Carbon Capture Scientific LLC, has performed lab-scale development and evaluation of novel additives that lower the viscosity of water-lean amine solvents for post-combustion carbon dioxide capture. This project focuses on developing additives that minimize the formation of long-range electrostatic and hydrogen bonding (HB) networks, decreasing the solvent viscosity, improving diffusion, and improving the process economics. The project objectives included: 1) performing computer simulation to understand the molecular interactions of the additive molecules in water-lean CO 2 capture solvents, 2) design and synthesis of HB disruptors additives, 3) performance testing with additive molecules on model amine solvents, and 4) demonstration of the effectiveness of the optimized additives in the presence of synthetic flue gas. To meet the abovementioned objectives, the project team utilized a holistic approach that combines molecular simulation, experimental testing, and economic analysis studies. The project team developed ab initio molecular model and then perform computer simulation to develop relationship between hydrogen bonding, viscosity, and performed quantitative analysis of additive on the viscosity of the solvent. The team completed computational comparative study on a range of organic functional groups such as ethers, esters, cyclic carbonates, alkanes, and ammonium salts for their effect on viscosity gaining key insights into molecular interactions and the impact of various functional groups and molecular shapes on viscosity. Assisted with molecular simulation insights, the project team conducted additive synthesis and testing, including a proof-of-concept study, additive screening, optimization, and synthetic flue gas testing. The experimental proof-of-concept study proved that the hydrogen bonding acceptors result in significant decrease of viscosities. Detailed additive screening (exploring various functionalities and molecular structures) has been performed. Several promising additives showed excellent reduction in viscosity (30-41%) at 5% additive loading, and over 50% viscosity reduction at 10% additive loading for the model solvents. The team also performed complex screening studies on additive loadings and mixing effect among additives using the design of experiments. Based on multiple screening experiments, one additive-solvent candidate was down-selected for synthetic flue gas testing. A 100-hour continuous absorption/desorption study was conducted under simulated flue gas using a lab-scale continuous capture and separation system. No degradation (for both solvent and additive) was observed based on the GC results of the solvent samples collected from the continuous study. The team conducted preliminary engineering analyses and cost-benefit analyses to quantify the potential economic benefits of the additive approach for solvent viscosity reduction. Based on the experimental data, CO 2 capture cost savings from the capital and operating cost savings are estimated at $\$$4.7/tonne and $\$$0.3/tonne CO 2 captured, respectively. Considering the additive cost, the net benefit is estimated to be between $\$$4.32~$\$$4.86/tonne CO 2 captured.

20 FOSSIL-FUELED POWER PLANTS↗

Avoiding pitfalls in simulating real-time computer systems

The software simulation of a computer target system on a computer host system, known as an interpretive computer simulator (ICS), functionally models and implements the action of the target hardware. For an ICS to function as efficiently as possible and to avoid certain pitfalls in designing an ICS, it is important that the details of the hardware architectural design of both the target and the host computers be known. This paper discusses both host selection considerations and ICS design features that, without proper consideration, could make the resulting ICS too slow to use or too costly to maintain and expand.

Smith, R. S.↗

Increasing Compressed Gas Energy Storage Density Using CO2–N2 Gas Mixture

This paper demonstrates a new method by which the energy storage density of compressed air systems is increased by 56.8% by changing the composition of the compressed gas to include a condensable component. A higher storage density of 7.33 MJ/m3 is possible using a mixture of 88% CO2 and 12% N2 compared to 4.67 MJ/m3 using pure N2. This ratio of gases representing an optimum mixture was determined through computer simulations that considered a variety of different proportions from pure CO2 to pure N2. The computer simulations are based on a thermodynamic equilibrium model that predicts the mixture composition as a function of volume and pressure under progressive compression to ultimately identify the optimal mixture composition (88% CO2 + 12% N2). The model and simulations predict that the optimal gas mixture attains a higher energy storage density than using either of the pure gases.

25 ENERGY STORAGE↗

An empirical analysis of the distribution of the duration of overshoots in a stationary gaussian stochastic process

This analysis utilizes computer simulation and statistical estimation. Realizations of stationary gaussian stochastic processes with selected autocorrelation functions are computer simulated. Analysis of the simulated data revealed that the mean and the variance of a process were functionally dependent upon the autocorrelation parameter and crossing level. Using predicted values for the mean and standard deviation, by the method of moments, the distribution parameters was estimated. Thus, given the autocorrelation parameter, crossing level, mean, and standard deviation of a process, the probability of exceeding the crossing level for a particular length of time was calculated.

Parrish, R. S.↗

A directional gamma-ray telescope using coded aperture techniques

A directional detector for gamma-ray astronomy has been developed to image sources in the energy range 0.1 to 5 MeV. An array of 35 gain stabilized bismuth germanate detectors, together with a coded aperture mask based on a uniformly redundant array allows imaging in 4 deg square sky bins over a 16 x 24 deg field-of-view. The position of a strong point source, such as the Crab Nebula, can be determined to within not more than about 1 deg. A complementary 'anti-mask' greatly reduces systematic effects arising from nonuniform background rates amongst the detectors. The telescope has an effective area of 190 sq cm and an energy resolution of 19.5 percent FWHM at 662 keV. Results of laboratory tests of the imaging system, including the ability to image multiple sources, uniformity of response over the field-of-view, and the effect of the 'anti-mask', are in good agreement with computer simulations. Features of the flight detector system are described and results of laboratory tests and computer simulations are reviewed. A balloon flight of the telescope is planned for the fall of 1982.

Mcconnell, M. L.↗

Unsteady RANS and Scale Resolving Simulations of Open Rotor Noise

Computational simulations using structured curvilinear overset grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) framework are presented for predicting the aerodynamic and aeroacoustic behavior of a contra-rotating open rotor (CROR) propulsion system using a newly developed and computationally efficient sliding mesh technique. A comprehensive evaluation of two turbulence closure treatments is carried out as part of the scope of this work: Unsteady Reynolds-Averaged Navier-Stokes Simulation (URANS), and Hybrid RANS/LES (HRLES). The computed results are then validated against existing experimental measurements at both the nominal take-off and cruise flight conditions. To propagate the noise generated by the near-field CROR system, aeroacoustic analysis based on a frequency-domain permeable surface Ffowcs Williams-Hawkings (FW-H) approach is applied.

ARMD↗