Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian modeling”

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

The evolution of voids in the adhesion approximation

We apply the adhesion approximation to study the formation and evolution of voids in the universe. Our simulations-carried out using 128(exp 3) particles in a cubical box with side 128 Mpc-indicate that the void spectrum evolves with time and that the mean void size in the standard Cosmic Background Explorer Satellite (COBE)-normalized cold dark matter (CDM) model with H(sub 50) = 1 scals approximately as bar D(z) = bar D(sub zero)/(1+2)(exp 1/2), where bar D(sub zero) approximately = 10.5 Mpc. Interestingly, we find a strong correlation between the sizes of voids and the value of the primordial gravitational potential at void centers. This observation could in principle, pave the way toward reconstructing the form of the primordialpotential from a knowledge of the observed void spectrum. Studying the void spectrum at different cosmological epochs, for spectra with a built in k-space cutoff we find that the number of voids in a representative volume evolves with time. The mean number of voids first increases until a maximum value is reached (indicating that the formation of cellular structure is complete), and then begins to decrease as clumps and filaments erge leading to hierarchical clustering and the subsequent elimination of small voids. The cosmological epoch characterizing the completion of cellular structure occurs when the length scale going nonlinear approaches the mean distance between peaks of the gravitaional potential. A central result of this paper is that voids can be populated by substructure such as mini-sheets and filaments, which run through voids. The number of such mini-pancakes that pass through a given void can be measured by the genus characteristic of an individual void which is an indicator of the topology of a given void in intial (Lagrangian) space. Large voids have on an average a larger measure than smaller voids indicating more substructure within larger voids relative to smaller ones. We find that the topology of individual voids is strongly epoch dependent, with void topologies generally simplifying with time. This means that as voids grow older they become progressively more empty and have less structure within them. We evaluate the genus measure both for individual voids as well as for the entire ensemble of voids predicted by CDM model. As a result we find that the topology of voids when taken together with the void spectrum is a very useful statistical indicator of the evolution of the structure of the universe on large scales.

Sahni, Varun↗
DOE OSTI2023

Relativistic Hartree–Fock chiral Lagrangians with confinement, nucleon finite size and short-range effects

A relativistic Hartree–Fock Lagrangian including a chiral potential and nucleon polarisation is investigated in hopes of providing a better description of dense nuclear matter. Here, we fully consider the contribution of the exchange Fock term to the energy and the self-energies, and in addition we investigate the nucleon’s compositeness and finite size effects (confinement and form factors) and short range correlations modeled by a Jastrow ansatz. These effects are added step by step, such that their impact on the dense matter properties can be analysed in details. The parameters of the model are adjusted to reproduce fundamental properties related to the QCD theory at low energy, such as the chiral symmetry breaking, nucleon’s quark substructure and Lattice-QCD predictions, as well as two empirical properties at saturation: the binding energy and the density. All other empirical parameters, e.g., symmetry energy and its slope, incompressibility modulus, effective mass, as well as spin–isospin Landau–Midgal parameter are predictions of the models and can be used to evaluate the gain of the different approximation schemes in describing nuclear properties. Bayesian statistics is employed in order to propagate parameter uncertainties into predictions for the nuclear matter properties. We show that the splitting of the effective Landau mass is largely influenced by the value of the ρ T coupling, and we show that the fit to the symmetry energy, which induces an increase of the coupling constant g ρ by about 20–25% compared to the case where it is fixed by the quark model, provides a very good EoS compatible with the present nuclear physics knowledge.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2000

A Mixed Finite Analysis of Axial Foreshortening

In previous studies, a mixed finite element was derived and shown to be effective in facilitating accurate modal reduction of rotor blades. This study builds upon that earlier work through the development of algorithms that allow the element's Lagrangian axial displacement degrees of freedom to be eliminated in favor of the axial force degrees of freedom. This reduction process, which had been demonstrated previously for the case of a single mixed finite element, may be viewed as a generalization of the UMARC blade analysis methodology to arbitrary topologies. A unique strength of the method. developed here is that the model may consist of two dimensional or even three dimensional elements.

