Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Variational energy principle for compressible, baroclinic flow. 1: First and second variations of total kinetic action

The case of a cold gas in the absence of external force fields is considered. Since the only energy involved is kinetic energy, the total kinetic action (i.e., the space-time integral of the kinetic energy density) should serve as the total free-energy functional in this case, and as such should be a local minimum for all possible fluctuations about stable flow. This conjecture is tested by calculating explicit, manifestly covariant expressions for the first and second variations of the total kinetic action in the context of Lagrangian kinematics. The general question of the correlation between physical stability and the convexity of any action integral that can be interpreted as the total free-energy functional of the flow is discussed and illustrated for the cases of rectillinear and rotating shearing flows.

Schmid, L. A.↗
DOE OSTI2021

Optimization-based algorithms for nonlinear mechanics and frictional contact

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

42 ENGINEERING↗
NASA NTRS1982

Environmental solid particle effects on compressor cascade performance

The effect of suspended solid particles on the performance of the compressor cascade was investigated experimentally in a specially built cascade tunnel, using quartz sand particles. The cascades were made of NACA 65(10)10 airfoils. Three cascades were tested, one accelerating cascade and two diffusing cascades. The theoretical analysis assumes inviscid and incompressible two dimensional flow. The momentum exchange between the fluid and the particle is accounted for by the interphase force terms in the fluid momentum equation. The modified fluid phase momentum equations and the continuity equation are reduced to the conventional stream function vorticity formulation. The method treats the fluid phase in the Eulerian system and the particle phase in Lagrangian system. The experimental results indicate a small increase in the blade surface static pressures, while the theoretical results indicate a small decrease. The theoretical analysis, also predicts the loss in total pressure associated with the particulate flow through the cascade.

Tabakoff, W.↗
NASA NTRS2000

An Overview of the NCC Spray/Monte-Carlo-PDF Computations

This paper advances the state-of-the-art in spray computations with some of our recent contributions involving scalar Monte Carlo PDF (Probability Density Function), unstructured grids and parallel computing. It provides a complete overview of the scalar Monte Carlo PDF and Lagrangian spray computer codes developed for application with unstructured grids and parallel computing. Detailed comparisons for the case of a reacting non-swirling spray clearly highlight the important role that chemistry/turbulence interactions play in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that some of the combustion occurs in a predominantly premixed-flame environment and the rest in a predominantly diffusion-flame environment. However, the non-PDF solution predicts wrongly for the combustion to occur in a vaporization-controlled regime. Near the premixed flame, the Monte Carlo particle temperature distribution shows two distinct peaks: one centered around the flame temperature and the other around the surrounding-gas temperature. Near the diffusion flame, the Monte Carlo particle temperature distribution shows a single peak. In both cases, the computed PDF's shape and strength are found to vary substantially depending upon the proximity to the flame surface. The results bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations. Finally, we end the paper by demonstrating the computational viability of the present solution procedure for its use in 3D combustor calculations by summarizing the results of a 3D test case with periodic boundary conditions. For the 3D case, the parallel performance of all the three solvers (CFD, PDF, and spray) has been found to be good when the computations were performed on a 24-processor SGI Origin work-station.

Raju, M. S.↗
DOE OSTI2024

Influence functions for the Delaware Oil and Gas Basin Tower Network

This data consists of influence functions created for the Delaware tower network for the purpose of calculating methane enhancements in the region. Two different footprint models were used. The Lagrangian Particle Dispersion model (LPDM, Uliasz 1993) was run from March 2020-May 2023. The model uses in-house WRF simulations for the meteorological transport as described in Barkley et al., (2023). Nine hundred (900) particles were released during each hourly period and are traced back 72 hours backwards in time. A static binning height of 110 m AGL was used to calculate surface influence. After May 2023, the LPDM simulations were phased out and replaced with influence functions created using HYSPLIT under the STILT configuration (Lin et al., 2003). The STILT influence functions use HRRR for the meteorological transport (Dowell et al., 2022). Nine hundred (900) particles were released during each hourly period and are traced back 72 hours backwards in time. A binning height of half the boundary layer height was used to calculate surface influence. Footprints are currently available from Sept 2020-December 2024 and will be updated periodically as more are generated in time. No STILT influence functions were created for tower site Maljimar, which was decommissioned in March 2022. Lat/Lon grid information for the LPDM influence functions is contained in lpdm-lat-lon-grid.nc. STILT grid information is contained within the influence functions themselves. Tower coordinates and heights can be found in tower-info.csv

Barkley, Z. R.↗
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 OSTI2022

CoSyR: A novel beam dynamics code for the modeling of synchrotron radiation effects

