Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 559 records · Page 31

A Predictive GNSS Performance Monitor for Autonomous Air Vehicles in Urban Environments

The emergence and development of advanced technologies and vehicle types has created a growing demand for the introduction of new forms of flight operations. These new and increasingly complex operational paradigms such as Advanced and Urban Air Mobility (AAM/UAM) present regulatory authorities and the aviation community with several design and implementation challenges – particularly for highly autonomous vehicles. An overarching and daunting task is finding methods to integrate these emerging operations without compromising safety or disrupting traditional airspace operations. A shift toward a more predictive risk mitigation capability becomes critical to meet this challenge. This paper focus on the development and testing of a prognostic service aimed at estimating the quality of Global Navigation Satellite System (GNSS) performance for an autonomous aircraft in complex environments. The intent of this function is to proactively reduce a flight operations risk of exposure to states that may induce poor or unacceptable navigation system performance by factoring in estimates of GNSS quality into pre-flight and/or in-flight route planning. Methodologies for producing quality estimates are specified and results are provided for selected simulation and flight test cases.

GNSS↗

Perturbation Modeling and Navigation Tuning for a Crewed Station in Near Rectilinear Halo Orbit

NASA's Gateway program will build a crew-tended station in an Earth-Moon Near Rectilinear Halo Orbit (NRHO) to support deep space missions to the lunar surface and heliocentric space. The station in the NRHO will be visited by crewed vehicles including lunar landers that are high mass relative to Gateway. Docking events with large vehicles induce significant thruster plume and docking contact velocity perturbations. Crewed vehicles induce significant venting perturbations that can stress the navigation filter's ability to remain converged about an estimated state. The torques that arise from venting and lunar gravity gradient require reaction control system (RCS) thruster intervention, which themselves induce a velocity perturbation. The best method to model RCS thruster perturbations is currently in development. This paper simulates different RCS perturbation models, schedules, and an increasing magnitude to investigate the navigation sensitivity to different RCS perturbation models. From the RCS perturbation models simulated, a worst-case scenario is chosen for further scrutiny. The navigation filter is further tuned to handle the most disruptive perturbation modeling and the resulting performance is analyzed.

Navigation↗

Integrated Process-Structure-Property Simulations for Additive Manufacturing Using the Open-Source Materialite Package

The microstructure and properties of additively manufactured (AM) metals are strongly dependent on process conditions. Therefore, process-structure-property (PSP) simulations are a useful tool for exploring process parameter space, studying process variations, and quantifying uncertainty in material properties. However, integrating process-structure and structure-property simulations often involves connecting multiple software packages. Each package may use unique data structures and require substantial domain knowledge. This presentation demonstrates PSP simulation capabilities of Materialite, an open-source package developed at NASA Langley Research Center. Materialite simplifies model linkages by using a common data structure and model interface, enabling straightforward simulation across a PSP model chain. Physics-based models, including kinetic Monte Carlo and crystal plasticity, are implemented within the package. The model interface is also intended to simplify implementation of new models and enable integration with external simulation tools. Example use cases include uncertainty quantification with PSP models and GPU-accelerated powder bed fusion AM process models.

additive manufacturing↗

Astronomical imaging by pupil plane interferometry

Comparing rotational shear interferometry to standard speckle interferometry, it is found that it is easier in the first case to separate the atmospheric phases from the object transform phases. Phase closure and blind deconvolution should be directly applicable. Laboratory simulations were conducted to verify theoretical predictions and computer simulations for the phase closure case, and preliminary results show promise.

Ribak, Erez↗

Parallel ALLSPD-3D: Speeding Up Combustor Analysis Via Parallel Processing

The ALLSPD-3D Computational Fluid Dynamics code for reacting flow simulation was run on a set of benchmark test cases to determine its parallel efficiency. These test cases included non-reacting and reacting flow simulations with varying numbers of processors. Also, the tests explored the effects of scaling the simulation with the number of processors in addition to distributing a constant size problem over an increasing number of processors. The test cases were run on a cluster of IBM RS/6000 Model 590 workstations with ethernet and ATM networking plus a shared memory SGI Power Challenge L workstation. The results indicate that the network capabilities significantly influence the parallel efficiency, i.e., a shared memory machine is fastest and ATM networking provides acceptable performance. The limitations of ethernet greatly hamper the rapid calculation of flows using ALLSPD-3D.

