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 415 records · Page 23

Antiferromagnetic Interlayer Exchange Coupling in All-Semiconducting EuS/PbS/EuS Trilayers

A comprehensive experimental study on the antiferromagnetic interlayer exchange coupling in high quality epitaxial all-semiconducting EuSPbSEuS trilayers is reported. The influence of substrates, the thickness of the non-magnetic PbS spacer layer, and of temperature, was investigated by means of SQUID magnetometry. In trilayers with a PbS thickness between 4 and 12 deg A the low temperature hysteresis loops showed the signature of antiferromagnetic coupling. The value of the interlayer exchange coupling energy was determined by simulating the data with a modified Stoner model, including Zeeman, anisotropy, and exchange coupling energies. An important observation was of a strong dependence of the interlayer exchange coupling energy on temperature, consistent with a power law dependence of the exchange coupling constant on the saturation magnetization of the EuS layers. While no theoretical description is readily available, we conjecture that the observed behavior is due to a dependence of the interlayer exchange coupling energy on the exchange splitting of the EuS conduction band.

Smits, C. J. P.↗

Improved Mars Upper Atmosphere Climatology

The detailed characterization of the Mars upper atmosphere is important for future Mars aerobraking activities. Solar cycle, seasonal, and dust trends (climate) as well as planetary wave activity (weather) are crucial to quantify in order to improve our ability to reasonably depict the state of the Mars upper atmosphere over time. To date, our best information is found in the Mars Global Surveyor (MGS) Accelerometer (ACC) database collected during Phase 1 (Ls = 184 - 300; F10.7 = 70 - 90) and Phase 2 (Ls = 30 - 90; F10.7 = 90 - 150) of aerobraking. This database (100 - 170 km) consists of thermospheric densities, temperatures, and scale heights, providing our best constraints for exercising the coupled Mars General Circulation Model (MGCM) and the Mars Thermospheric General Circulation Model (MTGCM). The Planetary Data System (PDS) contains level 0 and 2 MGS Accelerometer data, corresponding to atmospheric densities along the orbit track. Level 3 products (densities, temperatures, and scale heights at constant altitudes) are also available in the PDS. These datasets provide the primary model constraints for the new MGCM-MTGCM simulations summarized in this report. Our strategy for improving the characterization of the Mars upper atmospheres using these models has been three-fold : (a) to conduct data-model comparisons using the latest MGS data covering limited climatic and weather conditions at Mars, (b) to upgrade the 15-micron cooling and near-IR heating rates in the MGCM and MTGCM codes for ad- dressing climatic variations (solar cycle and seasonal) important in linking the lower and upper atmospheres (including migrating tides), and (c) to exercise the detailed coupled MGCM and MTGCM codes to capture and diagnose the planetary wave (migrating plus non-migrating tidal) features throughout the Mars year. Products from this new suite of MGCM-MTGCM coupled simulations are being used to improve our predictions of the structure of the Mars upper atmosphere for the upcoming MRO aerobraking exercises in 2006. A Michigan website, containing MTGCM output fields from previous climate simulations, is being expanded to include new MGCM-MTGCM simulations addressing planetary wave influences upon thermospheric aerobraking fields (densities and temperatures). In addition, similar MTGCM output fields have been supplied to the MSFC MARSGRAM - 200X empirical model, which will be used in mission operations for conducting aerobraking maneuvers.

Bougher, S. W.↗

Evaluating the Dominant Components of Warming in Pliocene Climate Simulations

The Pliocene Model Intercomparison Project (PlioMIP) is the first coordinated climate model comparison for a warmer palaeoclimate with atmospheric CO2 significantly higher than pre-industrial concentrations. The simulations of the mid-Pliocene warm period show global warming of between 1.8 and 3.6 C above pre-industrial surface air temperatures, with significant polar amplification. Here we perform energy balance calculations on all eight of the coupled ocean-atmosphere simulations within PlioMIP Experiment 2 to evaluate the causes of the increased temperatures and differences between the models. In the tropics simulated warming is dominated by greenhouse gas increases, with the cloud component of planetary albedo enhancing the warming in most of the models, but by widely varying amounts. The responses to mid-Pliocene climate forcing in the Northern Hemisphere midlatitudes are substantially different between the climate models, with the only consistent response being a warming due to increased greenhouse gases. In the high latitudes all the energy balance components become important, but the dominant warming influence comes from the clear sky albedo, only partially offset by the increases in the cooling impact of cloud albedo. This demonstrates the importance of specified ice sheet and high latitude vegetation boundary conditions and simulated sea ice and snow albedo feedbacks. The largest components in the overall uncertainty are associated with clouds in the tropics and polar clear sky albedo, particularly in sea ice regions. These simulations show that albedo feedbacks, particularly those of sea ice and ice sheets, provide the most significant enhancements to high latitude warming in the Pliocene.

