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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

The nature of Trojan asteroid 624 Hektor

Near-simultaneous visual and thermal IR (20-micron) photometry of the Trojan asteroid 624 Hektor is reported which was performed when the asteroid was observed nearly along its rotation axis. The results confirm and refine the low albedo and large size of this asteroid and confirm the general rotational-pole position and aspect angle predicted by Dunlap and Gehrels (1969). Hektor is found to be a truly extraordinary object in that it is larger and far more irregular in shape than other measured Trojans and far more irregular than other belt asteroids of comparable size. It is proposed that Hektor could be a partially coalesced pair of Trojan asteroids which collided with energy too low to cause complete fragmentation, thus forming a dumbbell-shaped object. A possible scenario is outlined according to which the two pre-Hektor objects were neighboring relatively large primitive spheroidal planetesimals trapped in Jupiter's Lagrangian cloud. Observational and theoretical tests of this model are suggested.

Hartmann, W. K.↗
NASA NTRS1982

Three-dimensional relativistic field-electron interaction in a multicavity high-power klystron. 1: Basic theory

A theoretical investigation of three dimensional relativistic klystron action is described. The relativistic axisymmetric equations of motion are derived from the time-dependent Lagrangian function for a charged particle in electromagnetic fields. An analytical expression of the fringing RF electric and magnetic fields within and in the vicinity of the interaction gap and the space-charge forces between axially and radially elastic deformable rings of charges are both included in the formulation. This makes an accurate computation of electron motion through the tunnel of the cavities and the drift tube spaces possible. Method of analysis is based on Lagrangian formulation. Bunching is computed using a disk model of electron stream in which the electron stream is divided into axisymmetric disks of equal charge and each disk is assumed to consist of a number of concentric rings of equal charges. The Individual representative groups of electrons are followed through the interaction gaps and drift tube spaces. Induced currents and voltages in interacting cavities are calculated by invoking the Shockley-Ramo theorem.

Kosmahl, H. G.↗
NASA NTRS1991

Centralized, decentralized, and independent control of a flexible manipulator on a flexible base

The dynamics and control of a flexible manipulator arm with payload mass on a flexible base in space are considered. The controllers are provided by one torquer at the center of the base and one torquer at the connection joint of the robot and the base. The nonlinear dynamics of the system is modeled by applying the finite element method and Lagrangian formula. Three control strategies are considered and compared, i.e., centralized control, decentralized control, and independent control. All these control designs are based on the linear quadratic regulator theory. A mathematical decomposition is used in the decentralization process so that the coupling between the subsystems is weak, while a physical decomposition is used in the independent control design process. For both the decentralized and the independent controls, the stability of the overall linear system is checked before a numerical simulations is initiated. Two numerical examples show that the response of the independent control system are close to those of the centralized control system, while the responses of the decentralized control system are not.

Li, Feiyue↗
NASA NTRS1994

Project JOVE

The goal of this project is to investigate new areas of research pertaining to free surface-interface fluids mechanics and/or microgravity which have potential commercial applications. This paper presents an introduction to ferrohydrodynamics (FHD), and discusses some applications. Also, computational methods for solving free surface flow problems are presented in detail. Both have diverse applications in industry and in microgravity fluids applications. Three different modeling schemes for FHD flows are addressed and the governing equations, including Maxwell's equations, are introduced. In the area of computational modeling of free surface flows, both Eulerian and Lagrangian schemes are discussed. The state of the art in computational methods applied to free surface flows is elucidated. In particular, adaptive grids and re-zoning methods are discussed. Additional research results are addressed and copies of the publications produced under the JOVE Project are included.

Lyell, M. J.↗
NASA NTRS1999

Transport of Passive Tracers in Baroclinic Wave Life Cycles

