Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled 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 55 records · Page 3

Simulation of a Controlled Airfoil with Jets

Numerical simulations of a two-dimensional airfoil, controlled by an applied moment in pitch and an airfoil controlled by jets, were investigated. These simulations couple the Reynolds-averaged Navier-Stokes equations and Euler's equations of rigid body motion, with an active control system. Controllers for both systems were designed to track altitude commands and were evaluated by simulating a closed-loop altitude step response using the coupled system. The airfoil controlled by a pitching moment used an optimal state feedback controller. A closed-loop simulation, of the airfoil with an applied moment, showed that the trajectories compared very well with quasi-steady aerodynamic theory, providing a measure of validation. The airfoil with jets used a controller designed by robust control methods. A linear plant model for this system was identified using open-loop data generated by the nonlinear coupled system. A closed-loop simulation of the airfoil with jets, showed good tracking of an altitude command. This simulation also showed oscillations in the control input as a result of dynamics not accounted for in the control design. This research work demonstrates how computational fluid dynamics, coupled with rigid body dynamics, and a control law can be used to prototype control systems in problematic nonlinear flight regimes.

Allan, Brian G.↗

The Effects of Aircraft Wake Dynamics on Contrail Development

Results of large-eddy simulations of the development of young persistent ice contrails are presented, concentrating on the interactions between the aircraft wake dynamics and the ice cloud evolution over ages front a few seconds to approx. 30 min. The 3D unsteady evolution of the dispersing engine exhausts, trailing vortex pair interaction and breakup, and subsequent Brunt-Vaisala oscillations of the older wake plume are modeled in detail in high-resolution simulations, coupled with it bulk microphysics model for the contrail ice development. The simulations confirm that the early wake dynamics can have a strong influence on the properties of persistent contrails even at late times. The vortex dynamics are the primary determinant of the vertical extent of the contrail (until precipitate ton becomes significant): and this together with the local wind shear largely determines the horizontal extent. The ice density, ice crystal number density, and a conserved exhaust tracer all develop and disperse in different fashions from each other. The total ice crystal number can be significantly reduced due to adiabatic compression resulting from the downward motion of the vortex system, even for ambient conditions that are substantially supersaturated with respect to ice. The fraction of the initial ice crystals surviving, their spatial distribution and the ice mass distribution are all sensitive to the aircraft type, ambient humidity, assumed initial ice crystal number, and ambient turbulence conditions. There is a significant range of conditions for which a smaller transport such as a B737 produces as significant a persistent contrail as a larger transport such as a B747, even though the latter consumes almost five times as much fuel. The difficulties involved in trying to minimize persistent contrail production are discussed.

Lewellen, D. C.↗

Overview of the Predictive Simulation Capability Element of the Plume Surface Interaction Project

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The GCD PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the PSC Element of the PSI Project as well as providing descriptions of recent accomplishments and remaining work. The overall structure of the PSC Element is broken down into four areas of focus. The first area is the Prediction of Plume Flow in low pressure environments. The development approach taken is the augmentation of the existing production-mode computational fluid dynamics (CFD) tool Loci/Chem, with targeted extensions necessary to accurately model rarefied conditions found in both Martian and Lunar applications. Production readiness and validation of predictive capability are the major objectives of this task. The second area is the effect of mixed continuum/rarefied flow on crater development and ejecta sheets. A new CFD application, Loci/GGFS (Gas Granular Flow Solver), is being developed which implements an Eulerian/Eulerian two-phase model of gas- and soil-phases in order to simulate the soil erosion, crater formation, and soil ejecta transport in a fully coupled simulation. This task seeks to verify Loci/GGFS is production-ready as well as perform validation studies to determine the degree of predictive capability achieved by Loci/GGFS. The third area is focused on the details and extension of particle phase modeling of soil. In this task, Discrete Element Modeling (DEM) techniques are used to perform direct simulations of complex soil particles under the action of forcing similar to that to be cause by PSI. The simulation results are then used to construct closures to the Eulerian model of the soil phase used by Loci/GGFS. The fourth area is gas-particle interaction modeling. In this task, experiments are being conducted as well as detailed simulation results are being studied to further understand the complexities of gas-particle interactions in dilute, intermediate, and high soil volume fraction regimes. Improved models of particle drag and the particle turbulent kinetic energy (PTKE) resulting from the interaction of gas flows within particle clouds are the objective of this task.

