Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian Particle Tracking”

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 37 records · Page 2

"Hypothetical" Heavy Particles Dynamics in LES of Turbulent Dispersed Two-Phase Channel Flow

The extensive experimental study of dispersed two-phase turbulent flow in a vertical channel has been performed in Eaton's research group in the Mechanical Engineering Department at Stanford University. In Wang & Squires (1996), this study motivated the validation of LES approach with Lagrangian tracking of round particles governed by drag forces. While the computed velocity of the flow have been predicted relatively well, the computed particle velocity differed strongly from the measured one. Using Monte Carlo simulation of inter-particle collisions, the computation of Yamamoto et al. (2001) was specifically performed to model Eaton's experiment. The results of Yamamoto et al. (2001) improved the particle velocity distribution. At the same time, Vance & Squires (2002) mentioned that the stochastic simualtion of inter-particle collisions is too expensive, requiring significantly more CPU resources than one needs for the gas flow computation. Therefore, the need comes to account for the inter-particle collisions in a simpler and still effective way. To present such a model in the framework of LES/Lagrangian particle approach, and to compare the calculated results with Eaton's measurement and modeling of Yamamoto is the main objective of the present paper.

Gorokhovski, M.↗

CFD Analysis of the 24-inch JIRAD Hybrid Rocket Motor

A series of multispecies, multiphase computational fluid dynamics (CFD) analyses of the 24-inch diameter joint government industry industrial research and development (JIRAD) hybrid rocket motor is described. The 24-inch JIRAD hybrid motor operates by injection of liquid oxygen (LOX) into a vaporization plenum chamber upstream of ports in the hydroxyl-terminated polybutadiene (HTPB) solid fuel. The injector spray pattern had a strong influence on combustion stability of the JIRAD motor so a CFD study was initiated to define the injector end flow field under different oxidizer spray patterns and operating conditions. By using CFD to gain a clear picture of the flow field and temperature distribution within the JIRAD motor, it is hoped that the fundamental mechanisms of hybrid combustion instability may be identified and then suppressed by simple alterations to the oxidizer injection parameters such as injection angle and velocity. The simulations in this study were carried out using the General Algorithm for Analysis of Combustion SYstems (GALACSY) multiphase combustion codes. GALACSY consists of a comprehensive set of droplet dynamic submodels (atomization, evaporation, etc.) and a computationally efficient hydrocarbon chemistry package built around a robust Navier-Stokes solver optimized for low Mach number flows. Lagrangian tracking of dispersed particles describes a closely coupled spray phase. The CFD cases described in this paper represent various levels of simplification of the problem. They include: (A) gaseous oxygen with combusting fuel vapor blowing off the walls at various oxidizer injection angles and velocities, (B) gaseous oxygen with combusting fuel vapor blowing off the walls, and (C) liquid oxygen with combusting fuel vapor blowing off the walls. The study used an axisymmetric model and the results indicate that the injector design significantly effects the flow field in the injector end of the motor. Markedly different recirculation patterns are observed in the vaporization chamber as the oxygen velocity and/or spray pattern is varied. The ability of these recirculation patterns to stabilize the diffusion flame above the surface of the solid fuel gives a plausible explanation for the experimentally determined combustion stability characteristics of the JIRAD motor, and suggests how combustion stability can be assured by modifications to the injector design.

Liang, Pak-Yan↗

A combined Eulerian-Lagrangian two-phase flow analysis of SSME HPOTP nozzle plug trajectories. II - Results

An analysis of Space Shuttle Main Engine high-pressure oxygen turbopump nozzle plug trajectories has been performed, using a Lagrangian method to track nozzle plug particles expelled from a turbine through a high Reynolds number flow in a turnaround duct with turning vanes. Axisymmetric and parametric analyses reveal that if nozzle plugs exited the turbine they would probably impact the LOX heat exchanger with impact velocities which are significantly less than the penetration velocity. The finding that only slight to moderate damage will result from nozzle plug failure in flight is supported by the results of a hot-fire engine test with induced nozzle plug failures.

Mcconnaughey, P. K.↗

