Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian 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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

The influence of dynamic inflow and torsional flexibility on rotor damping in forward flight from symbolically generated equations

The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The Helicopter Equations for Stability and Loads (HESL) program is extended to derive the governing equations of motion for the blade, and a Lagrangian formulation is used to obtain the equations in generalized coordinates. The program generates the steady-state and linearized perturbation equations in symbolic form and then codes them into FORTRAN subroutines. The coefficients for each equation and for each mode are identified through a numerical program; the latter can also be used to obtain the harmonic balance equations. The governing multiblade equations are derived explicitly using HESL. These equations can accommodate any number of elastic blade modes. Stability results are presented for several hingeless rotor blade structural models, and the influence of dynamic inflow in forward flight with an elastic hingeless rotor is investigated.

Reddy, T. S. R.↗
DOE OSTI2026

Lagrange-Remap strategy for multi-material fluid-solid simulations using compressive limiters

In the present work, the Lagrange-Remap strategy proposed in [1] is extended to multi-material fluid-solid simulations. Both hypo-elastic and hyper-elastic material models are considered to describe the mechanical behavior of the solids. In practice, the deviatoric stress tensor (for hypo-elastic materials) and the left Cauchy-Green tensor (for isotropic hyper-elastic materials) are remapped, while the use of compressive limiters effectively reduces numerical diffusion during the remapping step. The simplicity of this diffuse interface approach is emphasized in the context of multi-material fluid-solid simulations. A series of Lagrange-Remap test cases, involving both solids and fluids, are conducted and compared with reference Lagrangian simulations, demonstrating the robustness and accuracy of the overall numerical strategy.

Compressive limiters↗
NASA NTRS1994

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗
NASA NTRS1991

The formation of terrains antipodal to major impacts

The SALE, or Simplified Arbitrary Lagrangian Eulerian code is presently used to calculate the antipodal pressures generated by large impacts on icy satellites and planets of various compositions and core sizes. While in a high density core/lower density mantle planet there is a refraction of the incident pressure wave away from the antipode by the core, resulting in a low antipodal pressure, a planet with a low velocity core with a higher-wave speed mantle will focus the waveform to maximize antipodal pressure. Specific impact basins on the moon, Mercury, and several icy satellites are modeled to ascertain whether the impacts were capable of producing the observed antipodal features.

Watts, Alison W.↗
NASA NTRS2017

Estimating the Distribution of Colored Dissolved Organic Matter During the Southern Ocean Gas Exchange Experiment Using Four-Dimensional Variational Data Assimilation

We diagnosed the effect of various physical processes on the distribution of mixed-layer colored dissolved organic matter (CDOM) and a sulfur hexauoride (SF6) tracer during the Southern Ocean Gas Exchange Experiment (SO GasEx). The biochemical upper ocean state estimate uses in situ and satellite biochemical and physical data in the study region, including CDOM (absorption coefcient and spectral slope), SF6, hydrography, and sea level anomaly. Modules for photobleaching of CDOM and surface transport of SF6 were coupled with an ocean circulation model for this purpose. The observed spatial and temporal variations in CDOM were captured by the state estimate without including any new biological source term for CDOM, assuming it to be negligible over the 26 days of the state estimate. Thermocline entrainment and photobleaching acted to diminish the mixed-layer CDOM with time scales of 18 and 16 days, respectively. Lateral advection of CDOM played a dominant role and increased the mixed-layer CDOM with a time scale of 12 days, whereas lateral diffusion of CDOM was negligible. A Lagrangian view on the CDOM variability was demonstrated by using the SF6 as a weighting function to integrate the CDOM elds. This and similar data assimilation methods can be used to provide reasonable estimates of optical properties, and other physical parameters over the short-term duration of a research cruise, and help in the tracking of tracer releases in large-scale oceanographic experiments, and in oceanographic process studies.

Del Castillo, C. E.↗
NASA NTRS2002

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗
NASA NTRS1988

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

Statistical Distributions for Mesh Independent Solutions in ALEGRA

