Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrange”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Coupling equivalent plate and finite element formulations in multiple-method structural analyses

A coupled multiple-method analysis procedure for use late in conceptual design or early in preliminary design of aircraft structures is described. Using this method, aircraft wing structures are represented with equivalent plate models, and structural details such as engine/pylon structure, landing gear, or a 'stick' model of a fuselage are represented with beam finite element models. These two analysis methods are implemented in an integrated multiple-method formulation that involves the assembly and solution of a combined set of linear equations. The corresponding solution vector contains coefficients of the polynomials that describe the deflection of the wing and also the components of translations and rotations at the joints of the beam members. Two alternative approaches for coupling the methods are investigated; one using transition finite elements and the other using Lagrange multipliers. The coupled formulation is applied to the static analysis and vibration analysis of a conceptual design model of a fighter aircraft. The results from the coupled method are compared with corresponding results from an analysis in which the entire model is composed of finite elements.

VIBRATION ANALYSIS↗
NASA NTRS2004

L2 Earth atmosphere observatory : formation guidance, metrology, and control synthesis

This paper discusses the results of research sponsored by the NASA Revolutionary Aerospace Systems Concepts (RASC) program, and includes the synthesis and analysis of the guidance, metrology and control for a two-spacecraft formation in a unique continuously powered orbit near the Sun-Earth L2 Lagrange point observing the illuminated atmosphere of the Earth while it is continuously occulting the Sun.

observatory↗
NASA NTRS2004

L2 Earth Atmosphere Observatory: Formation Guidance, Metrology and Control Synthesis

The Earth Observatory Formation at L2, a Lagrange libration point, is a unique large aperture (25 m diameter) space telescope concept that will improve the knowledge and understanding of dynamic, chemical and radiative mechanisms that cause changes in the atmosphere, and can lead to the development of models and techniques to predict short and long-term climate changes. The results of this concept definition study show that the telescope concept is feasible, and can have technology readiness in the 2020 time frame. Further advanced development in several subsystems is needed, such as higher efficiency Xenon ion thrusters with throttling, and optical quality large membrane mirror with active shape control. It presents an analysis and solution of guidance, sensing, control, and propulsion problems for a formation of two spacecraft on the Sun-Earth line in the neighborhood of the Sun-Earth L2 point, that observes Earth s atmosphere during continuous solar occultation by the Earth. A system architecture is described for the observatory, and its components that include unique mission specific metrology. The formation must follow a powered trajectory with strictly limited fuel use to observe solar occultation. A configuration of ion thrusters and reaction wheels for translation and attitude control is designed along with algorithms for orbit following and formation control. Simulation results of the orbital and formation dynamics are presented that verify performance of the control systems.

Earth Atmosphere Observatory Formation↗
NASA NTRS2004

Earth Atmosphere Observatory Formation at L2

This paper is a product of research supported by NASA under RASC (the Revolutionary Aerospace Systems Concepts) program. It presents an overall system architecture, and covers issues of deployment, navigation, and control related to a formation of two spacecraft in the neighborhood of the Sun-Earth L2 Lagrange point (on the Sun-Earth line), that serves as an observatory of Earth's atmosphere. The observatory concept definition study was a multi-center NASA effort conducted in 2003, and covered a much wider scope than is presented in this focused paper.The Earth observatory at L2 is a unique design concept that can improve the knowledge and understanding of dynamic, chemical and radiative mechanisms that cause changes in the atmosphere, and can lead to the development of models and techniques to predict short and long-term climate changes.

observatory↗
NASA NTRS2012

Orbital Operations for Phobos and Deimos Exploration

One of the deep-space human exploration activities proposed for the post-Shuttle era is a mission to one of the moons of Mars, Phobos or Deimos. There are several options available to the mission architect for operations around these bodies. These options include distant retrograde orbits (DROs), Lagrange-point orbits such as halos and Lyapunov orbits, and fixed-point stationkeeping or "hovering." These three orbit options are discussed in the context of the idealized circular restricted three body problem, full-dynamics propagations, and a concept of operations. The discussion is focused on Phobos, but all results hold for Deimos

distant retrograde orbits (DROs)↗
NASA NTRS2022

On-Orbit Performance of the Spitzer Space Telescope: Science Meets Engineering

The Spitzer Space Telescope operated for over 16 years in an Earth-trailing solar orbit, returning not only a wealth of scientific data but also, as a by-product, spacecraft and instrument engineering data that will be of interest to future mission planners. These data will be particularly useful because Spitzer operated in an environment essentially identical to that at the L2 LaGrange point, where many future astrophysics missions will operate. In particular, the radiative cooling demonstrated by Spitzer has been adopted by other infrared space missions, from JWST to SPHEREx. We aim to facilitate the utility of the Spitzer engineering data by collecting the more unique and potentially useful portions into a single, readily accessible publication. We avoid discussion of less unique systems, such as the telecom, flight software, and electronics systems, and do not address the innovations in mission and science operations that the Spitzer team initiated. These and other items of potential interest are addressed in references supplied in an appendix to this paper.

