Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrange 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 91 records · Page 5

Contribution to the optimal shape design of two-dimensional internal flows with embedded shocks

We explore the practicability of optimal shape design for flows modeled by the Euler equations. We define a functional whose minimum represents the optimality condition. The gradient of the functional with respect to the geometry is calculated with the Lagrange multipliers, which are determined by solving a co-state equation. The optimization problem is then examined by comparing the performance of several gradient-based optimization algorithms. In this formulation, the flow field can be computed to an arbitrary order of accuracy. Finally, some results for internal flows with embedded shocks are presented, including a case for which the solution to the inverse problem does not belong to the design space.

Iollo, Angelo↗

On a Lagrange-Hamilton formalism describing position and momentum uncertainties

According to Heisenberg's uncertainty relation, in quantum mechanics it is not possible to determine, simultaneously, exact values for the position and the momentum of a material system. Calculating the mean value of the Hamiltonian operator with the aid of exact analytic Gaussian wave packet solutions, these uncertainties cause an energy contribution additional to the classical energy of the system. For the harmonic oscillator, e.g., this nonclassical energy represents the ground state energy. It will be shown that this additional energy contribution can be considered as a Hamiltonian function, if it is written in appropriate variables. With the help of the usual Lagrange-Hamilton formalism known from classical particle mechanics, but now considering this new Hamiltonian function, it is possible to obtain the equations of motion for position and momentum uncertainties.

Schuch, Dieter↗

New discretization of Caputo–Fabrizio derivative

We derive a numerical approximation, namely L1–2 formula, to the Caputo–Fabrizio derivative by using a quadratic interpolation. Quadratic and cubic convergence rates are achieved for L1 and L1–2 formulas using Lagrange interpolation, respectively. We compute Caputo–Fabrizio derivatives of some known functions both theoretically and numerically. In addition, we solve non/linear sub-diffusion equations to test theoretical findings. Numerical results confirm the theoretically observed convergence rates.

97 MATHEMATICS AND COMPUTING↗

Electrodynamic Shaker Capability Estimation Through Experimental Dynamic Substructuring [Thesis]

Electrodynamic shaker systems are an essential tool in shock and vibration testing of dynamic environments. However, the specific performance capability of these systems is difficult to characterize. The dynamics of the shaker itself, the device under test and the specific test configuration used all couple to create a dynamic response unique to each test. Poorly predicted limitations in shaker capability affect the ability to achieve test specifications, delay testing schedules, and create difficulties for choosing test equipment. To predict shaker capability for a specific test configuration prior to setup, a lumped parameter model of the shaker system and a modal model of a device under test was developed. These models were then analytically coupled using LaGrange multiplier frequency based substructuring to estimate their coupled frequency response functions. The coupled frequency response functions were used to predict electrical inputs required to meet a given test specification. These input requirements were finally compared to a validation test using the specification and setup. Input requirements estimated using the substructuring estimated frequency response functions showed significant error. However, results using an ideal frequency response function showed very little error. These results indicate that with a better method of experimental dynamic substructuring employed it would be possible to accurately predict shaker capability for a given test configuration prior to setup.

42 ENGINEERING↗

Analytic theory of orbit contraction

The motion of a satellite in orbit, subject to atmospheric force and the motion of a reentry vehicle are governed by gravitational and aerodynamic forces. This suggests the derivation of a uniform set of equations applicable to both cases. For the case of satellite motion, by a proper transformation and by the method of averaging, a technique appropriate for long duration flight, the classical nonlinear differential equation describing the contraction of the major axis is derived. A rigorous analytic solution is used to integrate this equation with a high degree of accuracy, using Poincare's method of small parameters and Lagrange's expansion to explicitly express the major axis as a function of the eccentricity. The solution is uniformly valid for moderate and small eccentricities. For highly eccentric orbits, the asymptotic equation is derived directly from the general equation. Numerical solutions were generated to display the accuracy of the analytic theory.

Vinh, N. X.↗

On solving the irregularly-shaped plate buckling problem

This paper develops methods for finding critical buckling loads of isotropic plates. In particular, this paper is motivated by the lack of tabular data for non-rectangular plates such as general trapezoids. High order displacement functions are used, and boundary conditions are imposed by utilizing the Lagrange multiplier method, which allows for a variety of constraint conditions. The element is derived with the use of minimum potential energy and subparametric mappings. An example is given, and several test cases are presented to show the validity of the approach in the one element case.

Smith, James P.↗

Optimization-based algorithms for nonlinear mechanics and frictional contact

