Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled process modeling”

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 127 records · Page 7

Hybrid propulsion for launch vehicle boosters: A program status update

Results obtained in studying the origin and suppression of large-amplitude pressure oscillations in a 24 in. diameter hybrid motor using a liquid oxygen/hydroxylterminated polybutadiene/polycyclopentadiene propellant system are discussed. Tests conducted with liquid oxygen flow rates varying from 10 to 40 lbm/sec were designed to gauge the effectiveness of various vaporization chamber flow fields, injector designs, and levels of heat addition in suppressing high-frequency longitudinal mode oscillations. Longitudinal acoustic modes did not arise in any tests. However, initial testing revealed the presence of high-amplitude, sinusoidal, nonacoustic oscillations persisting throughout the burn durations. Analysis showed this to be analogous to chug mode instability in liquid rocket engines brought about by a coupling of motor combustion processes and the liquid oxygen feed system. Analytical models were developed and verified by test data to predict the amplitude and frequency of feed-system-coupled combustion pressure oscillations. Subsequent testing showed that increasing the feed system impedance eliminated the bulk mode instability. This paper documents the work completed to date in performance of the Hybrid Propulsion Technology for Launch Vehicle Boosters Program (NAS8-39942) sponsored by NASA's George C. Marshall Space Flight Center.

Carpenter, R. L.↗

System Level Uncertainty Assessment for Collaborative RLV Design

A collaborative design process utilizing Probabilistic Data Assessment (PDA) is showcased. Given the limitation of financial resources by both the government and industry, strategic decision makers need more than just traditional point designs, they need to be aware of the likelihood of these future designs to meet their objectives. This uncertainty, an ever-present character in the design process, can be embraced through a probabilistic design environment. A conceptual design process is presented that encapsulates the major engineering disciplines for a Third Generation Reusable Launch Vehicle (RLV). Toolsets consist of aerospace industry standard tools in disciplines such as trajectory, propulsion, mass properties, cost, operations, safety, and economics. Variations of the design process are presented that use different fidelities of tools. The disciplinary engineering models are used in a collaborative engineering framework utilizing Phoenix Integration's ModelCenter and AnalysisServer environment. These tools allow the designer to join disparate models and simulations together in a unified environment wherein each discipline can interact with any other discipline. The design process also uses probabilistic methods to generate the system level output metrics of interest for a RLV conceptual design. The specific system being examined is the Advanced Concept Rocket Engine 92 (ACRE-92) RLV. Previous experience and knowledge (in terms of input uncertainty distributions from experts and modeling and simulation codes) can be coupled with Monte Carlo processes to best predict the chances of program success.

Charania, A. C.↗

Simulation of 1986 South China Sea Monsoon with a Regional Climate Model

A Regional Land-Atmosphere Climate Simulation System (RELACS) project is being developed at NASA Goddard Space Flight Center. One of the major goals of RELACS is to use a regional scale model with improved physical processes and in particular land-related processes, to understand the role of the land surface and its interaction with convection and radiation as well as the water/energy cycles in the IndoChina/South China Sea (SCS) region. The Penn State/NCAR MM5 atmospheric modeling system, a state of the art atmospheric numerical model designed to simulate regional weather and climate, has been successfully coupled to the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model. The original MM5 model (without PLACE) includes the option for either a simple slab soil model or a five-layer soil model (MRF) in which the soil moisture availability evolves over time. However, the MM5 soil models do not include the effects of vegetation, and thus important physical processes such as evapotranspiration and interception are precluded. The PLACE model incorporates vegetation type and has been shown in international comparisons to accurately predict evapotranspiration and runoff over a wide variety of land surfaces. The coupling of MM5 and PLACE creates a numerical modeling system with the potential to more realistically simulate atmosphere and land surface processes including land-sea interaction, regional circulations such as monsoons, and flash flood events. In addition, the Penn State/NCAR MM5 atmospheric modeling system has been: (1) coupled to the Goddard Ice Microphysical scheme; (2) coupled to a turbulent kinetic energy (TKE) scheme; (3) modified to ensure cloud budget balance; and (4) incorporated initialization with the Goddard EOS data sets at NASA/Goddard Laboratory for Atmospheres. The improved MM5 with two nested domains (60 and 20 km horizontal resolution) was used to simulate convective activity over IndoChina and the South China Sea, during the monsoon season, from May 6 to May 20, 1986. The model results captured several dominant observed features, such as twin cyclones, a depression system over the Bay of Bengal, strong south-westerly winds over IndoChina before and during the on-set of convection over the SCS, and a vortex over the SCS. Two additional MM5 runs with different land process models, Blackadar and MRF, were performed, and their results are compared to the run with PLACE. The preliminary results indicate that the MM5 results using PLACE and Blackadar are in very good agreement, but the results using MRF do not contain the south-westerly wind over IndoChina prior to the on-set of convection over the SCS.