A Lagrangian method for a two layer shallow water model

One of the biggest drawbacks of Eulerian methods is the presence of the nonlinear convective terms which leads to inaccurate representations of advection and fronts. An alternative and conceptually simpler approach is the Lagrangian approach in which the fluid particles themselves are tracked and equations are derived based on local spatial interaction. In this formulation the nonlinear convective terms do not appear. The main drawback with the Lagrangian formulation, however, is that the local spatial interactions are time dependent. That is, at each time level, one must know the current neighbors of a given particle to accurately compute the forces acting on the particle. The Voroni Mesh which is continuously deforming and always linking nearest neighbors is needed. The effectiveness of a simple Lagrangian method for simulating the onset of a Kelvin-Helmholtz instability in a 2 layer fluid is demonstrated.

Augenbaum, J. M.↗

Lagrangian Approach to Jet Mixing and Optimization of the Reactor for Production of Carbon Nanotubes

This study was motivated by an attempt to optimize the High Pressure carbon oxide (HiPco) process for the production of carbon nanotubes from gaseous carbon oxide, The goal is to achieve rapid and uniform heating of catalyst particles by an optimal arrangement of jets. A mixed Eulerian and Lagrangian approach is implemented to track the temperature of catalyst particles along their trajectories as a function of time. The FLUENT CFD software with second-order upwind approximation of convective terms and an algebraic multigrid-based solver is used. The poor performance of the original reactor configuration is explained in terms of features of particle trajectories. The trajectories most exposed to the hot jets appear to be the most problematic for heating because they either bend towards the cold jet interior or rotate upwind of the mixing zone. To reduce undesirable slow and/or oscillatory heating of catalyst particles, a reactor configuration with three central jets is proposed and the optimal location of the central and peripheral nozzles is determined.

Povitsky, Alex↗

Demonstration of Initial GlennICE Relative Frame Capability: Axial-Flow Propeller

Modifications have been implemented in the GlennICE software to accommodate a non-inertial reference frame. GlennICE accepts a flow solution from an external flow solver. It then introduces particles and tracks them through the flow field in a Lagrangian manner. Centrifugal and Coriolis terms were added to the GlennICE software to account for relative frame simulations. The objective of the present paper is twofold. First, to check that the new terms are implemented correctly and that the code still behaves as expected with respect to convergence. And second, to provide some initial insight into an upcoming propeller experiment in the NASA Icing Research Tunnel. The paper presents a description of the code modifications. In addition, results are presented for two operating conditions, and three particle sizes. Each case was simulated with four different grid densities to assess grid dependence.

Aircraft Icing↗

Demonstration of Initial GlennICE Relative Frame Capability: Axial-Flow Propeller

Modifications have been implemented in the GlennICE software to accommodate a non-inertial reference frame. GlennICE accepts a flow solution from an external flow solver. It then introduces particles and tracks them through the flow field in a Lagrangian manner. Centrifugal and Coriolis terms were added to the GlennICE software to account for relative frame simulations. The objective of the present paper is twofold. First, to check that the new terms are implemented correctly and that the code still behaves as expected with respect to convergence. And second, to provide some initial insight into an upcoming propeller experiment in the NASA Icing Research Tunnel. The paper presents a description of the code modifications. In addition, results are presented for two operating conditions, and three particle sizes. Each case was simulated with four different grid densities to assess grid dependence.

Aircraft Icing↗

A laser based computer aided non-intrusive technique for full field flow characterization in macroscopic curved channels

This paper presents the application of a laser based, computer aided image processing technique for the nonintrusive evaluation of field velocities and accelerations in a 180 deg curved channel geometry, characteristic of inlet passages of the cowls of airbreathing aircraft engines. The method based on the illumination of the test section with a laser generated planar sheet of light, allows microscopic or macroscopic surveillance of fluid flow across or along large cross sections for relatively long (greater than 1 sec) periods of time. The method is Lagrangian in concept and permits identification and tracking of the same particle, thus facilitating construction of more comprehensive trajectories and the calculation of velocities and accelerations.

Braun, M. J.↗

