Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian methods”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

An investigation of using an RQP based method to calculate parameter sensitivity derivatives

Estimation of the sensitivity of problem functions with respect to problem variables forms the basis for many of our modern day algorithms for engineering optimization. The most common application of problem sensitivities has been in the calculation of objective function and constraint partial derivatives for determining search directions and optimality conditions. A second form of sensitivity analysis, parameter sensitivity, has also become an important topic in recent years. By parameter sensitivity, researchers refer to the estimation of changes in the modeling functions and current design point due to small changes in the fixed parameters of the formulation. Methods for calculating these derivatives have been proposed by several authors (Armacost and Fiacco 1974, Sobieski et al 1981, Schmit and Chang 1984, and Vanderplaats and Yoshida 1985). Two drawbacks to estimating parameter sensitivities by current methods have been: (1) the need for second order information about the Lagrangian at the current point, and (2) the estimates assume no change in the active set of constraints. The first of these two problems is addressed here and a new algorithm is proposed that does not require explicit calculation of second order information.

Beltracchi, Todd J.↗
NASA NTRS1990

Elliptic grid generation with orthogonality and spacing control on an arbitrary number of boundaries

A procedure for the generation of two and quasi-three-dimensional grids with control of orthogonality and spacing with respect to any and/or all boundaries of the domain is described. The elliptic grid generation equations of Thompson are solved implicitly. Control of the grid behavior is achieved through the introduction of forcing functions terms in the manner of Steger and Sorenson or in a modification of the method of Hilgenstock. The forcing function terms are constructed on the boundaries and propagated into the domain using transfinite Lagrangian bivariate interpolation. An anisotropic transfinite stencil is introduced and is shown to produce excellent grid behavior particularly in the vicinity of corner singularities. Emphasis is placed on the generation of viscous grids and the method is shown to be suited for use in the generation of grids for internal as well as external flow geometries. A FORTRAN program named PISCES has been written to implement the algorithm. Examples of grids for internal and external flows are given that highlight the characteristics and behavior of the algorithm.

White, J. A.↗
NASA NTRS1987

Linear optimum trim solution for aircraft with three longitudinal control effectors

A simple and fast linear optimum method for trimming aircraft with a combination of three longitudinal aerodynamic or propulsive control effectors is developed, with application to preliminary design or in-flight scheduling. A Lagrangian formulation is used to determine the longitudinal load distribution which minimizes trim drag. Trim drag results obtained with the linear optimum method are similar to those obtained with a nonlinear optimum method, and significant trim drag reductions are found using either optimum method, in comparison with results obtained using various ad hoc methods.

Goodrich, Kenneth H.↗
NASA NTRS1972

Experimental determination of the turbulence in a liquid rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a liquid rocket combustion chamber were determined experimentally using the tracer gas diffusion method. The results indicate that the turbulent diffusion process can be adequately modeled by the one-dimensional Taylor theory; however, the numerical values show significant disagreement with previously accepted values. The intensity of turbulence is higher by a factor of about two, while the Lagrangian correlation coefficient which was assumed to be unity in the past is much less than unity.

Hara, J.↗
DOE OSTI2025

General field evaluation in high-order meshes on GPUs

Robust and scalable function evaluation at any arbitrary point in the finite/spectral element mesh is required for querying the partial differential equation solution at points of interest, comparison of solution between different meshes, and Lagrangian particle tracking. This is a challenging problem, particularly for high-order unstructured meshes partitioned in parallel with MPI, as it requires identifying the element that overlaps a given point and computing the corresponding reference space coordinates. Here, we present a robust and efficient technique for general field evaluation in large-scale high-order meshes with quadrilaterals and hexahedra. In the proposed method, a combination of globally partitioned and processor-local maps are used to first determine a list of candidate MPI ranks, and then locally candidate elements that could contain a given point. Next, element-wise bounding boxes further reduce the list of candidate elements. Finally, Newton’s method with trust region is used to determine the overlapping element and corresponding reference space coordinates. Since GPU-based architectures have become popular for accelerating computational analyses using meshes with tensor-product elements, specialized kernels have been developed to utilize the proposed methodology on GPUs. The method is also extended to enable general field evaluation on surface meshes. The paper concludes by demonstrating the use of the proposed method in various applications ranging from mesh-to-mesh transfer during r-adaptivity to Lagrangian particle tracking.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Improved perfect-fluid energy-momentum tensor with spin in Einstein-Cartan space-time

The description of the spin given here is classical in that it is intrinsic but not quantized. The approach in this matter is similar to, for example, the work of Bailey and Israel (1973, 1975, 1979), where the fluid particles, which have intrinsic spin, may be galaxies or clusters of galaxies. The elementary particles of these objects and the 'ferromagnetic alignment' of their quantum spins are not resorted to in order to describe a fluid with spin. Physically this means that the equation of motion for the spin tensor is a modified Fermi-Walker transport equation (Misner et al., 1973), arising as a direct result of the inclusion of spin as an intrinsic variable in the thermodynamic description of the internal energy. The variables in this description are classical variables throughout and are not microscopic fields. An improved perfect-fluid energy-momentum tensor that includes spin and torsion is presented. Use is made of a Lagrangian variational principle based on the tetrad formalism of Halbwach (1960) and the method od constraints of Ray (1972).