Jeff West↗

ENSO Mechanisms

This current leading theory of the El Nino-Southern Oscillation (ENSO) phenomenon involves reflection of equatorial waves on the western ocean boundary to shift the ENSO phase from El Nino to La Nina and vice versa. However, recent satellite altimetry and In situ observations indicate that this theory has several flaws. These include imperfect equatorial wave reflection on the western ocean boundary, a maxima of the simulated coupled wind-sea surface temperature (SST) interaction located too far into the eastern equatorial Pacific, and an excessive role of thermocline displacement on SST in this region. We show that the central equatorial Pacific SST is fundamental in driving the EIVSO ocean-atmosphere coupled system, through the east-west displacement of the eastern edge of the warm pool, perfectly related with ENSO. Within the equatorial wave-guide, the dominance of surface zonal advection in these displacements is demonstrated with four different satellite and in situ data sets and three ocean models. This demonstration is supported by the evidence of a convergence of water masses into the eastern edge of the warm pool, resulting in a well-defined salinity front. All these results lead us to propose a notable modification of the leading delayed-action oscillator theory for the oscillatory nature of ENSO. Simulations with a linearized coupled ocean-atmosphere model result in 3-6 year ENSO-like oscillations, with many Of the variable model parameters found very close to their observed values. This simple model suggests that ocean processes ignored or underestimated in the delayed action oscillator theory, such as zonal current convergence, zonal advection of sea surface temperature and equatorial wave reflection from the eastern-ocean boundary, are fundamental to the development of ENSO, in particular to its manifestations in the central equatorial Pacific.

Picaut, Joel↗

Influence of Idealized Heterogeneity on Wet and Dry Planetary Boundary Layers Coupled to the Land Surface

We examine the influence of surface heterogeneity on boundary layers using a large-eddy simulation coupled to a land-surface model. Heterogeneity, imposed in strips varying from 2-30 km (1 less than lambda/z(sub i) less than 18), is found to dramatically alter the structure of the free convective boundary layer by inducing significant organized circulations. A conditional sampling technique, based on the scale of the surface heterogeneity (phase averaging), is used to identify and quantify the organized surface fluxes and motions in the atmospheric boundary layer. The impact of the organized motions on turbulent transport depends critically on the scale of the heterogeneity lambda, the boundary layer height zi and the initial moisture state of the boundary layer. Dynamical and scalar fields respond differently as the scale of the heterogeneity varies. Surface heterogeneity of scale 4 less than lamba/z(sub i) less than 9 induces the strongest organized flow fields (up, wp) while heterogeneity with smaller or larger lambda/z(sub i) induces little organized motion. However, the organized components of the scalar fields (virtual potential temperature and mixing ratio) grow continuously in magnitude and horizontal scale, as lambda/z(sub i) increases. For some cases, the organized motions can contribute nearly 100% of the total vertical moisture flux. Patch-induced fluxes are shown to dramatically impact point measurements that assume the time-average vertical velocity to be zero. The magnitude and sign of this impact depends on the location of the measurement within the region of heterogeneity.

Houser, Paul↗

