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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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 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,↗
NASA NTRS1972

Spherical balloon response to three-dimensional time-dependent flows

The concept of the Lagrangian displacement of a balloon is introduced. It is shown that the general balloon response problem is extremely complicated because the wind-forcing functions in the balloon equations of motion are functions of the wind velocity vector and its Eulerian first derivatives evaluated at the location of the balloon. The linear perturbation equations for a spherical balloon are derived by perturbing the components of velocity of the balloon about a terminal velocity state which is in equilibrium with a space-time invariant mean horizontal flow. The atmospheric flow is also perturbed such that the resulting equations can be used to analyze the responses of spherical balloons to three-dimensional time-dependent flows. The wind field is represented in terms of a four-fold Fourier integral that involves three orthogonal wave numbers and a frequency, while the balloon components of velocity are represented as Fourier integrals involving a frequency which, in turn, is a function of the wind field wave numbers and frequency and the unperturbed flow components of velocity.

Fichtl, G. H.↗
NASA NTRS2013

Risk-Constrained Dynamic Programming for Optimal Mars Entry, Descent, and Landing

A chance-constrained dynamic programming algorithm was developed that is capable of making optimal sequential decisions within a user-specified risk bound. This work handles stochastic uncertainties over multiple stages in the CEMAT (Combined EDL-Mobility Analyses Tool) framework. It was demonstrated by a simulation of Mars entry, descent, and landing (EDL) using real landscape data obtained from the Mars Reconnaissance Orbiter. Although standard dynamic programming (DP) provides a general framework for optimal sequential decisionmaking under uncertainty, it typically achieves risk aversion by imposing an arbitrary penalty on failure states. Such a penalty-based approach cannot explicitly bound the probability of mission failure. A key idea behind the new approach is called risk allocation, which decomposes a joint chance constraint into a set of individual chance constraints and distributes risk over them. The joint chance constraint was reformulated into a constraint on an expectation over a sum of an indicator function, which can be incorporated into the cost function by dualizing the optimization problem. As a result, the chance-constraint optimization problem can be turned into an unconstrained optimization over a Lagrangian, which can be solved efficiently using a standard DP approach.

Ono, Masahiro↗
NASA NTRS1988

Predicting the velocity and azimuth of fragments generated by the range destruction or random failure of rocket casings and tankage

A calculational method is described which provides a powerful tool for predicting solid rocket motor (SRM) casing and liquid rocket tankage fragmentation response. The approach properly partitions the available impulse to each major system-mass component. It uses the Pisces code developed by Physics International to couple the forces generated by an Eulerian-modeled gas flow field to a Lagrangian-modeled fuel and casing system. The details of the predictive analytical modeling process and the development of normalized relations for momentum partition as a function of SRM burn time and initial geometry are discussed. Methods for applying similar modeling techniques to liquid-tankage-overpressure failures are also discussed. Good agreement between predictions and observations are obtained for five specific events.

Eck, Marshall↗
NASA NTRS1993

Solution of free-boundary problems using finite-element/Newton methods and locally refined grids - Application to analysis of solidification microstructure

A new method is presented for the solution of free-boundary problems using Lagrangian finite element approximations defined on locally refined grids. The formulation allows for direct transition from coarse to fine grids without introducing non-conforming basis functions. The calculation of elemental stiffness matrices and residual vectors are unaffected by changes in the refinement level, which are accounted for in the loading of elemental data to the global stiffness matrix and residual vector. This technique for local mesh refinement is combined with recently developed mapping methods and Newton's method to form an efficient algorithm for the solution of free-boundary problems, as demonstrated here by sample calculations of cellular interfacial microstructure during directional solidification of a binary alloy.

Tsiveriotis, K.↗
NASA NTRS2017

Osiris-REx Spacecraft Current Status and Forward Plans