Full field flow visualization and computer-aided velocity measurements in a bank of cylinders in a wind tunnel

The full field flow tracking (FFFT) method that is presented in this paper uses a laser-generated, mechanically strobed planar sheet of light, a low luminosity TV camera coupled with a long distance microscope, and a computer-controlled video recorder to study nonintrusively and qualitatively the flow structures in a bank of cylinders that are placed in a wind tunnel. This setup simulates an upscale version of the geometry of internal cooling passageways characteristic of small air-cooled radial turbines. The qualitative images supplied by the FFFT system are processed by means of a computer-integrated image quantification (CIIQ) method into quantitative information, trajectories and velocities, that describe the flow upstream of and within the bank of cylinders. The tracking method is Lagrangian in concept, and permits identification and tracking of the same particle, thus facilitating construction of time dependent trajectories and the calculation of true velocities and accelerations. The error analysis evaluates the accuracy with which the seed particles follow the flow and the errors incurred during the quantitative processing of the raw data derived from the FFFT/CIIQ method.

Braun, M. J.↗

Evaluation of an efficient statistical transport model for turbulent droplet dispersion in dilute combusting sprays

Evaluation of an efficient statistical transport model for turbulent droplet dispersion is made for a dilute spray of methanol droplets injected into a turbulent, axisymmnetric methane-fuelled diffusion flame burning in stagnant air. In the dispersion model, a computational parcel representing a group of physical particles (droplets) is considered to have a normal probability density function (pdf) in space. The mean is determined by Lagrangian tracking through a sequence of stochastically generated turbulent eddies and the variance is evaluated from a statistical formulation based on the linearized particle equations of motion. The basic validity of this model is established through comparison with available experimental data and with theoretical predictions using a conventional stochastic direct modeling approach. The conclusion of the evaluation is that the proposed dispersion model compares favorably with experimental data and provides a valid technique for simulating turbulent combusting sprays with significant computational savings over conventional methods.

Litchford, Ron J.↗

A computationally efficient model for turbulent droplet dispersion in spray combustion

A novel model for turbulent droplet dispersion is formulated having significantly improved computational efficiency in comparison to the conventional point source stochastic sampling methodology. In the proposed model, a computational parcel representing a group of physical particles is considered to have a normal (Gaussian) probability density function (PDF) in three-dimensional space. The mean of each PDF is determined by Lagrangian tracking of each computational parcel, either deterministically or stochastically. The variance is represented by a turbulence-induced mean squared dispersion which is based on statistical inferences from the linearized direct modeling formulation for particle/eddy interactions. Convolution of the computational parcel PDF's produces a single PDF for the physical particle distribution profile. The validity of the new model is established by comparison with the conventional stochastic sampling method, where in each parcel is represented by a delta function distribution, for non-evaporating particles injected into simple turbulent air flows.

Litchford, Ron J.↗

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 cooling processes in a large-scale liquid hydrogen storage tank under normal gravity. The self-pressurization model utilizes the kinetics-based Schrage equation with the Volume-Of-Fluid (VOF) model to account for the evaporative and condensing interfacial mass fluxes. It uses a laminar approach to model natural convection and compares the pressures and temperatures predicted by the CFD model with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) self-pressurization experiment. Additionally, the predicted interfacial mass flux is presented. A CFD model for simulating pressure control using a spray bar is also presented. An Eulerian-Lagrangian approach models the interactions between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by tracking particles in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. The T-sat model is used to calculate the droplet-ullage heat and mass transfer. In this model, droplets warm up to the saturation temperature corresponding to the ullage vapor pressure thereafter evaporating while remaining at the saturation temperature. The tank pressure, vapor, and liquid temperature evolutions are validated against data from the MHTB spray bar mixing experiment.

Self-Pressurization↗

Computation of two-phase shear-layer flow using an Eulerian-Lagrangian analysis

