Search NASA⌕ Search

SEARCH · Search NASA

Results for “lagrangian functions”

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 253 records · Page 14

Computational Methods for Frictional Contact With Applications to the Space Shuttle Orbiter Nose-Gear Tire

A computational procedure is presented for the solution of frictional contact problems for aircraft tires. A Space Shuttle nose-gear tire is modeled using a two-dimensional laminated anisotropic shell theory which includes the effects of variations in material and geometric parameters, transverse-shear deformation, and geometric nonlinearities. Contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with both contact and friction conditions. The contact-friction algorithm is based on a modified Coulomb friction law. A modified two-field, mixed-variational principle is used to obtain elemental arrays. This modification consists of augmenting the functional of that principle by two terms: the Lagrange multiplier vector associated with normal and tangential node contact-load intensities and a regularization term that is quadratic in the Lagrange multiplier vector. These capabilities and computational features are incorporated into an in-house computer code. Experimental measurements were taken to define the response of the Space Shuttle nose-gear tire to inflation-pressure loads and to inflation-pressure loads combined with normal static loads against a rigid flat plate. These experimental results describe the meridional growth of the tire cross section caused by inflation loading, the static load-deflection characteristics of the tire, the geometry of the tire footprint under static loading conditions, and the normal and tangential load-intensity distributions in the tire footprint for the various static vertical loading conditions. Numerical results were obtained for the Space Shuttle nose-gear tire subjected to inflation pressure loads and combined inflation pressure and contact loads against a rigid flat plate. The experimental measurements and the numerical results are compared.

Tanner, John A.↗

Computational methods for frictional contact with applications to the Space Shuttle orbiter nose-gear tire: Comparisons of experimental measurements and analytical predictions

A computational procedure is presented for the solution of frictional contact problems for aircraft tires. A Space Shuttle nose-gear tire is modeled using a two-dimensional laminated anisotropic shell theory which includes the effects of variations in material and geometric parameters, transverse-shear deformation, and geometric nonlinearities. Contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with both contact and friction conditions. The contact-friction algorithm is based on a modified Coulomb friction law. A modified two-field, mixed-variational principle is used to obtain elemental arrays. This modification consists of augmenting the functional of that principle by two terms: the Lagrange multiplier vector associated with normal and tangential node contact-load intensities and a regularization term that is quadratic in the Lagrange multiplier vector. These capabilities and computational features are incorporated into an in-house computer code. Experimental measurements were taken to define the response of the Space Shuttle nose-gear tire to inflation-pressure loads and to inflation-pressure loads combined with normal static loads against a rigid flat plate. These experimental results describe the meridional growth of the tire cross section caused by inflation loading, the static load-deflection characteristics of the tire, the geometry of the tire footprint under static loading conditions, and the normal and tangential load-intensity distributions in the tire footprint for the various static vertical-loading conditions. Numerical results were obtained for the Space Shuttle nose-gear tire subjected to inflation pressure loads and combined inflation pressure and contact loads against a rigid flat plate. The experimental measurements and the numerical results are compared.

Tanner, John A.↗

LSPRAY-V: A Lagrangian Spray Module

LSPRAY-V is a Lagrangian spray solver developed for application with unstructured grids and massively parallel computers. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray encountered over a wide range of operating conditions in modern aircraft engine development. It could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. 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-V, we have advanced the state-of-the-art in spray computations in several important ways.

CFD↗

Background-field Method and QCD Factorization

One method for deriving a factorization for QCD processes is to use successive integration over fields in the functional integral. In this approach, we separate the fields into two categories: dynamical fields with momenta above a relevant cutoff, and background fields with momenta below the cutoff. The dynamical fields are then integrated out in the background of the low-momentum background fields. This strategy works well at tree level, allowing us to quickly derive QCD factorization formulas at leading order. However, to extend the approach to higher loops, it is necessary to rigorously define the functional integral over dynamical fields in an arbitrary background field. This framework was carefully developed for the calculation of the effective action in a background field at the two-loop level in the classic paper by Abbott «The Background Field Method Beyond One Loop», Nucl. Phys. B 185 , 189 (1981). Building on this work, I specify the renormalized background-field Lagrangian and define the notion of the quantum average of an operator in a background field, consistent with the “separation of scales” scheme mentioned earlier. As examples, I discuss the evolution of the twist-2 gluon light-ray operator and the one-loop gluon propagator in a background field near the light cone.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mixed models and reduced/selective integration displacement models for nonlinear analysis of curved beams

