Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian Particles”

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 109 records · Page 6

On the Lagrangian Method for Steady and Unsteady Flow

A new and general Lagrangian formulation of fluid motion is given in which the independent variables are three material functions and a Lagrangian time, which differs for different fluid particles and is distinct from the Eulerian time. For steady flow it requires only three independent variables - the Lagrangian time and two stream functions - in contrast with the conventional Lagrangian formulation which apparently still requires four independent variables for describing a steady flow. This places the Lagrangian formulation for steady flow on the same footing as the Eulerian. For unsteady flow, the new formulation includes the conventional formulation as a special case when the Lagrangian time is identified with the Eulerian time and when the material functions are taken to be the fluid particle's position at some given time. The distinction between the Lagrangian and Eulerian time, however, is found useful in applications, e.g., to problems involving a free boundary.

Hui, W. H.↗

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.↗

Three-dimensional relativistic field-electron interaction in a multicavity high-power klystron. 1: Basic theory

A theoretical investigation of three dimensional relativistic klystron action is described. The relativistic axisymmetric equations of motion are derived from the time-dependent Lagrangian function for a charged particle in electromagnetic fields. An analytical expression of the fringing RF electric and magnetic fields within and in the vicinity of the interaction gap and the space-charge forces between axially and radially elastic deformable rings of charges are both included in the formulation. This makes an accurate computation of electron motion through the tunnel of the cavities and the drift tube spaces possible. Method of analysis is based on Lagrangian formulation. Bunching is computed using a disk model of electron stream in which the electron stream is divided into axisymmetric disks of equal charge and each disk is assumed to consist of a number of concentric rings of equal charges. The Individual representative groups of electrons are followed through the interaction gaps and drift tube spaces. Induced currents and voltages in interacting cavities are calculated by invoking the Shockley-Ramo theorem.

Kosmahl, H. G.↗

Modeling of Water Injection into a Vacuum

A loosely coupled two-phase vacuum water plume model has been developed. This model consists of a source flow model to describe the expansion of water vapor, and the Lagrangian equations of motion for particle trajectories. Gas/Particle interaction is modeled through the drag force induced by the relative velocities. Particles are assumed traveling along streamlines. The equations of motion are integrated to obtain particle velocity along the streamline. This model has been used to predict the mass flux in a 5 meter radius hemispherical domain resulting from the burst of a water jet of 1.5 mm in diameter, mass flow rate of 24.2 g/s, and stagnation pressure of 21.0 psia, which is the nominal Orbiter water dump condition. The result is compared with an empirical water plume model deduced from a video image of the STS-29 water dump. To further improve the model, work has begun to numerically simulate the bubble formation and bursting present in a liquid stream injected into a vacuum. The technique of smoothed particle hydrodynamics was used to formulate this simulation. A status and results of the on-going effort are presented and compared to results from the literature.

Alred, John W.↗

Comparison of ALE and SPH Simulations of Vertical Drop Tests of a Composite Fuselage Section into Water

Simulation of multi-terrain impact has been identified as an important research area for improved prediction of rotorcraft crashworthiness within the NASA Subsonic Rotary Wing Aeronautics Program on Rotorcraft Crashworthiness. As part of this effort, two vertical drop tests were conducted of a 5-ft-diameter composite fuselage section into water. For the first test, the fuselage section was impacted in a baseline configuration without energy absorbers. For the second test, the fuselage section was retrofitted with a composite honeycomb energy absorber. Both tests were conducted at a nominal velocity of 25-ft/s. A detailed finite element model was developed to represent each test article and water impact was simulated using both Arbitrary Lagrangian Eulerian (ALE) and Smooth Particle Hydrodynamics (SPH) approaches in LS-DYNA, a nonlinear, explicit transient dynamic finite element code. Analytical predictions were correlated with experimental data for both test configurations. In addition, studies were performed to evaluate the influence of mesh density on test-analysis correlation.

Jackson, Karen E.↗

The narrow rings of Jupiter, Saturn and Uranus

The origin of the newly discovered narrow ring systems around Jupiter, Saturn and Uranus is considered. It is pointed out that both the Uranian and Jovian ring systems have mean orbital radii of 1.8 planetary radii and lie within the Roche zones of their respective planets, and it is suggested that the Jovian ring is the product of the disintegration of a satellite that entered the Roche zone, and that large numbers of small particles are now in horseshoe orbits about the Lagrangian equilibrium points of the remnant chunks. Analysis of the path of a ring particle in a horseshoe orbit is shown to result in ring structures in agreement with those observed for the circular rings of Jupiter and the highly eccentric ring of Uranus. The stability of these ring systems is then considered, and it is suggested that the F ring of Saturn, which lies outside the Roche zone, represents primordial matter not yet accreted by small satellites just inside the Mimas first-order resonances.