The transport of passive tracers in idealized baroclinic wave life cycles is studied using output from the National Center for Atmospheric Research Community Climate Model (CCM2). Two life cycles, LCn and LCs, are simulated, starting with baroclinically unstable initial conditions similar to those used by Thorncroft et al. in their study of two life cycle paradigms. The two life cycles LCn and LCs have different initial horizontal wind shear structures that result in distinctive nonlinear development. In terms of potential vorticity-potential temperature (PV-theta) diagnostics, the LCn case is characterized by thinning troughs that are advected anti-cyclonically and equatorward, while the LCs case has broadening troughs that wrap up cyclonically and poleward. Four idealized passive tracers are included in the model to be advected by the semi-Lagrangian transport scheme of the CCM2, and their evolutions are investigated throughout the life cycles. Tracer budgets are analyzed in terms of the transformed Eulerian mean constituent transport formalism in pressure coordinates and also in isentropic coordinates. Results for both LCn and LCs show transport that is downgradient with respect to the background structure of the tracer field, but with a characteristic spatial structure that maximizes in the middle to high latitudes. For the idealized tropospheric tracers in this study, this represents a net upward and poleward transport that enhances concentrations at high latitudes. These results vary little with the initial distribution of the constituent field. The time tendency of the tracer is influenced most strongly by the eddy flux term. with the largest transport occurring during the nonlinear growth stage of the life cycle. The authors also study the transport of a lower-stratospheric tracer, to examine stratosphere-troposphere exchange for baroclinic waves.

Stone, Elizabeth M.↗
DOE OSTI2022

Optimization of artificial viscosity in production codes based on Gaussian Regression surrogate models

To accurately model flows with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity has to be introduced to convert kinetic energy into internal energy, thereby increasing the entropy across shocks. Determining the appropriate strength of the artificial viscosity is an art and strongly depends on the particular problem and experience of the researcher. The objective of this study is to pose the problem of finding the appropriate strength of artificial viscosity as an optimization problem and solve this problem using machine learning (ML) tools, specifically using surrogate models based on Gaussian Process regression and Bayesian analysis. We describe the optimization method and discuss various practical details of its implementation. The shock-containing problems for which we apply this method all have been implemented in the LANL code FLAG. First, we apply ML to find optimal values to isolated shock problems of different strengths. Second, we apply ML to optimize viscosity for a 1D propagating detonation problem based on Zel’dovich-von Neumann-Doring (ZND) detonation theory using a reactive burn model. We compare results for default (currently used values in FLAG) and optimized values of artificial viscosity for these problems demonstrating the potential for significant improvement in the accuracy of computations.

42 ENGINEERING↗
DOE OSTI2021

A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges

Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.

58 GEOSCIENCES↗
NASA NTRS1985

Computing three-dimensional incompressible flows with vortex elements

The techniques, capabilities and applicability of numerical models of three-dimensional, unsteady vortical flows with high Re are assessed. Vorticity is calculated only in appropriate regions and the velocity field is derived from the boundary conditions. Vorticity is assumed to take the shape of tubes with uniform core structures in the case of turbulence. The efforts being made to simplify equations for dense collections of vortex filaments in order to make them tractable to computer simulations are described. The effectiveness of vorticity arrow representations for accurately describing vorticity fields near surfaces is discussed, along with Lagrangian vortex elements, which may be of use in modelling the rotational part of flows around bluff bodies, nonuniform density flows and chemically reacting flows.

Leonard, A.↗
NASA NTRS2016

The Framework for 0-D Atmospheric Modeling (F0AM) v3.1

The Framework for 0-D Atmospheric Modeling(F0AM) is a flexible and user-friendly MATLAB-based platform for simulation of atmospheric chemistry systems. The F0AM interface incorporates front-end configuration of observational constraints and model setups, making it readily adaptable to simulation of photochemical chambers, Lagrangian plumes, and steady-state or time-evolving solar cycles. Six different chemical mechanisms and three options for calculation of photolysis frequencies are currently available. Example simulations are presented to illustrate model capabilities and, more generally, highlight some of the advantages and challenges of 0-D box modeling.

chemistry↗
DOE OSTI2022

Modelling the Lyman-α forest with Eulerian and SPH hydrodynamical methods