The self-consistent nonlinear dynamics of a relativistic charged particle beam interacting with its complete self-fields is a fundamental problem underpinning many of the accelerator design issues in high brightness beam applications, as well as the development of advanced accelerators. Particularly, synchrotron radiation induced effects in a magnetic dispersive beamline element can lead to collective beam instabilities and emittance growth. A novel beam dynamics code is developed based on a Lagrangian method for the calculation of the particles’ radiation near-fields using wavefront/wavelet meshes via the Green’s function of the Maxwell equations. These fields are then interpolated onto a moving mesh for dynamic update of the beam. This method allows radiation co-propagation and self-consistent interaction with the beam in 2D/3D simulations at greatly reduced numerical errors. Multiple levels of parallelisms are inherent in this method and implemented in our code CoSyR to enable at-scale simulations of nonlinear beam dynamics on modern computing platforms using MPI, multi-threading, and GPUs. Here, the current 2D implementation of CoSyR has been used to evaluate the transverse and longitudinal coherent radiation effects on the beam and to investigate beam optics designs proposed for mitigation of beam brightness degradation in a magnetic bunch compressor. In this paper, the design of CoSyR, as well as the benchmark with other coherent synchrotron radiation models, are described and discussed. Extension of the core algorithms to 3D is possible and planned.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

A new theory of gravity.

A new relativistic theory of gravity is presented. This theory agrees with all experiments to date. It is a metric theory; it is Lagrangian-based; and it possesses a preferred frame with conformally flat space slices. With an appropriate choice of certain adjustable functions and parameters and of the cosmological model, this theory possesses precisely the same post-Newtonian limit as general relativity.

Ni, W.-T.↗
NASA NTRS1983

Lagrange thermodynamic potential and intrinsic variables for He-3 He-4 dilute solutions

For a two-fluid model of dilute solutions of He-3 in liquid He-4, a thermodynamic potential is constructed that provides a Lagrangian for deriving equations of motion by a variational procedure. This Lagrangian is defined for uniform velocity fields as a (negative) Legendre transform of total internal energy, and its primary independent variables, together with their thermodynamic conjugates, are identified. Here, similarities between relations in classical physics and quantum statistical mechanics serve as a guide for developing an alternate expression for this function that reveals its character as the difference between apparent kinetic energy and intrinsic internal energy. When the He-3 concentration in the mixtures tends to zero, this expression reduces to Zilsel's formula for the Lagrangian for pure liquid He-4. An investigation of properties of the intrinsic internal energy leads to the introduction of intrinsic chemical potentials along with other intrinsic variables for the mixtures. Explicit formulas for these variables are derived for a noninteracting elementary excitation model of the fluid. Using these formulas and others also derived from quantum statistical mechanics, another equivalent expression for the Lagrangian is generated.

Jackson, H. W.↗
NASA NTRS1972

A new theory of gravity

A new relativistic theory of gravity is presented. This theory agrees with all experiments to date. It is a metric theory, it is Lagrangian-based, and it possesses a preferred frame with conformally-flat space slices. With an appropriate choice of certain adjustable functions and parameters, this theory possesses precisely the same post-Newtonian limit as general relativity.

Ni, W.↗
NASA NTRS1980

Solutions of contact problems by the assumed stress hybrid model

A method was developed for contact problems which may be either frictional or frictionless and may involve extensive sliding between deformable bodies. It was based on an assumed stress hybrid approach and on an incremental variational principle for which the Euler's equations of the functional include the equilibrium and compatibility conditions at the contact surface. The tractions at an assumed contact surface were introduced as Lagrangian multipliers in the formulation. It was concluded from the results of several example solutions that the extensive sliding contact between deformable bodies can be solved by the present method.

Kubomura, K.↗
NASA NTRS1982

Long Periodic Terms in the Solar System

The long period variations of the first eight planets in the solar system are studied. First, the Lagrangian solution is calculated and then the long period terms with fourth order eccentricities and inclinations are introduced into the perturbation function. A second approximation was made taking into account the short period terms' contribution, namely the perturbations of first order with respect to the masses. Special attention was paid to the determination of the integration constants. The relative importance of the different contributions is shown. It is useless, for example, to introduce the long period terms of fifth order if no account has been taken of the short period terms. Meanwhile, the terms that have been neglected would not introduce large changes in the integration constants. Even so, the calculation should be repeated with higher order short period terms and fifth order long periods.

Bretagnon, P.↗
NASA NTRS1994

Stably stratified shear turbulence: A new model for the energy dissipation length scale

