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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS2003

"Hypothetical" Heavy Particles Dynamics in LES of Turbulent Dispersed Two-Phase Channel Flow

The extensive experimental study of dispersed two-phase turbulent flow in a vertical channel has been performed in Eaton's research group in the Mechanical Engineering Department at Stanford University. In Wang & Squires (1996), this study motivated the validation of LES approach with Lagrangian tracking of round particles governed by drag forces. While the computed velocity of the flow have been predicted relatively well, the computed particle velocity differed strongly from the measured one. Using Monte Carlo simulation of inter-particle collisions, the computation of Yamamoto et al. (2001) was specifically performed to model Eaton's experiment. The results of Yamamoto et al. (2001) improved the particle velocity distribution. At the same time, Vance & Squires (2002) mentioned that the stochastic simualtion of inter-particle collisions is too expensive, requiring significantly more CPU resources than one needs for the gas flow computation. Therefore, the need comes to account for the inter-particle collisions in a simpler and still effective way. To present such a model in the framework of LES/Lagrangian particle approach, and to compare the calculated results with Eaton's measurement and modeling of Yamamoto is the main objective of the present paper.

Gorokhovski, M.↗

DOE OSTI2023

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

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

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1988

Advances in contact algorithms and their application to tires

Currently used techniques for tire contact analysis are reviewed. Discussion focuses on the different techniques used in modeling frictional forces and the treatment of contact conditions. A status report is presented on a new computational strategy for the modeling and analysis of tires, including the solution of the contact problem. The key elements of the proposed strategy are: (1) use of semianalytic mixed finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynomials in the meridional direction; (2) use of perturbed Lagrangian formulation for the determination of the contact area and pressure; and (3) application of multilevel iterative procedures and reduction techniques to generate the response of the tire. Numerical results are presented to demonstrate the effectiveness of a proposed procedure for generating the tire response associated with different Fourier harmonics.

Noor, Ahmed K.↗

NASA NTRS2005

Dynamic Deployment Simulations of Inflatable Space Structures

The feasibility of using Control Volume (CV) method and the Arbitrary Lagrangian Eulerian (ALE) method in LSDYNA to simulate the dynamic deployment of inflatable space structures is investigated. The CV and ALE methods were used to predict the inflation deployments of three folded tube configurations. The CV method was found to be a simple and computationally efficient method that may be adequate for modeling slow inflation deployment sine the inertia of the inflation gas can be neglected. The ALE method was found to be very computationally intensive since it involves the solving of three conservative equations of fluid as well as dealing with complex fluid structure interactions.

Wang, John T.↗

NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗

DOE OSTI2023

Classical integrability of root-$T\overline{T}$ flows

The root-$T\overline{T}$ flow was recently introduced as a universal and classically marginal deformation of any two-dimensional translation-invariant field theory. The flow commutes with the (irrelevant) $T\overline{T}$ flow, and it can be integrated explicitly for a large class of actions, leading to nonanalytic Lagrangians reminiscent of the four-dimensional modified-Maxwell theory (ModMax). It is not a priori obvious whether the root-$T\overline{T}$ flow preserves integrability, as is the case for the $T\overline{T}$ flow. In this paper we demonstrate that this is the case for a large class of classical models by explicitly constructing a deformed Lax connection. We discuss the principal chiral model and the nonlinear sigma models on symmetric and semisymmetric spaces, without or with the Wess-Zumino term. We also construct Lax connections for the two-parameter families of theories deformed by both root-$T\overline{T}$ and $T\overline{T}$ for all of these models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2012

The Kinematic and Microphysical Control of Storm Integrated Lightning Flash Extent