Tao, W. -K.↗

Performance of Coupled Physics Solvers for Multidisciplinary Hypersonic Flow Simulations on Several Classes of Computer Architectures

The application of hypersonic flow simulation tools to realistic flight scenarios will require the coupling of multiple physical effects to the baseline fluid dynamics. Such multiphysics effects can include the aerooelastic response of the airframe or engine components, dynamic transport of atmospheric particles, the deformation of solid-fluid interfaces that can ablate, pyrolyze, or erode, as well as a host of other processes, all of which are governed by unique sets of physical equations and models. Coupling multiple (and potentially disparate) physics solvers to a robust compressible flow solver poses additional challenges related to the stability, performance and scalability of the combined solver. The choices made during the software design process can therefore lead to a variation in simulation efficiency across different computer architectures. In this paper, we will consider two representative multiphysics hypersonic flow scenarios: the interaction of solid particulates with the flow field created by a hypersonic lifting body and the aerooelastic deformation of a model airframe under high-Mach-number flow conditions. For these simulations we explore the behavior of several hypersonic simulation tools, including Kestrel, FUN3D, US3D, and JENRE multiphysics framework, on several high performance computing systems containing various CPU and GPU architectures.

architecture↗

Low Stretch Solid-Fuel Flame Transient Response to a Step Change in Gravity

The effect of a step change in gravity level on the stability of low stretch diffusion flames over a solid fuel is studied both numerically and experimentally. Drop tower experiments have been conducted in NASA Glenn Research Center's 5.2 Zero Gravity Facility. In the experiments burning PMMA cylinders, a dynamic transition is observed when the steadily burning 1g flame is dropped and becomes a 0g flame. To understand the physics behind this dynamic transition, a transient stagnation point model has been developed which includes gas-phase radiation and solid phase coupling to describe this dynamic process. In this paper, the experimental results are compared with the model predictions. Both model and experiment show that the interior of the solid phase does not have time to change significantly in the few seconds of drop time, so the experimental results are pseudo-steady in the gas-phase, but the solid is inherently unsteady over long time scales. The model is also used to examine the importance of fractional heat losses on extinction, which clearly demonstrates that as the feedback from the flame decreases, the importance of the ongoing heat losses becomes greater, and extinction is observed when these losses represent 80% or more of the flame feedback.

Armstrong, J. B.↗

The Continual Intercomparison of Radiation Codes (CIRC) Assessing Anew the Quality of GCM Radiation Algorithms

The simulation of changes in the Earth's climate due to solar and thermal radiative processes with global climate models (GCMs) is highly complex, depending on the parameterization of a multitude of nonlinearly coupled physical processes. In contrast, the germ of global climate change, the radiative forcing from enhanced abundances of greenhouse gases, is relatively well understood. The impressive agreement between detailed radiation calculations and highly resolved spectral radiation measurements in the thermal infrared under cloudless conditions (see, for example, Fig. 1) instills confidence in our knowledge of the sources of gaseous absorption. That the agreement spans a broad range of temperature and humidity regimes using instruments mounted on surface, aircraft, and satellite platforms not only attests to our capability to accurately calculate radiative fluxes under present conditions, but also provides confidence in the spectroscopic basis for computation of fluxes under conditions that might characterize future global climate (e.g., radiative forcing). Alas, the computational costs of highly resolved spectral radiation calculations cannot be afforded presently in GCMs. Such calculations have instead been used as the foundation for approximations implemented in fast but generally less accurate algorithms performing the needed radiative transfer (RT) calculations in GCMs. Credible climate simulations by GCMs cannot be ensured without accurate solar and thermal radiative flux calculations under all types of sky conditions: pristine cloudless, aerosol-laden, and cloudy. The need for accuracy in RT calculations is not only important for greenhouse gas forcing scenarios, but is also profoundly needed for the robust simulation of many other atmospheric phenomena, such as convective processes.

Oreopoulos, Lazaros↗

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

combustion↗

The Sensitivity of Model Ozone to Advective and Photochemical Processes in the High Latitude Winter Lower Stratosphere