Effects of Real-Time NASA Vegetation Data on Model Forecasts of Severe Weather

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed a Greenness Vegetation Fraction (GVF) dataset, which is updated daily using swaths of Normalized Difference Vegetation Index data from the Moderate Resolution Imaging Spectroradiometer (MODIS) data aboard the NASA-EOS Aqua and Terra satellites. NASA SPoRT started generating daily real-time GVF composites at 1-km resolution over the Continental United States beginning 1 June 2010. A companion poster presentation (Bell et al.) primarily focuses on impact results in an offline configuration of the Noah land surface model (LSM) for the 2010 warm season, comparing the SPoRT/MODIS GVF dataset to the current operational monthly climatology GVF available within the National Centers for Environmental Prediction (NCEP) and Weather Research and Forecasting (WRF) models. This paper/presentation primarily focuses on individual case studies of severe weather events to determine the impacts and possible improvements by using the real-time, high-resolution SPoRT-MODIS GVFs in place of the coarser-resolution NCEP climatological GVFs in model simulations. The NASA-Unified WRF (NU-WRF) modeling system is employed to conduct the sensitivity simulations of individual events. The NU-WRF is an integrated modeling system based on the Advanced Research WRF dynamical core that is designed to represents aerosol, cloud, precipitation, and land processes at satellite-resolved scales in a coupled simulation environment. For this experiment, the coupling between the NASA Land Information System (LIS) and the WRF model is utilized to measure the impacts of the daily SPoRT/MODIS versus the monthly NCEP climatology GVFs. First, a spin-up run of the LIS is integrated for two years using the Noah LSM to ensure that the land surface fields reach an equilibrium state on the 4-km grid mesh used. Next, the spin-up LIS is run in two separate modes beginning on 1 June 2010, one continuing with the climatology GVFs while the other uses the daily SPoRT/MODIS GVFs. Finally, snapshots of the LIS land surface fields are used to initialize two different simulations of the NU-WRF, one running with climatology LIS and GVFs, and the other running with experimental LIS and NASA/SPoRT GVFs. In this paper/presentation, case study results will be highlighted in regions with significant differences in GVF between the NCEP climatology and SPoRT product during severe weather episodes.

Case, Jonathan L.↗

Dynamical Simulation of Probabilities

It has been demonstrated that classical probabilities, and in particular, probabilistic Turing machine, can be simulated by combining chaos and non-Lipschitz dynamics, without utilization of any man-made devices(such as random number generators). Self-orgainizing properties of systems coupling simulated and calculated probabilities and their link to quantum computations are discussed. Special attention was focused upon coupled stochastic processes, defined in terms of conditional probabilities, for which joint probability does not exist. Simulations of quantum probabilities are also discussed.

number generators quantum probabilities↗

System Simulation by Recursive Feedback: Coupling a Set of Stand-Alone Subsystem Simulations

Conventional construction of digital dynamic system simulations often involves collecting differential equations that model each subsystem, arran g them to a standard form, and obtaining their numerical gin solution as a single coupled, total-system simultaneous set. Simulation by numerical coupling of independent stand-alone subsimulations is a fundamentally different approach that is attractive because, among other things, the architecture naturally facilitates high fidelity, broad scope, and discipline independence. Recursive feedback is defined and discussed as a candidate approach to multidiscipline dynamic system simulation by numerical coupling of self-contained, single-discipline subsystem simulations. A satellite motion example containing three subsystems (orbit dynamics, attitude dynamics, and aerodynamics) has been defined and constructed using this approach. Conventional solution methods are used in the subsystem simulations. Distributed and centralized implementations of coupling have been considered. Numerical results are evaluated by direct comparison with a standard total-system, simultaneous-solution approach.

Nixon, D. D.↗

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation↗

3D Multistage Simulation of Each Component of the GE90 Turbofan Engine

A 3D multistage simulation of each component of the GE90 Turbofan engine has been made. This includes 49 blade rows. A coupled simulation of all blade rows will be made very soon. The simulation is running using two levels of parallelism. The first level is on a blade row basis with information shared using files. The second level is using a grid domain decomposition with information shared using MPI. Timings will be shown for running on the SP2, an SGI Origin and a distributed system of HP workstations. On the HP workstations, the CHIMP version of MPI is used, with queuing supplied by LSF (Load Sharing Facility). A script-based control system is used to ensure reliability. An MPEG movie illustrating the flow simulation of the engine has been created using PV3, a parallel visualization library created by Bob Haimes of MIT. PVM is used to create a virtual machine from 10 HP workstations and display on an SGI workstation. A representative component simulation will be compared to rig data to demonstrate its usefulness in turbomachinery design and analysis.

