Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATION”

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,243 records · Page 69

Computational Simulations of Convergent Nozzles for the AIAA 1st Propulsion Aerodynamics Workshop

Computational Fluid Dynamics (CFD) simulations were completed for a series of convergent nozzles in participation of the American Institute of Aeronautics and Astronautics (AIAA) 1st Propulsion Aerodynamics Workshop. The simulations were performed using the Wind-US flow solver. Discharge and thrust coefficients were computed for four axisymmetric nozzles with nozzle pressure ratios (NPR) ranging from 1.4 to 7.0. The computed discharge coefficients showed excellent agreement with available experimental data; the computed thrust coefficients captured trends observed in the experimental data, but over-predicted the thrust coefficient by 0.25 to 1.0 percent. Sonic lines were computed for cases with NPR >= 2.0 and agreed well with experimental data for NPR >= 2.5. Simulations were also performed for a 25 deg. conic nozzle bifurcated by a flat plate at NPR = 4.0. The jet plume shock structure was compared with and without the splitter plate to the experimental data. The Wind-US simulations predicted the shock structure well, though lack of grid resolution in the plume reduced the sharpness of the shock waves. Unsteady Reynolds-Averaged Navier-Stokes (URANS) simulations and Detached Eddy Simulations (DES) were performed at NPR = 1.6 for the 25 deg conic nozzle with splitter plate. The simulations predicted vortex shedding from the trailing edge of the splitter plate. However, the vortices of URANS and DES solutions appeared to dissipate earlier than observed experimentally. It is believed that a lack of grid resolution in the region of the vortex shedding may have caused the vortices to break down too soon

DES↗

Generating Inviscid and Viscous Fluid-Flow Simulations over an Aircraft Surface Using a Fluid-Flow Mesh

Fluid-flow simulation over a computer-generated aircraft surface is generated using inviscid and viscous simulations. A fluid-flow mesh of fluid cells is obtained. At least one inviscid fluid property for the fluid cells is determined using an inviscid fluid simulation that does not simulate fluid viscous effects. A set of intersecting fluid cells that intersects the aircraft surface are identified. One surface mesh polygon of the surface mesh is identified for each intersecting fluid cell. A boundary-layer prediction point for each identified surface mesh polygon is determined. At least one boundary-layer fluid property for each boundary-layer prediction point is determined using the at least one inviscid fluid property of the corresponding intersecting fluid cell and a boundary-layer simulation that simulates fluid viscous effects. At least one updated fluid property for at least one fluid cell is determined using the at least one boundary-layer fluid property and the inviscid fluid simulation.

Rodriguez, David L.↗

Generating Inviscid and Viscous Fluid Flow Simulations over a Surface Using a Quasi-simultaneous Technique

A fluid-flow simulation over a computer-generated surface is generated using a quasi-simultaneous technique. The simulation includes a fluid-flow mesh of inviscid and boundary-layer fluid cells. An initial fluid property for an inviscid fluid cell is determined using an inviscid fluid simulation that does not simulate fluid viscous effects. An initial boundary-layer fluid property a boundary-layer fluid cell is determined using the initial fluid property and a viscous fluid simulation that simulates fluid viscous effects. An updated boundary-layer fluid property is determined for the boundary-layer fluid cell using the initial fluid property, initial boundary-layer fluid property, and an interaction law. The interaction law approximates the inviscid fluid simulation using a matrix of aerodynamic influence coefficients computed using a two-dimensional surface panel technique and a fluid-property vector. An updated fluid property is determined for the inviscid fluid cell using the updated boundary-layer fluid property.

Sturdza, Peter↗

The Trick Simulation Toolkit: A NASA/Opensource Framework for Running Time Based Physics Models