Three dimensional chemistry and transport models (CTMs) contain a set of coupled continuity equations which describe the evolution of constituents such as ozone and other minor species which affect ozone. Both advection and photochemical processes contribute to constituent evolution, and a CTM provides a means to evaluate these contributions separately. Such evaluation is particularly useful when both terms are important to the modeled tendency. An example is the ozone tendency in the high latitude winter lower stratosphere, where advection tends to increase ozone, and catalytic processes involving chlorine radicals tend to decrease ozone. The Goddard three dimensional chemistry and transport model uses meteorological fields from the Goddard Earth Observing System Data Assimilation System, thus the modeled ozone evolution may reproduce the observed evolution and provide a test of the model representation of photochemical processes if the transport is shown to be modeled appropriately. We have investigated the model advection further using diabatic trajectory calculations. For long lived constituents such as N2O, the model field for a particular time on a potential temperature surface is compared with a field produced by calculating 15 day back trajectories for a fixed latitude longitude grid, and mapping model N2O at the terminus of the back trajectories onto the initial grid. This provides a quantitative means to evaluate two aspects of the CTM transport: one, the model horizontal gradient between middle latitudes and the polar vortex is compared with the gradient produced using the non-diffusive trajectory calculation; two, the model vertical advection, which is produced by the divergence of the horizontal winds, is compared with the vertical transport expected from diabatic cooling.

Douglass, A.↗

Numerical study of chemically reacting flows using an LU scheme

A new computational fluid dynamic code has been developed for the study of mixing and chemical reactions in the flow fields of ramjets and scramjets. The code employs an implicit finite volume, lower-upper symmetric successive overrelaxation scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a fully-coupled and very efficient manner. The combustion processes are modeled by an 8-species, 14-step finite rate chemistry model whereas turbulence is simulated by a Baldwin-Lomax algebraic model. The validity of the code is demonstrated by comparing the numerical calculations with both experimental data and previous calculations of a cold flow helium injection into a straight channel and premixed hydrogen-air reacting flows in a ramped duct. The code is then used to calculate the mixing and chemical reactions of a hydrogen jet transversely injected into a supersonic airstream. Results are presented describing the flow field, the recirculation regions in front and behind the injector, and the chemical reactions.

Shuen, Jian Shun↗

Numerical study of chemically reacting flows using an LU scheme

A new computational fluid dynamic code has been developed for the study of mixing and chemical reactions in the flow fields of ramjets and scramjets. The code employs an implicit finite volume, lower-upper symmetric successive overrelaxation scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a fully-coupled and very efficient manner. The combustion processes are modeled by an 8-species, 14-step finite rate chemistry model whereas turbulence is simulated by a Baldwin-Lomax algebraic model. The validity of the code is demonstrated by comparing the numerical calculations with both experimental data and previous calculations of a cold flow helium injection into a straight channel and premixed hydrogen-air reacting flows in a ramped duct. The code is then used to calculate the mixing and chemical reactions of a hydrogen jet transversely injected into a supersonic airstream. Results are presented describing the flow field, the recirculation regions in front and behind the injector, and the chemical reactions.

Shuen, Jian Shun↗

A numerical study of chemically reacting flow in nozzles

The space station uses small rocket motors, called thrusters, for orientation control. Because of the lack of viable design tools for small rockets, the initial thruster design was basically a very small version of a large rocket motor. Thrust measurements of the initial design were lower than predicted. To improve predictions it was decided to develop a verison of the RPLUS2D reacting flow code for thruster calculations. RPLUS2D employs an implicit finite volume, lower-upper symmetric successive overrelaxation (LU-SSOR) scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a coupled and very efficient manner. The combustion processes are modeled by a 9-species, 18 step finite-rate chemistry model, and the turbulence is simulated by a Baldwin-Lomax algebraic model. The code is extended to handle multiple subsonic inlet conditions where the total mass flow is governed by conditions calculated at the thruster-throat. Results are shown for a thruster design where the overall mixture ratio is hydrogen rich. A calculation of a large area ratio divergent nozzle is also presented.

Vanoverbeke, Thomas J.↗

A numerical study of chemically reacting flow in nozzles