Fricker, David M.↗

The Influences of Airmass Histories on Radical Species During POLARIS

The POLARIS mission focused on understanding the processes associated with the decrease of polar stratospheric ozone from spring to fall at high latitudes. This decrease is linked primarily to in situ photochemical destruction by reactive nitrogen species, NO and NO2, which also control other catalytic loss cycles. Steady state models have been used to test photochemistry and radical behavior but are not always adequate in simulating radical species observations. In some cases, air mass history can be important and trajectory models give an improved simulation of the radical species. Trajectory chemistry models, however, still consistently underestimate NO and NO2 abundances compared to measurements along the ER-2 flight track. The Goddard chemistry on trajectory model has been used to test updated rate constants for NO2 + OH, NO2 + O and OH + HNO3, key reactions that affect NO and NO2 abundances. We present comparisons between the modified Goddard chemistry on trajectory model, the JPL steady state model and observations from selected flights.

Pierson, James M.↗

Structure of the dayside reconnection layer in resistive MHD and hybrid models

Numerical simulations were performed to investigate the structure of the reconnection layer at the dayside magnetopause. Two typical cases are examined in detail; both are asymmetric in magnetic field and plasma density. In case 1, the guide fields in the magnetosheath and in the magnetosphere are set at zero and thus the tangential magnetic fields on the two sides of the initial current sheet are exactly antiparallel. In case 2, the angle between the tangential magnetic fields on the two sides of the initial current sheet is 145 deg. The results obtained from a resistive MHD model and from a hybrid model are found to be different. In the MHD simulation of case 1, a 2-4 intermediate shock is found to bound the reconnection layer on the magnetosheath side, while an Alfven wave pulse bounds the reconnection layer on the magnetospheric side. In case 2, it is found that a time-dependent intermediate shock (TDIS) bounds the reconnection layer on the magnetosheath side, with a slow expansion wave propagating behind. With the MHD simulations, in the general case in which the tangential magnetic fields on the two sides of the initial current sheet are not exactly antiparallel, a rotational discontinuity across which the tangential magnetic field rotates, a large angle is found to bound the reconnection layer on the magnetosheath side.

Lin, Y.↗

Some observations on mesh refinement schemes applied to shock wave phenomena

This workshop's double-wedge test problem is taken from one of a sequence of experiments which were performed in order to classify the various canonical interactions between a planar shock wave and a double wedge. Therefore to build up a reasonably broad picture of the performance of our mesh refinement algorithm we have simulated three of these experiments and not just the workshop case. Here, using the results from these simulations together with their experimental counterparts, we make some general observations concerning the development of mesh refinement schemes for shock wave phenomena.

Quirk, James J.↗

Distributed simulation, no special tools required

A tool kit of C language functions that can be linked with SIMSCRIPT programs to provide the data communication primitives necessary for distributed simulation is presented. A test case is discussed and some timing data are presented. Additionally some metrics, developed to determine the applicability of the server model decomposition for particular simulations, are discussed.

Paterra, Frank↗

Options and Uncertainties in Planetary Defense: Mission Planning and Vehicle Design for Flexible Response

This paper is part of an integrated study by NASA and the NNSA to quantitatively understand the response timeframe should a threatening Earth-impacting near-Earth object (NEO) be identified. The two realistic responses considered are the use of a spacecraft functioning as either a kinetic impactor or a nuclear explosive carrier to deflect the approaching NEO. The choice depends on the NEO size and mass, the available response time prior to Earth impact, and the various uncertainties. Whenever practical, the kinetic impactor is the preferred approach, but various factors, such as large uncertainties or short available response time, reduce the kinetic impactor's suitability and, ultimately, eliminate its sufficiency. Herein we examine response time and the activities that occur between the time when an NEO is recognized as being a sufficient threat to require a deflection and the time when the deflection impulse is applied to the NEO. To use a kinetic impactor for successful deflection of an NEO, it is essential to minimize the reaction time and maximize the time available for the impulse delivered to the NEO by the kinetic impactor to integrate forward in time to the eventual deflection of the NEO away from Earth impact. To shorten the response time, we develop tools to survey the profile of needed spacecraft launches and the possible mission payloads. We further present a vehicle design capable of either serving as a kinetic impactor, or, if the need arises, serving as a system to transport a nuclear explosive to the NEO. These results are generated by analyzing a specific case study in which the simulated Earth-impacting NEO is modeled very closely after the real NEO known as 101955 Bennu (1999 RQ36). Bennu was selected for our case study in part because it is the best-studied of the known NEOs. It is also the destination of NASA's OSIRIS-REx sample return mission, which is, at the time of this writing, enroute to Bennu following a September 2016 launch.

