Multi-material continuum topology optimization with arbitrary volume and mass constraints
Not Available
Sanders, Emily D.↗
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.
A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method
Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.
42 ENGINEERING↗
New Variational Formulations of Hybrid Stress Elements
In the variational formulations of finite elements by the Hu-Washizu and Hellinger-Reissner principles the stress equilibrium condition is maintained by the inclusion of internal displacements which function as the Lagrange multipliers for the constraints. These versions permit the use of natural coordinates and the relaxation of the equilibrium conditions and render considerable improvements in the assumed stress hybrid elements. These include the derivation of invariant hybrid elements which possess the ideal qualities such as minimum sensitivity to geometric distortions, minimum number of independent stress parameters, rank sufficient, and ability to represent constant strain states and bending moments. Another application is the formulation of semiLoof thin shell elements which can yield excellent results for many severe test cases because the rigid body nodes, the momentless membrane strains, and the inextensional bending modes are all represented.
Pian, T. H. H.↗
Theoretical transfer functions for the Saturn 5 Final report
Laplace transformations and Euler-Lagrange equations for theoretical transfer functions for Saturn 5 launch vehicle
Soifer, M. T.↗
The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation
We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.
42 ENGINEERING↗
On variation of action integral in Finsler gravity
In this paper, a generalized action integral of both gravity and matter is defined on the sphere bundle over Finsler space–time manifold M with a Lorentz–Finsler metric. The Euler–Lagrange equation of this functional, a generalization of the Riemann–Einstein gravity equation is obtained by using some divergence theorems. Fibres of the sphere bundle are unbounded according to the pseudo-Finsler metric. Moreover, solutions of vacuum Finsler gravity equation under the weakly Landsberg condition are discussed and some concrete examples are provided. At last, we raise some questions for further study.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Users manual: Dynamics of two bodies connected by an elastic tether, six degrees of freedom forebody and five degrees of freedom decelerator
The equations of motion and a computer program for the dynamics of a six degree of freedom body joined to a five degree of freedom body by a quasilinear elastic tether are presented. The forebody is assumed to be a completely general rigid body with six degrees of freedom; the decelerator is also assumed to be rigid, but with only five degrees of freedom (symmetric about its longitudinal axis). The tether is represented by a spring and dashpot in parallel, where the spring constant is a function of tether elongation. Lagrange's equation is used to derive the equations of motion with the Lagrange multiplier technique used to express the constraint provided by the tether. A computer program is included which provides a time history of the dynamics of both bodies and the tension in the tether.
Doyle, G. R., Jr.↗
Modeling Principles Using the Relation Between Lagrange Multipliers and Bond Graphs
Final document is attached. Modeling dynamic systems by bond graphs has become state of the art technology since hundreds of researchers around the world have incorporated the technology in many fields of engineering and science. The legacy of its invertor Prof. Henry Paynter at MIT in 1959 is now a fundamental and practical technique to understand reality by building computer models. This paper addresses a particular aspect of this technology when modeling of mechanical systems require relaxation of constraints by means of Lagrange principles. Lagrange's equations are a useful means of describing and solving systems with kinematic constraints. Lagrange multipliers are variables used in equations to find the extremes of multivariate functions. Here we explore the relation of Lagrange multipliers to solve modeling difficulties of a space vehicle with equations with dependent derivatives. Lagrange multipliers were used in conjunction with bond graphs to simulate a system where joints of kinematic linkages produce dependent derivatives. NASA's Morpheus Project lunar lander was used as a case study. The Morpheus Project is a terrestrial test vehicle designed to fly the terminal descent trajectory of a lunar lander to advance the Autonomous Landing Hazard Avoidance Technology (ALHAT). An objective of this study is to apply the modeling approach herein to capture the dynamic movement of the lander as the propellant is sloshed and consumed. This paper expands further the analysis presented by (Granda, J J. Nguyen, L, Carlson, T, Sahragard-Monfared, G., Fornalski, E., Brocker 2016). Using an automated approach bond graph models of state space equations were generated using the Computer Aided Modeling Program (CAMPG). Integral causality models and derivative causality models were considered in order to find the simpler solution for the mathematical dependencies produced in modeling this vehicle.
Granda, Jose J.↗
On stochastic extremum problems - calculus.
Lagrange multiplier technique for determining stationary points /of functions/ or stationary functions /of integrals/ of expected value of random functions with certain random constraints
LAGRANGE MULTIPLIER↗
Aerodynamic design and optimization in one shot
This paper describes an efficient numerical approach for the design and optimization of aerodynamic bodies. As in classical optimal control methods, the present approach introduces a cost function and a costate variable (Lagrange multiplier) in order to achieve a minimum. High efficiency is achieved by using a multigrid technique to solve for all the unknowns simultaneously, but restricting work on a design variable only to grids on which their changes produce nonsmooth perturbations. Thus, the effort required to evaluate design variables that have nonlocal effects on the solution is confined to the coarse grids. However, if a variable has a nonsmooth local effect on the solution in some neighborhood, it is relaxed in that neighborhood on finer grids. The cost of solving the optimal control problem is shown to be approximately two to three times the cost of the equivalent analysis problem. Examples are presented to illustrate the application of the method to aerodynamic design and constraint optimization.
Ta'asan, Shlomo↗
Shape optimization governed by the Euler equations using an adjoint method
A numerical approach for the treatment of optimal shape problems governed by the Euler equations is discussed. Focus is on flows with embedded shocks. A very simple problem is considered: the design of a quasi-one-dimensional Laval nozzle. A cost function and a set of Lagrange multipliers are introduced to achieve the minimum. The nature of the resulting costate equations is discussed. A theoretical difficulty that arises for cases with embedded shocks is pointed out and solved. Finally, some results are given to illustrate the effectiveness of the method.
Iollo, Angelo↗
Optimal Control Strategy With Efficiency and Reliability Improvement for Offshore DC Microgrids
Offshore microgrids, due to their remote location and lack of external energy support, face significant challenges in wide-range load operation and maintenance. Consequently, efficiency and reliability are critical concerns for converters in offshore dc microgrids. This article presents an optimal control strategy aimed at enhancing both efficiency and reliability. A normalized nonlinear relationship between power loss and thermal stress of a paralleled converter is first established. Based on this, a dual-objective optimization function with an active weight function as well as a system overall performance index is established. The active weight function dynamically adjusts the control priority based on converter efficiency and switching device thermal stress. Then, the optimal power-sharing strategy is derived by the Lagrange multiplier method with the proposed optimal function. Additionally, to accommodate a wide load range, an optimal selection strategy for operating converter combinations is proposed, requiring only low-bandwidth communication. Experiment verification is given to validate the effectiveness of the proposed control strategy. The experiment results demonstrate that the proposed control strategy can improve the overall performance of offshore microgrids by optimizing efficiency and reliability.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Nonlinear programming extensions to rational function approximations of unsteady aerodynamics
This paper deals with approximating unsteady generalized aerodynamic forces in the equations of motion of a flexible aircraft. Two methods of formulating these approximations are extended to include both the same flexibility in constraining them and the same methodology in optimizing nonlinear parameters as another currently used 'extended least-squares' method. Optimal selection of 'nonlinear' parameters is made in each of the three methods by use of the same nonlinear (nongradient) optimizer. The objective of the nonlinear optimization is to obtain rational approximations to the unsteady aerodynamics whose state-space realization is of lower order than that required when no optimization of the nonlinear terms is performed. The free 'linear' parameters are determined using least-squares matrix techniques on a Lagrange multiplier formulation of an objective function which incorporates selected linear equality constraints. State-space mathematical models resulting from the different approaches are described, and results are presented which show comparative evaluations from application of each of the extended methods to a numerical example. The results obtained for the example problem show a significant (up to 63 percent) reduction in the number of differential equations used to represent the unsteady aerodynamic forces in linear time-invariant equations of motion as compared to a conventional method in which nonlinear terms are not optimized.
Tiffany, Sherwood H.↗
Nonlinear programming extensions to rational function approximation methods for unsteady aerodynamic forces
The approximation of unsteady generalized aerodynamic forces in the equations of motion of a flexible aircraft are discussed. Two methods of formulating these approximations are extended to include the same flexibility in constraining the approximations and the same methodology in optimizing nonlinear parameters as another currently used extended least-squares method. Optimal selection of nonlinear parameters is made in each of the three methods by use of the same nonlinear, nongradient optimizer. The objective of the nonlinear optimization is to obtain rational approximations to the unsteady aerodynamics whose state-space realization is lower order than that required when no optimization of the nonlinear terms is performed. The free linear parameters are determined using the least-squares matrix techniques of a Lagrange multiplier formulation of an objective function which incorporates selected linear equality constraints. State-space mathematical models resulting from different approaches are described and results are presented that show comparative evaluations from application of each of the extended methods to a numerical example.
Tiffany, Sherwood H.↗
Healthy degenerate theories with higher derivatives
In the context of classical mechanics, we study the conditions under which higher-order derivative theories can evade the so-called Ostrogradsky instability. More precisely, we consider general Lagrangians with second order time derivatives, of the form L(ϕ{sup ¨a}, ϕ-dot {sup a},ϕ{sup a}; q-dot {sup i},q{sup i}) with a=1,⋯,n and i=1,⋯,m. For n=1, assuming that the q{sup i}’s form a nondegenerate subsystem, we confirm that the degeneracy of the kinetic matrix eliminates the Ostrogradsky instability. The degeneracy implies, in the Hamiltonian formulation of the theory, the existence of a primary constraint, which generates a secondary constraint, thus eliminating the Ostrogradsky ghost. For n>1, we show that, in addition to the degeneracy of the kinetic matrix, one needs to impose extra conditions to ensure the presence of a sufficient number of secondary constraints that can eliminate all the Ostrogradsky ghosts. When these conditions that ensure the disappearance of the Ostrogradsky instability are satisfied, we show that the Euler-Lagrange equations, which involve a priori higher order derivatives, can be reduced to a second order system.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Crystal-field calculations for transition-metal ions by application of an opposing potential
We propose a fully ab initio method, the opposing crystal potential (OCP), to calculate the crystal-field parameters of transition-metal impurities in insulator hosts. Through constrained density functional calculations, OCP obtains the constraining Lagrange multipliers, which act as a cancellation potential against the crystal field and lead to spherical d-electron distribution. Furthermore, the method is applied to several insulators doped with Mn 4+ and Mn 2+ ions and shown to be in good agreement with experiment.
36 MATERIALS SCIENCE↗
Optimal shielding design for minimum materials cost or mass
The mathematical underpinnings of cost optimal radiation shielding designs based on an extension of optimal control theory are presented, a heuristic algorithm to iteratively solve the resulting optimal design equations is suggested, and computational results for a simple test case are discussed. A typical radiation shielding design problem can have infinitely many solutions, all satisfying the problem's specified set of radiation attenuation requirements. Each such design has its own total materials cost. For a design to be optimal, no admissible change in its deployment of shielding materials can result in a lower cost. This applies in particular to very small changes, which can be restated using the calculus of variations as the Euler-Lagrange equations. Furthermore, the associated Hamiltonian function and application of Pontryagin's theorem lead to conditions for a shield to be optimal.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Comment of fitting ab initio intermolecular potentials for scattering calculations
The angular dependence of the interaction between two rigid linear rotors is considered, and one of the two expansions in current use is shown to be inappropriate. Interpolation of tabulated radial interactions is also discussed, and numerical tests are presented which demonstrate the accuracy of Lagrange and spline algorithms as a function of the spacing between points.
Green, S.↗