The Space Station uses small rocket motors, called thrusters, for orientation control. Because of the lack of viable design tools for small rockets, the initial thruster design was basically a very small version of a large rocket motor. Thrust measurements of the initial design were lower than predicted. To improve predictions it was decided to develop a version of the RPLUS2D reacting flow code for thruster calculations. RPLUS2D employs an implicit finite volume, lower-upper symmetric successive overrelaxation (LU-SSOR) scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a coupled and very efficient manner. The combustion processes are modeled by a 9-species, 18 step finite-rate chemistry model, and the turbulence is simulated by a Baldwin-Lomax algebraic model. The code is extended to handle multiple subsonic inlet conditions where the total mass flow is governed by conditions calculated at the thruster-throat. Results are shown for a thruster design where the overall mixture ratio is hydrogen rich. A calculation of a large area ratio divergent nozzle is also presented.

Vanoverbeke, Thomas J.↗

I-deas TMG to NX Space Systems Thermal Model Conversion and Computational Performance Comparison

CAD/CAE packages change on a continuous basis as the power of the tools increase to meet demands. End -users must adapt to new products as they come to market and replace legacy packages. CAE modeling has continued to evolve and is constantly becoming more detailed and complex. Though this comes at the cost of increased computing requirements Parallel processing coupled with appropriate hardware can minimize computation time. Users of Maya Thermal Model Generator (TMG) are faced with transitioning from NX I -deas to NX Space Systems Thermal (SST). It is important to understand what differences there are when changing software packages We are looking for consistency in results.

Thermal Model Generator (TMG)↗

An Extensible Perturbed Parameter Ensemble for the Community Atmosphere Model Version 6

This paper documents the methodology and preliminary results from a Perturbed Parameter Ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2-CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations are analyzed and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks. The use of various emulators is explored in the multi- dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that using PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2-CAM6 PPE is publicly available, and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g. coupling to the land surface).

Machine learning↗

Conceptual Feasibility Study of the Hyperloop Vehicle for Next-Generation Transport

The Hyperloop concept is proposed as a faster, cheaper alternative to high-speed rail and traditional short-haul aircraft. It consists of a passenger pod traveling through a tube under light vacuum while being propelled and levitated by a combination of permanent and electro-magnets. The concept addresses NASA's research thrusts for growth in demand, sustainability, and technology convergence for high-speed transport. Hyperloop is a radical departure from other advanced aviation concepts, however it remains an aeronautics concept that tackles the same strategic goals of low-carbon propulsion and ultra-effcient vehicles. System feasibility was investigated by building a multidisciplinary vehicle sizing model that takes into account aerodynamic, thermodynamic, structures, electromagnetic, weight, and mission analyses. The sizing process emphasized the strong coupling between the two largest systems: the tube and the passenger pod. The model was then exercised to examine Hyperloop from a technical and cost perspective. The structural sizing analysis of the travel tube demonstrates potential for signi cant capital cost reductions by considering an underwater route. Examination of varying passenger capacity indicates that the system can be operated with a wide range of passenger loads without significant change in operating expenses. Lastly, a high-level sizing study simulated variations in tube area, pressure, pod speed, and passenger capacity showing that there is a tube pressure that minimizes operating energy usage. The value of this optimal tube pressure is highly sensitive to numerous design details. These combined estimates of energy consumption, passenger throughput, and mission analyses all support Hyperloop as a faster and cheaper alternative to short-haul flights. The tools and expertise used to quantify these results also demonstrate how traditional aerospace design methods can be leveraged to handle the complex and coupled design process. Much of the technology development required for the Hyperloop is shared with next-generation aircraft. Furthermore, the substantial public interest and active commercial development make it an ideal candidate as an aircraft technology driver and test bed.

Hyperloop Vehicle↗

An Ensemble Investigation of the Causes for Regional Air-Quality Model Critical Load Exceedances Prediction Variability in European and North American Domains Using Diagnostics From Phase 4 of the Air Quality Model Evaluation International Initiative