Simple mixed models are developed for use in the geometrically nonlinear analysis of deep arches. A total Lagrangian description of the arch deformation is used, the analytical formulation being based on a form of the nonlinear deep arch theory with the effects of transverse shear deformation included. The fundamental unknowns comprise the six internal forces and generalized displacements of the arch, and the element characteristic arrays are obtained by using Hellinger-Reissner mixed variational principle. The polynomial interpolation functions employed in approximating the forces are one degree lower than those used in approximating the displacements, and the forces are discontinuous at the interelement boundaries. Attention is given to the equivalence between the mixed models developed herein and displacement models based on reduced integration of both the transverse shear and extensional energy terms. The advantages of mixed models over equivalent displacement models are summarized. Numerical results are presented to demonstrate the high accuracy and effectiveness of the mixed models developed and to permit a comparison of their performance with that of other mixed models reported in the literature.

Noor, A. K.↗

Generalized Roche potential for misaligned binary systems - Properties of the critical lobe

The paper considers the Roche potential for binary systems where the stellar rotation axis is not aligned with the orbital revolution axis. It is shown that, as the degree of misalignment varies, internal Lagrangian points and external Lagrangian points may switch their roles. A systematic method to identify the internal Lagrangian point and to calculate the volume of the critical lobe is developed, and numerical results for a wide range of parameters of binary systems with circular orbits are presented. For binary systems with large enough misalignment, discrete changes occur in the topological structure of the equipotential surfaces as the orbital phase varies. The volume of the critical lobe has minima, as a function of orbital phase, at the two instances when the secondary crosses the equatorial plane of the primary. In semidetached systems, mass transfer may be confined to the vicinity of these two instances.

Avni, Y.↗

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

Breakup dynamics in a pressure-swirl injector for urea-water solution applications: A computational study

The co-optimization of in-cylinder combustion and after-treatment technology has become a major aspect in engine design and development, with the goal of meeting the increasingly restrictive emission regulations in the transportation industry. Selective Catalytic Reduction is a robust technology to control the emission of NO x , and the injection of urea in water solution is the exhaust tailpipe is a key aspect of its operation. The proposed work uses high-fidelity Computational Fluid Dynamics to characterize the atomization dynamics of the liquid jet in relevant cross-flow conditions. The study focuses on a commercial low-pressure (9 bar) pressure-swirl injector which is characterized in its internal geometry through high-resolution X-ray micro-computational tomography. The internal two-phase flow has been modeled according to the volume-of-fluid approach in a large eddy simulation framework and validated against near-nozzle X-ray radiography measurement. Moreover, characterizing the breakup dynamics for the swirling hollow cone formation, and assessing the influence of the cross-flow in the breakup dynamics was completed. The results have been reported proposing Re-Oh maps and probability density functions of the spray kinematics. Higher cross-flow momentum generates an increase in the jet intact length and a reduction of the liquid droplet diameters. The axial momentum of the jet is affected by the cross-flow already in the near-nozzle region, determining a relevant deviation of the spray velocities. In conclusion, this work aims to inform the initialization of Eulerian-Lagrangian spray models through the assignment of droplet kinematics and static one-way coupling between volume-of-fluid results and Lagrangian spray parcels, to be used for system-size domain simulations.

33 ADVANCED PROPULSION SYSTEMS↗

Theory of local thermal instability in spherical systems

The gasdynamical properties of local thermal instability in optically thin astrophysical plasmas as it occurs in spherical accretion and winds is investigated. In a medium characterized by both thermal and hydrostatic equilibrium, if the cooling function is not an explicit function of position and does not display isentropic thermal instability, then isobaric thermal instability by the Field criterion is present if and only if convective instability is present by the Schwarzschild criterion. In this case, thermal overstability cannot occur. Convective instability by the Schwarzschild criterion will also occur in accretion flows locally dominated by external heating or in marginally unbound, radiatively cooling outflows. A very general Lagrangian equation for the development of nonradial thermal instability in flows with spherical symmetry is derived and is solved analytically in certain regimes. The results are applied to cluster X-ray cooling flows.

Balbus, Steven A.↗

Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models

Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.