An optimization-based strategy for solving nonlinear mechanics problems is proposed. In contrast to typical nonlinear equation solver algorithms that aim to find zeros in the residual force function, we minimize an energy (or energy-like) function to encourage solutions which are locally stable equilibria. These smooth and potentially non-convex objective functions are minimized using a preconditioned conjugate-gradient trust-region algorithm. Contact is formulated as an inequality constrained minimization problem, and is solved with an augmented Lagrangian algorithm. Friction is included in the approach via a regularized quasi-potential energy, and other dissipative behavior is included through the use of variational constitutive updates. Finally, to accelerate convergence rates for the Lagrange multipliers, we propose a novel multiplier update algorithm utilizing the Fischer-Burmeister function, and demonstrate super-linear solver convergence for some applications.

42 ENGINEERING↗

Use of a minimum rate of change formalism to quantify variability of extragalactic X-ray sources

A method to obtain rigorous, quantitative constraints on the rate of variability is suggested. The method is motivated by the case when the existence of time variability is unequivocal, but the statistical uncertainties are too large to apply more direct methods such as power spectrum analysis. Conceptually the data are fitted to a function of time and a set of free parameters. The method of Lagrange multipliers is used to solve for that parameter set which minimizes the value of the derivative at time t(m). The result can be related physically to the minimum efficiency with which rest mass must be converted to radiation energy. Slightly different formulations of the general principal are used to estimate maximum source size scales associated with the variability.

Schwartz, Daniel A.↗

Application of Newton modified barrier method (NMBM) to structural optimization

This paper presents the application of the NMBM to obtain a minimum weight structure with constraints on displacements and minimum sizes. The solution to the problem is obtained via minimizing the Modified Barrier Function (MBF) at each step by using the Newton Method and updating Lagrange multipliers. The Lagrange multipliers are updated by using the value of the constraints at the minimum of the MBF. Three truss problems with a different number of design variables are solved. The convergence to the minimum weight design was found to be monotonic and the algorithm is potentially robust for solving problems with a large number of design variables.

Khot, N. S.↗

Further Improvement in 3DGRAPE

3DGRAPE/AL:V2 denotes version 2 of the Three-Dimensional Grids About Anything by Poisson's Equation with Upgrades from Ames and Langley computer program. The preceding version, 3DGRAPE/AL, was described in Improved 3DGRAPE (ARC-14069) NASA Tech Briefs, Vol. 21, No. 5 (May 1997), page 66. These programs are so named because they generate volume grids by iteratively solving Poisson's Equation in three dimensions. The grids generated by the various versions of 3DGRAPE have been used in computational fluid dynamics (CFD). The main novel feature of 3DGRAPE/AL:V2 is the incorporation of an optional scheme in which anisotropic Lagrange-based trans-finite interpolation (ALBTFI) is coupled with exponential decay functions to compute and blend interior source terms. In the input to 3DGRAPE/AL:V2 the user can specify whether or not to invoke ALBTFI in combination with exponential-decay controls, angles, and cell size for controlling the character of grid lines. Of the known programs that solve elliptic partial differential equations for generating grids, 3DGRAPE/AL:V2 is the only code that offers a combination of speed and versatility with most options for controlling the densities and other characteristics of grids for CFD.

Alter, Stephen↗

Effects of Proton Irradiation on the Performance of Skipper CCDs

Skipper CCDs are a mature detector technology that has been suggested for future space telescope instruments requiring sub-electron readout noise in the near-ultraviolet to the near-infrared. While modern skipper CCDs inherit from the radiation-tolerant p-channel detectors developed by LBNL, the effects of high doses of ionizing radiation on skipper CCDs (such as those expected in space) remains largely unmeasured. We report preliminary results on the performance of p-channel skipper CCDs following irradiation with 217-MeV protons at the Northwestern Medicine Proton Center. The total nonionizing energy loss (NIEL) experienced by the detectors exceeds 6 years at the Sun-Earth Lagrange Point 2 (L2). We demonstrate that the skipper amplifier continues to function as expected following this irradiation. Owing to the low readout noise of these detectors, controlled irradiation tests can be used to sensitively characterize the charge transfer inefficiency, dark current, and the density and time constants of charge traps as a function of proton fluence. We conclude with a brief outlook toward future tests of these detectors at other proton and gamma-ray facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Including the vacuum energy in stellarator coil design

Being three-dimensional, stellarators have the advantage that plasma currents are not essential for creating rotational-transform; however, the external current-carrying coils in stellarators can have strong geometrical shaping, which can complicate the construction. Reducing the inter-coil electromagnetic forces acting on strongly shaped 3D coils and the stress on the support structure while preserving the favorable properties of the magnetic field is a design challenge. In this work, we recognize that the inter-coil forces are the gradient of the vacuum magnetic energy. We introduce an objective functional built on the usual quadratic flux on a prescribed target surface together with a weighed penalty on the vacuum energy. The Euler–Lagrange equation for stationary states is derived, and numerical illustrations are computed using a modern stellarator optimization framework. A study of the effect of the energy functional on the inter-coil forces is conducted and the energy is shown to be a promising quantity in producing coils with low forces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gaussian variational equations for osculating elements of an arbitrary separable reference orbit.