Dermott, S. F.↗

Assessment of Particle Drag Models Using GBK Experimental Data

Developing a simulation capability to model dust-laden flows requires the ability to accurately predict the velocity of a particle as it travels through a shock layer, nozzle plume, or the flow of an experimental facility. The primary mechanism that determines the particle velocity is the drag force acting on the particle. The drag force is typically expressed in terms of the frontal area of the particle, its velocity relative to the surrounding fluid, the density of the surrounding fluid, and a non-dimensional drag coefficient. As part of the ESM project, an integrated CFD-particle trajectory code named US3D-DUST is being developed]. This code uses a Lagrangian-based framework that computes the trajectories of individual particles in the flow. Several of the most commonly-used spherical drag models have been implemented in US3D-DUST. In this study, the US3D-DUST code will be used to simulate particle flows in the GBK facility. Differences in particle velocities predicted using the different drag models will be compared against recent GBK data

Dusty flow↗

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

Hydrometeor development in cold clouds in FIRE

The microphysical characteristics of cold clouds were examined using aircraft. The clouds investigated ranged in temperature from minus 20 to minus 46 C. Liquid water was detected in these clouds at minus 35 C and may exist at even colder temperatures. Aircraft patterns in the form of Lagrangian spirals were used to interpret particle growth processes. Significant broadening of the particle size spectra is observed with minor changes in the spectra at small sizes. Virtually all of the broadening observed is attributable to ice particle aggregation which occurs at all temperatures. The crystals comprising the aggregates are of comparable size, joined either at tips or edges, and are usually 2 in number at temperatures lower than minus 25 C, increasing to 3 or more at warmer temperatures. The data strongly suggest that sintering is the mechanism through which the crystals aggregate. Aggregation appears to be important in the transfer of water mass from upper to lower levels in clouds.

Heymsfield, Andrew J.↗

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.↗

Dynamics of the baryonic component in hierarchical clustering universes

I present self-consistent 3-D simulations of the formation of virialized systems containing both gas and dark matter in a flat universe. A fully Lagrangian code based on the Smoothed Particle Hydrodynamics technique and a tree data structure has been used to evolve regions of comoving radius 2-3 Mpc. Tidal effects are included by coarse-sampling the density of the outer regions up to a radius approx. 20 Mpc. Initial conditions are set at high redshift (z greater than 7) using a standard Cold Dark Matter perturbation spectrum and a baryon mass fraction of 10 percent (omega(sub b) = 0.1). Simulations in which the gas evolves either adiabatically or radiates energy at a rate determined locally by its cooling function were performed. This allows us to investigate with the same set of simulations the importance of radiative losses in the formation of galaxies and the equilibrium structure of virialized systems where cooling is very inefficient. In the absence of radiative losses, the simulations can be rescaled to the density and radius typical of galaxy clusters. A summary of the main results is presented.

Navarro, Julio↗

A combined Eulerian-Lagrangian two-phase analysis of the SSME HPOTP nozzle plug trajectories

As a result of high cycle fatigue, hydrogen embrittlement, and extended engine use, it was observed in testing that the trailing edge on the first stage nozzle plug in the High Pressure Oxygen Turbopump (HPOTP) could detach. The objective was to predict the trajectories followed by particles exiting the turbine. Experiments had shown that the heat exchanger soils, which lie downstream of the turbine, would be ruptured by particles traveling in the order of 360 ft/sec. An axisymmetric solution of the flow was obtained from the work of Lin et. al., who used INS3D to obtain the solution. The particle trajectories were obtained using the method of de Jong et. al., which employs Lagrangian tracking of the particle through the Eulerian flow field. The collision parameters were obtained from experiments conducted by Rocketdyne using problem specific alloys, speeds, and projectile geometries. A complete 3-D analysis using the most likely collision parameters shows maximum particle velocities of 200 ft/sec. in the heat exchanger region. Subsequent to this analysis, an engine level test was conducted in which seven particles passed through the turbine but no damage was observed on the heat exchanger coils.

Garcia, Robert↗