The Trick Simulation Toolkit is a simulation development environment used to create high fidelity training and engineering simulations at the NASA Johnson Space Center and many other NASA facilities. Its purpose is to generate a simulation executable from a collection of user-supplied models and a simulation definition file. For each Trick-based simulation, Trick automatically provides job scheduling, numerical integration, the ability to write and restore human readable checkpoints, data recording, interactive variable manipulation, a run-time interpreter, and many other commonly needed capabilities. This allows simulation developers to concentrate on their domain expertise and the algorithms and equations of their models. Also included in Trick are tools for plotting recorded data and various other supporting utilities and libraries. Trick is written in C/C++ and Java and supports both Linux and MacOSX computer operating systems. This paper describes Trick's design and use at NASA Johnson Space Center.

Penn, John M.↗

ISS Double-Gimbaled CMG Subsystem Simulation Using the Agile Development Method

This paper presents an evolutionary approach in simulating a cluster of 4 Control Moment Gyros (CMG) on the International Space Station (ISS) using a common sense approach (the agile development method) for concurrent mathematical modeling and simulation of the CMG subsystem. This simulation is part of Training systems for the 21st Century simulator which will provide training for crew members, instructors, and flight controllers. The basic idea of how the CMGs on the space station are used for its non-propulsive attitude control is briefly explained to set up the context for simulating a CMG subsystem. Next different reference frames and the detailed equations of motion (EOM) for multiple double-gimbal variable-speed control moment gyroscopes (DGVs) are presented. Fixing some of the terms in the EOM becomes the special case EOM for ISS's double-gimbaled fixed speed CMGs. CMG simulation development using the agile development method is presented in which customer's requirements and solutions evolve through iterative analysis, design, coding, unit testing and acceptance testing. At the end of the iteration a set of features implemented in that iteration are demonstrated to the flight controllers thus creating a short feedback loop and helping in creating adaptive development cycles. The unified modeling language (UML) tool is used in illustrating the user stories, class designs and sequence diagrams. This incremental development approach of mathematical modeling and simulating the CMG subsystem involved the development team and the customer early on, thus improving the quality of the working CMG system in each iteration and helping the team to accurately predict the cost, schedule and delivery of the software.

Inampudi, Ravi↗

Assessment of the Impacts of ACLS on the ISS Life Support System Using Dynamic Simulations in V-HAB

The Advanced Closed Loop System (ACLS) is currently under development by Airbus Defense and Space and is slated for launch to the International Space Station (ISS) in 2017. The addition of new hardware into an already complex system such as the ISS life support system (LSS) always poses operational risks. It is therefore important to understand the impacts ACLS will have on the existing systems to ensure smooth operations for the ISS. This analysis can be done by using dynamic computer simulations and one possible tool for such a simulation is the Virtual Habitat (V-HAB). Based on MATLAB, V-HAB has been under development at the Institute of Astronautics of the Technical University of Munich (TUM) since 2004 and in the past has been successfully used to simulate the ISS life support systems. The existing V-HAB ISS simulation model treated the interior volume of the space station as one large, ideally-stirred container. This model was improved to allow the calculation of the atmospheric composition inside individual modules of the ISS by splitting it into twelve distinct volumes. The virtual volumes are connected by a simulation of the inter-module ventilation flows. This allows for a combined simulation of the LSS hardware and the atmospheric composition aboard the ISS. A dynamic model of ACLS is added to the ISS Simulation and several different operating modes for both ACLS and the existing ISS life support systems are studied and the impacts of ACLS on the rest of the system are determined. The results suggest that the US, Russian and ACLS CO2 systems can operate at the same time without impeding each other. Furthermore, based on the results of this analysis, the US and ACLS Sabatier systems can be operated in parallel as well to a achieve a very low CO2 concentration in the cabin atmosphere.

Putz, Daniel↗

Assessment of the Impacts of ACLS on the ISS Life Support System using Dynamic Simulations in V-HAB