A model is presented to compute the turbulent kinetic energy dissipation length scale l(sub epsilon) in a stably stratified shear flow. The expression for l(sub epsilon) is derived from solving the spectral balance equation for the turbulent kinetic energy. The buoyancy spectrum entering such equation is constructed using a Lagrangian timescale with modifications due to stratification. The final result for l(sub epsilon) is given in algebraic form as a function of the Froude number Fr and the flux Richardson number R(sub f), l(sub epsilon) = l(sub epsilon)(Fr, R(sub f). The model predicts that for R(sub f) less than R(sub fc), l(sub epsilon) decreases with stratification. An attractive feature of the present model is that it encompasses, as special cases, some seemingly different models for l(sub epsilon) that have been proposed in the past by Deardorff, Hunt et al., Weinstock, and Canuto and Minotti. An alternative form for the dissipation rate epsilon is also discussed that may be useful when one uses a prognostic equation for the heat flux. The present model is applicable to subgrid-scale models, which are needed in large eddy simulations (LES), as well as to ensemble average models. The model is applied to predict the variation of l(sub epsilon) with height z in the planetary boundary layer. The resulting l(sub epsilon) versus z profile reproduces very closely the nonmonotonic profile of l(sub epsilon) exhibited by many LES calculations, beginning with the one by Deardorff in 1974.

Cheng, Y.↗
DOE OSTI2022

Modelling nearest neighbour distributions of biased tracers using hybrid effective field theory

ABSTRACT We investigate the application of hybrid effective field theory (HEFT) – which combines a Lagrangian bias expansion with subsequent particle dynamics from N-body simulations – to the modelling of k-nearest neighbour cumulative distribution functions (kNN-CDFs) of biased tracers of the cosmological matter field. The kNN-CDFs are sensitive to all higher order connected N-point functions in the data, but are computationally cheap to compute. We develop the formalism to predict the kNN-CDFs of discrete tracers of a continuous field from the statistics of the continuous field itself. Using this formalism, we demonstrate how kNN-CDF statistics of a set of biased tracers, such as haloes or galaxies, of the cosmological matter field can be modelled given a set of low-redshift HEFT component fields and bias parameter values. These are the same ingredients needed to predict the two-point clustering. For a specific sample of haloes, we show that both the two-point clustering and the kNN-CDFs can be well-fit on quasi-linear scales (≳ 20h−1Mpc) by the second-order HEFT formalism with the same values of the bias parameters, implying that joint modelling of the two is possible. Finally, using a Fisher matrix analysis, we show that including kNN-CDF measurements over the range of allowed scales in the HEFT framework can improve the constraints on σ8 by roughly a factor of 3, compared to the case where only two-point measurements are considered. Combining the statistical power of kNN measurements with the modelling power of HEFT, therefore, represents an exciting prospect for extracting greater information from small-scale cosmological clustering.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

THE STRUCTURE FUNCTION OF THE FREE NEUTRON AT HIGH X-BJORKEN

Understanding the internal structure of nucleons is one of the primary goal of nuclear physicists. As protons and neutrons are only the bound state solution of the QCD lagrangian (at least inside atomic nuclei), studying protons and neutrons helps uncover nuclear struc ture. Due to its easy availability, many studies on protons have been done on a wide range of kinematics. However, free neutron targets are not readily achievable. So, any information on neutrons has to be extracted from neutron-rich nuclei, and some nuclear models have to be used to subtract the contributions from other nucleons to extract the information on neutrons. So, the Barely Off-shell Nucleon Structure (BONuS12) experiment at Jefferson Lab was conducted to overcome these challenges by using spectator tagging. The experiment effectively created a quasi-free neutron target by scattering electrons off a deuterium target and detecting low-momentum, backward-moving protons using a custom-built Radial Time Projection Chamber (RTPC). Selecting the low momentum and backward-moving spectators would enable us to minimize the model-dependent effects due to final state interactions and target fragmentation. The RTPC was a 40 cm-long cylindrical detector that works on the principle of gaseous ionization. It had three layers of Gas Electron Multipliers (GEMs) for charge amplification and a surrounding readout pad. The scattered electrons were measured using the CLAS12 detector, and data were collected using a 10.4 GeV electron beam dur ing Spring and Summer 2020. Using spectator tagging, we extracted the structure function ratio Fn 2 of the quasi-free neutron in the deep inelastic scattering at high x, upto x ~ 0.8. The result was extracted in the region with the invariant mass W > 1.8 GeV/c2, and Q2 in the range 1.3 to 11 GeV2. This dissertation presents the methodology, event selection criteria and refinements, estimation and subtraction of backgrounds, and complete analysis of extraction of Fn 2/Fp 2 in a model-independent way. Also, systematic uncertainties in our final analysis will be discussed in detail.

Pokhrel, Madhusudhan [Old Dominion Univ., Norfolk,↗