Spitzer↗
DOE OSTI2020

Optimal Operation of a Hybrid Hydraulic Electric Architecture (HHEA) for Off-Road Vehicles Over Discrete Operating Decisions

Many off-highway machines including construction and agriculture equipment use hydraulics for power transmission and throttling as a means for control. Trends towards better efficiency and electrification have led to the creation of a novel Hybrid Hydraulic-Electric Architecture (HHEA) which could significantly reduce energy consumption and maintain control performance, even in machines that are too large to be directly electrified. This is achieved by using a set of common pressure rails to transmit the majority of power via power dense hydraulics and modulating the power with small electric motor-drives to achieve precise control. This paper proposes a computationally efficient method for computing the optimal sequence of pressure rail selections for the HHEA over finite drive cycles. This is useful for fairly comparing the novel architecture’s energy performance to existing architectures and for use in iterative optimal design of the architecture. The optimal control technique makes use of a static model of the architecture and losses. Constraints are enforced such that the drive cycle is repeatable. The constrained optimal operation is solved using a Lagrange multiplier technique that transforms the optimization into a small set of sub-problems by considering combinations of active constraints. Each of these sub-problems can be solved efficiently because loss calculations for all time steps can be computed in parallel. A case study of an off-road construction machine demon-strates that the HHEA reduces energy consumption by 2/3 compared to the baseline load sensing architecture.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗
DOE OSTI2023

Simultaneous CTEQ-TEA extraction of PDFs and SMEFT parameters from jet and $t\bar{t}$ data

Recasting phenomenological Lagrangians in terms of SM effective field theory (SMEFT) provides a valuable means of connecting potential BSM physics at momenta well above the electroweak scale to experimental signatures at lower energies. In this work we jointly fit the Wilson coefficients of SMEFT operators as well as the PDFs in an extension of the CT18 global analysis framework, obtaining self-consistent constraints to possible BSM physics effects. Global fits are boosted with machine-learning techniques in the form of neural networks to ensure efficient scans of the full PDF+SMEFT parameter space. We focus on several operators relevant for top-quark pair and jet production at hadron colliders and obtain constraints on the Wilson coefficients with Lagrange Multiplier scans. We find mild correlations between the extracted Wilson coefficients, PDFs, and other QCD parameters, and see indications that these correlations may become more prominent in future analyses based on data of higher precision. This work serves as a new platform for joint analyses of SM and BSM physics based on the CTEQ-TEA framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Covariant bit threads

We derive several new reformulations of the Hubeny-Rangamani-Takayanagi covariant holographic entanglement entropy formula. These include: (1) a minimax formula, which involves finding a maximal-area achronal surface on a timelike hypersurface homologous to D(A) (the boundary causal domain of the region A whose entropy we are calculating) and minimizing over the hypersurface; (2) a max V-flow formula, in which we maximize the flux through D(A) of a divergenceless bulk 1-form V subject to an upper bound on its norm that is non-local in time; and (3) a min U-flow formula, in which we minimize the flux over a bulk Cauchy slice of a divergenceless timelike 1-form U subject to a lower bound on its norm that is non-local in space. The two flow formulas define convex programs and are related to each other by Lagrange duality. For each program, the optimal configurations dynamically find the HRT surface and the entanglement wedges of A and its complement. The V-flow formula is the covariant version of the Freedman-Headrick bit thread reformulation of the Ryu-Takayanagi formula. We also introduce a measure-theoretic concept of a “thread distribution”, and explain how Riemannian flows, V-flows, and U-flows can be expressed in terms of thread distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Algebraic Structures in the Coupling of Gravity to Gauge Theories

Highlights: • Geometric foundation for Quantum General Relativity coupled to Quantum Electrodynamics. • Rigorous introduction to Hopf algebraic renormalization. • Appropriate modifications in case renormalization Hopf algebra is ill-defined. • Generalization of Furry’s Theorem to include gravitons and ghosts. • Discussion of obstructions to multiplicative renormalization. This article is an extension of the author’s second master thesis, Prinz (2017). It aims to introduce to the theory of perturbatively quantized General Relativity coupled to Spinor Electrodynamics, provide the results thereof and set the notation to serve as a starting point for further research in this direction. It includes the differential geometric and Hopf algebraic background, as well as the corresponding Lagrange density and some renormalization theory. Then, a particular problem in the renormalization of Quantum General Relativity coupled to Quantum Electrodynamics is addressed and solved by a generalization of Furry’s Theorem. Next, the restricted combinatorial Green’s functions for all two-loop propagators and all one-loop divergent subgraphs thereof are presented. Finally, relations between these one-loop restricted combinatorial Green’s functions necessary for multiplicative renormalization are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗
DOE OSTI2022