climate models↗

Electro-Mechanical Simulation of a Large Aperture MOEMS Fabry-Perot Tunable Filter

We are developing a micro-machined electrostatically actuated Fabry-Perot tunable filter with a large clear aperture for application in high through-put wide-field imaging spectroscopy and lidar systems. In the first phase of this effort, we are developing key components based on coupled electro-mechanical simulations. In particular, the movable etalon plate design leverages high coating stresses to yield a flat surface in drum-head tension over a large diameter (12.5 mm). In this approach, the cylindrical silicon movable plate is back etched, resulting in an optically coated membrane that is suspended from a thick silicon support ring. Understanding the interaction between the support ring, suspended membrane, and coating is critical to developing surfaces that are flat to within stringent etalon requirements. In this work, we present the simulations used to develop the movable plate, spring suspension system, and electrostatic actuation mechanism. We also present results from tests of fabricated proof of concept components.

Kuhn, Jonathan L.↗

Development of the GEOS-5 Atmospheric General Circulation Model: Evolution from MERRA to MERRA2.

The Modern-Era Retrospective Analysis for Research and Applications-2 (MERRA2) version of the GEOS-5 (Goddard Earth Observing System Model - 5) Atmospheric General Circulation Model (AGCM) is currently in use in the NASA Global Modeling and Assimilation Office (GMAO) at a wide range of resolutions for a variety of applications. Details of the changes in parameterizations subsequent to the version in the original MERRA reanalysis are presented here. Results of a series of atmosphere-only sensitivity studies are shown to demonstrate changes in simulated climate associated with specific changes in physical parameterizations, and the impact of the newly implemented resolution-aware behavior on simulations at different resolutions is demonstrated. The GEOS-5 AGCM presented here is the model used as part of the GMAO's MERRA2 reanalysis, the global mesoscale "nature run", the real-time numerical weather prediction system, and for atmosphere-only, coupled ocean-atmosphere and coupled atmosphere-chemistry simulations. The seasonal mean climate of the MERRA2 version of the GEOS-5 AGCM represents a substantial improvement over the simulated climate of the MERRA version at all resolutions and for all applications. Fundamental improvements in simulated climate are associated with the increased re-evaporation of frozen precipitation and cloud condensate, resulting in a wetter atmosphere. Improvements in simulated climate are also shown to be attributable to changes in the background gravity wave drag, and to upgrades in the relationship between the ocean surface stress and the ocean roughness. The series of "resolution aware" parameters related to the moist physics were shown to result in improvements at higher resolutions, and result in AGCM simulations that exhibit seamless behavior across different resolutions and applications.

AGCM↗

Approach to Modeling Boundary Layer Ingestion Using a Fully Coupled Propulsion-RANS Model

Airframe-propulsion integration concepts that use boundary layer ingestion have the potential to reduce aircraft fuel burn. One concept that has been recently explored is NASA's Starc-ABL aircraft configuration, which offers the potential for 12% mission fuel burn reduction by using a turbo-electric propulsion system with an aft-mounted electrically driven boundary layer ingestion propulsor. This large potential for improved performance motivates a more detailed study of the boundary layer ingestion propulsor design, but to date, analyses of boundary layer ingestion have used uncoupled methods. These methods account for only aerodynamic effects on the propulsion system or propulsion system effects on the aerodynamics, but not both simultaneously. This work presents a new approach for building fully coupled propulsive-aerodynamic models of boundary layer ingestion propulsion systems. A 1D thermodynamic cycle analysis is coupled to a RANS simulation to model the Starc-ABL aft propulsor at a cruise condition and the effects variation in propulsor design on performance are examined. The results indicates that both propulsion and aerodynamic effects contribute equally toward the overall performance and that the fully coupled model yields substantially different results compared to uncoupled. The most significant finding is that boundary layer ingestion, while offering substantial fuel burn savings, introduces throttle dependent aerodynamics effects that need to be accounted for. This work represents a first step toward the multidisciplinary design optimization of boundary layer ingestion propulsion systems.