Ruzicka, Gene↗
NASA NTRS2011

A New Limit on Planck Scale Lorentz Violation from Gamma-ray Burst Polarization

Constraints on possible Lorentz invariance violation (UV) to first order in E/M(sub Plank) for photons in the framework of effective field theory (EFT) are discussed, taking cosmological factors into account. Then. using the reported detection of polarized soft gamma-ray emission from the gamma-ray burst GRB041219a that is indicative' of an absence of vacuum birefringence, together with a very recent improved method for estimating the redshift of the burst, we derive constraints on the dimension 5 Lorentz violating modification to the Lagrangian of an effective local QFT for QED. Our new constraints are more than five orders of magnitude better than recent constraints from observations of the Crab Nebula.. We obtain the upper limit on the Lorentz violating dimension 5 EFT parameter absolute value of zeta of 2.4 x 10(exp -15), corresponding to a constraint on the dimension 5 standard model extension parameter. Kappa (sup 5) (sub (v)oo) much less than 4.2 X 10(exp -3)4 / GeV.

Stecker, Floyd W.↗
DOE OSTI2024

Chiral limit of 2d QCD revisited with lightcone conformal truncation

We study the chiral limit of 2d QCD with a single quark flavor at finite Nc using LCT. By modifying the LCT basis according to the quark mass in a manner motivated by ’t Hooft’s analysis, we are able to restore convergence for quark masses much smaller than the QCD strong coupling scale. For such small quark masses, the IR of the theory is expected to be well described by the Sine-Gordon model. We verify that LCT numerics are able to capture in detail the spectrum and correlation functions of the Sine-Gordon model. This opens up the possibility for studying deformations of various integrable CFTs using LCT by considering the chiral limit of QCD like theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1989

On the stream-accretion disk interaction - Response to increased mass transfer rate

The time-dependent interaction between the stream of mass from the inner Lagrangian point and the accretion disk, resulting from an increasing mass transfer rate is calculated. The calculation is fully three-dimensional, using a pseudoparticle description of the hydrodynamics. It is demonstrated that the results of such calculations, when combined with specific observations, have the potential of both determining essential parameters, such as the viscosity parameter alpha, and can distinguish between different models of dwarf nova eruptions.

Dgani, Ruth↗
DOE OSTI2020

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Position-dependent Voronoi probability distribution functions for matter and halos

Here, we measure the Voronoi density probability distribution function (PDF) for both dark matter and halos in N-body simulations. For the dark matter, Voronoi densities represent the matter density field smoothed on a uniform mass scale, which approximates the Lagrangian density field. For halos, the Voronoi densities contain information about the local environment of each halo. We measure the halo virial masses, the total amount of dark matter within each halo Voronoi cell, and the halo Voronoi cell volumes, and we show how halo abundances depend on these three quantities. We then study the position-dependent Voronoi density PDF, measured within finite subregions of the Universe, using separate universe simulations. We demonstrate that the spatial variation of the position-dependent PDF is due to large-scale density fluctuations, indicating that the position-dependent PDF is a biased tracer of large-scale structure. We measure this bias for the dark matter, and interpret it as the bias of regions of the Lagrangian density field that are selected based on density. For the halos, this bias can be interpreted as a form of assembly bias. We present the mapping from late-time to early-time Voronoi density for each simulation dark matter particle, which is highly stochastic. We compare the median of this stochastic map with spherical collapse calculations and discuss challenges involved in modeling the evolution of the density field on these scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

A hydrodynamic study of a slow nova outburst

