Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation modelling”

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,387 records · Page 77

Acceleration techniques for dependability simulation

As computer systems increase in complexity, the need to project system performance from the earliest design and development stages increases. We have to employ simulation for detailed dependability studies of large systems. However, as the complexity of the simulation model increases, the time required to obtain statistically significant results also increases. This paper discusses an approach that is application independent and can be readily applied to any process-based simulation model. Topics include background on classical discrete event simulation and techniques for random variate generation and statistics gathering to support simulation.

Barnette, James David↗

Second Generation Crop Yield Models Review

Second generation yield models, including crop growth simulation models and plant process models, may be suitable for large area crop yield forecasting in the yield model development project. Subjective and objective criteria for model selection are defined and models which might be selected are reviewed. Models may be selected to provide submodels as input to other models; for further development and testing; or for immediate testing as forecasting tools. A plant process model may range in complexity from several dozen submodels simulating (1) energy, carbohydrates, and minerals; (2) change in biomass of various organs; and (3) initiation and development of plant organs, to a few submodels simulating key physiological processes. The most complex models cannot be used directly in large area forecasting but may provide submodels which can be simplified for inclusion into simpler plant process models. Both published and unpublished models which may be used for development or testing are reviewed. Several other models, currently under development, may become available at a later date.

Hodges, T.↗

High-Performance Modeling and Simulation of Anchoring in Granular Media for NEO Applications

NASA is interested in designing a spacecraft capable of visiting a near-Earth object (NEO), performing experiments, and then returning safely. Certain periods of this mission would require the spacecraft to remain stationary relative to the NEO, in an environment characterized by very low gravity levels; such situations require an anchoring mechanism that is compact, easy to deploy, and upon mission completion, easy to remove. The design philosophy used in this task relies on the simulation capability of a high-performance multibody dynamics physics engine. On Earth, it is difficult to create low-gravity conditions, and testing in low-gravity environments, whether artificial or in space, can be costly and very difficult to achieve. Through simulation, the effect of gravity can be controlled with great accuracy, making it ideally suited to analyze the problem at hand. Using Chrono::Engine, a simulation pack age capable of utilizing massively parallel Graphic Processing Unit (GPU) hardware, several validation experiments were performed. Modeling of the regolith interaction has been carried out, after which the anchor penetration tests were performed and analyzed. The regolith was modeled by a granular medium composed of very large numbers of convex three-dimensional rigid bodies, subject to microgravity levels and interacting with each other with contact, friction, and cohesional forces. The multibody dynamics simulation approach used for simulating anchors penetrating a soil uses a differential variational inequality (DVI) methodology to solve the contact problem posed as a linear complementarity method (LCP). Implemented within a GPU processing environment, collision detection is greatly accelerated compared to traditional CPU (central processing unit)- based collision detection. Hence, systems of millions of particles interacting with complex dynamic systems can be efficiently analyzed, and design recommendations can be made in a much shorter time. The figure shows an example of this capability where the Brazil Nut problem is simulated: as the container full of granular material is vibrated, the large ball slowly moves upwards. This capability was expanded to account for anchors of different shapes and penetration velocities, interacting with granular soils.

Quadrelli, Marco B.↗

Application and Evaluation of a Snowmelt Runoff Model in the Tamor River Basin, Eastern Himalaya Using a Markov Chain Monte Carlo (MCMC) Data Assimilation Approach

Previous studies have drawn attention to substantial hydrological changes taking place in mountainous watersheds where hydrology is dominated by cryospheric processes. Modelling is an important tool for understanding these changes but is particularly challenging in mountainous terrain owing to scarcity of ground observations and uncertainty of model parameters across space and time. This study utilizes a Markov Chain Monte Carlo data assimilation approach to examine and evaluate the performance of a conceptual, degree-day snowmelt runoff model applied in the Tamor River basin in the eastern Nepalese Himalaya. The snowmelt runoff model is calibrated using daily streamflow from 2002 to 2006 with fairly high accuracy (average Nash-Sutcliffe metric approx. 0.84, annual volume bias <3%). The Markov Chain Monte Carlo approach constrains the parameters to which the model is most sensitive (e.g. lapse rate and recession coefficient) and maximizes model fit and performance. Model simulated streamflow using an interpolated precipitation data set decreases the fractional contribution from rainfall compared with simulations using observed station precipitation. The average snowmelt contribution to total runoff in the Tamor River basin for the 2002-2006 period is estimated to be 29.7+/-2.9% (which includes 4.2+/-0.9% from snowfall that promptly melts), whereas 70.3+/-2.6% is attributed to contributions from rainfall. On average, the elevation zone in the 4000-5500m range contributes the most to basin runoff, averaging 56.9+/-3.6% of all snowmelt input and 28.9+/-1.1% of all rainfall input to runoff. Model simulated streamflow using an interpolated precipitation data set decreases the fractional contribution from rainfall versus snowmelt compared with simulations using observed station precipitation. Model experiments indicate that the hydrograph itself does not constrain estimates of snowmelt versus rainfall contributions to total outflow but that this derives from the degree-day melting model. Lastly, we demonstrate that the data assimilation approach is useful for quantifying and reducing uncertainty related to model parameters and thus provides uncertainty bounds on snowmelt and rainfall contributions in such mountainous watersheds.

