Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eulerian model”

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

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, where droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Computational Fluid Dynamics↗

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, wherein droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Self-Pressurization↗

Modeling Primary Breakup: A Three-Dimensional Eulerian Level Set/Vortex Sheet Method for Two-Phase Interface Dynamics

This paper is divided into four parts. First, the level set/vortex sheet method for three-dimensional two-phase interface dynamics is presented. Second, the LSS model for the primary breakup of turbulent liquid jets and sheets is outlined and all terms requiring subgrid modeling are identified. Then, preliminary three-dimensional results of the level set/vortex sheet method are presented and discussed. Finally, conclusions are drawn and an outlook to future work is given.

Herrmann, M.↗

Calculation of confined gas-particle two-phase turbulent flows

Numerical calculations have been carried out for gas-particle confined turbulent flows using a recently developed two-phase turbulence closure model. The present modeling scheme utilizes Eulerian formulations of the transport equations and accounts for the combined effects of interphase slip and turbulent dispersion of particles. A multiple-scale turbulence model is used for the turbulent field modeling of the underlying fluid flow. For the particle size and particle loading considered in this study, the fluid turbulence transport equations must be modified to include the damping effects of particles. Predictions and comparisons are made in the fully developed gas-solid pipe flow and the confined particle-laden jet. Numerical results are in reasonably good agreement with the published experimental information.

Chen, C. P.↗

A model of stratospheric chemistry and transport on an isentropic surface

This paper presents a new photochemical transport model designed to simulate the behavior of stratospheric trace species in the middle stratosphere. The model has an Eulerian grid with the latitude and longitude coordinates on a single isentropic surface (hemispheric or global), in which both the dynamical and the photochemical processes can be accurately represented. The model is intgegrated for 12 days with winds and temperatures supplied by three-dimensional integration of an idealized wavenumber-one disturbance. The results for the long-lived tracers such as N2O showed excellent correlation with the potential vorticity distribution, validating the transport scheme. Calculations with zonally averaged wind and temperature fields showed that discrepancies in the calculation of the zonal mean were less than 10 percent for O3 and HNO3, compared with the zonal mean of the previous results.

Austin, John↗

Multiphase Modeling of Water Injection on Flame Deflector

This paper describes the use of an Eulerian Dispersed Phase (EDP) model to simulate the water injected from the flame deflector and its interaction with supersonic rocket exhaust from a proposed Space Launch System (SLS) vehicle. The Eulerian formulation, as part of the multi-phase framework, is described. The simulations show that water cooling is only effective over the region under the liquid engines. Likewise, the water injection provides only minor effects over the surface area under the solid engines.

Vu, Bruce T.↗

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗

Sympathetic reaction tests and analyses

The critical separation distances for explosive reactions of a solid rocket propellant were measured. Explosive reactions included low order explosion, low order detonation, and high order detonation. The effects of sample size, shape, damage and temperature on sympathetic reaction were determined experimentally. The sympathetic detonation of small cubes of solid rocket propellant was modelled numerically, using the Eulerian reactive hydrodynamic code 2DE with Forest Fire burn rates. The model was applied to cubes of 2.54 - 7.62 cm (1 - 3 in.), with agreement between calculated and experimental results.

Ricardson, D. E.↗

Modeling of combustion processes of stick propellants via combined Eulerian-Lagrangian approach

This research is motivated by the improved ballistic performance of large-caliber guns using stick propellant charges. A comprehensive theoretical model for predicting the flame spreading, combustion, and grain deformation phenomena of long, unslotted stick propellants is presented. The formulation is based upon a combined Eulerian-Lagrangian approach to simulate special characteristics of the two phase combustion process in a cartridge loaded with a bundle of sticks. The model considers five separate regions consisting of the internal perforation, the solid phase, the external interstitial gas phase, and two lumped parameter regions at either end of the stick bundle. For the external gas phase region, a set of transient one-dimensional fluid-dynamic equations using the Eulerian approach is obtained; governing equations for the stick propellants are formulated using the Lagrangian approach. The motion of a representative stick is derived by considering the forces acting on the entire propellant stick. The instantaneous temperature and stress fields in the stick propellant are modeled by considering the transient axisymmetric heat conduction equation and dynamic structural analysis.

Kuo, K. K.↗

Impact ejecta on the moon

The partitioning of energy and the distribution of the resultant ejecta on the moon is numerically modeled using a Eulerian finite difference grid. The impact of an iron meteoroid at 15 km/sec on a gabbroic anorthosite lunar crust is examined. The high speed impact induced flow is described over the entire hydrodynamic regime from a time where the peak pressures are 6 Mbar until the stresses everywhere in the flow are linearly elastic, and less than 5 kbar. Shock-induced polymorphic phase changes, (plagioclase and pyroxene to hollandite and perovskite), and the subsequent reversion to low pressure phases are demonstrated to enhance shock wave attenuation. A rate-dependent equation of state is used for describing the hysteretic effect of the phase change. Ballistic equations for a spherical planet are then applied to material with net velocity away from the moon.