The paper reports use of a Lagrangian implicit hydrodynamics computer code incorporating a full nuclear-reaction network to follow a thermonuclear runaway in the hydrogen-rich envelope of a 1.25 solar-mass white dwarf. In this evolutionary sequence the envelope was assumed to be of normal (solar) composition and the resulting outburst closely resembles that of the slow nova HR Del. In contrast, previous CNO-enhanced models resemble fast nova outbursts. The slow-nova model ejects material by radiation pressure when the high luminosity of the rekindled hydrogen shell source exceeds the local Eddington luminosity of the outer layers. This is in contrast to the fast nova outburst where ejection is caused by the decay of the beta(+)-unstable nuclei. Nevertheless, radiation pressure probably plays a major role in ejecting material from the fast nova remnants. Therefore, the sequence from slow to fast novae can be interpreted as a sequence of white dwarfs with increasing amounts of enhanced CNO nuclei in their hydrogen envelopes, although other parameters such as the white-dwarf mass and accretion rate probably contribute to the observed variation between novae.

Sparks, W. M.↗
NASA NTRS1999

A Priori Subgrid Analysis of Temporal Mixing Layers with Evaporating Droplets

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5. In the DNS, the gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. The Large Eddy Simulation (LES) equations corresponding to the DNS are first derived, and key assumptions in deriving them are first confirmed by computing the terms using the DNS database. Since LES of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be the sum of the filtered variables and a correction based on the filtered standard deviation; this correction is then computed from the Subgrid Scale (SGS) standard deviation. This model predicts the unfiltered variables at droplet locations considerably better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: the Smagorinsky approach, the Gradient model and the Scale-Similarity formulation. When the proportionality constant inherent in the SGS models is properly calculated, the Gradient and Scale-Similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗
NASA NTRS1989

An investigation of using an RQP based method to calculate parameter sensitivity derivatives

Estimation of the sensitivity of problem functions with respect to problem variables forms the basis for many of our modern day algorithms for engineering optimization. The most common application of problem sensitivities has been in the calculation of objective function and constraint partial derivatives for determining search directions and optimality conditions. A second form of sensitivity analysis, parameter sensitivity, has also become an important topic in recent years. By parameter sensitivity, researchers refer to the estimation of changes in the modeling functions and current design point due to small changes in the fixed parameters of the formulation. Methods for calculating these derivatives have been proposed by several authors (Armacost and Fiacco 1974, Sobieski et al 1981, Schmit and Chang 1984, and Vanderplaats and Yoshida 1985). Two drawbacks to estimating parameter sensitivities by current methods have been: (1) the need for second order information about the Lagrangian at the current point, and (2) the estimates assume no change in the active set of constraints. The first of these two problems is addressed here and a new algorithm is proposed that does not require explicit calculation of second order information.

Beltracchi, Todd J.↗
NASA NTRS2015

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗
DOE OSTI2021

A Lagrangian perspective on tropical anvil cloud lifecycle in present and future climate

The evolution of tropical anvil clouds from their origin in deep convective cores to their slow decay determines the climatic effects of clouds in tropical convective regions. Despite the relevance of anvil clouds for climate and responses of clouds to global warming, processes dominating their evolution are not well understood. Currently available observational data reveal instantaneous snapshots of anvil cloud properties, but cannot provide a process-based perspective on anvil evolution. We therefore conduct simulations with the high resolution version of the Exascale Earth System Model in which we track mesoscale convective systems over the Tropical Western Pacific and compute trajectories that follow air parcels detrained from peaks of convective activity. With this approach we gain new insight into the anvil cloud evolution both in present day and future climate. Comparison with geostationary satellite data shows that the model is able to simulate maritime mesoscale convective systems reasonably well. Trajectory results indicate that anvil cloud lifetime is about 15 hours with no significant change in a warmer climate. The anvil ice mixing ratio is larger in a warmer climate due to a larger source of ice by detrainment and larger depositional growth leading to a more negative net cloud radiative effect along detrained trajectories. However, the increases in sources are counteracted by increases in sinks of ice, particularly snow formation and sedimentation. Furthermore, we find that the mean anvil cloud feedback along trajectories is positive and consistent with results from more traditional cloud feedback calculation methods.

climate, weather, mesocale cellular convection↗
DOE OSTI2022

Lagrangian particle simulation of hydrogen pellets and SPI into runaway electron beam in ITER