runoff↗

Precipitation Characteristics and Future Changes Over the Southern Slope of Tibetan Plateau Simulated by a High-Resolution Global Nonhydrostatic Model

Current climate model resolution cannot accurately describe the complex topography over the Tibetan Plateau, which limits our understanding of past and future precipitation over this region. This study investigates the daily precipitation characteristics and its future changes over the Tibetan Plateau, especially for the southern slope of Tibetan Plateau (SSTP), by 14-km Nonhydrostatic ICosahedral Atmospheric Model (NICAM) with explicitly calculated convection for historical and future 30-years period. By comparing with the satellite Global Precipitation Measurement (GPM), NICAM well reproduces the historical precipitation spatial pattern, seasonal cycle, and the extreme precipitation belt over SSTP, but overestimates precipitation amount by~35%. It is found that heavy precipitation probability decreases as elevation becomes higher while the light precipitation generally shows the opposite. For the precipitation changes during June to September from 1979–2008 to 2075–2104, NICAM predicts that mean precipitation will decrease over low-level SSTP but increase over high-level SSTP. Nonprecipitation and heavy precipitation probability will increase while light precipitation probability will decrease in the future over SSTP. The extreme precipitation probability and intensity will increase ~50%/°C and ~8%/°C over SSTP, and this increase is more obvious as elevation becomes higher. The robust increase of extreme precipitation along the SSTP topography is unique and has not been identified by the climate model simulations before. The strong meridional gradient of specific humidity over SSTP is found to be further enhanced under global warming, and this gradient enhancement is suggested to be responsible for the increase in the extreme precipitation over SSTP.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rotor systems research aircraft simulation mathematical model

An analytical model developed for evaluating and verifying advanced rotor concepts is discussed. The model was used during in both open loop and real time man-in-the-loop simulation during the rotor systems research aircraft design. Future applications include: pilot training, preflight of test programs, and the evaluation of promising concepts before their implementation on the flight vehicle.

Houck, J. A.↗

Efficient and Accurate Terrain Integration & Application to Modeling and Simulation of Planetary Landing Systems

As planetary spacecraft landing system complexity increases, pre-flight performance assessments include correspondingly more detailed and complex models that represent the various environmental elements and vehicle subsystems. Specifically, the relevance and criticality of planetary terrain modeling in the context of precision landing and hazard avoidance (PL&HA) has become evident as both commercial and national agencies return to the Lunar surface and beyond.

Terrain↗

NAME Modeling and Data Assimilation: A Strategic Overview

In this presentation a strategic overview of modeling and related data analysis and assimilation components of the North American Monsoon Experiment (NAME) are given. Building on the NAME science plan, a strategy is outlined for accelerating progress on the fundamental modeling issues pertaining to NAME science goals. The strategy takes advantage of NAME enhanced observations, and should simultaneously provide model-based guidance to the evolving multi-tiered NAME observing program. The NAME modeling strategy recognizes three distinct, but related, roles that observations play in model development and assessment. These are: (1) to guide model development by providing constraints on model simulations at the process level (e.g. convection, land/atmosphere and ocean/atmosphere interactions); (2) to help assess the veracity of model simulations of the various key NAMS phenomena (e.g. low level jets, land sea breezes, tropical storms), and the linkages to regional and larger-scale climate variability; and (3) to provide initial and boundary conditions, and verification data for model predictions.

Schubert, Siegfried↗

Users guide: The LaRC human-operator-simulator-based pilot model

A Human Operator Simulator (HOS) based pilot model has been developed for use at NASA LaRC for analysis of flight management problems. The model is currently configured to simulate piloted flight of an advanced transport airplane. The generic HOS operator and machine model was originally developed under U.S. Navy sponsorship by Analytics, Inc. and through a contract with LaRC was configured to represent a pilot flying a transport airplane. A version of the HOS program runs in batch mode on LaRC's (60-bit-word) central computer system. This document provides a guide for using the program and describes in some detail the assortment of files used during its operation.

Bogart, E. H.↗

Best Practices for Crash Modeling and Simulation

Aviation safety can be greatly enhanced by the expeditious use of computer simulations of crash impact. Unlike automotive impact testing, which is now routine, experimental crash tests of even small aircraft are expensive and complex due to the high cost of the aircraft and the myriad of crash impact conditions that must be considered. Ultimately, the goal is to utilize full-scale crash simulations of aircraft for design evaluation and certification. The objective of this publication is to describe "best practices" for modeling aircraft impact using explicit nonlinear dynamic finite element codes such as LS-DYNA, DYNA3D, and MSC.Dytran. Although "best practices" is somewhat relative, it is hoped that the authors' experience will help others to avoid some of the common pitfalls in modeling that are not documented in one single publication. In addition, a discussion of experimental data analysis, digital filtering, and test-analysis correlation is provided. Finally, some examples of aircraft crash simulations are described in several appendices following the main report.

