Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LAGRANGIAN FUNCTION”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Modeling and Analysis of Power Processing Systems (MAPPS). Volume 1: Technical report

Computer aided design and analysis techniques were applied to power processing equipment. Topics covered include: (1) discrete time domain analysis of switching regulators for performance analysis; (2) design optimization of power converters using augmented Lagrangian penalty function technique; (3) investigation of current-injected multiloop controlled switching regulators; and (4) application of optimization for Navy VSTOL energy power system. The generation of the mathematical models and the development and application of computer aided design techniques to solve the different mathematical models are discussed. Recommendations are made for future work that would enhance the application of the computer aided design techniques for power processing systems.

Lee, F. C.↗
NASA NTRS1985

Three-dimensional hybrid-stress finite element analysis of composite laminates with cracks and cutouts

A three-dimensional hybrid-stress finite element analysis of composite laminates containing cutouts and cracks is presented. Fully three-dimensional, hexahedral isoparametric elements of the hybrid-stress model are formulated on the basis of the Hellinger-Reissner variational principle. Traction-free edges, cutouts, and crack surfaces are modeled by imposition of exact traction boundary conditions along element surfaces. Special boundary and surface elements are constructed by introducing proper constraints on assumed stress functions. The Lagrangian multiplier technique is used to enforce ply-interface continuity conditions in hybrid bimaterial composite elements for modeling the interface region in a composite laminate. Two examples are given to illustrate the capability of the present method of approach: (1) the well-known delamination problem in an angle-ply laminate, and (2) the important problem of a composite laminate containing a circular hole. Results are presented in detail for each case. Implications of interlaminar and intralaminar crack initiation, growth and fracture in composites containing cracks and cutouts are discussed.

Wang, S. S.↗
NASA NTRS1989

Grid generation and inviscid flow computation about aircraft geometries

Grid generation and Euler flow about fighter aircraft are described. A fighter aircraft geometry is specified by an area ruled fuselage with an internal duct, cranked delta wing or strake/wing combinations, canard and/or horizontal tail surfaces, and vertical tail surfaces. The initial step before grid generation and flow computation is the determination of a suitable grid topology. The external grid topology that has been applied is called a dual-block topology which is a patched C (exp 1) continuous multiple-block system where inner blocks cover the highly-swept part of a cranked wing or strake, rearward inner-part of the wing, and tail components. Outer-blocks cover the remainder of the fuselage, outer-part of the wing, canards and extend to the far field boundaries. The grid generation is based on transfinite interpolation with Lagrangian blending functions. This procedure has been applied to the Langley experimental fighter configuration and a modified F-18 configuration. Supersonic flow between Mach 1.3 and 2.5 and angles of attack between 0 degrees and 10 degrees have been computed with associated Euler solvers based on the finite-volume approach. When coupling geometric details such as boundary layer diverter regions, duct regions with inlets and outlets, or slots with the general external grid, imposing C (exp 1) continuity can be extremely tedious. The approach taken here is to patch blocks together at common interfaces where there is no grid continuity, but enforce conservation in the finite-volume solution. The key to this technique is how to obtain the information required for a conservative interface. The Ramshaw technique which automates the computation of proportional areas of two overlapping grids on a planar surface and is suitable for coding was used. Researchers generated internal duct grids for the Langley experimental fighter configuration independent of the external grid topology, with a conservative interface at the inlet and outlet.

Smith, Robert E.↗
NASA NTRS1991

Grid generation and compressible flow computations about a high-speed civil transport configuration

Techniques and software are discussed for generating grids about a high-speed civil transport configuration. The configuration is defined by a computer-aided design system in wing, fuselage, tail and engine-nacelle components. Grid topology and the surfaces outlining the blocks of the topology are computed with interactive software. The volume grid is computed using software based on transfinite interpolation and Lagrangian blending functions. Several volume grids for inviscid and viscous flow have been generated using this system of codes. Demonstration flowfields around this vehicle are described.

