Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian approach”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Load Shedding for Voltage Regulation With Probabilistic Agent Compliance

With the increased observability and controllability of distribution systems, the share of behind-the-meter systems is trending upwards rapidly. As a consequence, the impact of human behaviors on system performance can no longer be ignored and should be reflected in the energy management system models. In this paper, we discuss the problem of distribution system voltage control by active power curtailment where the agent compliance of the load curtailment signal is probabilistic. We discuss the modeling of the optimal voltage control problem with probabilistic agent compliance as a chance-constrained optimization problem, its tractable safe approximation using convex restriction, and a scenario-based mixed-integer reformulation as well as the associated solution method based on augmented Lagrangian method. The numerical simulation on IEEE test system validates the effectiveness of the proposed approach in obtaining high-quality feasible load curtailment signal with low computational cost, which makes it a viable tool for real time decision making.

augmented Lagrangian method↗
DOE OSTI2021

Fractons and Exotic Symmetries from Branes

The emerging study of fractons, a new type of quasi-particle with restricted mobility, has motivated the construction of several classes of interesting continuum quantum field theories with novel properties. One such class consists of foliated field theories which, roughly, are built by coupling together fields supported on the leaves of foliations of spacetime. Another approach, which we refer to as exotic field theory, focuses on constructing Lagrangians consistent with special symmetries (like subsystem symmetries) that are adjacent to fracton physics. A third framework is that of infinite-component Chern-Simons theories, which attempts to generalize the role of conventional Chern-Simons theory in describing (2+1)D Abelian topological order to fractonic order in (3+1)D. The study of these theories is ongoing, and many of their properties remain to be understood. Historically, it has been fruitful to study QFTs by embedding them into string theory. One way this can be done is via D-branes, extended objects whose dynamics can, at low energies, be described in terms of conventional quantum field theory. Further, QFTs that can be realized in this way can then be analyzed using the rich mathematical and physical structure of string theory. In this paper, we show that foliated field theories, exotic field theories, and infinite-component Chern-Simons theories can all be realized on the world-volumes of branes. We hope that these constructions will ultimately yield valuable insights into the physics of these interesting field theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Computational study of the ECN spray C via one-way coupling of internal nozzle flow and ensuing spray

Here, in this study, computational fluid dynamics simulations of the Engine Combustion Network Spray C injector were performed to investigate the effect of the internal nozzle flow on the development of the ensuing spray in automotive fuel injectors. A set of simulation best practices was first developed to achieve an accurate representation of the internal nozzle flow. Detailed, spatiotemporally defined information was extracted at the injector’s orifice exit and stored in maps that were subsequently employed to initialize a Lagrangian spray by means of a static, one-way coupling (OWC) approach. To highlight the impact of the in-nozzle flow on the spray, the OWC approach was compared with the conventional rate-of-injection (ROI) approach, for which the spray parcels are commonly initialized using a blob injection model. The results indicated that, using the ROI approach, considerable evaporation was observed within the plume in the near-nozzle region. On the other hand, the OWC approach showed evidence of evaporation farther downstream and at the plume periphery. The different behavior observed with the two approaches was mostly due to substantial differences in fuel mass distribution. The liquid fuel concentration predicted by the OWC approach matched well with X-ray tomography measurements of projected fuel mass per unit volume, in terms of both magnitude and distribution. In contrast, an unphysical distribution of fuel mass characterized by an extremely low liquid concentration was predicted by the ROI approach. In the far-field, the effect of the in-nozzle flow was found to be weaker as the spray penetrated and interacted with the ambient gas. In this region, the OWC and ROI approaches provided similar predictions of vapor penetration, spray morphology, and mixing-related characteristics.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effective field theories for collective excitations of atomic nuclei

Collective modes emerge as the relevant degrees of freedom that govern low-energy excitations of atomic nuclei. These modes—rotations, pairing rotations, and vibrations—are separated in energy from non-collective excitations, making it possible to describe them in the framework of effective field theory. Rotations and pairing rotations are the remnants of Nambu–Goldstone modes from the emergent breaking of rotational symmetry and phase symmetries in finite deformed and finite superfluid nuclei, respectively. The symmetry breaking severely constrains the structure of low-energy Lagrangians and thereby clarifies what is essential and simplifies the description. The approach via effective field theories exposes the essence of nuclear collective excitations and is defined with a breakdown scale in mind. This permits one to make systematic improvements and to estimate and quantify uncertainties. Effective field theories of collective excitations have been used to compute spectra, transition rates, and other matrix elements of interest. In particular, predictions of the nuclear matrix element for neutrinoless double beta decay then come with quantified uncertainties. This review summarizes these results and also compares the approach via effective field theories to well-known models and ab initio computations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Emergence of the polydeterminant in QCD