characterization↗

Options and Uncertainties in Planetary Defense: Mission Planning and Vehicle Design for Flexible Response

This paper is part of an integrated study by NASA and the NNSA to quantitatively understand the response timeframe should a threatening Earth-impacting near-Earth object (NEO) be identified. The two realistic responses considered are the use of a spacecraft functioning as either a kinetic impactor or a nuclear explosive carrier to deflect the approaching NEO. The choice depends on the NEO size and mass, the available response time prior to Earth impact, and the various uncertainties. Whenever practical, the kinetic impactor is the preferred approach, but various factors, such as large uncertainties or short available response time, reduce the kinetic impactor's suitability and, ultimately, eliminate its sufficiency. Herein we examine response time and the activities that occur between the time when an NEO is recognized as being a sufficient threat to require a deflection and the time when the deflection impulse is applied to the NEO. To use a kinetic impactor for successful deflection of an NEO, it is essential to minimize the reaction time and maximize the time available for the impulse delivered to the NEO by the kinetic impactor to integrate forward in time to the eventual deflection of the NEO away from Earth impact. To shorten the response time, we develop tools to survey the profile of needed spacecraft launches and the possible mission payloads. We further present a vehicle design capable of either serving as a kinetic impactor, or, if the need arises, serving as a system to transport a nuclear explosive to the NEO. These results are generated by analyzing a specific case study in which the simulated Earth-impacting NEO is modeled very closely after the real NEO known as 101955 Bennu (1999 RQ36). Bennu was selected for our case study in part because it is the best-studied of the known NEOs. It is also the destination of NASA's OSIRIS-REx sample return mission, which is, at the time of this writing, enroute to Bennu following a September 2016 launch.

Brent W. Barbee↗

GCSS Idealized Cirrus Model Comparison Project

The GCSS Working Group on Cirrus Cloud Systems (WG2) is conducting a systematic comparison and evaluation of cirrus cloud models. This fundamental activity seeks to support the improvement of models used for climate simulation and numerical weather prediction through assessment and improvement of the "process" models underlying parametric treatments of cirrus cloud processes in large-scale models. The WG2 Idealized Cirrus Model Comparison Project is an initial comparison of cirrus cloud simulations by a variety of cloud models for a series of idealized situations with relatively simple initial conditions and forcing. The models (16) represent the state-of-the-art and include 3-dimensional large eddy simulation (LES) models, two-dimensional cloud resolving models (CRMs), and single column model (SCM) versions of GCMs. The model microphysical components are similarly varied, ranging from single-moment bulk (relative humidity) schemes to fully size-resolved (bin) treatments where ice crystal growth is explicitly calculated. Radiative processes are included in the physics package of each model. The baseline simulations include "warm" and "cold" cirrus cases where cloud top initially occurs at about -47C and -66C, respectively. All simulations are for nighttime conditions (no solar radiation) where the cloud is generated in an ice supersaturated layer, about 1 km in depth, with an ice pseudoadiabatic thermal stratification (neutral). Continuing cloud formation is forced via an imposed diabatic cooling representing a 3 cm/s uplift over a 4-hour time span followed by a 2-hour dissipation stage with no cooling. Variations of these baseline cases include no-radiation and stable-thermal-stratification cases. Preliminary results indicated the great importance of ice crystal fallout in determining even the gross cloud characteristics, such as average vertically-integrated ice water path (IWP). Significant inter-model differences were found. Ice water fall speed is directly related to the shape of the particle size distribution and the habits of the ice crystal population, whether assumed or explicitly calculated. In order to isolate the fall speed effect from that of the associated ice crystal population, simulations were also performed where ice water fall speed was set to the same constant value everywhere in each model. Values of 20 and 60 cm/s were assumed. Current results of the project will be described and implications will be drawn. In particular, this exercise is found to strongly focus the definition of issues resulting in observed inter-model differences and to suggest possible strategies for observational validation of the models. The next step in this project is to perform similar comparisons for well observed case studies with sufficient high quality data to adequately define model initiation and forcing specifications and to support quantitative validation of the results.