Turner, Mark↗

Combustion Instabilities Modeled

NASA Lewis Research Center's Advanced Controls and Dynamics Technology Branch is investigating active control strategies to mitigate or eliminate the combustion instabilities prevalent in lean-burning, low-emission combustors. These instabilities result from coupling between the heat-release mechanisms of the burning process and the acoustic flow field of the combustor. Control design and implementation require a simulation capability that is both fast and accurate. It must capture the essential physics of the system, yet be as simple as possible. A quasi-one-dimensional, computational fluid dynamics (CFD) based simulation has been developed which may meet these requirements. The Euler equations of mass, momentum, and energy have been used, along with a single reactive species transport equation to simulate coupled thermoacoustic oscillations. A very simple numerical integration scheme was chosen to reduce computing time. Robust boundary condition procedures were incorporated to simulate various flow conditions (e.g., valves, open ends, and choked inflow) as well as to accommodate flow reversals that may arise during large flow-field oscillations. The accompanying figure shows a sample simulation result. A combustor with an open inlet, a choked outlet, and a large constriction approximately two thirds of the way down the length is shown. The middle plot shows normalized, time-averaged distributions of the relevant flow quantities, and the bottom plot illustrates the acoustic mode shape of the resulting thermoacoustic oscillation. For this simulation, the limit cycle peak-to-peak pressure fluctuations were 13 percent of the mean. The simulation used 100 numerical cells. The total normalized simulation time was 50 units (approximately 15 oscillations), which took 26 sec on a Sun Ultra2.

Paxson, Daniel E.↗

Probability Simulations by Non-Lipschitz Chaos

It has been demonstrated that classical probabilities, and in particular, probabilistic Turing machine, can be simulated by combining chaos and non-Lipschitz dynamics, without utilization of any man-made devices. Self-organizing properties of systems coupling simulated and calculated probabilities and their link to quantum computations are discussed.

classical dynamics classical probabilities probabi↗

System Simulation by Recursive Feedback: Coupling A Set of Stand-Alone Subsystem Simulations

Recursive feedback is defined and discussed as a framework for development of specific algorithms and procedures that propagate the time-domain solution for a dynamical system simulation consisting of multiple numerically coupled self-contained stand-alone subsystem simulations. A satellite motion example containing three subsystems (other dynamics, attitude dynamics, and aerodynamics) has been defined and constructed using this approach. Conventional solution methods are used in the subsystem simulations. Centralized and distributed versions of coupling structure have been addressed. Numerical results are evaluated by direct comparison with a standard total-system simultaneous-solution approach.

Nixon, Douglas D.↗

System Simulation by Recursive Feedback: Coupling a Set of Stand-Alone Subsystem Simulations

Recursive feedback is defined and discussed as a framework for development of specific algorithms and procedures that propagate the time-domain solution for a dynamical system simulation consisting of multiple numerically coupled, self-contained, stand-alone subsystem simulations. A satellite motion example containing three subsystems (orbit dynamics, attitude dynamics, and aerodynamics) has been defined and constructed using this approach. Conventional solution methods are used in the subsystem simulations. Centralized and distributed versions of coupling structure have been addressed. Numerical results are evaluated by direct comparison with a standard total-system, simultaneous-solution approach.

Nixon, Douglas D.↗

Simulations of Probabilities for Quantum Computing

It has been demonstrated that classical probabilities, and in particular, probabilistic Turing machine, can be simulated by combining chaos and non-LIpschitz dynamics, without utilization of any man-made devices (such as random number generators). Self-organizing properties of systems coupling simulated and calculated probabilities and their link to quantum computations are discussed.

Turing machines↗

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration - dissociation - recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration-dissociation recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