Sensitivity of puff dynamics and airborne droplet nuclei distribution to variations in violent expiration events

Here, we use large eddy simulations to investigate the puff and droplet dynamics from violent expiratory events such as coughs and sneezes in the first few seconds following an ejection. For each of the eleven simulations considered, over 60,000 droplets are ejected and individually tracked using the point-particle Euler–Lagrange approach. We test the sensitivity of the puff and droplet dynamics to various parameters including the ejection volume, momentum, and orientation. We also explore the effect of the mouth shape on the aforementioned dynamics by considering elliptical and circular inlet cross-sections. The results from the simulations compare favorably with a recent theoretical framework put forth by Balachandar et al. in terms of the puff size and propagation velocity. More importantly however, the theory is able to accurately predict the number and size spectra of the potentially virus-laden droplet nuclei that remain airborne within the puff. We observe that the ejection angle and mouth shape do not significantly affect the puff and droplet dynamics. Additionally, we quantify the carrying capacity of the detached puff portions in terms of the number and size spectra of droplets/droplet nuclei suspended within.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2022

Multi-area parameter error identification for large power systems

Power grid model parameters may contain errors due to various reasons. Detecting and correcting parameter errors typically requires significant computational effort due to the size and complexity of the parameter database. While the normalized Lagrange multiplier (NLM) method can effectively detect, identify and correct parameter errors, its computational burden could rapidly grow with increasing system size. This paper addresses this issue by proposing a multi-area parameter error identification method. Each area has its own outlier detection tool for detecting the incorrect parameters and measurements within the area. On the other hand, due to the reduced redundancy at area boundaries, parameter errors on branches incident to boundary buses may not be detected. Such errors are subsequently detected by a coordination level estimator completing the system-wide parameter detection procedure. In conclusion, performance of the developed method is demonstrated using the IEEE 118-bus and 2000-bus Texas synthetic systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A consistent and conservative volume distribution algorithm and its applications to multiphase flows using Phase-Field models

In the present study, the multiphase volume distribution problem, where there can be an arbitrary number of phases, is addressed using a consistent and conservative volume distribution algorithm. The proposed algorithm satisfies the summation constraint, the conservation constraint, and the consistency of reduction. The first application of the volume distribution algorithm is to determine the Lagrange multipliers in multiphase Phase-Field models that enforce the mass conservation, and a multiphase conservative Allen-Cahn model that satisfies the consistency of reduction is developed. A corresponding consistent and conservative numerical scheme is developed for the model. The multiphase conservative Allen-Cahn model has a better ability than the multiphase Cahn-Hilliard model to preserve under-resolved structures. The second application is to develop a numerical procedure, called the boundedness mapping, to map the order parameters, obtained numerically from a multiphase model, into their physical interval, and at the same time to preserve the physical properties of the order parameters. Along with the consistent and conservative schemes for the multiphase Phase-Field models, the numerical solutions of the order parameters are reduction consistent, conservative, and bounded, which are theoretically analyzed and numerically validated. Then, the multiphase Phase-Field models are coupled with the momentum equation by satisfying the consistency of mass conservation and the consistency of mass and momentum transport, thanks to the consistent formulation. Finally, it is demonstrated that the proposed model and scheme converge to the sharp-interface solution and are capable of capturing the complicated multiphase dynamics even when there is a large density and/or viscosity ratio.

42 ENGINEERING↗
DOE OSTI2020

Overset meshes for incompressible flows: On preserving accuracy of underlying discretizations

This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Specifically, using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical performance where the underlying discretization is second-order accurate. The first aspect of information exchange surrounds interpolation of solution where we examine Lagrange and point-cloud-based interpolation for creating constraint equations between overset meshes. To maintain overall second-order accuracy, higher-order interpolation is required for elliptic problems, but linear interpolation is sufficient for hyperbolic problems in first-order form. Higher-order point-cloud-based interpolation provides a pathway to maintaining accuracy in unstructured meshes, but at higher complexity. The second aspect of information exchange focuses on comparing the approaches of overset single system (OSS) and overset Additive Schwarz (OAS) for coupling the linear systems of the overlapping meshes. While the former involves a single linear system, in the latter the discrete linear systems are solved separately, and solving the global system is accomplished through outer iterations and sequential information exchange in a Jacobi fashion. For the test cases studied, accuracy for hyperbolic problems is maintained by performing two outer iterations, whereas many outer iterations are required for elliptic systems. The order-of-accuracy studies explored here are critical for verifying the overset-mesh coupling algorithms used in engineering simulations. Accuracy of these simulations themselves is, however, quantified using engineering quantities of interest such as drag, power, etc. Consequently, we conclude with numerical experiments using NS equations for incompressible flows where we show that linear interpolation and few outer iterations are sufficient for achieving asymptotic convergence of engineering quantities of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows

In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.

97 MATHEMATICS AND COMPUTING↗