ABSTRACT We compare two state-of-the-art numerical codes to study the overall accuracy in modelling the intergalactic medium and reproducing Lyman-α forest observables for DESI and high-resolution data sets. The codes employ different approaches to solving both gravity and modelling the gas hydrodynamics. The first code, Nyx, solves the Poisson equation using the Particle-Mesh (PM) method and the Euler equations using a finite-volume method. The second code, CRK-HACC , uses a Tree-PM method to solve for gravity, and an improved Lagrangian smoothed particle hydrodynamics (SPH) technique, where fluid elements are modelled with particles, to treat the intergalactic gas. We compare the convergence behaviour of the codes in flux statistics as well as the degree to which the codes agree in the converged limit. We find good agreement overall with differences being less than observational uncertainties, and a particularly notable ≲1 per cent agreement in the 1D flux power spectrum. This agreement was achieved by applying a tessellation methodology for reconstructing the density in CRK-HACC instead of using an SPH kernel as is standard practice. We show that use of the SPH kernel can lead to significant and unnecessary biases in flux statistics; this is especially prominent at high redshifts, z ∼ 5, as the Lyman-α forest mostly comes from lower-density regions that are intrinsically poorly sampled by SPH particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Rocking response of liquid-filled cylindrical tanks

The hydrodynamic response of liquid-filled storage tanks subjected to translational earthquake shaking has been the focus of studies for more than 75 years. In service, tanks experience six components of seismic input, with rocking due to support flexibility and/or rotational components of ground motion. Published theory for predicting hydrodynamic responses in rigid cylindrical tanks due to rotational seismic input is examined. Analytical solutions for pressure in the fluid and on the tank wall, and base moment are modified to account for the effects of angular displacements at the base. A finite element model of a water-filled cylindrical tank using the Arbitrary Lagrangian and Eulerian (ALE) solver in LS-DYNA is used to demonstrate the efficacy of the proposed expressions for three base rocking motions. Results from the modified analytical solutions and the finite element model are in excellent agreement.

58 GEOSCIENCES↗
NASA NTRS1980

Instrumented aircraft verification of clear-air radar detection of low-level wind shear

A technique is described in which the longitudinal wind component at airports, measured by Doppler radar, is put through a Fourier transformation to produce an energy density spectrum. The result is multiplied by the aircraft response characteristics to produce quantitative data on aircraft performance in a wind shear situation. The Doppler radar measurements are made along the intended flight path. Two case studies, along a 3 deg approach path, are described from measurements in two modes: Lagrangian, where predictions were computed one range gate ahead of the plane, and Eulerian, where samples were measured instantaneously along the flight path. Wind data obtained every second were interpolated to fit a numerical model by using a cubic spline. Eulerian data were superior to Lagrangian data in terms of ranging, and nearly equal for wind speeds. Improvements in current 75% accuracy are noted to be possible with shrouded antennas to reduce sidelobes and use of Doppler radar with higher resolution to eliminate effects of small scale disturbances within the pulse volume.

Mccarthy, J.↗
NASA NTRS2005

Simulations of Evaporating Multicomponent Fuel Drops

A paper presents additional information on the subject matter of Model of Mixing Layer With Multicomponent Evaporating Drops (NPO-30505), NASA Tech Briefs, Vol. 28, No. 3 (March 2004), page 55. To recapitulate: A mathematical model of a three-dimensional mixing layer laden with evaporating fuel drops composed of many chemical species has been derived. The model is used to perform direct numerical simulations in continuing studies directed toward understanding the behaviors of sprays of liquid petroleum fuels in furnaces, industrial combustors, and engines. The model includes governing equations formulated in an Eulerian and a Lagrangian reference frame for the gas and drops, respectively, and incorporates a concept of continuous thermodynamics, according to which the chemical composition of a fuel is described by use of a distribution function. In this investigation, the distribution function depends solely on the species molar weight. The present paper reiterates the description of the model and discusses further in-depth analysis of the previous results as well as results of additional numerical simulations assessing the effect of the mass loading. The paper reiterates the conclusions reported in the cited previous article, and states some new conclusions. Some new conclusions are: 1. The slower evaporation and the evaporation/ condensation process for multicomponent-fuel drops resulted in a reduced drop-size polydispersity compared to their single-component counterpart. 2. The inhomogeneity in the spatial distribution of the species in the layer increases with the initial mass loading. 3. As evaporation becomes faster, the assumed invariant form of the molecular- weight distribution during evaporation becomes inaccurate.