A combined Eulerian-Lagrangian analysis which combines a linearized block implicit Navier-Stokes analysis for the continuous phase with a Lagrangian analysis for the discrete phase has been utilized to simulate a two-phase shear-layer flow in a two-dimensional channel flow configuration. The motion of the particles is tracked in computational coordinate space, resulting in improved computational efficiency, and interphase coupling terms for the Eulerian analysis are computed from the instantaneous distribution of the particles rather than the trajectory information. The use of the implicit Navier-Stokes analysis for the continuous phase has made it possible to efficiently use a highly stretched mesh. A low Reynolds number form of the k-epsilon turbulence model was used with near-wall resolution, thus eliminating the need to use the 'wall-function' approach. The streamwise velocity profiles for the continuous and the discrete phases have been compared with experimental data for two test cases.

Sabnis, Jayant S.↗

Modeling Droplet Heat and Mass Transfer during Spray Bar Pressure Control of the Multipurpose Hydrogen Test Bed (MHTB) Tank in Normal Gravity

A CFD model for simulating pressure control in cryogenic storage tanks through the injection of a subcooled liquid into the ullage is presented and applied to the 1g MHTB spray bar cooling experiments. An Eulerian-Lagrangian approach is utilized to track the spray droplets and capture the interaction between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. A new model for calculating the droplet-ullage heat and mass transfer is developed. In this model, a droplet is allowed to warm up to the saturation temperature corresponding to the ullage vapor pressure, after which it evaporates while remaining at the saturation temperature. The droplet model is validated against the results of the MHTB spray-bar cooling experiments with 50% and 90% tank fill ratios. The predictions of the present T-sat based model are compared with those of a previously developed kinetic-based droplet mass transfer model. The predictions of the two models regarding the evolving tank pressure and temperature distributions, as well as the droplets' trajectories and temperatures, are examined and compared in detail. Finally, the ullage pressure and local vapor and liquid temperature evolutions are validated against the corresponding data provided by the MHTB spray bar mixing experiment.

mass transfer↗

Stratospheric aircraft exhaust plume and wake chemistry studies

This report documents progress to date in an ongoing study to analyze and model emissions leaving a proposed High Speed Civil Transport (HSCT) from when the exhaust gases leave the engine until they are deposited at atmospheric scales in the stratosphere. Estimates are given for the emissions, summarizing relevant earlier work (CIAP) and reviewing current propulsion research efforts. The chemical evolution and the mixing and vortical motion of the exhaust are analyzed to track the exhaust and its speciation as the emissions are mixed to atmospheric scales. The species tracked include those that could be heterogeneously reactive on the surfaces of the condensed solid water (ice) particles and on exhaust soot particle surfaces. Dispersion and reaction of chemical constituents in the far wake are studied with a Lagrangian air parcel model, in conjunction with a radiation code to calculate the net heating/cooling. Laboratory measurements of heterogeneous chemistry of aqueous sulfuric acid and nitric acid hydrates are also described. Results include the solubility of HCl in sulfuric acid which is a key parameter for modeling stratospheric processing. We also report initial results for condensation of nitric acid trihydrate from gas phase H2O and HNO3.

Miake-Lye, R. C.↗

Ice Cloud Formation and Dehydration in the Tropical Tropopause Layer

Stratospheric water vapor is important not only for its greenhouse forcing, but also because it plays a significant role in stratospheric chemistry. several recent studies have focused on the potential for dehydration due to ice cloud formation in air rising slowly through the tropical tropopause layer. Holton and Gettelman showed that temperature variations associated with horizontal transport of air in the tropopause layer can drive ice cloud formation and dehydration, and Gettelman et al. recently examined the cloud formation and dehydration along kinematic trajectories using simple assumptions about the cloud properties. In this study, we use a Lagrangian, one-dimensional cloud model to further investigate cloud formation and dehydration as air is transported horizontally and vertically through the tropical tropopause layer. Time-height curtains of temperature are extracted from meteorological analyses. The model tracks the growth and sedimentation of individual cloud particles. The regional distribution of clouds simulated in the model is comparable to the subvisible cirrus distribution indicated by SAGE II. The simulated cloud properties depend strongly on the assumed ice supersaturation threshold for ice nucleation. with effective nuclei present (low supersaturation threshold), ice number densities are high (0.1--10 cm(circumflex)-3), and ice crystals do not grow large enough to fall very far, resulting in limited dehydration. With higher supersaturation thresholds, ice number densities are much lower (less than 0.01 cm(circumflex)-3), and ice crystals grow large enough to fall substantially; however, supersaturated air often crosses the tropopause without cloud formation. The clouds typically do not dehydrate the air along trajectories down to the temperature minimum saturation mixing ratio. Rather the water vapor mixing ratio crossing the tropopause along trajectories is typically 10-50% larger than the saturation mixing ratio.