The Advanced Closed Loop System (ACLS) is currently under development by Airbus Defense and Space and is slated for launch to the International Space Station (ISS) in 2017. The addition of new hardware into an already complex system such as the ISS life support system (LSS) always poses operational risks. It is therefore important to understand the impacts ACLS will have on the existing systems to ensure smooth operations for the ISS. This analysis can be done by using dynamic computer simulations and one possible tool for such a simulation is Virtual Habitat (V-HAB). Based on Matlab (Registered Trademark) V-HAB has been under development at the Institute of Astronautics of the Technical University Munich (TUM) since 2006 and in the past has been successfully used to simulate the ISS life support systems. The existing V-HAB ISS simulation model treated the interior volume of the space station as one large ideally-stirred container. This model was improved to allow the calculation of the atmospheric composition inside the individual modules of the ISS by splitting it into ten distinct volumes. The virtual volumes are connected by a simulation of the inter-module ventilation flows. This allows for a combined simulation of the LSS hardware and the atmospheric composition aboard the ISS. A dynamic model of ACLS is added to the ISS simulation and different operating modes for both ACLS and the existing ISS life support systems are studied to determine the impacts of ACLS on the rest of the system. The results suggest that the US, Russian and ACLS CO2 systems can operate at the same time without impeding each other. Furthermore, based on the results of this analysis, the US and ACLS Sabatier systems can be operated in parallel as well to achieve the highest possible CO2 recycling together with a low CO2 concentration.

Puetz, Daniel↗

Simulations of Cloud-Radiation Interaction Using Large-Scale Forcing Derived from the CINDY/DYNAMO Northern Sounding Array

The recently completed CINDY/DYNAMO field campaign observed two Madden-Julian oscillation (MJO) events in the equatorial Indian Ocean from October to December 2011. Prior work has indicated that the moist static energy anomalies in these events grew and were sustained to a significant extent by radiative feedbacks. We present here a study of radiative fluxes and clouds in a set of cloud-resolving simulations of these MJO events. The simulations are driven by the large-scale forcing data set derived from the DYNAMO northern sounding array observations, and carried out in a doubly periodic domain using the Weather Research and Forecasting (WRF) model. Simulated cloud properties and radiative fluxes are compared to those derived from the S-PolKa radar and satellite observations. To accommodate the uncertainty in simulated cloud microphysics, a number of single-moment (1M) and double-moment (2M) microphysical schemes in the WRF model are tested. The 1M schemes tend to underestimate radiative flux anomalies in the active phases of the MJO events, while the 2M schemes perform better, but can overestimate radiative flux anomalies. All the tested microphysics schemes exhibit biases in the shapes of the histograms of radiative fluxes and radar reflectivity. Histograms of radiative fluxes and brightness temperature indicate that radiative biases are not evenly distributed; the most significant bias occurs in rainy areas with OLR less than 150 W/ cu sq in the 2M schemes. Analysis of simulated radar reflectivities indicates that this radiative flux uncertainty is closely related to the simulated stratiform cloud coverage. Single-moment schemes underestimate stratiform cloudiness by a factor of 2, whereas 2M schemes simulate much more stratiform cloud.

Indian Ocean↗

Statistical Emulator for Expensive Classification Simulators

Expensive simulators prevent any kind of meaningful analysis to be performed on the phenomena they model. To get around this problem the concept of using a statistical emulator as a surrogate representation of the simulator was introduced in the 1980's. Presently, simulators have become more and more complex and as a result running a single example on these simulators is very expensive and can take days to weeks or even months. Many new techniques have been introduced, termed criteria, which sequentially select the next best (most informative to the emulator) point that should be run on the simulator. These criteria methods allow for the creation of an emulator with only a small number of simulator runs. We follow and extend this framework to expensive classification simulators.

Ross, Jerret↗

Simulation of Sweep-Jet Flow Control, Single Jet and Full Vertical Tail