The objective of this preliminary study is to investigate the kinematic and microphysical control of lightning properties, particularly those that may govern the production of nitrogen oxides (NOx) in thunderstorms, such as flash rate, type and extent. The mixed-phase region is where the noninductive charging (NIC) process is thought to generate most storm electrification during rebounding collisions between ice particles in the presence of supercooled water. As a result, prior radar-based studies have demonstrated that lightning flash rate is well correlated to kinematic and microphysical properties in the mixed-phase region of thunderstorms such as updraft volume, graupel mass, or ice mass flux. There is also some evidence that lightning type is associated with the convective state. Intracloud (IC) lightning tends to dominate during the updraft accumulation of precipitation ice mass while cloud-to-ground (CG) lightning is more numerous during the downdraft-driven descent of radar echo associated with graupel and hail. More study is required to generalize these relationships, especially regarding lightning type, in a wide variety of storm modes and meteorological conditions. Less is known about the co-evolving relationship between storm kinematics, microphysics, morphology and three-dimensional flash extent, despite its importance for lightning NOx production. To address this conceptual gap, the NASA MSFC Lightning Nitrogen Oxides Model (LNOM) is applied to North Alabama Lightning Mapping Array (NALMA) and Vaisala National Lightning Detection NetworkTM (NLDN) observations following ordinary convective cells through their lifecycle. LNOM provides estimates of flash type, channel length distributions, lightning segment altitude distributions (SADs) and lightning NOx production profiles. For this study, LNOM is applied in a Lagrangian sense to well isolated convective cells on 3 April 2007 (single cell and multi-cell hailstorm, non-severe multicell) and 6 July 2007 (non-severe multi-cell) over Northern Alabama. The LNOM lightning characteristics are compared to the evolution of updraft and precipitation properties inferred from dual-Doppler and polarimetric radar analyses applied to observations from a nearby Doppler radar network, including the UA Huntsville Advanced Radar for Meteorological and Operational Research (ARMOR, C-band, polarimetric). The LNOM estimated SAD and lightning NOx production profiles are placed in the context of radar derived profiles of vertical motion, precipitation types and amounts. Finally, these analyses are used to determine if storm integrated flash channel extent is as well correlated to volumetric updraft and precipitation ice characteristics in the mixed phase region as flash rate for these individual convective cells.

Carey, Lawrence↗

NASA NTRS1995

Lagrangian Transport Calculations Using UARS Data: Ozone - Part 2

Trajectory calculations are used to examine ozone transport in the polar winter stratosphere during periods of the Upper Atmosphere Research Satellite (UARS) observations. The value of these calculations for determining mass transport was demonstrated previously using UARS observations of long-lived tracers, In the middle stratosphere, the overall ozone behavior observed by the Microwave Limb Sounder in the polar vortex is reproduced by this purely dynamical model. Calculations show the evolution of ozone in the lower stratosphere during early winter to be dominated by dynamics in December 1992 in the Arctic. Calculations for June 1992 in the Antarctic show evidence of chemical ozone destruction and indicate that approx. 50% of the chemical destruction may be masked by dynamical effects, mainly diabatic descent, which bring higher ozone into the lower-stratospheric vortex. Estimating differences between calculated and observed fields suggests that dynamical changes masked approx. 20% - 35% of chemical ozone loss during late February and early March 1993 in the Arctic. In the Antarctic late winter, in late August and early September 1992, below approx. 520 K, the evolution of vortex-averaged ozone is entirely dominated by chemical effects; above this level, however, chemical ozone depletion can be partially or completely masked by dynamical effects. Our calculations for 1992 showed that chemical loss was nearly completely compensated by increases due to diabatic descent at 655 K.

Manney, Gloria L.↗

DOE OSTI2024

Wilsonian effective potentials and closed string field theory