Numerical studies of the ablation of pellets and shattered pellet injection (SPI) fragments into a runaway electron beam in ITER have been performed using a time-dependent pellet ablation code [Samulyak et al., Nucl. Fusion, 61(4), 046007 (2021)]. The code resolves detailed ablation physics near pellet fragments and large-scale expansion of ablated clouds. The study of a single-fragment ablation quantifies the influence of various factors, in particular, the impact ionization by runaway electrons and cross-field transport models, on the dynamics of ablated plasma and its penetration into the runaway beam. Simulations of SPI performed using different numbers of pellet fragments study the formation and evolution of the ablation clouds and their large-scale dynamics in ITER. In conclusion, the penetration depth of the ablation clouds is found to be of the order of 50 cm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Planetary rings - 2-2/3 centuries of nearly total ignorance, 4 years of information explosion

Saturn lies at nearly twice Jupiter's distance from the sun and nearly all parts of its system are characterized by much smaller scales than those which are important in the case of Jupiter. This appears in the structures of the planet's atmosphere, in the sizes of classical satellites other than Titan vis-a-vis those of the Galilean satellites, in the plethora of small Saturnian satellites, especially Lagrangian co-orbiters, in the structure of Saturn's F-Ring as contrasted with that of Jupiter's Ring and finally in the highly structured detail in Saturn's Rings, much finer than seriously considered in past theoretical discussions. Uranus' Rings were unknown until five years ago. The discovery and observation of these rings have revived contributions to theory originally intended for application to Saturn's Rings. Models have also been generated for eccentric rings for application to Uranus' Rings which also apply to those of Saturn. These two classes of model are reviewed in the present paper along with the first tentative steps made down the road to unravelling the complexity of Saturn's Rings.

Cook, A. F.↗
DOE OSTI2022

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

A Mixed Integer Linear Programming-basedDistributed Energy Management for Three-phaseUnbalanced Active Distribution Network

A mixed integer linear programming (MILP)–baseddistributed energy management for three-phase unbalancedactive distribution network is proposed. Modern distributionnetworks have becoming more and more active with increasingdeployment of microgrids, distributed energy resources (DERs)as well as controllable loads. Considering various ownership andcontrol models of microgrids, DERs and controllable loads, adistributed energy management was formulated using the alternatingdirection method of multipliers (ADMM) algorithm. ByADMM, the distribution management system (DMS) and theseactive components are coordinated through price signals, whichare adjusted according to the generation-load mismatch per nodeper phase. To enable resolution of the ADMM-based distributedoptimization using more accessible and popular MILP solver,different linearization techniques were proposed to linearize theaugmented Lagrangian terms and other nonlinear terms. Resultsof case studies on a three-phase active distribution network withthree microgrids and several DERs and controllable loads validatedthe effectiveness of proposed MILP-based distributed energymanagement. In addition, the capability of proposed method inmitigating phase power unbalance has been demonstrated.

Liu, Guodong↗
DOE OSTI2022

An investigation of shock-induced phase transition in soda-lime glass

There exists a large body of evidence from experiments and molecular dynamics simulations to suggest the occurrence of phase transitions in soda-lime glass (SLG) and other silica glasses subject to shock compression to pressures above 3 GPa. In light of these findings, the current work investigated the existence of phase transition in SLG using shock and release experiments. The experiments employed symmetric SLG–SLG impact to achieve complete unloading to zero stress after shock compression to stresses in the range of 3–7 GPa. The stress–strain response and the Lagrangian release wave speed behavior of SLG obtained from these experiments are seen to reveal a mismatch between the loading and unloading paths of the pressure–strain curve for the material, which serves as compelling evidence for the occurrence of a shock-induced phase transition in the material at relatively low pressures. Furthermore, the release wave speed vs strain data obtained from experiments were used to construct a methodology for modeling the shock and release behavior of SLG. Lastly, this scheme implemented in numerical simulations was able to capture the release behavior of shock compressed SLG, for which a robust and satisfactory model was previously unavailable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