Starr, David OC.↗

Influence of Ice-phase of Hydrometeors on Moist-Convection

Climate models often ignore the influence of ice-phase physics (IPP) of hydrometeors as a second order effect. This has also been true for McRAS (Microphysics of clouds with Relaxed Arakawa Schubert Scheme) developed by the authors. Recognizing that the temperature sounding is critical for moist-convection, and, that IPP would modify it, we investigated the influence of introducing IPP into McRAS coupled to FvGCM (finite volume General Circulation Model with NCAR physics). We analyzed three 3-yr long simulations; the first called Control Case, CC and had no IPP; the other two called Experiments El and E2 had IPP introduced with two different in-cloud freezing assumptions. Simulation El assumed that all hydrometeors remain liquid in the updraft and freeze upon detrainment. Simulation E2 invoked the in-cloud freezing of new condensate generated at subfreezing temperatures in the updraft while old cloud water continued to ascend as liquid. Upon detrainment, this cloud water also froze like in E1. With these assumptions, about 50% of hydrometeors froze in the tower and the rest froze in the anvil. However, in both El and E2, the frozen hydrometeors melted during fall at the first encounter of above freezing ambient temperature. Comparative analysis revealed that El simulated far more mid-level and far less deep clouds while E2 had modified deep and more mid-level clouds as compared to CC along with some major changes around the melt-level. We infer that IPP produced a more realistic response in E2. At the basic level, the results show that ice-phase processes influence convective detrainment at mid- and deep levels in accord with TOGAGOARE observations. The results suggest that IPP can help to mitigate less-than-observed mid-level and over-abundance of deep convective clouds in McRAS.

Sud, Y. C.↗

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗

Sensitivity of a Cloud-Resolving Model to the Bulk and Explicit Bin Microphysical Schemes: Validations with a PRE-STORM Case - Part 1

A cloud-resolving model is used to study sensitivities of two different microphysical schemes, one is the bulk type, and the other is an explicit bin scheme, in simulating a mid-latitude squall line case (PRE-STORM, June 10-11, 1985). Simulations using different microphysical schemes are compared with each other and also with the observations. Both the bulk and bin models reproduce the general features during the developing and mature stage of the system. The leading convective zone, the trailing stratiform region, the horizontal wind flow patterns, pressure perturbation associated with the storm dynamics, and the cool pool in front of the system all agree well with the observations. Both the observations and the bulk scheme simulation serve as validations for the newly incorporated bin scheme. However, it is also shown that, the bulk and bin simulations have distinct differences, most notably in the stratiform region. Weak convective cells exist in the stratiform region in the bulk simulation, but not in the bin simulation. These weak convective cells in the stratiform region are remnants of the previous stronger convections at the leading edge of the system. The bin simulation, on the other hand, has a horizontally homogeneous stratiform cloud structure, which agrees better with the observations. Preliminary examinations of the downdraft core strength, the potential temperature perturbation, and the evaporative cooling rate show that the differences between the bulk and bin models are due mainly to the stronger low-level evaporative cooling in convective zone simulated in the bulk model. Further quantitative analysis and sensitivity tests for this case using both the bulk and bin models will be presented in a companion paper.

Li, Xiao-Wen↗

Hierarchical nonlinear behavior of hot composite structures

Hierarchical computational procedures are described to simulate the multiple scale thermal/mechanical behavior of high temperature metal matrix composites (HT-MMC) in the following three broad areas: (1) behavior of HT-MMC's from micromechanics to laminate via METCAN (Metal Matrix Composite Analyzer), (2) tailoring of HT-MMC behavior for optimum specific performance via MMLT (Metal Matrix Laminate Tailoring), and (3) HT-MMC structural response for hot structural components via HITCAN (High Temperature Composite Analyzer). Representative results from each area are presented to illustrate the effectiveness of computational simulation procedures and accompanying computer codes. The sample case results show that METCAN can be used to simulate material behavior such as the entire creep span; MMLT can be used to concurrently tailor the fabrication process and the interphase layer for optimum performance such as minimum residual stresses; and HITCAN can be used to predict the structural behavior such as the deformed shape due to component fabrication. These codes constitute virtual portable desk-top test laboratories for characterizing HT-MMC laminates, tailoring the fabrication process, and qualifying structural components made from them.

Murthy, P. L. N.↗

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Bui, Trong T.↗