Abolhassani, J. S.↗
DOE OSTI2024

Nonlocal chiral contributions to generalized parton distributions of the proton at nonzero skewness

We compute the one-loop contributions to spin-averaged generalized parton distributions (GPDs) in the proton from pseudoscalar mesons with intermediate octet and decuplet baryon states at nonzero skewness. Our framework is based on nonlocal covariant chiral effective theory, with ultraviolet divergences regularized by introducing a relativistic regulator derived consistently from the nonlocal Lagrangian. Using the splitting functions calculated from the nonlocal Lagrangian, we find the nonzero skewness GPDs from meson loops by convoluting with the phenomenological pion GPD and the generalized distribution amplitude, and verify that these satisfy the correct polynomiality properties. We also compute the lowest two moments of GPDs to quantify the meson loop effects on the Dirac, Pauli, and gravitational form factors of the proton.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Eulerian and Lagrangian second-order statistics of superfluid 4 He grid turbulence

Here, we use particle-tracking velocimetry to study Eulerian and Lagrangian second-order statistics of superfluid 4 He grid turbulence. The Eulerian energy spectra at scales larger than the mean distance between quantum vortex lines behave classically with close to Kolmogorov-1941 scaling and are almost isotropic. The Lagrangian second-order structure functions and frequency power spectra, measured at scales comparable with the intervortex distance, demonstrate a sharp transition from nearly classical behavior to a regime dominated by the motion of quantum vortex lines. Employing the homogeneity of the flow, we verify a set of relations that connect various second-order statistical objects that stress different aspects of turbulent behavior, allowing a multifaceted analysis. We use the two-way bridge relations between Eulerian energy spectra and second-order structure functions to reconstruct the energy spectrum from the known second-order velocity structure function and vice versa. The Lagrangian frequency spectrum reconstructed from the measured Eulerian spectrum using the Eulerian-Lagrangian bridge differs from the measured Lagrangian spectrum in the quasiclassical range, which calls for further investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Comparison of Eulerian Bin and Lagrangian Particle-Based Microphysics in Simulations of Nonprecipitating Cumulus

A single nonprecipitating cumulus congestus setup is applied to compare droplet spectra grown by the diffusion of water vapor in Eulerian bin and particle-based Lagrangian microphysics schemes. Bin microphysics represent droplet spectral evolution applying the spectral density function. In the Lagrangian microphysics, computational particles referred to as superdroplets are followed in time and space with each superdroplet representing a multiplicity of natural cloud droplets. The same cloud condensation nuclei (CCN) activation and identical representation of the droplet diffusional growth allow the comparison. The piggybacking method is used with the two schemes operating in a single simulation, one scheme driving the dynamics and the other one piggybacking the simulated flow. Piggybacking allows point-by-point comparison of droplet spectra predicted by the two schemes. Overall, the results show the impact of inherent limitations of the two microphysics simulation methods, numerical diffusion in the Eulerian scheme and a limited number of superdroplets in the Lagrangian scheme. Numerical diffusion in the Eulerian scheme results in a more dilution of the cloud upper half and thus smaller cloud droplet mean radius. The Lagrangian scheme typically has larger spatial fluctuations of droplet spectral properties. A significantly larger mean spectral width in the bin microphysics across the entire cloud depth is the largest difference between the two schemes. A fourfold increase of the number of superdroplets per grid volume and a twofold increase of the spectral resolution and thus the number of bins have small impact on the results and provide only minor changes to the comparison between simulated cloud properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Comparison of Eulerian Bin and Lagrangian Particle-Based Schemes in Simulations of Pi Chamber Dynamics and Microphysics

