Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational modeling and 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 775 records · Page 43

CFD Simulations of the Space Launch System Ascent Aerodynamics and Booster Separation

This paper presents details of Computational Fluid Dynamic modeling of the Space Launch System during ascent. The primary focus of the paper is the flow simulation of the vehicle during ascent using the Overflow Navier-Stokes code. Computations of 739 first-stage flight conditions covering a range of Mach numbers, angles of attack, and roll angles were computed. The overset grid system contained 375 million grid points, and over 28 mil- lion CPU hours were used in the simulations. The simulations were run on the Pleiades supercomputer at the NASA Advanced Supercomputer Center at Ames Research Center. The data products from this work include integrated line-loads, surface pressure coefficients, venting pressures, and protuberance air-loads. Detailed comparisons were made of the aerodynamic performance predicted by Overflow and the wind-tunnel derived aero- dynamic database. A small number of the cases were run with two different turbulence models and with two differencing schemes. These results were used to quantify the sensitivity to the choice of the turbulence model and to the differencing scheme. The paper also introduces an effort to use the inviscid, unstructured Cartesian solver Cart3D to compute the aerodynamics during booster separation. Adaptive mesh refinement is being used to enable accurate simulations of sixteen booster-separation-motor plumes. The use of this tool is explored in preparation for building a booster-separation aerodynamic database.

Ascent Aerodynamics↗

ETARA PC version 3.3 user's guide: Reliability, availability, maintainability simulation model

A user's manual describing an interactive, menu-driven, personal computer based Monte Carlo reliability, availability, and maintainability simulation program called event time availability reliability (ETARA) is discussed. Given a reliability block diagram representation of a system, ETARA simulates the behavior of the system over a specified period of time using Monte Carlo methods to generate block failure and repair intervals as a function of exponential and/or Weibull distributions. Availability parameters such as equivalent availability, state availability (percentage of time as a particular output state capability), continuous state duration and number of state occurrences can be calculated. Initial spares allotment and spares replenishment on a resupply cycle can be simulated. The number of block failures are tabulated both individually and by block type, as well as total downtime, repair time, and time waiting for spares. Also, maintenance man-hours per year and system reliability, with or without repair, at or above a particular output capability can be calculated over a cumulative period of time or at specific points in time.

Hoffman, David J.↗

Approximate simulation model for analysis and optimization in engineering system design

Computational support of the engineering design process routinely requires mathematical models of behavior to inform designers of the system response to external stimuli. However, designers also need to know the effect of the changes in design variable values on the system behavior. For large engineering systems, the conventional way of evaluating these effects by repetitive simulation of behavior for perturbed variables is impractical because of excessive cost and inadequate accuracy. An alternative is described based on recently developed system sensitivity analysis that is combined with extrapolation to form a model of design. This design model is complementary to the model of behavior and capable of direct simulation of the effects of design variable changes.

Sobieszczanski-Sobieski, Jaroslaw↗

Computational Non-Destructive Evaluation Improving Ultrasonic Interrogation of Complex Geometry Composite Parts

Finite element simulation was employed in modeling the ultrasound (UT) pressure pulse propagation through a coupled liquid-composite medium to reproduce experimental data. From the simulation point of view, the proposed approach is challenging when there is a large simulation domain. For example, it is shown that a sub-micron wavelength of an ultrasonic wave requires a mesh size of several microns and this in turn requires significant computational resources, as well as special care in modeling. Some of the simulation results are presented considering that such modeling should reproduce experimental data for a healthy and faulty composite structure with complex geometry. Many possible experimental setups are simulated to demonstrate the non-destructive testing technique. This setup includes the generation of pressure pulse propagating through the tested composite plate and possible scattering by discontinuities (area of different impedance) that may be present in the panel. This scattered pulse together with the baseline pressure pulse generates a signature on the probe element which can be used to locate the position of defects in the structures.

modeling↗

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗

Accelerating Climate Simulations Through Hybrid Computing

Unconventional multi-core processors (e.g., IBM Cell B/E and NYIDIDA GPU) have emerged as accelerators in climate simulation. However, climate models typically run on parallel computers with conventional processors (e.g., Intel and AMD) using MPI. Connecting accelerators to this architecture efficiently and easily becomes a critical issue. When using MPI for connection, we identified two challenges: (1) identical MPI implementation is required in both systems, and; (2) existing MPI code must be modified to accommodate the accelerators. In response, we have extended and deployed IBM Dynamic Application Virtualization (DAV) in a hybrid computing prototype system (one blade with two Intel quad-core processors, two IBM QS22 Cell blades, connected with Infiniband), allowing for seamlessly offloading compute-intensive functions to remote, heterogeneous accelerators in a scalable, load-balanced manner. Currently, a climate solar radiation model running with multiple MPI processes has been offloaded to multiple Cell blades with approx.10% network overhead.