Ray, J. R.↗
DOE OSTI2024

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

Life Cycle of Midlatitude Deep Convective Systems in a Lagrangian Framework

Deep Convective Systems (DCSs) consist of intense convective cores (CC), large stratiform rain (SR) regions, and extensive non-precipitating anvil clouds (AC). This study focuses on the evolution of these three components and the factors that affect convective AC production. An automated satellite tracking method is used in conjunction with a recently developed multi-sensor hybrid classification to analyze the evolution of DCS structure in a Lagrangian framework over the central United States. Composite analysis from 4221 tracked DCSs during two warm seasons (May-August, 2010-2011) shows that maximum system size correlates with lifetime, and longer-lived DCSs have more extensive SR and AC. Maximum SR and AC area lag behind peak convective intensity and the lag increases linearly from approximately 1-hour for short-lived systems to more than 3-hours for long-lived ones. The increased lag, which depends on the convective environment, suggests that changes in the overall diabatic heating structure associated with the transition from CC to SR and AC could prolong the system lifetime by sustaining stratiform cloud development. Longer-lasting systems are associated with up to 60% higher mid-tropospheric relative humidity and up to 40% stronger middle to upper tropospheric wind shear. Regression analysis shows that the areal coverage of thick AC is strongly correlated with the size of CC, updraft strength, and SR area. Ambient upper tropospheric wind speed and wind shear also play an important role for convective AC production where for systems with large AC (radius greater than 120-km) they are 24% and 20% higher, respectively, than those with small AC (radius=20 km).

Feng, Zhe↗
NASA NTRS1989

Recovering the full velocity and density fields from large-scale redshift-distance samples

A new method for extracting the large-scale three-dimensional velocity and mass density fields from measurements of the radial peculiar velocities is presented. Galaxies are assumed to trace the velocity field rather than the mass. The key assumption made is that the Lagrangian velocity field has negligible vorticity, as might be expected from perturbations that grew by gravitational instability. By applying the method to cosmological N-body simulations, it is demonstrated that it accurately reconstructs the velocity field. This technique promises a direct determination of the mass density field and the initial conditions for the formation of large-scale structure from galaxy peculiar velocity surveys.

Bertschinger, Edmund↗
NASA NTRS1993

Clouding tracing: Visualization of the mixing of fluid elements in convection-diffusion systems

This paper describes a highly interactive method for computer visualization of the basic physical process of dispersion and mixing of fluid elements in convection-diffusion systems. It is based on transforming the vector field from a traditionally Eulerian reference frame into a Lagrangian reference frame. Fluid elements are traced through the vector field for the mean path as well as the statistical dispersion of the fluid elements about the mean position by using added scalar information about the root mean square value of the vector field and its Lagrangian time scale. In this way, clouds of fluid elements are traced and are not just mean paths. We have used this method to visualize the simulation of an industrial incinerator to help identify mechanisms for poor mixing.

Ma, Kwan-Liu↗
NASA NTRS1984

Lagrangian least-squares prediction of solar flux (F-bar 10.7)

The result of a study on the application of an improved statistical prediction method for estimating the intermediate-term (months) and long-term (years) behavior of solar flux is discussed. The study indicates that better predictions, in a chi square sense, are possible by selecting sets of the solar flux data such that each set (cycle) starts and ends at the maxima (or minima) for the data base and initialization point of the procedure. Then one applies a Lagrangian least-squares statistical technique. Evidence is also presented to support the existence of an aperiodic variation in the periods as well as the amplitudes.

Holland, R. L.↗
DOE OSTI2024

Source Levels of In‐Cloud Air in Shallow Cumulus: Consistency Between Paluch Diagram and Lagrangian Particle Tracking