Boundary Layer Ingestion↗

A Framework for Optimization-Based ISRU Tool Design Using Discrete Element Modeling

Novel robotic excavation technologies are needed to perform in-situ resource utilization (ISRU) tasks at levels required to sustain a long-term presence on the lunar surface. Developing and testing multiple iterations of functional hardware is time and cost prohibitive, thus slowing down the pace of progress and delaying humanity’s settlement of the Moon. High-fidelity, physics-based simulation can reduce the time and effort required to develop and deploy robotic systems [1]. We have adopted this approach to create high-fidelity models of robotic test hardware to enable rapid virtual design and optimization of excavation technologies [2]. Such models can leverage modern computational tools like Discrete Element Method (DEM) simulations that can be coupled with automated design approaches like topology optimization to reduce the amount of prototyping and physical testing needed to realize useful tools.

ISRU↗

Convection Effects in Three-dimensional Dendritic Growth

A phase-field model is developed to simulate free dendritic growth coupled with fluid flow for a pure material in three dimensions. The preliminary results presented here illustrate the strong influence of convection on the three-dimensional (3D) dendrite growth morphology. The detailed knowledge of the flow and temperature fields in the melt around the dendrite from the simulations allows for a detailed understanding of the convection effects on dendritic growth.

Lu, Yili↗

Effects of Gyroscopic Cross Coupling Between Pitch and Yaw on the Handling Qualities of VTOL Aircraft

Gyroscopic cross coupling between pitch and yaw was simulated with a variable-stability helicopter while hovering in a yawing maneuver to ascertain the effect of cross coupling on handling qualities. Pilot ratings of the controllability of cross coupling were obtained for various combinations of longitudinal control power, angular velocity about the yaw axis, and simulated engine angular momentum. A theoretical investigation, supplemented by simulator data, was undertaken to determine the effect of longitudinal damping on the coupling controllability. Also, a comparison was made between flight and simulator data. The results indicated that for an aircraft with otherwise satisfactory longitudinal handling qualities, the level of cross coupling is satisfactory when less than 30 percent of the available longitudinal control will trim out the largest gyroscopic coupling moment which might be encountered. Increased longitudinal damping resulted in a significant increase in the controllability of pitch-yaw gyroscopic cross coupling.

Garren, John F., Jr.↗

Impact of Land Model Calibration on Coupled Land-Atmosphere Prediction

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface heat and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry and wet land surface conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through calibration of the Noah land surface model using the new optimization and uncertainty estimation subsystem in NASA's Land Information System (LIS-OPT/UE). The impact of the calibration on the a) spinup of the land surface used as initial conditions, and b) the simulated heat and moisture states and fluxes of the coupled WRF simulations is then assessed. Changes in ambient weather and land-atmosphere coupling are evaluated along with measures of uncertainty propagation into the forecasts. In addition, the sensitivity of this approach to the period of calibration (dry, wet, average) is investigated. Results indicate that the offline calibration leads to systematic improvements in land-PBL fluxes and near-surface temperature and humidity, and in the process provide guidance on the questions of what, how, and when to calibrate land surface models for coupled model prediction.

Santanello, Joseph A., Jr.↗

Application of Rotor Disk Model in the Overflow CFD Code

To reduce the complexity and to reduce the cost of turnaround time of simulating unsteady flows of rotating blades with moving overset grids, amid-fidelity option is provided in the OVERFLOW computational fluid dynamics analysis to model rotors and propellers as infinitely thin disks. The effect of the flow through the disk is represented by the momentum source terms in the governing Navier-Stokes equations. The model is loosely coupled with the rotorcraft comprehensive analysis CAMRADII. A brief description of the coupling algorithm is provided. Simulations of a single rotor and multiple rotors are shown, including comparisons with measurements and with the full rotor where all blades are modeled, and unsteady simulation performed with moving grids for rotating geometry,

ARMD↗

Quantifying the Effect of Mixing on the Mean Age of Air in CCMVal-2 and CCMI-1 Models