Zhou, Shujia↗

Propulsion Simulations with the Unstructured-Grid CFD Tool TetrUSS

A computational investigation has been completed to assess the capability of the NASA Tetrahedral Unstructured Software System (TetrUSS) for simulation of exhaust nozzle flows. Three configurations were chosen for this study: (1) a fluidic jet effects model, (2) an isolated nacelle with a supersonic cruise nozzle, and (3) a fluidic pitchthrust- vectoring nozzle. These configurations were chosen because existing data provided a means for measuring the ability of the TetrUSS flow solver USM3D for simulating complex nozzle flows. Fluidic jet effects model simulations were compared with structured-grid CFD (computational fluid dynamics) data at Mach numbers from 0.3 to 1.2 at nozzle pressure ratios up to 7.2. Simulations of an isolated nacelle with a supersonic cruise nozzle were compared with wind tunnel experimental data and structured-grid CFD data at Mach numbers of 0.9 and 1.2, with a nozzle pressure ratio of 5. Fluidic pitch-thrust-vectoring nozzle simulations were compared with static experimental data and structured-grid CFD data at static freestream conditions and nozzle pressure ratios from 3 to 10. A fluidic injection case was computed with the third configuration at a nozzle pressure ratio of 4.6 and a secondary pressure ratio of 0.7. Results indicate that USM3D with the S-A turbulence model provides accurate exhaust nozzle simulations at on-design conditions, but does not predict internal shock location at overexpanded conditions or pressure recovery along a boattail at transonic conditions.

Deere, Karen A.↗

Uncertainty Reduction With Multi-Model Monte Carlo for Crystal Plasticity Simulations of Additively Manufactured Metals

In this work, multi-model Monte Carlo estimators are developed to reduce uncertainty in quantities of interest (QoIs) extracted from crystal plasticity simulations of additively manufactured (AM) metals. A significant concern in AM parts is uncertainty in mechanical properties caused in part by complex microstructures that arise from the AM process. Quantifying uncertainty in microstructure-sensitive behavior using experiments alone is costly, especially when mechanical allowables must be established. Quantitative relationships among microstructure, micromechanical metrics like slip accumulation, crack initiation, and failure are also difficult to capture with limited experiments. Crystal plasticity material models instead enable computational prediction of micromechanical stress and strain fields given a discretized microstructure. However, high-fidelity finely discretized crystal plasticity simulations are computationally expensive, while lower-fidelity models are less accurate and generally biased, making uncertainty quantification and reduction computationally difficult as well. Multi-model Monte Carlo methods leverage correlations between high- and low-fidelity models to produce unbiased estimators for QoIs with reduced uncertainty relative to standard Monte Carlo. Crystal plasticity QoIs considered in this work include yield strength and the mean and extreme values of micromechanical fields that are relevant to crack initiation. Multi-model Monte Carlo estimators are developed for each individual QoI and several groups of QoIs. The results of this work establish relationships among model correlations, sample allocation, and uncertainty reduction for different combinations of QoIs and demonstrate a trend of less uncertainty reduction as QoIs become more sensitive to local microstructure. Limitations from using pilot samples to estimate model covariances and train low-fidelity models are also addressed. The uncertainty reduction achieved by multi-model Monte Carlo is an important step toward using computational mechanics models to predict microstructure-sensitive crack initiation and failure in AM parts.

uncertainty quantification↗

Users guide for guidance and control Launch and Abort Simulation for Spacecraft (LASS), volume 1

The mathematical models and computer program which are used to implement LASS are described. The computer program provides for a simulation of boost to orbit and abort capability from boost trajectories to a prescribed target. The abort target provides a decision point for engine shutdown from which the vehicle coasts to the vicinity of the selected abort recovery site. The simulation is a six degree of freedom simulation describing a rigid body. The vehicle is influenced by forces and moments from nondistributed aerodynamics. An adaptive autopilot is provided to control vehicle attitudes during powered and unpowered flight. A conventional autopilot is provided for study of vehicle during powered flight.