The NASA New Frontiers OSIRIS-REx spacecraft executed a flawless launch on September 8, 2016 to begin its 23-month journey to near-Earth asteroid (101955). The primary objective of the OSIRIS-REx mission is to collect and return to Earth a pristine sample of regolith from the asteroid surface. The sampling event will occur after a two-year period of remote sensing that will ensure a high probability of successful sampling of a region on the asteroid surface having high science value and within well-defined geological context. The OSIRIS-REx instrument payload includes three high-resolution cameras (OCAMS), a visible and near-infrared spectrometer (OVIRS), a thermal imaging spectrometer (OTES), an X-ray imaging spectrometer (REXIS), and a laser altimeter (OLA). As the spacecraft follows its nominal outbound-cruise trajectory, the propulsion, power, communications, and science instruments have undergone basic functional tests, with no major issues. Outbound cruise science investigations include a search for Earth Trojan asteroids as the spacecraft approaches the Sun-Earth L4 Lagrangian point in February 2017. Additional instrument checkouts and calibrations will be carried out during the Earth gravity assist maneuver in September 2017. During the Earth-moon flyby, visual and spectral images will be acquired to validate instrument command sequences planned for Bennu remote sensing. The asteroid Bennu remote sensing campaign will yield high resolution maps of the temperature and thermal inertia, distributions of major minerals and concentrations of organic matter across the asteroid surface. A high resolution 3d shape model including local surface slopes and a high-resolution gravity field will also be determined. Together, these data will be used to generate four separate maps that will be used to select the sampling site(s). The Safety map will identify hazardous and safe operational regions on the asteroid surface. The Deliverability map will quantify the accuracy with which the navigation team can deliver the spacecraft to and from specific sites on the asteroid surface. The Sampleability map quantifies the regolith properties, providing an estimation of how much material would be sampled at different points on the surface. The final Science Value map synthesizes the chemical, mineralogical, and geological, observations to identify the areas of the asteroid surface with the highest science value. Here, priority is given to organic, water-rich regions that have been minimally altered by surface processes. Asteroid surface samples will be acquired with a touch-and-go sample acquisition system (TAGSAM) that uses high purity pressurized N2 gas to mobilize regolith into a stainless steel canister. Although the mission requirement is to collect at least 60 g of material, tests of the TAGSAM routinely exceeded 300 g of simulant in micro-gravity tests. After acquiring the sample, the spacecraft will depart Bennu in 2021 to begin its return journey, with the sample return capsule landing at the Utah Test and Training Range on September 23, 2023. The OSIRIS-REx science team will carry out a series of detailed chemical, mineralogical, isotopic, and spectral studies that will be used to determine the origin and history of Bennu and to relate high spatial resolution sample studies to the global geological context from remote sensing. The outline of the sample analysis plan is described in a companion abstract.

Messenger, Scott↗
NASA NTRS2021

Air Pollution Scenario over Pakistan: Characterization and Ranking of Extremely Polluted Cities using Long-Term Concentrations of Aerosols and Trace Gases

Pakistan ranks third in the world in terms of mortality attributable to air pollution, with aerosol mass concentrations (PM2.5) consistently well above WHO (World Health Organization) air quality guidelines (AQG). However, regulation is dependent on a sparse network of air quality monitoring stations and insufficient ground data. This study utilizes long-term observations of aerosols and trace gases to characterize and rank the air pollution scenarios and pollution characteristics of 80 selected cities in Pakistan. Datasets used include (1) the Aqua and Terra (AquaTerra) MODIS (Moderate Resolution Imaging Spectroradiometer) Level 2 Collection 6.1 merged Dark Target and Deep Blue (DTB) aerosol optical depth (AOD) retrieval products; (2) the CAMS (Copernicus Atmosphere Monitoring Service) reanalysis PM1, PM2.5, and PM10 data; (3) the MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, Version 2) reanalysis PM2.5 data, (4) the OMI (Ozone Monitoring Instrument) tropospheric vertical column density (TVCD) of nitrogen dioxide (NO2), and VCD of sulfur dioxide (SO2) in the Planetary Boundary Layer (PBL), (5) the VIIRS (Visible Infrared Imaging Radiometer Suite) Nighttime Lights data, (6) MODIS Collection 6 Version 2 global monthly fire location data (MCD14ML), (7) population density, (8) MODIS Level 3 Collection 6 land cover types, (9) AERONET (AErosol RObotic NETwork) Version 3 Level 2.0 data, and (10) ground-based PM2.5 concentrations from air quality monitoring stations. Potential Source Contribution Function (PSCF) analyses were performed by integrating with ground-based PM2.5 concentrations and the NOAA (National Oceanic and Atmospheric Administration) HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) air parcel back trajectories to identify potential pollution source areas which are responsible for extreme air pollution in Pakistan. Results show that the ranking of the top polluted cities depends on the type of pollutant considered and the metric used. For example, Jhang, Multan, and Vehari were characterized as the top three polluted cities in Pakistan when considering AquaTerra DTB AOD products; for PM1, PM2.5, and PM10 Lahore, Gujranwala, and Okara were the top three; for tropospheric NO2 VCD Lahore, Rawalpindi, and Islamabad and for PBL SO2 VCD Lahore, Mirpur, and Gujranwala. The results demonstrate that Pakistan’s entire population has been exposed to high PM2.5 concentrations for many years, with a mean annual value of 54.7 μg/cu. m, over all Pakistan from 2003 to 2020. This value exceeds Pakistan’s National Environmental Quality Standards (Pak-NEQS, i.e., <15 μg/cu. m annual mean) for ambient air defined by the Pakistan Environmental Protection Agency (Pak-EPA) as well as the WHO Interim Target-1 (i.e., mean annual PM2.5 <35 μg/cu. m). The spatial analyses of the concentrations of aerosols and trace gases in terms of population density, nighttime lights, land cover types, and fire location data, and the PSCF analysis indicate that Pakistan’s air quality is strongly affected by anthropogenic sources inside of Pakistan, with contributions from surrounding countries. Statistically significant positive (increasing) trends in PM1, PM2.5, PM10, tropospheric NO2 VCD, and SO2 VCD were observed in ~89%, ~67%, ~48%, 91%, and ~88% of the Pakistani cities (80 cities), respectively. This comprehensive analysis of aerosol and trace gas levels, their characteristics in spatio-temporal domains, and their trends over Pakistan, is the first of its kind. Results will be helpful to the Ministry of Climate Change (Government of Pakistan), Pak-EPA, SUPARCO (Pakistan Space and Upper Atmosphere Research Commission), policymakers, and the local research community to mitigate air pollution and its effects on human health.