The stratospheric age of air (AoA) is a useful measure of the overall capabilities of a general circulation model (GCM) to simulate stratospheric transport. Previous studies have reported a large spread in the simulation of AoA by GCMs and coupled chemistry-climate models (CCMs). Compared to observational estimates, simulated AoA is mostly too low. Here we attempt to untangle the processes that lead to the AoA differences between the models and between models and observations. AoA is influenced by both mean transport by the residual circulation and two-way mixing; we quantify the effects of these processes using data from the CCM inter-comparison projects CCMVal-2 (Chemistry-Climate Model Validation Activity 2) and CCMI-1 (Chemistry-Climate Model Initiative, phase 1). Transport along the residual circulation is measured by the residual circulation transit time (RCTT). We interpret the difference between AoA and RCTT as additional aging by mixing. Aging by mixing thus includes mixing on both the resolved and subgrid scale. We find that the spread in AoA between the models is primarily caused by differences in the effects of mixing and only to some extent by differences in residual circulation strength. These effects are quantified by the mixing efficiency, a measure of the relative increase in AoA by mixing. The mixing efficiency varies strongly between the models from 0.24 to 1.02. We show that the mixing efficiency is not only controlled by horizontal mixing, but by vertical mixing and vertical diffusion as well. Possible causes for the differences in the models' mixing efficiencies are discussed. Differences in subgrid-scale mixing (including differences in advection schemes and model resolutions) likely contribute to the differences in mixing efficiency. However, differences in the relative contribution of resolved versus parameterized wave forcing do not appear to be related to differences in mixing efficiency or AoA.

CCMI↗

Influence of Coupled Radiation and Ablation on the Aerothermodynamic Environment of Planetary Entry Vehicles

A review of recently published coupled radiation and ablation capabilities involving the simulation of hypersonic flowfields relevant to Earth, Mars, or Venus entry is presented. The three fundamental mechanisms of radiation coupling are identified as radiative cooling, precursor photochemistry, and ablation-radiation interaction. The impact of these mechanisms are shown to be significant for a 3 m radius sphere entering Earth at hypothetical Mars return conditions (approximately 15 km/s). To estimate the influence precursor absorption on the radiative flux for a wide range of conditions, a simplified approach is developed that requires only the non-precursor solution. Details of a developed coupled ablation approach, which is capable of treating both massively ablating flowfields in the sublimation regime and weakly ablating diffusion Climited oxidation cases, are presented. A review of the two primary uncoupled ablation approximations, identified as the blowing correction and film coefficient approximations, is made and their impact for Earth and Mars entries is shown to be significant for recession and convective heating predictions. Fully coupled ablation and radiation simulations are presented for the Mars return sphere throughout its entire trajectory. Applying to the Mars return sphere the Pioneer- Venus heritage carbon phenolic heatshield, which has properties available in the open literature, the differences between steady state ablation and coupling to a material response code are shown to be significant.

Johnston, Christopher O.↗

Relative Contribution of Greenhouse Gases and Ozone Change to Temperature Trends in the Stratosphere: A Chemistry/Climate Model Study

Long-term changes in greenhouse gases, primarily carbon dioxide, are expected to lead to a warming of the troposphere and a cooling of the stratosphere. We examine the cooling of the stratosphere and compare the contributions greenhouse gases and ozone change for the decades between 1980 and 2000. We use 150 years of simulation done with our coupled chemistry/climate model (GEOS 4 GCM with GSFC CTM chemistry) to calculate temperatures and constituents fiom,1950 through 2100. The contributions of greenhouse gases and ozone to temperature change are separated by a time-series analysis using a linear trend term throughout the period to represent the effects of greenhouse gases and an equivalent effective stratospheric chlorine (EESC) term to represent the effects of ozone change. The temperature changes over the 150 years of the simulation are dominated by the changes in greenhouse gases. Over the relatively short period (approx. 20 years) of ozone decline between 1980 and 2000 changes in ozone are competitive with changes in greenhouse gases. The changes in temperature induced by the ozone change are comparable to, but smaller than, those of greenhouse gases in the upper stratosphere (1-3 hPa) at mid latitudes. The ozone term dominates the temperature change near both poles with a negative temperature change below about 3-5 hPa and a positive change above. At mid latitudes in the upper stratosphere and mesosphere (above about 1 hPa) and in the middle stratosphere (3 to 70 ma), the greenhouse has term dominates. From about 70 hPa down to the tropopause at mid latitudes, cooling due to ozone changes is the largest influence on temperature. Over the 150 years of the simulation, the change in greenhouse gases is the most important contributor to temperature change. Ozone caused a perturbation that is expected to reverse over the coming decades. We show a model simulation of the expected temperature change over the next two decades (2006-2026). The simulation shows a crossover between lower atmospheric heating and upper atmospheric cooling that is located at about 90 hPa in the tropics and 30-40 hPa in the polar regions. This results from the combination of continuing increases in greehouse gases and recovery from ozone depletion.