Lagrange-type equations are often used in planetary theory and sometimes in satellite theory. The equations express the variation of osculating Keplerian elements in terms of derivatives of a disturbing function or a perturbing potential. When the perturbing force is derivable from a potential, it is possible to convert the Lagrange-type equations to another form. This form, usually, attributed to Gauss, contains the perturbing forces instead of the derivatives of the potential. The case of an arbitrary separable reference orbit is discussed together with a lemma, a Keplerian check, and questions of the applicability of the equations to the spheroidal method.

Vinti, J. P.↗

On restricting first order form of gauge theories to one-loop order

Highlights: • We prove the equivalence of the first and second-order formalism. • We use Lagrange multiplier fields to limit radiative effects to one-loop order. • We prove that the quantized theory is both renormalizable and unitary. The first order form of the Yang–Mills and Einstein–Hilbert actions are quantized, and it is shown how Green’s functions computed using the first and the second order form of these theories are related. Next we show how by use of Lagrange multiplier fields (LM), radiative effects beyond one-loop order can be eliminated. This allows one to compute Green’s functions exactly without loss of unitarity. The consequences of this restriction on radiative effects are examined for the Yang–Mills and Einstein–Hilbert actions. In these two gauge theories, we find that the quantized theory is both renormalizable and unitary once the LM field is used to eliminate effects beyond one-loop order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A new computational approach for the solutions of generalized pantograph-delay differential equations

In this study, a new computational approach is presented to solve the generalized pantograph-delay differential equations (PDDEs). The solutions obtained by our scheme represent by a linear combination of a special kind of basis functions, and can be deduced in a straightforward manner. Firstly, using the least squares approximation method and the Lagrange-multiplier method, the given PDDE is converted to a linear system of algebraic equations, and those unknown coefficients of the solution of the problem are determined by solving this linear system. Secondly, a PDDE related to the error function of the approximate solution is constructed based on the residual error function technique, and error estimation is presented for the suggested method. The convergence of the approximate solution is proved. Several numerical examples are given to demonstrate the accuracy and efficiency. Comparisons are made between the proposed method and other existing methods.

97 MATHEMATICS AND COMPUTING↗

Finite element analysis of hyperelastic structures

The use of the penalty function, to account for near incompressibility is discussed and compared to that of Lagrange multiplier. A scheme to use Lagrange multiplier, without having to treat it as unknown, is also presented.

Tabaddor, F.↗

A survey of parametrized variational principles and applications to computational mechanics

This survey paper describes recent developments in the area of parametrized variational principles (PVP's) and selected applications to finite-element computational mechanics. A PVP is a variational principle containing free parameters that have no effect on the Euler-Lagrange equations. The theory of single-field PVP's based on gauge functions (also known as null Lagrangians) is a subset of the inverse problem of variational calculus that has limited value. On the other hand, multifield PVP's are more interesting from theoretical and practical standpoints. Following a tutorial introduction, the paper describes the recent construction of multifield PVP's in several areas of elasticity and electromagnetics. It then discusses three applications to finite-element computational mechanics: the derivation of high-performance finite elements, the development of element-level error indicators, and the constructions of finite element templates. The paper concludes with an overview of open research areas.

Felippa, Carlos A.↗

Analytic error analysis of cross section interpolation methods in nodal diffusion codes - I : Theory

This paper discusses two cross section interpolation methods commonly found in popular nodal codes; the partial derivatives and multiple tables models. The motivation for choosing a model, and thus a case matrix structure, is a trade off between accuracy and computational cost. Due to decades of experience, there are default structures that are sufficient for current light water reactor analysis. However, this is not necessarily the case for advanced reactor designs. Therefore, it is advantageous to understand the sources of error in cross section interpolation models so that the quality of a case matrix may be improved. A mathematical framework for these models is presented in this work that provides a more rigorous connection between the nuclear engineering field's cross section interpolation methods and the broader mathematical field of function approximation. The two cross section models examined in this paper were found to utilize Lagrange interpolation and are a subset of Lagrange tensor products. Classical results of Lagrange polynomial error analysis were then applied to the partial derivative and multiple tables models to derive expressions for the total point-wise error. The analytical results classify the total error into two parts: the model form error and interpolation error. Finally, based on our observations, a better foundation for improving the quality of a case matrix is proposed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