A generalization of the determinant appears in particle physics in effective Lagrangian interaction terms that model the chiral anomaly in quantum chromodynamics (Giacosa et al. in Phys Rev D 97(9):091901, 2018, Phys Rev D 109(7):L071502, 2024), in particular in connection to mesons. This polydeterminant function, known in the mathematical literature as a mixed discriminant, associates N distinct N x N complex matrices into a complex number and reduces to the usual determinant when all matrices are taken as equal. Here, we explore the main properties of the polydeterminant applied to (quantum) fields by using a formalism and a language close to high-energy physics approaches. We discuss its use as a tool to write down novel chiral anomalous Lagrangian terms and present an explicit illustrative model for mesons. Finally, the extension of the polydeterminant as a function of tensors is shown.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Nonlinear, real-time optimization for actuator management in tokamaks

Experiments in DIII-D have been carried out to test a novel actuator management approach in tokamaks. Here, the actuator management scheme is posed as a nonlinear-optimization problem in which the actuator commands are calculated in real time according to the changing control priorities, plasma state, and actuator availability. Such optimization problem is solved using the augmented Lagrangian method, combined with a gradient projection method and a conjugate-gradient iteration algorithm. The algorithmic approach followed in this work does not depend on the particular control objectives or actuators considered, which facilitates its integration with other independently-designed control components within a plasma-control system. In addition, the actuator-management algorithm is able to handle the optimization problem in a computationally efficient manner, making it suitable for real-time implementations. Initial DIII-D results in the steady-state high-q min scenario have demonstrated the capabilities of the actuator manager to perform both simultaneous multiple mission and repurposing sharing, which will be required in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A computationally efficient algorithm for computing convex hull prices

Electricity markets worldwide allow participants to bid non-convex production offers. While non-convex offers can more accurately reflect a resource's capabilities, they create challenges for market clearing processes. For example, system operators may be required to execute side payments to participants whose costs are not covered through energy sales as determined via traditional locational marginal pricing schemes. Convex hull pricing minimizes this and other types of side payments while providing uniform (i.e., locationally and temporally consistent) prices. Computing convex hull prices involves solving either a large-scale linear program or the Lagrangian dual of the corresponding non-convex scheduling problem. Further, the former approach requires explicit descriptions of market participants' convex hulls. While linear programs for computing convex hull prices are large, their structure is naturally decomposable by generators. Here, in this work, we propose and empirically analyze a Benders decomposition approach to computing convex hull prices that leverages recent advances in convex hull formulations for thermal generating units. We demonstrate across a large set of test instances that our decomposition approach only requires modest computational effort, obtaining solutions at least an order of magnitude faster than the equivalent large-scale linear programming approach. Overall, we provide a computationally feasible method for computing convex hull prices for industrial scale market clearing problems, enabling the possibility of practical adoption of this advanced pricing mechanism.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Bayesian Entropy Neural Networks for physics-aware prediction

This article addresses the need for deep learning models to integrate well-defined constraints into their outputs, driven by their application in surrogate models, learning with limited data and partial information, and scenarios requiring flexible model behavior to incorporate non-data sample information. We introduce Bayesian Entropy Neural Networks (BENN), a framework grounded in Maximum Entropy (MaxEnt) principles, designed to impose constraints on Bayesian Neural Network (BNN) predictions. BENN is capable of constraining not only the predicted values but also their derivatives and variances, ensuring a more robust and reliable model output. To achieve simultaneous uncertainty quantification and constraint satisfaction, we employ the method of multipliers approach. This allows for the concurrent estimation of neural network parameters and the Lagrangian multipliers associated with the constraints. Our experiments, spanning diverse applications such as beam deflection modeling and microstructure generation, demonstrate the effectiveness of BENN. The results highlight significant improvements over traditional BNNs and showcase competitive performance relative to contemporary constrained deep learning methods.

14 SOLAR ENERGY↗
NASA NTRS1985