Muhammad Bilal↗
NASA NTRS1992

Adaptation of multidimensional group particle tracking and particle wall-boundary condition model to the FDNS code

A particulate two-phase flow CFD model was developed based on the FDNS code which is a pressure based predictor plus multi-corrector Navier-Stokes flow solver. Turbulence models with compressibility correction and the wall function models were employed as submodels. A finite-rate chemistry model was used for reacting flow simulation. For particulate two-phase flow simulations, a Eulerian-Lagrangian solution method using an efficient implicit particle trajectory integration scheme was developed in this study. Effects of particle-gas reaction and particle size change to agglomeration or fragmentation were not considered in this investigation. At the onset of the present study, a two-dimensional version of FDNS which had been modified to treat Lagrangian tracking of particles (FDNS-2DEL) had already been written and was operational. The FDNS-2DEL code was too slow for practical use, mainly because it had not been written in a form amenable to vectorization on the Cray, nor was the full three-dimensional form of FDNS utilized. The specific objective of this study was to reorder to calculations into long single arrays for automatic vectorization on the Cray and to implement the full three-dimensional version of FDNS to produce the FDNS-3DEL code. Since the FDNS-2DEL code was slow, a very limited number of test cases had been run with it. This study was also intended to increase the number of cases simulated to verify and improve, as necessary, the particle tracking methodology coded in FDNS.

Chen, Y. S.↗
DOE OSTI2024

A magnetic analog of pressure–strain interaction

We study the evolution equation for magnetic energy density for a non-relativistic magnetized plasma in the (Lagrangian) reference frame comoving with the electron bulk velocity. Analyzing the terms that arise due to the ideal electric field, namely, perpendicular electron compression and magnetic field line bending, we recast them to reveal a quantity with a functional form analogous to the often-studied pressure–strain interaction term that describes one piece of internal energy density evolution of the species in a plasma, except with the species pressure tensor replaced by the magnetic stress tensor. We dub it the “magnetic stress–strain interaction.” We discuss decompositions of the magnetic stress–strain interaction analogous to those used for pressure–strain interaction. These analogies facilitate the interpretation of the evolution of the various forms of energy in magnetized plasmas and should be useful for a wide array of applications, including magnetic reconnection, turbulence, collisionless shocks, and wave–particle interactions. We display and analyze all the terms that can change magnetic energy density in the Lagrangian reference frame of the electrons using a particle-in-cell simulation of magnetic reconnection.

Barbhuiya, M. Hasan (ORCID:0000000163301650)↗
NASA NTRS1988

An investigation of new methods for estimating parameter sensitivities

Parameter sensitivity is defined as the estimation of changes in the modeling functions and the design variables due to small changes in the fixed parameters of the formulation. There are currently several methods for estimating parameter sensitivities requiring either difficult to obtain second order information, or do not return reliable estimates for the derivatives. Additionally, all the methods assume that the set of active constraints does not change in a neighborhood of the estimation point. If the active set does in fact change, than any extrapolations based on these derivatives may be in error. The objective here is to investigate more efficient new methods for estimating parameter sensitivities when the active set changes. The new method is based on the recursive quadratic programming (RQP) method and in conjunction a differencing formula to produce estimates of the sensitivities. This is compared to existing methods and is shown to be very competitive in terms of the number of function evaluations required. In terms of accuracy, the method is shown to be equivalent to a modified version of the Kuhn-Tucker method, where the Hessian of the Lagrangian is estimated using the BFS method employed by the RPQ algorithm. Inital testing on a test set with known sensitivities demonstrates that the method can accurately calculate the parameter sensitivity. To handle changes in the active set, a deflection algorithm is proposed for those cases where the new set of active constraints remains linearly independent. For those cases where dependencies occur, a directional derivative is proposed. A few simple examples are included for the algorithm, but extensive testing has not yet been performed.

Beltracchi, Todd J.↗
NASA NTRS1989

An investigation of new methods for estimating parameter sensitivities

The method proposed for estimating sensitivity derivatives is based on the Recursive Quadratic Programming (RQP) method and in conjunction a differencing formula to produce estimates of the sensitivities. This method is compared to existing methods and is shown to be very competitive in terms of the number of function evaluations required. In terms of accuracy, the method is shown to be equivalent to a modified version of the Kuhn-Tucker method, where the Hessian of the Lagrangian is estimated using the BFS method employed by the RQP algorithm. Initial testing on a test set with known sensitivities demonstrates that the method can accurately calculate the parameter sensitivity.

Beltracchi, Todd J.↗