This paper presents a comparison of simulations applying either a traditional Eulerian bin microphysics or a novel particle-based Lagrangian approach to represent CCN activation and cloud droplet growth. The Eulerian microphysics solve the evolution equation for the spectral density function, whereas the Lagrangian approach follows computational particles referred to as superdroplets. Each superdroplet represents a multiplicity of natural droplets that makes the Lagrangian approach computationally feasible. The two schemes apply identical representation of CCN activation and use the same droplet growth equation; these make direct comparison between the two schemes practical. The comparison, the first of its kind, applies an idealized simulation setup motivated by laboratory experiments with the Pi Chamber and previous model simulations of the Pi Chamber dynamics and microphysics. The Pi Chamber laboratory apparatus considers interactions between turbulence, CCN activation, and cloud droplet growth in moist Rayleigh–Bénard convection. Simulated steady-state droplet spectra averaged over the entire chamber are similar, with the mean droplet concentration, mean radius, and spectral width close in Eulerian and Lagrangian simulations. Small differences that do exist are explained by the inherent differences between the two schemes and their numerical implementation. The local droplet spectra differ substantially, again in agreement with the inherent limitations of the theoretical foundation behind each approach. There is a general agreement between simulations and Pi Chamber observations, with simplifications of the CCN activation and droplet growth equation used in the simulations likely explaining specific differences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulation-driven optimization of high-order meshes in ALE hydrodynamics

Here we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP). The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

The Newton Modified Barrier Method for QP Problems

The Modified Barrier Functions (MBF) have elements of both Classical Lagrangians (CL) and Classical Barrier Functions (CBF). The MBF methods find an unconstrained minimizer of some smooth barrier function in primal space and then update the Lagrange multipliers, while the barrier parameter either remains fixed or can be updated at each step. The numerical realization of the MBF method leads to the Newton MBF method, where the primal minimizer is found by using Newton's method. This minimizer is then used to update the Lagrange multipliers. In this paper, we examine the Newton MBF method for the Quadratic Programming (QP) problem. It will be shown that under standard second-order optimality conditions, there is a ball around the primal solution and a cut cone in the dual space such that for a set of Lagrange multipliers in this cut cone, the method converges quadratically to the primal minimizer from any point in the aforementioned ball, and continues to do so after each Lagrange multiplier update. The Lagrange multipliers remain within the cut cone and converge linearly to their optimal values. Any point in this ball will be called a "hot start". Starting at such a "hot start", at most Omicron(1n 1n epsilon(exp -1)) Newton steps are sufficient to perform the primal minimization which is necessary for the Lagrange multiplier update. Here, epsilon > 0 is the desired accuracy. Because of the linear convergence of the Lagrange multipliers, this means that only Omicron(1n epsilon(exp -1))omicron(ln 1n epsilon(exp-1)) Newton steps are required to reach an epsilon-approximation to the solution from any "hot start". In order to reach the "hot start", one has to perform Omicron(square root(m) 1n C) Newton steps, where m characterizes the size of the problem and C > 0 is the condition number of the QP problem. This condition number will be characterized explicitly in terms of key parameters of the QP problem, which in turn depend on the input data and the size of the problem.

Melman, A.↗
DOE OSTI2023

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

An adaptive sampling augmented Lagrangian method for stochastic optimization with deterministic constraints

The primary goal of this paper is to provide an efficient solution algorithm based on the augmented Lagrangian framework for optimization problems with a stochastic objective function and deterministic constraints. Our main contribution is combining the augmented Lagrangian framework with adaptive sampling, resulting in an efficient optimization methodology validated with practical examples. To achieve the presented efficiency, here we consider inexact solutions for the augmented Lagrangian subproblems, and through an adaptive sampling mechanism, we control the variance in the gradient estimates. Furthermore, we analyze the theoretical performance of the proposed scheme by showing equivalence to a gradient descent algorithm on a Moreau envelope function, and we prove sublinear convergence for convex objectives and linear convergence for strongly convex objectives with affine equality constraints. The worst-case sample complexity of the resulting algorithm, for an arbitrary choice of penalty parameter in the augmented Lagrangian function, is $\mathscr{O}$(ϵ -3-δ ) , where ϵ > 0 is the expected error of the solution and δ > 0 is a user-defined parameter. If the penalty parameter is chosen to be $\mathscr{O}$(ϵ -1 ), we demonstrate that the result can be improved to $\mathscr{O}$(ϵ -2 ) , which is competitive with the other methods employed in the literature. Moreover, if the objective function is strongly convex with affine equality constraints, we obtain $\mathscr{O}$(ϵ -1 log(1/ϵ)) complexity. Finally, we empirically verify the performance of our adaptive sampling augmented Lagrangian framework in machine learning optimization and engineering design problems, including topology optimization of a heat sink with environmental uncertainty.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Interactive grid generation for fighter aircraft geometries