Combined aerodynamic and structural dynamic problem emulating routines (CASPER): Theory and implementation

The Combined Aerodynamic and Structural Dynamic Problem Emulating Routines (CASPER) is a collection of data-base modification computer routines that can be used to simulate Navier-Stokes flow through realistic, time-varying internal flow fields. The Navier-Stokes equation used involves calculations in all three dimensions and retains all viscous terms. The only term neglected in the current implementation is gravitation. The solution approach is of an interative, time-marching nature. Calculations are based on Lagrangian aerodynamic elements (aeroelements). It is assumed that the relationships between a particular aeroelement and its five nearest neighbor aeroelements are sufficient to make a valid simulation of Navier-Stokes flow on a small scale and that the collection of all small-scale simulations makes a valid simulation of a large-scale flow. In keeping with these assumptions, it must be noted that CASPER produces an imitation or simulation of Navier-Stokes flow rather than a strict numerical solution of the Navier-Stokes equation. CASPER is written to operate under the Parallel, Asynchronous Executive (PAX), which is described in a separate report.

Jones, William H.↗
NASA NTRS2013

Large-Deformation Displacement Transfer Functions for Shape Predictions of Highly Flexible Slender Aerospace Structures

Large deformation displacement transfer functions were formulated for deformed shape predictions of highly flexible slender structures like aircraft wings. In the formulation, the embedded beam (depth wise cross section of structure along the surface strain sensing line) was first evenly discretized into multiple small domains, with surface strain sensing stations located at the domain junctures. Thus, the surface strain (bending strains) variation within each domain could be expressed with linear of nonlinear function. Such piecewise approach enabled piecewise integrations of the embedded beam curvature equations [classical (Eulerian), physical (Lagrangian), and shifted curvature equations] to yield closed form slope and deflection equations in recursive forms.

Deflection equations↗
DOE OSTI2022

New opportunities for Insight into Emergence of Hadron Mass from Studies of Nucleon Resonance Electroexcitation

Understanding the emergence of hadron mass (EHM) represents one of the most challenging and still open problems in the Standard Model. Results of EHM studies from the nucleon resonance electroexcitation amplitudes (gvpN* electrocouplings) determined from exclusive meson electroproduction data measured with the CLAS detector will be presented in this talk. These CLAS data have provided the first and only available results on the evolution of the ?vpN* electrocouplings with photon virtuality Q2 up to 5 GeV2 for most nucleon resonances in the mass range up to 1.8 GeV. A successful description of the CLAS results on the Q2-evolution of the ?(1232)3/2+, N(1440)1/2+, and ?(1600)3/2+ electrocouplings achieved within the CSM approach by employing the same dressed quark mass function inferred from the QCD Lagrangian, which also provided a successful description of the pion elastic electromagnetic form factor and PDF, and nucleon elastic form factor, conclusively demonstrated the capability of gaining insight into EHM. The expected results from the CLAS12 detector on the Q2-evolution of the gvpN* electrocouplings for all prominent resonances in the mass range up to 2 GeV at the still almost unexplored range of Q2>5 GeV2 will allow us to place tight constraints on the dressed quark mass function in the range of quark momenta where the substantial part of hadron mass ~50% is expected to be generated. Consistent results on the momentum dependence of the dressed quark mass obtained from analyses of the ?vpN* electrocouplings of nucleon resonances of different structure will validate insight into EHM in a nearly model-independent way, allowing us to address key open problems in hadron physics on the emergence of hadron mass and nucleon resonance structure from QCD.

Mokeev, Viktor↗
NASA NTRS1985

Translational and extensional energy release rates (the J- and M-integrals) for a crack layer in thermoelasticity

A number of papers have been presented on the evaluation of energy release rate for thermoelasticity and corresponding J integral. Two main approaches were developed to treat energy release rate in elasticity. The first is based on direct calculation of the potential energy rate with respect to crack length. The second makes use of Lagrangian formalism. The translational and expansional energy release rates in thermoelasticity are studied by employing the formalism of irreversible thermodynamics and the Crack Layer Approach.

Chudnovsky, A.↗
DOE OSTI2026

A Decomposition-Based Learn-To-Optimize Approach with Feasibility Layer Assistance for Sub-Hourly Unit Commitment