Okeefe, J. D.↗

A model-based diagnostic study of the rapid development phase of the Presidents' Day cyclone

A model simulation of the rapid development phaser of the Presidents' Day cyclone on February 19, 1979, is studied. Model trajectories and Eulerian analyses show that, during this phase, three airstreams converge into the cyclogenetic region. The evolution of the PJ-trough system and its associated tropopause fold is described, including the simulated surface low and the vertical profile of mass divergence which contributes to the decreasing sea-level pressure. The various terms of the vorticity equation are evaluated to identify the processes contributing to the increase in absolute vorticity. The nature and origin of airstreams entering the low-level circulation associated with rapid cyclogenesis are examined from a Lagrangian perspective.

Whitaker, Jeffrey S.↗

A diagnostic model for studying daytime urban air quality trends

A single cell Eulerian photochemical air quality simulation model was developed and validated for selected days of the 1976 St. Louis Regional Air Pollution Study (RAPS) data sets; parameterizations of variables in the model and validation studies using the model are discussed. Good agreement was obtained between measured and modeled concentrations of NO, CO, and NO2 for all days simulated. The maximum concentration of O3 was also predicted well. Predicted species concentrations were relatively insensitive to small variations in CO and NOx emissions and to the concentrations of species which are entrained as the mixed layer rises.

Brewer, D. A.↗

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗

Very Large Eddy Simulations of a Jet-A Spray Reacting Flow in a Single Element LDI Injector With and Without Invoking an Eulerian Scalar DWFDF Method

This paper presents the very large eddy simulations (VLES) of a Jet-A spray reacting flow in a single element lean direct injection (LDI) injector by using the National Combustion Code (NCC) with and without invoking the Eulerian scalar DWFDF method, in which DWFDF is defined as the density weighted time filtered fine grained probability density function. The flow field is calculated by using the time filtered compressible Navier-Stokes equations (TFNS) with nonlinear subscale turbulence models, and when the Eulerian scalar DWFDF method is invoked, the energy and species mass fractions are calculated by solving the equation of DWFDF. A nonlinear subscale model for closing the convection term of the Eulerian scalar DWFDF equation is used and will be briefly described in this paper. Detailed comparisons between the results and available experimental data are carried out. Some positive findings of invoking the Eulerian scalar DWFDF method in both improving the simulation quality and maintaining economic computing cost are observed.

Shih, Tsan-Hsing↗

Evaporating Spray in Supersonic Streams Including Turbulence Effects

Evaporating spray plays an important role in spray combustion processes. This paper describes the development of a new finite-conductivity evaporation model, based on the two-temperature film theory, for two-phase numerical simulation using Eulerian-Lagrangian method. The model is a natural extension of the T-blob/T-TAB atomization/spray model which supplies the turbulence characteristics for estimating effective thermal diffusivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating sprays. High speed cross flow spray results indicate the effectiveness of the T-blob/T-TAB model and point to the needed improvements in high speed evaporating spray modeling.

Balasubramanyam, M. S.↗

A two-phase restricted equilibrium model for combustion of metalized solid propellants

An Eulerian-Lagrangian two-phase approach was adopted to model the multi-phase reacting internal flow in a solid rocket with a metalized propellant. An Eulerian description was used to analyze the motion of the continuous phase which includes the gas as well as the small (micron-sized) particulates, while a Lagrangian description is used for the analysis of the discrete phase which consists of the larger particulates in the motor chamber. The particulates consist of Al and Al2O3 such that the particulate composition is 100 percent Al at injection from the propellant surface with Al2O3 fraction increasing due to combustion along the particle trajectory. An empirical model is used to compute the combustion rate for agglomerates while the continuous phase chemistry is treated using chemical equilibrium. The computer code was used to simulate the reacting flow in a solid rocket motor with an AP/HTPB/Al propellant. The computed results show the existence of an extended combustion zone in the chamber rather than a thin reaction region. The presence of the extended combustion zone results in the chamber flow field and chemical being far from isothermal (as would be predicted by a surface combustion assumption). The temperature in the chamber increases from about 2600 K at the propellant surface to about 3350 K in the core. Similarly the chemical composition and the density of the propellant gas also show spatially non-uniform distribution in the chamber. The analysis developed under the present effort provides a more sophisticated tool for solid rocket internal flow predictions than is presently available, and can be useful in studying apparent anomalies and improving the simple correlations currently in use. The code can be used in the analysis of combustion efficiency, thermal load in the internal insulation, plume radiation, etc.

Sabnis, J. S.↗

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

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 Predictive Simulation Capability (PSC) Element is 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 predictive simulation capability, 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 will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