Abstract The Paluch diagram is a widely used tool for interpreting aircraft measurements of shallow cumulus clouds. A prior study conducted by Heus et al. (2008,https://doi.org/10.1175/2008jas2572.1) concluded that the source levels of in‐cloud air inferred from the Paluch diagram exhibit biases, sometimes of several hundred meters, in comparison to those derived from Lagrangian particle tracking. In this short study we revisit this comparison. The results indicate that the upper source levels of in‐cloud air determined from the Lagrangian Particle Tracking and the Paluch diagram are consistent, and the choice of statistical methods is crucial. The significance of this research lies in confirming the reliability of the Paluch analysis, enabling its confident application to aircraft data.

Meteorology & Atmospheric Sciences↗
NASA NTRS1982

Environmental solid particle effects on compressor cascade performance

The effect of suspended solid particles on the performance of the compressor cascade was investigated experimentally in a specially built cascade tunnel, using quartz sand particles. The cascades were made of NACA 65(10)10 airfoils. Three cascades were tested, one accelerating cascade and two diffusing cascades. The theoretical analysis assumes inviscid and incompressible two dimensional flow. The momentum exchange between the fluid and the particle is accounted for by the interphase force terms in the fluid momentum equation. The modified fluid phase momentum equations and the continuity equation are reduced to the conventional stream function vorticity formulation. The method treats the fluid phase in the Eulerian system and the particle phase in Lagrangian system. The experimental results indicate a small increase in the blade surface static pressures, while the theoretical results indicate a small decrease. The theoretical analysis, also predicts the loss in total pressure associated with the particulate flow through the cascade.

Tabakoff, W.↗
NASA NTRS1983

Optimal mistuning for enhanced aeroelastic stability of transonic fans

An inverse design procedure was developed for the design of a mistuned rotor. The design requirements are that the stability margin of the eigenvalues of the aeroelastic system be greater than or equal to some minimum stability margin, and that the mass added to each blade be positive. The objective was to achieve these requirements with a minimal amount of mistuning. Hence, the problem was posed as a constrained optimization problem. The constrained minimization problem was solved by the technique of mathematical programming via augmented Lagrangians. The unconstrained minimization phase of this technique was solved by the variable metric method. The bladed disk was modelled as being composed of a rigid disk mounted on a rigid shaft. Each of the blades were modelled with a single tosional degree of freedom.

Hall, K. C.↗
NASA NTRS1993

Supersonic gasdispersional jets and jet noise

This paper examines the potential for controlling jet noise radiation using methods developed for modifying jet infra-red thermal radiation. The control of jet noise may be possible by properly adding different solid and liquid particles into the jet flow and by using special nozzle shapes to change the jet exhaust flow structure. The numerical methods used to achieve these objectives are outlined in this paper. A combined Lagrangian-Eulerian approach is used to numerically simulate a Jet flow with particle addition. The unsteady behavior of jet impingement is examined. The techniques for grid and boundary condition definition are discussed as related to the accuracy of the calculations. Preliminary comparisons to experimental data are presented.

Gilinsky, M. M.↗
NASA NTRS2003

Where was the Iron Synthesized in Cassiopeia A?

We investigate the properties of Fe-rich knots on the east limb of the Cassiopeia A supernova remnant observed with Chandra/AXAF CCD Imaging Spectrometer (ACIS). Using analysis methods developed in a companion paper, we constrain the ejecta density profile and the Lagrangian mass coordinates of the knots from their fitted ionization age and electron temperature. Fe-rich knots which also have strong emission from Si, S, Ar, and Ca are clustered around mass coordinates q approx. equal to 0.35 - 0.4 in the shocked ejecta of 2 solar masses; this places them 0.7 - 0.8 solar masses out from the center (or 2 - 2.1 solar masses, allowing for the mass of a compact object). We also find an Fe clump that is evidently devoid of line emission from lower mass elements, as would be expected for a region that had undergone alpha-rich freeze out. This clump has a similar mass coordinate to the other Fe knots.

Hwang, Una↗
NASA NTRS1992

Radiative equilibrium in rapidly expanding shells

The equation of radiative equilibrium for rapidly expanding spherical shells including the special relativistic terms is discussed, and a temperature-correction procedure for the solution of this equation is given. The temperature-correction scheme combines a partial linearization approach with the approximate Lambda-operator iteration method. For some test problems, the temperature structures of atmospheres in radiative equilibrium in the Lagrangian frame are calculated. The effects of the velocity field on the temperature stratification and the luminosity are also discussed. We find that the temperature structure of a moving photosphere in radiative equilibrium in the Lagrangian frame differs significantly from that of a static photosphere and that, at least, all first-order terms, in particular the advection and aberration terms, have to be included consistently both in the radiative transfer and the energy equation in studies of supernovae atmospheres and other high-speed radiating flows.

Hauschildt, P. H.↗
NASA NTRS1991

General flow field analysis methods for helicopter rotor aeroacoustics

Previous work in the analysis of rotor flow fields for aeroacoustic applications involved the preliminary development of an efficient and accurate Lagrangian simulation of the unsteady vorticity field in the vicinity of helicopter main rotor that could analyze a limited class of rotor/wake interactions. The capabilities of this analysis have subsequently been considerably enhanced to allow it to serve as the foundation for a general analysis of the rotor/wake interaction noise. This paper presents the details of these enhancements, which focus on the expansion of the reconstruction approach developed previously to handle arbitrary vortex wake interactions within three-dimensional regions located near or within the rotor disk. Also, the development of nearfield velocity corrections appropriate for the analysis of such interactions is described, as is a preliminary study of methods for using the new high-resolution flow field analysis for noise predictions. The results show that by employing this novel flow field reconstruction technique it is possible to employ full-span free wake analyses with temporal and spatial resolution suitable for acoustic applications while reducing the computation time required by one to two orders of magnitude relative to traditional methods.

Quackenbush, Todd R.↗