An interactive procedure to construct grids about fighter aircraft configurations using a 'state-of-the-art' workstation is described. A multiple-block grid topology is used to minimize grid skewness. Grid points on the surface of a configuration are precomputed and stored on the faces of appropriate blocks. Grid points at intermediate and far field boundary surfaces are computed with analytical functions and stored on appropriate faces. Grid points in the interior of blocks are computed using transfinite interpolation with Lagrangian and exponential blending functions. Interactive software has been written for the IRIS 3030 graphics workstation to control the shape and spacing of grid points on intermediate and far field boundaries, and to control the interior interpolation.

Smith, Robert E.↗
NASA NTRS1989

Pdf - Transport equations for chemically reacting flows

The closure problem for the transport equations for pdf and the characteristic functions of turbulent, chemically reacting flows is addressed. The properties of the linear and closed equations for the characteristic functional for Eulerian and Lagrangian variables are established, and the closure problem for the finite-dimensional case is discussed for pdf and characteristic functions. It is shown that the closure for the scalar dissipation term in the pdf equation developed by Dopazo (1979) and Kollmann et al. (1982) results in a single integral, in contrast to the pdf, where double integration is required. Some recent results using pdf methods obtained for turbulent flows with combustion, including effects of chemical nonequilibrium, are discussed.

Kollmann, W.↗
DOE OSTI2023

A Fast Temporal Decomposition Procedure for Long-Horizon Nonlinear Dynamic Programming

We propose a fast temporal decomposition procedure for solving long-horizon nonlinear dynamic programs. The core of the procedure is sequential quadratic programming (SQP) that utilizes a differentiable exact augmented Lagrangian as the merit function. Within each SQP iteration, we approximately solve the Newton system using an overlapping temporal decomposition strategy. We show that the approximate search direction is still a descent direction of the augmented Lagrangian provided the overlap size and penalty parameters are suitably chosen, which allows us to establish the global convergence. Moreover, we show that a unit step size is accepted locally for the approximate search direction and further establish a uniform, local linear convergence over stages. This local convergence rate matches the rate of the recent Schwarz scheme (Na et al. 2022). However, the Schwarz scheme has to solve nonlinear subproblems to optimality in each iteration, whereas we only perform a single Newton step instead. Numerical experiments validate our theories and demonstrate the superiority of our method.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

A generalization of Petterssen's frontogenesis function and its relation to the forcing of vertical motion

The classic Peterssen frontogenesis function, defined as the Lagrangian rate of change of the magnitude of the potential temperature gradient, is generalized to apply to the vector potential temperature gradient. The derivation of vector frontogenesis in natural coordinates is shown, and analytical examples of frontal evolution in nondivergent horizontal velocity fields are given which suggest that both the frontogenetic and rotational components of the vector frontogenesis function F may be comparable in developing frontal zones observed in nature. The relative importance of the magnitude and direction contributions to F is quantitatively investigated, and it is found that the Lagrangian rates of change of the magnitude and direction of the potential temperature gradient are comparable. The frontal circulation is found to be related to the magnitude component of the Q vector, whereas the background circulation is related to the direction component.

Keyser, Daniel↗