We investigate Wilsonian effective field theory as a model for the construction of the tachyon potential and nonperturbative vacua in closed string field theory. In a number of cases we are able to find the effective potential exactly, and observe what appear to be universal features. We find that the effective field theory contains the same nonperturbative vacuum structure as the bare Lagrangian, though this information is encoded less efficiently as the distance scale of the effective field theory is increased. The implication is that closed string field theory plausibly contains information about the nonperturbative vacuum structure of string theory, in spite of its similarities to effective field theory. We also truncate the effective potential at a fixed power of the field and investigate how the global structure of the effective potential may be approximated via Padé resummation. Qualitative comparisons suggest that computation of the eighth to sixteenth order closed string vertex should be enough to obtain reliable results for the closed string field theory action evaluated on the tachyon field.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Sierra/SolidMechanics 5.0 Goodyear Specific

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 5.0 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Sierra/SolidMechanics 4.56 Goodyear User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 4.56 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

42 ENGINEERING↗

NASA NTRS1983

The effect of a coronal shock wave on the solar wind ionization state

In connection with studies of solar wind ionization state freezing, it is shown that, by using a Lagrangian approach of following individual fluid parcels, the techniques used previously for calculating ionization state variations in a steady state case can be extended straightforwardly to time-varying flows. The specific ionization state calculations presented are for a relatively simple picture of time-dependent coronal flow, based on a well-known model of a self-similar shock wave propagating through the corona. Time-dependent ionization effects for the sudden transition between two otherwise steady flows are likely to be limited to a narrow range of gas parcels which, having been shocked within the coronal freezing-in radius, pass a fixed interplanetary observer in an interval of a few tens of minutes. The amplitude of any rise in interplanetary ionization temperature associated with the coronal shock is likely to be considerably smaller than the jump in electron temperature that actually occurs in the corona.

Owocki, S. P.↗

NASA NTRSDate not supplied

Computational Icing Analysis on NASA’s SIDRM Geometry to Investigate Collection Efficiency

Computational icing analysis results were compared to experimental icing tunnel data including aerothermal (e.g. dry air) and supercooled water droplet rime-ice conditions from tests conducted in early 2022 at the NASA Icing Research Tunnel (IRT). The Simulated Intercompressor Duct Research Model (SIDRM) test article was used in this study, and its geometry represents the inter-compressor duct region of a turbofan engine. The test article’s purpose is to study the physics of supercooled water icing and ice crystal icing. This study compared three different icing codes: FENSAP-ICE (Eulerian approach), LEWICE3D (Lagrangian approach), and GlennICE (Lagrangian approach). All three icing codes were conducted on SIDRM’s complex body flow-field and compared to different experimental supercooled water rime runs. The test article instrumentation (pressure taps, thermocouples, etc.) and 3D laser scans of final ice shapes were used to compare against the different icing code simulations. The overall objectives are to understand how the icing codes handle capturing collection efficiency on the complex test article’s unheated surfaces. In the aerothermal cases, pressure tap readings matched the CFD results, but dry air CFD underpredicted thermocouple readings. Collection efficiency results from all three icing codes matched well together on the main body leading edge, main body slope, and the strut leading edges. All three icing codes underpredicted the final ice shape using a single-shot constant ice density approach, with more difficulty coming from the strut leading edge ice shape due to the swept wing like flow field. NASA’s overall goal for this effort is to develop computational icing tools to assist in the design and certification of engines for flight in icing conditions.

Supercooled icing, SIDRM, aerodynamics, computatio↗

DOE OSTI2025

The non-relativistic effective field theory of dark matter-electron interactions

Electronic excitations in atomic, molecular, and crystal targets are at the forefront of the ongoing search for light, sub-GeV dark matter (DM). In many light DM-electron interactions the energy and momentum deposited is much smaller than the electron mass, motivating a non-relativistic (NR) description of the electron. Thus, for any target, light DM-electron phenomenology relies on understanding the interactions between the DM and electron in the NR limit. In this work we derive the NR effective field theory (EFT) of general DM-electron interactions from a top-down perspective, starting from general high-energy DM-electron interaction Lagrangians. This provides an explicit connection between high-energy theories and their low-energy phenomenology in electron excitation based experiments. Furthermore, we derive Feynman rules for the DM-electron NR EFT, allowing observables to be computed diagrammatically, which can systematically explain the presence of in-medium screening effects in general DM models. We use these Feynman rules to compute absorption, scattering, and dark Thomson scattering rates for a wide variety of high-energy DM models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1983

