The Similarity Structure Function Lagrangian CFD and the Finite Scale Theory
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
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.
Electromagnetic and gravitational interactions from Lagrangian mechanics
Highlights: • Background Maxwell and Einstein fields are universal features of Lagrangian systems. • The Lagrangian function is promoted to a dynamic field. • The Einstein–Maxwell theory is derived from Lagrangian mechanics. Background fields of electromagnetic and gravitational type emerge in the low kinetic energy limit of any regular Lagrangian system and, in particular, in the corresponding limit of any spacetime theory in which the free motion of test particles is described by an unspecified regular Lagrangian. Electromagnetic and gravitational type interactions are therefore a universal feature of low kinetic energy Lagrangian systems. These background interactions can be consistently turned into dynamic Einstein–Maxwell fields by promoting the Lagrangian function to a dynamic scalar field on the tangent bundle of the configuration space. Accordingly, Einstein–Maxwell theory can be deduced from the assumption that the motion of elementary test particles in spacetime is described by Lagrangian mechanics. For higher kinetic energy-type values, identified with the square of the invariant mass of the particle, the Lagrangian induces higher rank interactions that seem however too weak to have been detected in spacetime physics, but which might prove relevant at the Planck scale.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Arbitrary-order Bernstein basis functions for Lagrangian Hydrodynamics [Slides]
Abstract not provided.
97 MATHEMATICS AND COMPUTING↗
An Immersed Finite Element Lagrangian-Eulerian Code-Coupling Framework
This report presents an assessment of immersed Eulerian-Lagrangian code-coupling techniques suitable for use in a broad range of mechanics applications. The coupling algorithm is based on an immersed finite element method that considers the Lagrangian and Eulerian overlap regions in the overall variational formulation. In this report the basic formulation details are presented followed by various aspects of the code-coupling algorithm using OpenIFEM as the Lagrangian/coupling framework. A series of representative test cases that illustrate the code-coupling algorithm are discussed. The current work provides an in-depth investigation into the immersed finite element method for the purposes of providing a rigorous coupling technique that is minimally invasive in the respective Eulerian and Lagrangian codes. A number of extensions to the base immersed finite element method have been examined. These extension include nodal and quadrature-based indicator functions, a Lagrangian volume-fraction calculation in regions of overlap, and the use of penalty constraints between the Lagrangian and Eulerian domains. A unique MPI-based coupling strategy that retains the independent MPI structure of each code has been demonstrated.
97 MATHEMATICS AND COMPUTING↗
Ergodic Lagrangian dynamics in a superhero universe
We present a fictional scenario that, while undeniably whimsical, provides the foundation for a unique exercise in extended problem solving, physics analysis, and quantitative model development. Starting with the foundational premise of the Wild Cards shared-world superhero universe, we demonstrate how a variety of concepts appropriate to the advanced undergraduate level—ergodicity, functional analysis, Lagrangian mechanics, and the ever-important simplifying approximation—can be combined into a rich, coherent mathematical model. The goal of this case study is to develop a useful pedagogical exercise in exploring an open-ended research question that presents, at first glance, no clear path forward. Being both eclectic and lengthy, this exercise offers a unique way for students to apply their core physics and mathematics education. It is perhaps best used within a senior honors seminar or within a brief (e.g., January term) elective class.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Some aspects of the five-dimensional Lovelock black hole spacetime: Strong homotopy retract, perihelion precession and quasistationary levels
Highlights: • Some topological aspects are explored. • The perihelion precession is calculated. • The general Heun functions are applied to obtain analytical solutions. • The quasibound states and the radial wave eigenfunctions are obtained. In this work we explore some mathematical physics aspects of the spherically symmetric Lovelock black hole in high dimensions. Intended for this aim, we thoroughly consider the metric corresponding to the five-dimensional Lovelock black hole spacetime. We construct the strong retractions by the geodesic equations on the background under consideration. As a result, from the topological point of view, we construct the theory of strong homotopy retract, which will allow us, in principle, to better understand some of its suitable applications on astrophysics and cosmology, in particular, in the analysis of the spacetime singularities. We find the solutions of the equation of motion for both radial and angular coordinates, and then we describe the outer (“exterior”) and lower (“interior”) apparent horizons. Indeed, the outer apparent horizon is the last surface from which the light waves could still escape from the black hole. Thus, it is meaningful to analyze some physical phenomena related to quantum particles propagating outside the exterior apparent horizon, in particular, we discuss the quasistationary levels of scalar fields and their radial wave functions, which are given in terms of the general Heun functions. We also calculate the perihelion precession in this background.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Differentiable programming for online training of a neural artificial viscosity function within a staggered grid Lagrangian hydrodynamics scheme
Lagrangian methods to solve the inviscid Euler equations produce numerical oscillations near shock waves. A common approach to reducing these oscillations is to add artificial viscosity (AV) to the discrete equations. The AV term acts as a dissipative mechanism that attenuates oscillations by smearing the shock across a finite number of computational cells. However, AV introduces several control parameters that are not determined by the underlying physical model, and hence, in practice are tuned to the characteristics of a given problem. We seek to improve the standard quadratic-linear AV form by replacing it with a learned neural function that reduces oscillations relative to exact solutions of the Euler equations, resulting in a hybrid numerical-neural hydrodynamic solver. Because AV is an artificial construct that exists solely to improve the numerical properties of a hydrodynamic code, there is no offline ‘viscosity data’ against which a neural network can be trained before inserting into a numerical simulation, thus requiring online training. We achieve this via differentiable programming, i.e. end-to-end backpropagation or adjoint solution through both the neural and differential equation code, using automatic differentiation of the hybrid code in the Julia programming language to calculate the necessary loss function gradients. A novel offline pre-training step accelerates training by initializing the neural network to the default numerical AV scheme, which can be learned rapidly by space-filling sampling over the AV input space. We find that online training over early time steps of simulation is sufficient to learn a neural AV function that reduces numerical oscillations in long-term hydrodynamic shock simulations. These results offer an early proof-of-principle that online differentiable training of hybrid numerical schemes with novel neural network components can improve certain performance aspects existing in purely numerical schemes.
97 MATHEMATICS AND COMPUTING↗
Multi-parametric analysis for mixed integer linear programming: An application to transmission upgrade and congestion management
Upgrading the capacity of existing transmission lines is essential for meeting the growing energy demands, facilitating the integration of renewable energy, and ensuring the security of the transmission system. This study focuses on the selection of lines whose capacities and by how much should be expanded from the perspective of the Independent System Operators (ISOs) to minimize the total system cost. We employ advanced multi-parametric programming and an enhanced branch-and-bound algorithm to address complex mixed-integer linear programming (MILP) problems, considering multi-period time constraints and physical limitations of generators and transmission lines. To characterize the various decisions in transmission expansion, we model the increased capacity of existing lines as parameters within a specified range. This study first relaxes the binary variables to continuous variables and applies the Lagrange method and Karush-Kuhn-Tucker (KKT) conditions to obtain optimal solutions and identify critical regions associated with active and inactive constraints. Moreover, we extend the traditional branch-and-bound (B&B) method by determining the problem’s upper and lower bounds at each node of the B&B decision tree, helping to manage computational challenges in large-scale MILP problems. Here, we compare the difference between the upper and lower bounds to obtain an approximate optimal solution within the decision-makers’ tolerable error range. In addition, the first derivative of the objective function on the parameters of each line is used to inform the selection of lines for easing congestion and maximizing social welfare. Finally, the capacity upgrades are selected by weighing the reductions in system costs against the expense of upgrading line capacities. The findings are supported by numerical simulations and provide transmission-line planners with decision-making guidance.
24 POWER TRANSMISSION AND DISTRIBUTION↗
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↗
Quantum field theory with the generalized uncertainty principle II: Quantum Electrodynamics
Highlights: • Relativistic Generalized Uncertainty Principle gives Frame independent minimum length. • Quantum Gravity modified Dirac equation from Modified Klein–Gordon Equations. • There exists a Quantum Gravity modified Lagrangian for a spinor field theory. • Feynman vertices contain 2 fermions & up to 5 gauge bosons are allowed. • Minimum length modifies the amplitude of Electrodynamic electron–muon scattering. Continuing our earlier work on the application of the Relativistic Generalized Uncertainty Principle (RGUP) to quantum field theories, in this paper we study Quantum Electrodynamics (QED) with minimum length. We obtain expressions for the Lagrangian, Feynman rules and scattering amplitudes of the theory, and discuss their consequences for current and future high energy physics experiments. We hope this will provide an improved window for testing Quantum Gravity effects in the laboratory.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Helicity, spin, and infra-zilch of light: A Lorentz covariant formulation
Highlights: • Spin and orbital parts of angular momentum. • The helicity tensor contains the helicity, spin and infra-zilch currents. • The helicity tensor is a Noether current for a duality-symmetric action. • Symmetry generator for the helicity tensor. • The helicity tensor is conserved in Lorenz gauge. In this paper, a novel conserved Lorentz covariant tensor, termed the helicity tensor, is introduced in Maxwell theory. The conservation of the helicity tensor expresses the conservation laws contained in the helicity array, introduced by Cameron et al. (2012), including helicity, spin, as well as the spin-flux or infra-zilch. The Lorentz covariance of the helicity tensor is in contrast to previous formulations of the helicity hierarchy of conservation laws, which required the non-Lorentz covariant transverse gauge. The helicity tensor is shown to arise as a Noether current for a variational symmetry of a duality-symmetric Lagrangian for Maxwell theory. This symmetry transformation generalizes the duality symmetry and includes the symmetry underlying the conservation of the spin part of the angular momentum.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Quasinormal modes for non-minimally coupled scalar fields in regular black hole spacetimes: Grey-body factors, area spectrum and shadow radius
Black hole perturbation theory is a useful approach to study interactions between black holes and fundamental fields. A particular class of black hole solutions arising out of modification of Einstein’s general theory of relativity are regular black holes (RBHs) which can be constructed using a nonlinear electrodynamic Lagrangian. Because of their importance, we are interested in studying the behavior of three kinds of such RBHs under perturbations generated by an external field. Indeed, we investigate the quasinormal modes (QNMs) of a massive scalar field propagating near the RBHs which is non-minimally coupled to the Ricci scalar-tensor of background geometry. We will attempt to find the low-lying quasinormal frequencies of the perturbations by using WKB approximation. We shall also study the relationship between the QNMs and some characteristic properties of black holes such as grey-body factors, area quantization, and shadow radius.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
The repulsive Casimir-type forces in binary mixture of weakly interacting Bose–Einstein condensates gases
Highlights: • The equations of the state of a binary mixture weakly interacting Bose gases are established based on the Lagrangian of the system. • The formulae of the Casimir-type energy and the resulting Casimir-type force are produced from the quantum part of the grand canonical energy. • The results for the Casimir-type effect in the one-field case can be re-obtained by using the results of the two-field case. Within the framework of the quantum field theory (QFT) in the one-loop approximation, we investigate the Casimir-type effect in a mixture of two-component weakly interacting Bose-Einstein condensates (BECs) gases confined between two parallel plates at zero temperature. The Casimir-type forces in the miscible state of BECs for different boundary conditions (BCs) are discovered and compared to each other. The repulsive nature of the Casimir-type force in BECs is defined for the first time. Our calculations are finished in both grand canonical ensemble (GCE) and canonical ensemble (CE).
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Hierarchical Speed Planner for Automated Vehicles: A Framework for Lagrangian Variable Speed Limit in Mixed-Autonomy Traffic
Here, this article presents a novel hierarchical speed planning framework for variable speed limits in mixed-autonomy traffic environments, leveraging server-side macroscopic control and vehicle-side microscopic execution. The framework integrates real-time traffic state estimation (TSE) and reinforcement learning (RL)-based control to mitigate congestion and improve traffic flow. A TSE enhancement module combines macroscopic data from sources like INRIX with high-resolution observations from connected autonomous vehicles (CAVs), enabling predictive modeling to address latency and noise. The target speed design module employs kernel smoothing and a buffer zone strategy to optimize traffic density and flow around bottlenecks. The proposed system was validated in the largest open-road test to date with 100 CAVs, demonstrating an overall 8% traffic density decrease, with a specific decrease of 7% upstream, 10% downstream, and a 52% decrease during the congestion formation phase at bottlenecks.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Design, Preparation, and Execution of the 100-AV Field Test for the CIRCLES Consortium: Methodology and Implementation of the Largest Mobile Traffic Control Experiment to Date
This article presents the comprehensive design, setup, execution, and evaluation of the MegaVanderTest (MVT) experiment conducted by the Congestion Impacts Reduction via CAV-in-the-Loop Lagrangian Energy Smoothing (CIRCLES) Consortium, which aimed to mitigate traffic congestion using partially autonomous vehicles (AVs) (see “Summary”). The experiment involved 100 vehicles on Nashville’s Interstate 24 (I-24) highway, utilizing various control algorithms to smooth stop-and-go traffic waves. The execution of the MVT experiment required a coordinated effort from multiple teams. This article details the meticulous planning process, the coordinated efforts of multiple teams, and the innovative use of a dynamic agent-based simulation framework for traffic evaluation. Here, the contributions of this work include demonstrating and providing a detailed roadmap for large-scale live traffic experiments, illustrating the lessons learned from the MVT experiment, and introducing the other articles in this issue and their complementary relationship in the MVT experiment.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Energy-Efficient Driving in Connected Corridors via Minimum Principle Control: Vehicle-in-the-Loop Experimental Verification in Mixed Fleets
Connected and automated vehicles (CAVs) can plan and actuate control that explicitly considers performance, system safety, and actuation constraints in a manner more efficient than their human-driven counterparts. In particular, eco-driving is enabled through connected exchange of information from signalized corridors that share their upcoming signal phase and timing (SPaT). This is accomplished in the proposed control approach, which follows first principles to plan a free-flow acceleration-optimal trajectory through green traffic light intervals by Pontryagin's Minimum Principle in a feedback manner. Urban conditions are then imposed from exogeneous traffic comprised of a mixture of human-driven vehicles (HVs) - as well as other CAVs. As such, safe disturbance compensation is achieved by implementing a model predictive controller (MPC) to anticipate and avoid collisions by issuing braking commands as necessary. The control strategy is experimentally vetted through vehicle-in-the-loop (VIL) of a prototype CAV that is embedded into a virtual traffic corridor realized through microsimulation. Up to 36% fuel savings are measured with the proposed control approach over a human-modelled driver, and it was found connectivity in the automation approach improved fuel economy by up to 26% over automation without. Additionally, the passive energy benefits realizable for human drivers when driving behind downstream CAVs are measured, showing up to 22% fuel savings in a HV when driving behind a small penetration of connectivity-enabled automated vehicles.
33 ADVANCED PROPULSION SYSTEMS↗
Analytic error analysis of cross section interpolation methods in nodal diffusion codes - II: Numerical results
This paper is the second part of a two-part paper that documents the numerical results for the partial derivatives model presented in part I. In this paper, we derive the error bounds for the analytical point-wise error expression and verify our bounds with numerical experiments. The point-wise error expressions make available, and bound, the sources that contribute to the total error of the interpolated cross section in terms of the Lagrange interpolation errors and the model form error. MPACT is used to generate two-group homogenized cross sections for Westinghouse's AP1000 Region 4 lattice to evaluate the accuracy of the bounds. Error bounds calculated over a grid are compared to numerical data for uni-variate and multi-variate interpolation. The point-wise error bounds of a typical case matrix - two branches in each state variable - are displayed for bi-variate interpolation in the state variables: moderator density, fuel temperature, and boron concentration. The error bounds are shown to be highly accurate compared to numerical results, and in accordance with the underlying physics. We then discuss and show how the sources of error contribute to the total error, and consider the improvement of each error source. Finally, we mention future work such as propagating our cross section error bounds through a reactivity calculation. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
On the Detection of Habitable Trojan Planets in the Kepler Circumbinary Systems
We present the results of a study of the prospect of detecting habitable Trojan planets in the Kepler Habitable Zone circumbinary planetary systems (Kepler-16, -47, -453, -1647, and -1661). We integrated the orbits of 10,000 separate N-body systems (N = 4, 6), each with a one Earth-mass body in a randomly selected orbit near the L {sub 4} and L {sub 5} Lagrangian points of the host HZ circumbinary planet. We find that stable Trojan planets are restricted to a narrow range of semimajor axes in all five systems and limited to small eccentricities in Kepler-16, -47, and -1661. To assess the prospect of the detection of these habitable Trojan planets, we calculated the amplitudes of the variations they cause in the transit timing of their host bodies. Results show that the mean amplitudes of the transit timing variations (TTVs) correlate with the mass of the transiting planet and range from 70 minutes for Kepler-16b to 390 minutes for Kepler-47c. Our analysis indicates that the TTVs of the circumbinary planets caused by these Trojan bodies fall within the detectable range of timing precision obtained from the Kepler telescope’s long-cadence data. The latter points to Kepler data as a viable source to search for habitable Trojan planets.
47 OTHER INSTRUMENTATION↗