54 ENVIRONMENTAL SCIENCES↗

Quadrature Moments Method for the Simulation of Turbulent Reactive Flows

A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.

Raman, Venkatramanan↗

Chiral limit of 2d QCD revisited with lightcone conformal truncation

We study the chiral limit of 2d QCD with a single quark flavor at finite Nc using LCT. By modifying the LCT basis according to the quark mass in a manner motivated by ’t Hooft’s analysis, we are able to restore convergence for quark masses much smaller than the QCD strong coupling scale. For such small quark masses, the IR of the theory is expected to be well described by the Sine-Gordon model. We verify that LCT numerics are able to capture in detail the spectrum and correlation functions of the Sine-Gordon model. This opens up the possibility for studying deformations of various integrable CFTs using LCT by considering the chiral limit of QCD like theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Survey and development of finite elements for nonlinear structural analysis. Volume 2: Nonlinear shell finite elements

The development of two new shell finite elements for applications to large deflection problems is considered. The elements in question are doubly curved and of triangular and quadrilateral planform. They are restricted to small strains of elastic materials, and can accommodate large rotations. The elements described, which are based on relatively simple linear elements, make use of a new displacement function approach specifically designed for strongly nonlinear problems. The displacement function development for nonlinear applications is based on certain beam element formulations, and the strain-displacement equations are of a shallow shell type. Additional terms were included in these equations in an attempt to avoid the large errors characteristic of shallow shell elements in certain types of problems. An incremental nonlinear solution procedure specifically adopted to the element formulation was developed. The solution procedure is of combined incremental and total Lagrangian type, and uses a new updating scheme. A computer program was written to evaluate the developed formulations. This program can accommodate small element groups in arbitrary arrangements. Two simple programs were successfully solved. The results indicate that this new type of element has definite promise and should be a fruitful area for further research.

Source record↗

Evaluating the gradients of localized diabatic state energies and couplings at minimum cost

We calculate the gradients of Boys diabatic state energies and couplings when the electronic vector space is generated by configuration interaction singles. Our approach follows the Lagrangian approach of Paz and Glover (rather than direct differentiation of the adiabatic-to-diabatic approaches that have been published previously). The result is that we achieve a dramatic increase in savings over previous approaches, and the present approach should be immediately useful to scientists focused on electronic relaxation, especially chemists studying electron transfer who wish to go beyond the Condon approximation. Here, a future extension to time-dependent density functional theory in the Tamm–Dancoff approximation is clear.

Chemical physics↗

Nonlinear program based optimization of boost and buck-boost converter designs

The facility of an Augmented Lagrangian (ALAG) multiplier based nonlinear programming technique is demonstrated for minimum-weight design optimizations of boost and buck-boost power converters. Certain important features of ALAG are presented in the framework of a comprehensive design example for buck-boost power converter design optimization. The study provides refreshing design insight of power converters and presents such information as weight and loss profiles of various semiconductor components and magnetics as a function of the switching frequency.

Rahman, S.↗

PHILLS Vicarious Calibration: A WFS Experiment and Guidance for LEO-15 Calibration Efforts

On April 24, 2001, Florida Shelf Lagrangian Experiment (FSLE) cruises were accompanied by PHILLS hyperspectral overflights at 30,000 ft altitude. Robert Steward, and Jim Ivey collected remote sensing reflectance and water samples for absorption during two FSLE experiments. The PHILLS was flown on several transects over the area with vicarious calibration measurements conducted by USF from the R/V LINK. The USF slow-drop package was also deployed to collect inherent and apparent optical properties as a function of depth for an evaluation of effects of vertical structure on remote sensing spectra. The PHILLS imagery has been successfully calibrated. We show the procedures, results, derived bathymetry, and some noise effects with a 3 s period.

Carder, Kendall L.↗

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗

Hadamard products and BPS networks

We study examples of fourth-order Picard-Fuchs operators that are Hadamard products of two second-order Picard-Fuchs operators. Each second-order Picard-Fuchs operator is associated with a family of elliptic curves, and the Hadamard product computes period integrals on the fibred product of the two elliptic surfaces. We construct 3-cycles on this geometry as the union of 2-cycles in the fibre over contours on the base. We then use the special Lagrangian condition to constrain the contours on the base. This leads to a construction that is reminiscent of spectral networks and exponential networks that have previously appeared in string theory literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