Sub-hourly unit commitment (UC) with 15-min intervals is gaining significant attention as a way to respond rapidly to the fluctuations in electricity supply and demand introduced by renewable resources. However, the increased temporal resolution and complex inter-temporal dependencies pose substantial computational challenges for traditional optimization methods. To this end, this paper explores a decomposition-based learn-to-optimize approach. Building on recent advances in machine learning, our method revisits the long- overlooked Lagrangian relaxation framework, which is a classical decomposition technique that enables tractable subproblem solving. These smaller subproblems are inherently well-suited for machine learning, as their reduced dimensionality and structural regularity allow predictive models to efficiently learn and generalize solution patterns. We thus propose a generic predictive model, which embeds Gated Recurrent Units (GRUs) and Attention in the encoder-decoder structure, and integrate a rule-based feasibility layer to capture temporal dependencies, reduce training effort, and improve feasibility w.r.t. unit-level constraints. Our method has been validated on the IEEE 118-bus system, demonstrating promising performance in solving sub-hourly UC problems efficiently and feasibly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Toward Polydisperse Flows With MFIX-EXA

In the presence of large size disparities, single-grid neighbor search algorithms lead to inflated neighbor lists that significantly degrade the performance of Lagrangian particle solvers. If Eulerian–Lagrangian (EL) frameworks are to remain performant when simulating realistic systems, improved neighbor detection approaches must be adopted. To this end, we consider the application of a multigrid neighbor search (MGNS) algorithm in the mfix-exa software package, an exascale EL solver built upon the AMReX library. Here, details regarding the implementation and verification of MGNS are provided along with speedup curves for a bidisperse mixing layer. MGNS is shown to yield up to 15$\times$ speedup on CPU and 6$\times$ speedup on GPU for the problems considered here. The mfix-exa software is then validated for a variety of polydisperse flows. Finally, a brief discussion is given for how dynamic MGNS may be completed, with application to spatially varying particle size distributions.

42 ENGINEERING↗
NASA NTRS1992

Truss topology optimization with simultaneous analysis and design

Strategies for topology optimization of trusses for minimum weight subject to stress and displacement constraints by Simultaneous Analysis and Design (SAND) are considered. The ground structure approach is used. A penalty function formulation of SAND is compared with an augmented Lagrangian formulation. The efficiency of SAND in handling combinations of general constraints is tested. A strategy for obtaining an optimal topology by minimizing the compliance of the truss is compared with a direct weight minimization solution to satisfy stress and displacement constraints. It is shown that for some problems, starting from the ground structure and using SAND is better than starting from a minimum compliance topology design and optimizing only the cross sections for minimum weight under stress and displacement constraints. A member elimination strategy to save CPU time is discussed.

Sankaranarayanan, S.↗
DOE OSTI2026

CoLoRe-2LPT: Lyman-$α$ mock catalogues for the validation of DESI cosmological analyses

The Lyman-$α$ (Ly$α$) forest has become a crucial probe for studying the large-scale structure of the universe at high redshift ($z > 2$), providing powerful constraints on Baryon Acoustic Oscillations (BAO) and the full-shape (FS) clustering of matter. As a key ingredient for upcoming BAO and FS analyses, we present a new generation of fast cosmological Ly$α$ mocks based on second-order Lagrangian perturbation theory (2LPT). These new mocks significantly improve upon previous log-normal approaches, both at accurately capturing small scale clustering and at recovering the non-linear broadening of the BAO peak. They are able to reproduce Ly$α$ statistics within $10\%$ of the latest DESI measurement; including the Ly$α$ bias and the redshift-space distortion $β$ parameter, mean transmitted flux, and 1D power spectrum. The corresponding quasar (QSO) clustering is also improved with respect to previous approaches, calibrated against high-resolution Abacus simulations, recovering the observational QSO linear bias to less than $5\%$ and improving redshift-space distortions via 2LPT velocities and the addition of Fingers-of-God effects. Furthermore, these mocks incorporate high column density systems and metal lines, allowing us to explore the effects and systematics induced by these astrophysical contaminants. This new set of mocks has been key for enhancing the modeling and validation of the DESI DR2 Ly$α$ full shape cosmological analysis. This work provides a physically motivated and computationally efficient tool for simulating current and next-generation Ly$α$ surveys and validating FS and BAO analysis.

Bernal, M.F. Ruiz-Herrera [Madrid, CIEMAT] (ORCID:↗