Havig, T. F.↗

Photoactivation of the manganese catalyst of O2 evolution. II A two-quantum mechanism.

A kinetic analysis is presented of the data obtained in previously described kinetic experiments relating to the photoactivation of the Mn complex required for photosynthetic O2 evolution. The results obtained from the computer simulation of the derived kinetic model compare favorably with the experimentally derived values. They also suggest that photoactivation can be rationalized in terms of a two-quantum process.

Radmer, R.↗

Dynamic modeling of spacecraft in a collisionless plasma

A new computational model is described which can simulate the charging of complex geometrical objects in three dimensions. Two sample calculations are presented. In the first problem, the capacitance to infinity of a complex object similar to a satellite with solar array paddles is calculated. The second problem concerns the dynamical charging of a conducting cube partially covered with a thin dielectric film. In this calculation, the photoemission results in differential charging of the object.

Katz, I.↗

The inversion of stratospheric aerosol and ozone vertical profiles from spacecraft solar extinction measurements

The inversions of multi-channel solar extinction measurements have been analyzed for the 0.35-1.0 micron wavelength region to retrieve stratospheric aerosol and ozone vertical profiles using both the constrained linear inversion scheme and the iterative scheme. The inversions of the multi-wavelength solar extinction data obtained from spacecraft have been analyzed based on the inversion of computer simulated data using various atmospheric models with differing amounts of aerosol and ozone in the stratosphere. The sensitivities of the inversion schemes to different experimental errors are discussed in terms of accuracy and resolution of the retrieved profiles.

Chu, W. P.↗

Pilot decision making in a computer-aided flight management situation

An experimental representation of a computer-aided multi-task flight management situation has been developed. A computer aiding program was implemented to serve as a back-up decision maker. An experiment was conducted with a balanced design of several subject runs for different workload levels. This was achieved using three levels of subsystem event arrival rates, three levels of control task involvement, and three levels of availability of computer aiding. Experimental results compared quite favorably with those from a computer simulation which employed a queueing model. It was shown that the aiding had enhanced system performance as well as subjective ratings, and that the adaptive aiding policy further reduced subsystem delay.

Chu, Y. Y.↗

Two- and Three-Dimensional Galaxy Simulations

Dynamics of galaxies with up to 100,000 stars are simulated in set of computer programs available from COSMIC. In these models, large number of simulation stars move in combined gravitational field produced by stars themselves (self-consistent field) or by external matter (fixed field).

Hohl, F.↗

Computer simulation analysis of the behavior of renal-regulating hormones during hypogravic stress

A computer simulation of a mathematical circulation model is used to study the alterations of body fluids and their electrolyte composition that occur in weightlessness. The behavior of the renal-regulating hormones which control these alterations is compared in simulations of several one-g analogs of weightlessness and space flight. It is shown that the renal-regulating hormones represent a tightly coupled system that responds acutely to volume disturbances and chronically to electrolyte disturbances. During hypogravic conditions these responses lead to an initial suppression of hormone levels and a long-term effect which varies depending on metabolic factors that can alter the plasma electrolytes. In addition, it is found that if pressure effects normalize rapidly, a transition phase may exist which leads to a dynamic multiphasic endocrine response.

Leonard, J. I.↗

AC motor and generator requirements for isolated WECS

After surveying electrically driven loads used on productive farms, the investigators chose three pumps for testing at voltages and frequencies far outside the normal operating range. These loads extract and circulate water and move heat via air, and all are critical to farm productivity. The object was to determine the envelope of supply voltage and frequency over which these loads would operate stably for time intervals under 1 hour. This information is among that needed to determine the feasibility of supplying critical loads, in case of a utility outage, from a wind driven alternator whose output voltage and frequency will vary dramatically in most continental wind regimes. Other related work is surveyed. The salient features and limitations of the test configurations used and the data reduction are described. The development of simulation models suitable for a small computer are outlined. The results are primarily displayed on the voltage frequency plane with the general conclusion that the particular pump models considered will operate over the range of 50 to 90 Hz and a voltage band which starts below rated, decreases as frequency decreases, and is limited on the high side by excessive motor heating. For example, centrifugal pump operating voltage ranges as extensive .4 to 1.4 appear possible. Particular problems with starting, stalling due to lack of motor torque, high speed cavitation, and likely overheating are addressed in a listing of required properties for wind driven alternators and their controllers needed for use in the isolated or stand alone configuration considered.

Park, G. L.↗