Stolarski, Richard S.↗

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. However, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of ~8–12°C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model–data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

Pliocene Model Intercomparison Project Phase 2 (Pl↗

System Verification of MSL Skycrane Using an Integrated ADAMS Simulation

Mars Science Laboratory (MSL) will use the Skycrane architecture to execute final descent and landing maneuvers. The Skycrane phase uses closed-loop feedback control throughout the entire phase, starting with rover separation, through mobility deploy, and through touchdown, ending only when the bridles have completely slacked. The integrated ADAMS simulation described in this paper couples complex dynamical models created by the mechanical subsystem with actual GNC flight software algorithms that have been compiled and linked into ADAMS. These integrated simulations provide the project with the best means to verify key Skycrane requirements which have a tightly coupled GNC-Mechanical aspect to them. It also provides the best opportunity to validate the design of the algorithm that determines when to cut the bridles. The results of the simulations show the excellent performance of the Skycrane system.

White, Christopher↗

Improving Fidelity of Launch Vehicle Liftoff Acoustic Simulations

Launch vehicles experience high acoustic loads during ignition and liftoff affected by the interaction of rocket plume generated acoustic waves with launch pad structures. Application of highly parallelized Computational Fluid Dynamics (CFD) analysis tools optimized for application on the NAS computer systems such as the Loci/CHEM program now enable simulation of time-accurate, turbulent, multi-species plume formation and interaction with launch pad geometry and capture the generation of acoustic noise at the source regions in the plume shear layers and impingement regions. These CFD solvers are robust in capturing the acoustic fluctuations, but they are too dissipative to accurately resolve the propagation of the acoustic waves throughout the launch environment domain along the vehicle. A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed to improve such liftoff acoustic environment predictions. The framework combines the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin (DG) solver, Loci/THRUST, developed in the same computational framework. Loci/THRUST employs a low dissipation, high-order, unstructured DG method to accurately propagate acoustic waves away from the source regions across large distances. The DG solver is currently capable of solving up to 4th order solutions for non-linear, conservative acoustic field propagation. Higher order boundary conditions are implemented to accurately model the reflection and refraction of acoustic waves on launch pad components. The DG solver accepts generalized unstructured meshes, enabling efficient application of common mesh generation tools for CHEM and THRUST simulations. The DG solution is coupled with the CFD solution at interface boundaries placed near the CFD acoustic source regions. Both simulations are executed simultaneously with coordinated boundary condition data exchange.

Liever, Peter↗

Machine-Learned Committor Functions for Reactive Molecular Dynamics

Reactive molecular dynamics (MD) is a powerful tool for atomistic-scale modeling of a diverse range of chemical processes. However, scaling these simulations to large systems and long times scales remains a challenge because of the complexity of the potential energy function required. The authors previously developed a heuristic approach, called REACTER, that incorporates reactivity in MD simulations in a less general but much more computationally efficient manner. REACTER uses standard, fixed valence force fields as the underlying potentialenergy surface for describing all interatomic interactions but adds a procedure for enforcing user-defined reactions that occur when certain geometric constraints on relative atomic positions are satisfied. Further, these bonding changes can be accepted or rejected with a probability related tothe local thermal energy. This work seeks to generalize this approach by replacing the set of user defined geometric constraints and energetic criteria with a committor function that specifies the probability of a reaction occurring on the basis of the local atomic configuration. The committor function is a useful mathematical tool for modeling rare events but, unfortunately, is very difficult to compute for realistic systems in a general way. This work describes a method for approximating the committor function using a machine learning approach, specifically a deep neural network trained with data from reactive MD and DFT-based dynamics simulations. This network is coupled to the existing REACTER protocol, as implemented in the LAMMPS MD package, and used to make on-the-fly predictions of reaction probabilities without the more extensive user input previously required. The new method is demonstrated using the polymerization of polystyrene as a case study. Although very dependent on the quality and quantity of training data, machine-learned committor functions show promise as a method for incorporating reaction probability from higher level calculations into highly scalable MD simulations.

polymer simulations↗