This work is a simulation technology demonstrator, of sweep jet flow control used to suppress boundary layer separation and increase the maximum achievable load coefficients. A sweep jet is a discrete Coanda jet that oscillates in the plane parallel to an aerodynamic surface. It injects mass and momentum in the approximate streamwise direction. It also generates turbulent eddies at the oscillation frequency, which are typically large relative to the scales of boundary layer turbulence, and which augment mixing across the boundary layer to attack flow separation. Simulations of a fluidic oscillator, the sweep jet emerging from a nozzle downstream of the oscillator, and an array of sweep jets which suppresses boundary layer separation are performed. Simulation results are compared to data from a dedicated validation experiment of a single oscillator and its sweep jet, and from a wind tunnel test of a full-scale Boeing 757 vertical tail augmented with an array of sweep jets. A critical step in the work is the development of realistic time-dependent sweep jet inflow boundary conditions, derived from the results of the single-oscillator simulations, which create the sweep jets in the full-tail simulations. Simulations were performed using the computational fluid dynamics (CFD) solver Overow, with high-order spatial discretization and a range of turbulence modeling. Good results were obtained for all flows simulated, when suitable turbulence modeling was used.

flow control↗

A Reference Field for GCR Simulation and an LET-Based Implementation at NSRL

Exposure to galactic cosmic rays (GCR) on long duration deep space missions presents a serious health risk to astronauts, with large uncertainties connected to the biological response. In order to reduce the uncertainties and gain understanding about the basic mechanisms through which space radiation initiates cancer and other endpoints, radiobiology experiments are performed. Some of the accelerator facilities supporting such experiments have matured to a point where simulating the broad range of particles and energies characteristic of the GCR environment in a single experiment is feasible from a technology, usage, and cost perspective. In this work, several aspects of simulating the GCR environment in the laboratory are discussed. First, comparisons are made between direct simulation of the external, free space GCR field and simulation of the induced tissue field behind shielding. It is found that upper energy constraints at the NASA Space Radiation Laboratory (NSRL) limit the ability to simulate the external, free space field directly (i.e. shielding placed in the beam line in front of a biological target and exposed to a free space spectrum). Second, variation in the induced tissue field associated with shielding configuration and solar activity is addressed. It is found that the observed variation is within physical uncertainties, allowing a single reference field for deep space missions to be defined. Third, an approach for simulating the reference field at NSRL is presented. The approach allows for the linear energy transfer (LET) spectrum of the reference field to be approximately represented with discrete ion and energy beams and implicitly maintains a reasonably accurate charge spectrum (or, average quality factor). Drawbacks of the proposed methodology are discussed and weighed against alternative simulation strategies. The neutron component and track structure characteristics of the proposed strategy are discussed in this context.

Slaba, Tony C.↗

Manufactured Porous Ambient Surface Simulants

The planetary science decadal survey for 2013-2022 (Vision and Voyages, NRC 2011) has promoted mission concepts for sample acquisition from small solar system bodies. Numerous comet-sampling tools are in development to meet this standard. Manufactured Porous Ambient Surface Simulants (MPASS) materials provide an opportunity to simulate variable features at ambient temperatures and pressures to appropriately test potential sample acquisition systems for comets, asteroids, and planetary surfaces. The original "flavor" of MPASS materials is known as Manufactured Porous Ambient Comet Simulants (MPACS), which was developed in parallel with the development of the Biblade Comet Sampling System (Backes et al., in review). The current suite of MPACS materials was developed through research of the physical and mechanical properties of comets from past comet missions results and modeling efforts, coordination with the science community at the Jet Propulsion Laboratory and testing of a wide range of materials and formulations. These simulants were required to represent the physical and mechanical properties of cometary nuclei, based on the current understanding of the science community. Working with cryogenic simulants can be tedious and costly; thus MPACS is a suite of ambient simulants that yields a brittle failure mode similar to that of cryogenic icy materials. Here we describe our suite of comet simulants known as MPACS that will be used to test and validate the Biblade Comet Sampling System (Backes et al., in review).

Carey, Elizabeth M.↗

The Critical Role of the Routing Scheme in Simulating Peak River Discharge in Global Hydrological Models