Fasanella, Edwin L.↗

Deciphering the R Ψ m correlator in search for the chiral magnetic effect in relativistic heavy ion collisions

The chiral magnetic effect (CME) is extensively studied in heavy-ion collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). An azimuthal correlator called R Ψ m was proposed to measure the CME. By observing the same R Ψ 2 and R Ψ 3 (convex) distributions from A Multi-Phase Transport (ampt) model, by contrasting data and model as well as large and small systems and by event shape engineering (ESE), a recent preprint (arXiv:2006.04251v1) from STAR suggests that the R Ψ m observable is sensitive to the CME signal and relatively insensitive to backgrounds, and their Au + Au data are inconsistent with known background contributions. Our ampt results contradict those from a recent preprint by STAR. Our toy model simulations demonstrate that R Ψ 2 is sensitive to both the CME signal and physics backgrounds. Toy model simulations of isobar collisions show similar centrality dependence and magnitudes for the relative R Ψ 2 strengths as well as the relative Δ γ strengths. We conclude that R Ψ 2 and the inclusive Δ γ are essentially the same.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electrolytic hydrogen production: An analysis and review

The thermodynamics of water electrolysis cells is presented, followed by a review of current and future technology of commercial cells. The irreversibilities involved are analyzed and the resulting equations assembled into a computer simulation model of electrolysis cell efficiency. The model is tested by comparing predictions based on the model to actual commercial cell performance, and a parametric investigation of operating conditions is performed. Finally, the simulation model is applied to a study of electrolysis cell dynamics through consideration of an ideal pulsed electrolyzer.

Evangelista, J.↗

An integrative approach to space-flight physiology using systems analysis and mathematical simulation

An approach was developed to aid in the integration of many of the biomedical findings of space flight, using systems analysis. The mathematical tools used in accomplishing this task include an automated data base, a biostatistical and data analysis system, and a wide variety of mathematical simulation models of physiological systems. A keystone of this effort was the evaluation of physiological hypotheses using the simulation models and the prediction of the consequences of these hypotheses on many physiological quantities, some of which were not amenable to direct measurement. This approach led to improvements in the model, refinements of the hypotheses, a tentative integrated hypothesis for adaptation to weightlessness, and specific recommendations for new flight experiments.

Leonard, J. I.↗

Current Status and Application of Hazard Definition Technology

A research is performed: to define wake non-encounter & hazard, to provide requirements for sensors, and to obtain input from the user community. This research includes: validating wake encounter simulation models, establishing a metric to quantify the upset potential of a wake encounter, applying hazard metric and simulation models to the commercial fleet for development of candidate acceptable encounter limits, and applying technology to near term problems to evaluate current status of technology. The following lessons are learned from this project: technology is not adequate to determine absolute spacing requirements; time, not distance, determines the duration of the wake hazard; Optimum standards depend on the traffic; Wing span is an important parameter for characterizing both generator and follower; and Short span "biz jets" are easily rolled.

Greene, George C.↗

EFQPSK Versus CERN: A Comparative Study

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

Borah, Deva K.↗

Ground and Range Operations for a Heavy-Lift Vehicle: Preliminary Thoughts

This paper discusses the ground and range operations for a Shuttle derived Heavy-Lift Vehicle being launched from the Kennedy Space Center on the Eastern range. Comparisons will be made between the Shuttle and a heavy lift configuration (SLS-ETF MPCV April 2011) by contrasting their subsystems. The analysis will also describe a simulation configuration with the potential to be utilized for heavy lift vehicle processing/range simulation modeling and the development of decision-making systems utilized by the range. In addition, a simple simulation model is used to provide the required critical thinking foundations for this preliminary analysis.

Rabelo, Luis↗

Current Status of Post-combustor Trace Chemistry Modeling and Simulation at NASA Glenn Research Center

The overall objective of the current effort at NASA GRC is to evaluate, develop, and apply methodologies suitable for modeling intra-engine trace chemical changes over post combustor flow path relevant to the pollutant emissions from aircraft engines. At the present time, the focus is the high pressure turbine environment. At first, the trace chemistry model of CNEWT were implemented into GLENN-HT as well as NCC. Then, CNEWT, CGLENN-HT, and NCC were applied to the trace species evolution in a cascade of Cambridge University's No. 2 rotor and in a turbine vane passage. In general, the results from these different codes provide similar features. However, the details of some of the quantities of interest can be sensitive to the differences of these codes. This report summaries the implementation effort and presents the comparison of the No. 2 rotor results obtained from these different codes. The comparison of the turbine vane passage results is reported elsewhere. In addition to the implementation of trace chemistry model into existing CFD codes, several pre/post-processing tools that can handle the manipulations of the geometry, the unstructured and structured grids as well as the CFD solutions also have been enhanced and seamlessly tied with NCC, CGLENN-HT, and CNEWT. Thus, a complete CFD package consisting of pre/post-processing tools and flow solvers suitable for post-combustor intra-engine trace chemistry study is assembled.

Wey, Thomas↗

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