The representation of material heterogeneity (also referred to as "spatial variation") plays a key role in the material failure simulation method used in ALEGRA. ALEGRA is an arbitrary Lagrangian-Eulerian shock and multiphysics code developed at Sandia National Laboratories and contains several methods for incorporating spatial variation into simulations. A desirable property of a spatial variation method is that it should produce consistent stochastic behavior regardless of the mesh used (a property referred to as "mesh independence"). However, mesh dependence has been reported using the Weibull distribution with ALEGRA's spatial variation method. This report describes efforts towards providing additional insight into both the theory and numerical experiments investigating such mesh dependence. In particular, we have implemented a discrete minimum order statistic model with properties that are theoretically mesh independent.

36 MATERIALS SCIENCE↗
NASA NTRS1981

An objective method for determining the generalized transport tensor for two-dimensional Eulerian models

An objective method for deriving the components of a generalized transport tensor for a two-dimensional model is presented. Representative meridional and vertical velocities and thermodynamic scalars at a uniform grid are used to reduce the problem to the solution of two flux equations for two unknowns. One unknown is the stream-function, coefficient of an antisymmetric tensor, which corrects the Eulerian mean motions for Stokes drift; the other is a time constant, which converts the deviatory velocity tensor to a symmetric transport tensor. The complete asymmetric tensor, called a transport tensor, has a divergence which yields both advection and diffusion by the deviatory velocities. Advantages and disadvantages of Lagrangian and Eulerian averages are discussed, and meridional-vertical velocity correlations are provided.

Danielsen, E. F.↗
NASA NTRS1988

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented in an elementary fashion and includes the relationship with traditional point-vortex studies, the convergence to smooth solutions of the Euler equations, and the essential differences between two- and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. The overlap with the excellent review articles available is kept to a minimum and more emphasis is placed on the area of expertise, namely two-dimensional flows around bluff bodies. When solid walls are present, complete mathematical models are not available and a more heuristic attitude must be adopted. The imposition of inviscid and viscous boundary conditions without conformal mappings or image vortices and the creation of vorticity along solid walls are examined in detail. Methods for boundary-layer treatment and the question of the Kutta condition are discussed. Practical aspects and tips helpful in creating a method that really works are explained. The topics include the robustness of the method and the assessment of accuracy, vortex-core profiles, timemarching schemes, numerical dissipation, and efficient programming. Calculations of flows past streamlined or bluff bodies are used as examples when appropriate.

Spalart, Philippe R.↗
NASA NTRS1992

Thermal mechanical analysis of sprag clutches

Work done at Case Western Reserve University on the Thermal Mechanical analysis of sprag helicopter clutches is reported. The report is presented in two parts. The first part is a description of a test rig for the measurement of the heat generated by high speed sprag clutch assemblies during cyclic torsional loading. The second part describes a finite element modeling procedure for sliding contact. The test rig provides a cyclic torsional load of 756 inch-pounds at 5000 rpm using a four-square arrangement. The sprag clutch test unit was placed between the high speed pinions of the circulating power loop. The test unit was designed to have replaceable inner ad outer races, which contain the instrumentation to monitor the sprag clutch. The torque loading device was chosen to be a water cooled magnetic clutch, which is controlled either manually or through a computer. In the second part, a Generalized Eulerian-Lagrangian formulation for non-linear dynamic problems is developed for solid materials. This formulation is derived from the basic laws and axioms of continuum mechanics. The novel aspect of this method is that we are able to investigate the physics in the spatial region of interest as material flows through it without having to follow material points. A finite element approximation to the governing equations is developed. Iterative Methods for the solution of the discrete finite element equations are explored. A FORTRAN program to implement this formulation is developed and a number of solutions to problems of sliding contact are presented.

Mullen, Robert L.↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

The interaction of radiative and dynamical processes during a simulated sudden stratospheric warming