Bellan, Josette↗
DOE OSTI2022

Dark matter effective field theory and an application to vector dark matter

The Standard Model Effective Field Theory (SMEFT) and the Low Energy Effective Field Theory (LEFT) can be extended by adding additional spin 0, 1/2 and 1 dark matter particles which are singlets under the Standard Model (SM) gauge group. We classify all gauge invariant interactions in the Lagrangian up to terms of dimension six, and present the tree-level matching conditions between the two theories at the electroweak scale. The most widely studied dark matter models, such as those based on the Higgs portal or on kinetic mixing between the photon and a dark photon, are based on dimension-four interactions with the SM sector. We consider a model with dark vector particles with a $\mathbb{Z}_2$ symmetry, so that the lightest dark matter particle is stable. The leading interaction with the SM is through dimension-six operators involving two dark vector field-strength tensors and the electromagnetic field-strength tensor. This model is a viable dark matter model in the freeze-in scenario for a wide range of parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2019

Turbulent Transport in the Gray Zone: A Large Eddy Model Intercomparison Study of the CONSTRAIN Cold Air Outbreak Case

To quantify the turbulent transport at gray zone length scales between 1 and 10 km the Lagrangian evolution of the CONSTRAIN cold air outbreak (CAO) case was simulated with seven large eddy models. The case is characterized by rather large latent and sensible heat fluxes, and a rapid deepening rate of the boundary layer. In some models the entrainment velocity exceeds 4 cm/s. A significant fraction of this growth is attributed to a strong longwave radiative cooling of the inversion layer. The evolution and the timing of the breakup of the stratocumulus cloud deck differ significantly among the models. Sensitivity experiments demonstrate that a decrease in the prescribed cloud droplet number concentration, and the inclusion of ice microphysics, both act to speed up the thinning of the stratocumulus by enhancing the production of precipitation. In all models the formation of mesoscale fluctuations is clearly evident in the cloud fields but also in the horizontal wind velocity. Resolved vertical fluxes remain important for scales up to 10 km. The simulation results show that the resolved vertical velocity variance gradually diminishes with a coarsening of the horizontal mesh, but the total vertical fluxes of heat, moisture, and momentum are only weakly affected. This is a promising result as it demonstrates the potential use of a mesh size dependent turbulent length scale for convective boundary layers at gray zone model resolutions.

gray zone↗
NASA NTRS2016

What If We Had A Magnetograph at Lagrangian L5?

Synoptic Carrington charts of magnetic field are routinely used as an input for modelings of solar wind and other aspects of space weather forecast. However, these maps are constructed using only the observations from the solar hemisphere facing Earth. The evolution of magnetic flux on the "farside" of the Sun, which may affect the topology of coronal field in the "nearside," is largely ignored. It is commonly accepted that placing a magnetograph in Lagrangian L5 point would improve the space weather forecast. However, the quantitative estimates of anticipated improvements have been lacking. We use longitudinal magnetograms from the Synoptic Optical Long-term Investigations of the Sun (SOLIS) to investigate how adding data from L5 point would affect the outcome of two major models used in space weather forecast.

Pevtsov, Alexei A.↗
NASA NTRS1985

Transport characteristics in the middle atmosphere and two-dimensional modeling

Transport mechanisms of tracers in the atmosphere are briefly discussed from a viewpoint of Lagrangian motions of air parcels. The tracers are retained by air parcels themselves. In this regard, both momentum and energy are quite different from chemical tracers. The present problem is what interpretation is best for understanding transport mechanisms of tracers. A method is reviewed and proposed, although conceptually. The numerical model used was a whole hemispheric global climate model. Annual- and zonal-mean diabatic forcing was assumed to drive the general circulation. Time integrations were carried out for about 100 days. The model quasi-equilibrium states were analyzed.

Kida, H.↗