Global hydrological models (GHMs) have been applied to assess global flood hazards, but their capacity to capture the timing and amplitude of peak river discharge which is crucial in flood simulations has traditionally not been the focus of examination. Here we evaluate to what degree the choice of river routing scheme affects simulations of peak discharge and may help to provide better agreement with observations. To this end we use runoff and discharge simulations of nine GHMs forced by observational climate data (1971-2010) within the ISIMIP2a (Inter-Sectoral Impact Model Intercomparison Project phase 2a) project. The runoff simulations were used as input for the global river routing model CaMa-Flood (Catchment-based Macro-scale Floodplain). The simulated daily discharge was compared to the discharge generated by each GHM using its native river routing scheme. For each GHM both versions of simulated discharge were compared to monthly and daily discharge observations from 1701 GRDC (Global Runoff Data Centre) stations as a benchmark. CaMa-Flood routing shows a general reduction of peak river discharge and a delay of about two to three weeks in its occurrence, likely induced by the buffering capacity of floodplain reservoirs. For a majority of river basins, discharge produced by CaMa-Flood resulted in a better agreement with observations. In particular, maximum daily discharge was adjusted, with a multi-model averaged reduction in bias over about two-thirds of the analysed basin area. The increase in agreement was obtained in both managed and near-natural basins. Overall, this study demonstrates the importance of routing scheme choice in peak discharge simulation, where CaMa-Flood routing accounts for floodplain storage and backwater effects that are not represented in most GHMs. Our study provides important hints that an explicit parameterisation of these processes may be essential in future impact studies.

peak river discharge↗

Gas-Granular Simulation Framework for Spacecraft Landing Plume-Surface Interaction and Debris Transport Analysis

The Gas-Granular Flow Solver (GGFS) multi-phase flow computational framework has been developed to enable simulations of particle flows complex extra-terrestrial regolith materials. Particle flows of interest include the damage of unprepared landing sites from rocket plume impingement on Moon, Mars, and asteroids. The flow solver implements an Eulerian-Eulerian two-fluid model with fluid representation of the gas phase and granular phase to avoid the need to model billions of particle interactions. The granular phase is modeled as an Eulerian fluid with constituent physics closure models derived from first-principle Discrete Element Model (DEM) particle interaction simulations that capture the complex, non-linear granular particle interaction effects. Granular phase constituent models have been developed and integrated that address the complex, non-linear granular material mechanics complexities resulting from both: the irregular, jagged particle shapes and poly-disperse mixture effects encountered in extra-terrestrial regolith, with lunar regolith as the extreme. The GGFS capabilities are being integrated into a proven NASA plume-surface interaction and debris transport simulation framework featuring the Loci/CHEM CFD program and Debris Transport Analysis (DTA) post-processing tools for applications in robotic and human Moon and Mars lander development. Integration of the three simulation tool components. Loci/CHEM, GGFS, and DTA, into a coordinated simulation framework will enable time-accurate spacecraft landing simulations that account for the alteration of the landing surface through plume-induced cratering and the resulting redirection of plume impingement flow and debris transport. Initial implementation of this simulation framework and application examples will be presented.

Liever, Peter A.↗

Air-Sea Interactions in a High-Resolution Ocean-Atmosphere Simulation

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interaction between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. We present here results from a coupled GEOS-MIT "nature run" simulation, wherein we have coupled a cubed-sphere-720 (~ 1/8) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere simulation to equivalent atmosphere-only and ocean-only simulations. A particular focus of the comparisons is the coupled versus uncoupled differences in interactions between Sea Surface Temperature (SST) and ocean surface wind. We discuss, in particular, a several-day mode of temporal variability in the SST-wind cycle and how it is represented in the different model simulations and in observationally-based products. A mechanism for the cycle, which is driven by SST-wind feedback, is proposed.

Strobach, Ehud↗

Simulation of Sweep-Jet Flow Control, Single Jet and Full Vertical Tail