An analysis of a spontaneous sudden stratospheric warming that occurred during a 2-year integration of the Langley Research Center (LaRC) Atmospheric Simulation Model is presented. The simulated warming resembles observed 'wave 1' warmings in the Northern Hemisphere stratosphere and provides an opportunity to investigate the radiative and dynamical processes occurring during the warming event. Isentropic analysis of potential vorticity sources and sinks indicates that dynamically induced departures from radiative equilibrium play an important role in the warming event. Enhanced radiative cooling associated with a series of upper stratospheric warm pools leads to radiative dampening within the polar vortex. Within the 'surf zone' large-scale radiative cooling leads to diabatic advection of high potential vorticity air from aloft. Lagrangian area diagnostics of the simulated warming agree well with Limb Infrared Monitor of the Stratosphere (LIMS) analyses. Dynamical mixing is shown to account for the majority of the decrease in the size of the polar vortex during the simulated warming. An investigation of the nonlinear deformation of material lines that are initially coincident with diagnosed potential vorticity isopleths is conducted to clarify the relationship between the Lagrangian area diagnostics and potential vorticity advection during wave breaking events.

Pierce, R. B.↗
NASA NTRS1992

The appliation of potential CFD methods to helicopter hover flows

Fixed-wing code development is now aimed primarily at the solution of problems dominated by separation--based on the assumptions that the ability to solve such problems implies the ability to solve all other problems and that present inviscid method are already adequate for most other problems. Neither of the above assumptions are correct for rotary wing problems. This is because of the unique and overriding importance of wake modeling to rotor problems and also due to the well-known numerical diffusion problems which convectional Eulerian Computational fluid dynamics (CFD) method encounter when called on to convect strong vortical regions for long distances. The need for accurate wake analyses is probably the most fundamental difference between rotory and fixed-wing aerodynamics. In addition, rotary wing complexity requires a much more intimate relationship between test and analysis than is common in fixed-wing work. With these issues in mind, this paper will review some of our recent experience in using a unique-Eulerian-Lagrangian Computational fluid dynamics (CFC) method for the solution of a critical rotor-wake problem--the prediction of hover performance.

Caradonna, F. X.↗
NASA NTRS1981

The secondary flow near a baroclinic planetary wave critical line

A critical line (CL) is the surface where the phase speed of a wave in a fluid is equal to the speed of the background flow. The considered investigation is concerned with one aspect of the simplest model of a CL in which the CL is assumed to totally absorb energy from steady, stationary, planetary waves. The aspect of interest is the secondary mean circulation near a CL in a baroclinic atmosphere. The motivation for this study is the observation of the nearly vertical CL by O'Neill and Taylor (1979) which appeared during the sudden warming of 1976/77. Even though the treatment of the CL is highly idealized in the investigation, there is evidence which indicates a very large rate of change in the zonally averaged temperature along a CL may occur. The Lagrangian-mean properties of an idealized baroclinic CL are also examined. It is found that the Lagrangian jets may provide an important transport process for exchange of stratospheric and tropospheric air.

Schoeberl, M. R.↗
NASA NTRS2004

LSPRAY-II: A Lagrangian Spray Module

LSPRAY-II is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray because of its importance in aerospace application. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-II, we have advanced the state-of-the-art in spray computations in several important ways.

Raju, M. S.↗
NASA NTRS2008

LSPRAY-III: A Lagrangian Spray Module

LSPRAY-III is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray because of its importance in aerospace application. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-III, we have advanced the state-of-the-art in spray computations in several important ways.

Raju, M. S.↗
DOE OSTI2026

Evaluating the Collision‐Coalescence Process in Idealized Cloud Convection Using Large‐Eddy Simulations With Lagrangian Microphysics

Drizzle initiation through the collision and coalescence of cloud droplets plays a crucial role in warm cloud precipitation. Recent theoretical studies suggest that the influence of collisional growth on the droplet size distribution can be quantified by a non-dimensional drizzle number (Dz). Here, large-eddy simulations with Lagrangian microphysics are employed to evaluate the theory by simulating a tall convection-cloud chamber under various conditions. Results show that the smaller the Dz, the larger the impact of collisions on the right tail of the droplet size distribution, consistent with the theory. The simulations confirm that the collision rate can be estimated from the droplet size distribution interacting only with cloud droplets of the same size at the mode radius. This suggests that the idealized theory can be a useful tool to design a cloud chamber for drizzle investigation, as well as to represent drizzle formation in models of real atmospheric clouds.

54 ENVIRONMENTAL SCIENCES↗