We summarize tentative findings from multi air quality model ensembles for the years 2009 and 2010 in Europe (EU), and 2010 and 2016 in North America (NA), under AQMEII-4. The model predictions of sulphur and nitrogen deposition were used to estimate exceedances of critical loads for acidification and eutrophication, to show the extent to which the ensemble members agree in the magnitude and the trend of ecologically meaningful impacts. Model exceedance variability was analyzed using AQMEII-4 diagnostics. Evaluation against concentration and wet deposition observations, coupled with these diagnostics, identified specific process representations as the causes for variability between model predictions and for reduced model performance. All models predicted reductions in ecosystem acidification impacts in North America between the years 2010 and 2016, in accord with SO2 emissions reduction legislation which started in 2010 (SO2 SIP) However, all models in EU and NA domains had net negative biases for wet deposition of sulphur and nitrogen relative to observations. The wet S deposition average mean bias for the NA ensemble was -0.17 eq ha-1 d-1, and for the EU ensemble -1.15 eq ha-1 d-1. The NA daily wet deposition average mean bias for NH4+ was -0.37 eq ha-1d-1; EU -1.19 eq ha-1 d-1. The daily NA wet NO3- deposition average mean bias was -0.24 eq ha-1d-1; EU -0.69 eq ha-1 d-1. The members of the ensemble diverged (factor of 10) in their North American predictions for Ndep and consequently their eutrophication exceedances. The models with the highest eutrophication predictions also predicted the highest levels of gas-phase ammonia dry deposition (standard deviation of ammonia dry deposition flux across ensemble members was larger than the ensemble average). These models also had negative biases of predicted ammonia concentrations; average mean biases of -0.63 (satellite NH3) and -0.85 ppbv (surface NH3) compared to ensemble averages of -0.30 and -0.34 ppbv. Diagnostics showed that these differences resulted from the manner in which bidirectional ammonia fluxes were parameterized within these models. The second largest source of NA eutrophication prediction variability were models with positive biases in particulate ammonium and nitrate concentrations, and higher particle nitrogen deposition levels ( particle ammonium concentration bias +0.35 ug m-3; ensemble bias +0.15 ug m-3). We believe two factors may have led to these latter overestimates: higher levels of fine mode particle nitrate formation compared to other models (due to the use of an inorganic heterogeneous chemistry algorithm which did not take base cation chemistry into account), and updates to particle dry deposition velocities carried out in the absence of concurrent updates to wet scavenging algorithms. The relative importance of dry gas, dry particulate, and wet deposition towards total sulphur and nitrogen deposition totals differed between EU and North American domains, though all models had negative biases in wet deposition as noted above. Parallel and subsequent work suggests that multiphase hydrometeor scavenging may improve model wet deposition performance. An increased research focus is recommended for four model processes: multiphase hydrometeor scavenging, ammonia bidirectional fluxes, base cation chemistry and emissions, and particle dry deposition.

regional air-quality model↗

Modal Testing of a Flexible Wing on a Dynamically Active Test Fixture Using Fixed Base Correction Method - IFASD 2019

In modal testing and finite element model correlation, analysts desire modal results using free-free or rigid boundary conditions to ease comparisons of test versus analytical data. It is often expensive both in cost and schedule to build and test with boundary conditions that replicate the free-free or rigid boundaries. Static test fixtures for load testing are often large, heavy, and unyielding, but do not provide adequate boundaries for modal tests because they are dynamically too flexible and often contain natural frequencies within the test article frequency range of interest. The dynamic coupling between the test article and test fixture complicates the model updating process because significant effort needs to be spent on modeling the test fixture and boundary conditions in addition to the test article. If the modal results could be corrected for fixture coupling, then setups used for other structural testing could be adequate for modal testing and would allow significant schedule and cost savings by eliminating a unique setup for only modal testing. To simplify future modal tests, this paper reports on a fixed base correction method that was investigated during modal testing of a full-scale, half-span, flexible wing cantilevered from a static test fixture. The results of this fixed base correction approach look very promising. The method aided in producing similar wing modal characteristics for two different physical boundary configurations of a dynamically active test fixture.

modal testing↗

Modal Testing of a Flexible Wing on a Dynamically Active Test Fixture Using the Fixed Base Correction Method

In modal testing and finite element model correlation, analysts desire modal results using free-free or rigid boundary conditions to ease comparisons of test versus analytical data. It is often expensive both in cost and schedule to build and test with boundary conditions that replicate the free-free or rigid boundaries. Static test fixtures for load testing are often large, heavy, and unyielding, and not provide adequate boundaries for modal tests because they are dynamically too flexible and often contain natural frequencies within the test article frequency range of interest. The dynamic coupling between the test article and test fixture complicates the model updating process because significant effort needs to be spent on modeling the test fixture and boundary conditions in addition to the test article. If the modal results could be corrected for fixture coupling, then setups used for other structural testing could be adequate for modal testing and would allow significant schedule and cost savings by eliminating a unique setup for only modal testing. To simplify future modal tests, this report describes a Fixed Base Correction method that was investigated during modal testing of a full-scale, half-span, flexible wing cantilevered from a static test fixture. The results of this Fixed Base Correction approach look very promising. The method aided in producing similar wing modal characteristics for two different physical boundary configurations of a dynamically active test fixture.

Natalie D Spivey↗