Jensen, Eric↗

Ice Cloud Formation and Dehydration in the Tropical Tropopause Layer

Stratospheric water vapor is important not only for its greenhouse forcing, but also because it plays a significant role in stratospheric chemistry. Several recent studies have focused on the potential for dehydration due to ice cloud formation in air rising slowly through the tropical tropopause layer (TTL). Holton and Gettelman showed that temperature variations associated with horizontal transport of air in the TTL can drive ice cloud formation and dehydration, and Gettelman et al. recently examined the cloud formation and dehydration along kinematic trajectories using simple assumptions about the cloud properties. In this study, a Lagrangian, one-dimensional cloud model has been used to further investigate cloud formation and dehydration as air is transported horizontally and vertically through the TTL. Time-height curtains of temperature are extracted from meteorological analyses. The model tracks the growth, advection, and sedimentation of individual cloud particles. The regional distribution of clouds simulated in the model is comparable to the subvisible cirrus distribution indicated by SAGE II. The simulated cloud properties and cloud frequencies depend strongly on the assumed supersaturation threshold for ice nucleation. The clouds typically do not dehydrate the air along trajectories down to the temperature minimum saturation mixing ratio. Rather the water vapor mixing ratio crossing the tropopause along trajectories is 10-50% larger than the saturation mixing ratio. I will also discuss the impacts of Kelvin waves and gravity waves on cloud properties and dehydration efficiency. These simulations can be used to determine whether observed lower stratospheric water vapor mixing ratios can be explained by dehydration associated with in situ TTL cloud formation alone.

Jensen, Eric↗

Assessment of Oceanographic Conditions during the North Atlantic EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) Field Campaign

This manuscript presents an overview of NASA’s EXport Processes in the Ocean from Remote Sensing 2021 Field Campaign in the North Atlantic (EXPORTS NA) and provides quantitative and dynamical descriptions of the physical processes modulating water transformations during the study. A major programmatic goal was to conduct the sampling in a Lagrangian mode so that ocean ecological and biogeochemical changes can be observed independent from physical advective processes. To accomplish this goal, EXPORTS NA conducted a multi-ship, multi-asset field sampling program within a retentive, anticyclonic mode water eddy. Beneath depths of ∼100 m, Lagrangian sampling assets remained within the eddy core waters (ECWs) throughout the experiment, demonstrating that the ECWs within the mode water eddy were retentive. However, strong westerly winds from four storm events deepened the mixed layer (ML) of the surface core waters (SCWs) above the eddy’s mode water core by 25-40 m and exchanged some of the SCWs with surface waters outside of the eddy via Ekman transport. Estimates of flushing times ranged from 5-8 days, with surface exchange fractions ranging from 20-75% and were consistent with particle tracking advected by combined geostrophic and Ekman velocities. The relative contributions of horizontal and vertical advection on changes in ECW tracers depended on the horizontal and vertical gradients of that tracer. For example, horizontal advection played a large role in ECW salinity fluxes, while vertical entrainment played a larger role in the fluxes of nutrients into SCW ML. Each storm injected nutrients and low oxygen waters into the ML, after which the surface ocean ecosystem responded by reducing nutrient concentrations and increasing %O2 saturation levels. Overall, ECW values of chlorophyll and POC were the largest at the onset of the field program and decreased throughout the campaign. The analysis presented provides a physical oceanographic context for the many measurements made during the EXPORTS NA field campaign while illustrating the many challenges of conducting a production-flux experiment, even in a Lagrangian frame, and the inherent uncertainties of interpreting biological carbon pump observations that were collected in a Eulerian frame of reference.

oceanography↗