Finite element methods for nonlinear elastostatic problems in rubber elasticity

A number of finite element methods for the analysis of nonlinear problems in rubber elasticity are outlined. Several different finite element schemes are discussed. These include the augmented Lagrangian method, continuation or incremental loading methods, and associated Riks-type methods which have the capability of incorporating limit point behavior and bifurcations. Algorithms for the analysis of limit point behavior and bifurcations are described and the results of several numerical experiments are presented. In addition, a brief survey of some recent work on modelling contact and friction in elasticity problems is given. These results pertain to the use of new nonlocal and nonlinear friction laws.

Oden, J. T.↗

NASA NTRS1985

A Lagrangian method for the shallow water equations based on a Voronoi mesh - Flows on a rotating sphere

A Lagrangian scheme using the Voronoi mesh is applied to study shallow water flow on a sphere. Discrete approximations to the shallow water equations are obtained for the surfaces of a nonrotating and a rotating sphere, and discrete differential operators are defined for the gradient and the divergence on the sphere. Dissipation is put into the model, when needed, by merging fluid points when they get too close to each other. The full numerical scheme is described and results of numerical computations on various test cases are given, including zonal flow and the Riemann problem.

Augenbaum, J. M.↗

NASA NTRS2013

An Investigation of the Kinematic and Microphysical Control of Lightning Rate, Extent and NOX Production using DC3 Observations and the NASA Lightning Nitrogen Oxides Model (LNOM)

The Deep Convective Clouds and Chemistry (DC3) experiment seeks to quantify the relationship between storm physics, lightning characteristics and the production of nitrogen oxides via lightning (LNOx). The focus of this study is to investigate the kinematic and microphysical control of lightning properties, particularly those that may govern LNOx production, such as flash rate, type and extent across Alabama during DC3. Prior studies have demonstrated that lightning flash rate and type is correlated to kinematic and microphysical properties in the mixed‐phase region of thunderstorms such as updraft volume and graupel mass. More study is required to generalize these relationships in a wide variety of storm modes and meteorological conditions. Less is known about the co‐evolving relationship between storm physics, morphology and three‐dimensional flash extent, despite its importance for LNOx production. To address this conceptual gap, the NASA Lightning Nitrogen Oxides Model (LNOM) is applied to North Alabama Lightning Mapping Array (NALMA) and Vaisala National Lightning Detection NetworkTM (NLDN) observations following ordinary convective cells through their lifecycle. LNOM provides estimates of flash rate, flash type, channel length distributions, lightning segment altitude distributions (SADs) and lightning NOx production profiles. For this study, LNOM is applied in a Lagrangian sense to multicell thunderstorms over Northern Alabama on two days during DC3 (21 May and 11 June 2012) in which aircraft observations of NOx are available for comparison. The LNOM lightning characteristics and LNOX production estimates are compared to the evolution of updraft and precipitation properties inferred from dual‐Doppler and polarimetric radar analyses applied to observations from a nearby radar network, including the UAH Advanced Radar for Meteorological and Operational Research (ARMOR). Given complex multicell evolution, particular attention is paid to storm morphology, cell mergers and possible dynamical, microphysical and electrical interaction of individual cells when testing various hypotheses.

Carey, Lawrence↗

DOE OSTI2020

“Lagrangian disks” in M-theory

While the study of bordered (pseudo-)holomorphic curves with boundary on Lagrangian submanifolds has a long history, a similar problem that involves (special) Lagrangian submanifolds with boundary on complex surfaces appears to be largely overlooked in both physics and math literature. We relate this problem to geometry of coassociative submanifolds in G 2 holonomy spaces and to Spin(7) metrics on 8-manifolds with T 2 fibrations. As an application to physics, we propose a large class of brane models in type IIA string theory that generalize brane brick models on the one hand and 2d theories T[M 4 ] on the other.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