This work is a simulation technology demonstrator, of sweep jets used to suppress boundary layer separation and increase maximum achievable load coefficients. A sweep jet is a discrete Coanda jet that oscillates in the plane parallel to an aerodynamic surface. It injects mass and momentum in the approximate stream wise direction. It also generate turbulent eddies at the oscillation frequency, which are typically large relative to boundary layer turbulence, and which augmenting mixing across the boundary layer to attack flow separation. Simulations of a fluidic oscillator, the sweep jet emerging from the oscillator, and the suppression of boundary layer separation by an array of sweep jets are performed. Simulation results are compared to data from a dedicated CFD validation experiment of a single oscillator and its sweep jet, and from a study of a full-scale Boeing 757 vertical tail augmented with an array of sweep jets.2, 20 A critical step in the work is the development of realistic time-dependent sweep-jet in flow boundary conditions, derived from the results of the single-oscillator simulations, which create the sweep jets in the full-tail simulations. Simulations were performed using the Over flow CFD solver, with high-order spatial discretization and a range of turbulence modeling. Good results were obtained for all flows simulated, when suitable turbulence modeling was used.

flow control↗

Evaluation of MERRA-2-Based Ozone Profile Simulations with the Global Ozonesonde Network

Chemical transport model (CTM) hindcasts of ozone (O3) are useful for filling in observational gaps and providing context for observed O3 variability and trends. We use global networks of ozonesonde stations to evaluate the O3 profiles in two simulations running versions of the NASA Global Modeling Initiative (GMI) chemical mechanism. Both simulations are tied to the NASA Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) meteorological reanalysis: 1) The GMI CTM, and 2) The MERRA-2 GMI Replay (M2 GMI). Both simulations start in 1980, and are compared against >50,000 ozonesonde profiles from 37 global stations from the tropics to the poles. The comparisons allow us to evaluate how the Replay technique affects modeled O3 distribution, how an updated chemical mechanism in the GMI CTM affects simulated tropospheric O3 amounts, and how observed O3 distributions compare to the full set of model output. In general, M2 GMI O3 is ~10% higher than in the GMI CTM, and shows global near-surface and tropical upper troposphere/lower stratosphere (UT/LS) high biases. The updated chemical mechanism in the GMI CTM reduces these high biases. Both simulations show similar negative biases in tropical free-tropospheric O3, especially during typical biomass burning seasons. The simulations are highly-correlated with ozonesonde measurements, particularly in the UT/LS (r > 0.8), showing the ability of MERRA-2 to capture tropopause height variations. Both simulations show improved correlations with ozonesonde data and smaller O3 biases in recent years. We expect to use the sonde/model comparisons to diagnose causes of disagreement and to gauge the feasibility of calculating multidecadal O3 trends from the model output.

Stauffer, Ryan↗

Evaluation of MERRA-2-Based Ozone Profile Simulations with the Global Ozonesonde Network

Chemical transport model (CTM) hindcasts of ozone (O3) are useful for filling in observational gaps and providing context for observed O3 variability and trends. We use global networks of ozonesonde stations to evaluate the O3 profiles in two simulations running versions of the NASA Global Modeling Initiative (GMI) chemical mechanism. Both simulations are tied to the NASA Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) meteorological reanalysis: 1) The GMI CTM, and 2) The MERRA-2 GMI Replay (M2 GMI). Both simulations start in 1980, and are compared against >50,000 ozonesonde profiles from 37 global stations from the tropics to the poles. The comparisons allow us to evaluate how the Replay technique affects modeled O3 distribution, how an updated chemical mechanism in the GMI CTM affects simulated tropospheric O3 amounts, and how observed O3 distributions compare to the full set of model output. In general, M2 GMI O3 is ~10% higher than in the GMI CTM, and shows global near-surface and tropical upper troposphere/lower stratosphere (UT/LS) high biases. The updated chemical mechanism in the GMI CTM reduces these high biases. Both simulations show similar negative biases in tropical free-tropospheric O3, especially during typical biomass burning seasons. The simulations are highly-correlated with ozonesonde measurements, particularly in the UT/LS (r > 0.8), showing the ability of MERRA-2 to capture tropopause height variations. Both simulations show improved correlations with ozonesonde data and smaller O3 biases in recent years. We expect to use the sonde/model comparisons to diagnose causes of disagreement and to gauge the feasibility of calculating multidecadal O3 trends from the model output.

Stauffer, Ryan M.↗