Gas-Particle Interaction Model Development in 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 GCD PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling portion of the GCD PSI project. The objectives of this task are to investigate and construct models of particle drag and gas-particle cloud interactions leading to what is known as particle turbulent kinetic energy (PTKE). While drag on a lone particle suspended in a flow gas of gas is well-known, the drag and drag-induced dynamics of a cloud of particles in a gas flow are not. The PSC element of the PSI project requires models of gas-particle cloud interactions to implement in the Gas Granular Flow Solver (GGFS) in order to arrive at a predictive simulation capability for PSI-induced soil erosion, cratering and ejecta transport. Experiments of dilute mixtures of soil particles in supersonic gas flow are being conducted at Johns Hopkins University in order to provide a dataset describing gas-particle cloud interactions. Augmented by a separate funding source are efforts to perform small-scale particle-resolved direct numerical simulations (PR-DNS) and larger scale Eulerian-Lagrangian simulations. Together, these experiments and simulations are being used to devise models of particle cloud drag and particle cloud drag-induced dynamics. The final product of these efforts will be particle drag models suitable for implementation into the GGFS application and a PTKE model suitable for the same purpose.

Jeff West↗

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↗

Dynamics of narrow rings

The ring models described here were developed to account for the dynamical problems posed by the narrow rings of Uranus. Some of these rings are now known to be eccentric, inclined, nonuniform in width, optically thick, and narrow, with very sharp edges. The eccentric rings have common pericenters and large, positive eccentricity gradients. The theory of shepherding satellites successfully accounts for most of these features and can also account for some features of the narrow Saturnian rings, in particular, waves, kinks, and periodic variations in brightness. Outstanding problems include the putative relation between eccentricity and inclination displayed by eight of the nine Uranian rings, and the magnitudes of the tidal torques acting on the shepherding satellites. The horseshoe-orbit model, although viable, probably has more application to the narrow rings from which the Saturnian coorbital satellites formed. The angular momentum flow rate due to particle collisions is a minimum at the Lagrangian equilibrium points L(4) and L(5), and one can expect accretion to be rapid at these points.

Dermott, S. F.↗

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗

Decay of a coherent scalar disturbance in a turbulent flow

The time evolution of an initially coherent, sinusoidal passive-scalar disturbance is considered when the wavelength q is less than the length scale of the surrounding isotropic turbulent flow. In 64 sup 3 direct numerical simulations a Gaussian prescription for the average scalar amplitude breaks down after a timescale associated with the wavenumber of the disturbance and there is a transition to a new characteristic decay. The Gaussian prescription is given by exp(-(1/2) q-squared w(t)), where a form for w(t), the Lagrangian mean square displacement of a single fluid particle, is proposed. After the transition the decay is given by exp(-t/tau), where tau is the new characteristic timescale. If q k(sub K), then 1/tau = 1/tau(sub D) + 1/tau(sub K), where k(sub K) is the Kolmogorov wavenumber, tau(sub D) is the diffusive timescale and tau(sub K) is the Kolmogorov timescale. An experiment originally proposed by de Gennesis considered in which the evolution of a coherent laser-induced pattern is read by a diffracting laser. The theory of this experiment involves the dispersion of particle pairs, but it is shown that in a certain limit it reduces to the single Fourier-mode problem and can be described in terms of single particle diffusion. The decay of a single mode after the transition in the simulation best describes the experiment.

Kerr, R. M.↗

Planetary engineering

Assuming commercial fusion power, heavy lift vehicles and major advances in genetic engineering, the authors survey possible late-21st century methods of working major transformations in planetary environments. Much more Earthlike climates may be produced on Mars by generating low freezing point greenhouse gases from indigenous materials; on Venus by biological conversion of CO2 to graphite, by canceling the greenhouse effect with high-altitude absorbing fine particles, or by a sunshield at the first Lagrangian point; and on Titan by greenhouses and/or fusion warming. However, in our present state of ignorance we cannot guarantee a stable endstate or exclude unanticipated climatic feedbacks or other unintended consequences. Moreover, as the authors illustrate by several examples, many conceivable modes of planetary engineering are so wasteful of scarce solar system resources and so destructive of important scientific information as to raise profound ethical issues, even if they were economically feasible, which they are not. Global warming on Earth may lead to calls for mitigation by planetary engineering, e.g., emplacement and replenishment of anti-greenhouse layers at high altitudes, or sunshields in space. But here especially we must be concerned about precision, stability, and inadvertent side-effects. The safest and most cost-effective means of countering global warming - beyond, e.g., improved energy efficiency, CFC bans and alternative energy sources - is the continuing reforestation of approximately 2 times 107 sq km of the Earth's surface. This can be accomplished with present technology and probably at the least cost.

Pollack, James B.↗