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 685 records · Page 38

Investigation of turbulent inflow specification in Euler–Lagrange simulations of mid-field spray

The process of atomization of a liquid jet by a parallel high-speed gas stream results in a spray, whose downstream development is of considerable interest to several applications. The round jet spray can be spatially divided into (i) a near-field (near-nozzle) region of liquid atomization and (ii) a downstream mid-field region of fully-dispersed droplets. In order to accurately model mid-field droplet dispersion, this work aims at developing a rigorous and robust injection model for Euler–Lagrange spray simulations. Results from experiments are used to obtain the relevant droplet number density, size distribution, and mean and standard deviation velocity distributions of the injection model, systematically in a step-by-step process. Two-phase large eddy simulations are performed by stochastically generating the Lagrangian droplets at the inlet of the mid-field region. Number flux, diameter distribution, mean velocity, and other time-averaged statistics at several downstream locations are shown to agree well with the corresponding experimental data.

42 ENGINEERING↗

FIERRO V.X

FIERRO is a modern C++ code intended to simulate quasi-static solid mechanics problems and transient, compressible material dynamic problems with Lagrangian methods, which have meshes with constant mass elements that move with the material. FIERRO is designed to aid a) material model research that has historically been done using commercial implicit and explicit finite element codes, b) numerical methods research, and c) computer science research. FIERRO supports user developed material models that adhere to several industry standard formats by using a C++ to Fortran interface to couple the model to the numerical solvers. FIERRO is built on the ELEMENTS library that supports a diverse suite of element types, including high-order elements, and quadrature rules. The mesh class within the ELEMENTS library is designed for efficient calculations on unstructured meshes and to minimize memory usage. FIERRO is designed to readily accommodate a range of numerical methods including continuous finite element, finite volume, and discontinuous Galerkin methods. FIERRO is designed to support explicit and implicit time integration methods. FIERRO only supports a single material in an element. No physical data exists within the code.

Morgan, Nathaniel↗

Homogeneous and heterogeneous chemistry along air parcel trajectories

The study of coupled heterogeneous and homogeneous chemistry due to polar stratospheric clouds (PSC's) using Lagrangian parcel trajectories for interpretation of the Airborne Arctic Stratosphere Experiment (AASE) is discussed. This approach represents an attempt to quantitatively model the physical and chemical perturbation to stratospheric composition due to formation of PSC's using the fullest possible representation of the relevant processes. Further, the meteorological fields from the United Kingdom Meteorological office global model were used to deduce potential vorticity and inferred regions of PSC's as an input to flight planning during AASE.

Jones, R. L.↗

Hover performance analysis of advanced rotor blades

This is an effort aimed at validating recent hover prediction methods. The experimental basis for this validation work is an extensive set of loads, wake and performance data, which were obtained from a pressure instrumented model UH-60 rotor tested at the Sikorsky hover test facility and at Duits-Nederlandse Windtunnel (DNW). This model was equipped with replaceable tips - including a tapered and a BERP-type tip - which permitted studies of the effects of rotor geometry. The central prediction method studied is a free-wake, vortex embedded, full-potential CFD method - called HELIX-I. It is found that the HELIX-I code produces very good comparisons with the data including wake, surface pressure and performance. Comparisons with the measured radial load distributions have permitted an improved understanding of the wake resolution modelling requirements of CFD methods. Since HELIX-I is a combined Eulerian/Lagrangian method, limited comparisons are also made with a Lagrangian boundary element code (called EHPIC) and an Eulerian Navier-Stokes code (called TURNS). In most cases all methods produce good comparisons with the data. It is found that the HELIX-I code provides a good compromise between the speed of boundary integral methods and the comprehensive nature of Navier-Stokes methods.

Tung, C.↗

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

The Dynamics of Megafire Smoke Plumes in Climate Models: Why a Converged Solution Matters for Physical Interpretations

As the climate system warms, megafires have become more frequent with devastating effects. A byproduct of these events is the creation of smoke plumes that can rise into the stratosphere and spread across the globe where they reside for many months. To gain a deeper understanding of the plume dynamics, global climate simulations of a megafire were performed at a wide range of grid spacings from 2.0° down to 7 km, including a 7 km nonhydrostatic experiment. The analysis focuses on how the resolved dynamics affects the specification of the plume characteristics such as injection height and black carbon (BC) mass. Prior studies initialize the smoke plume at one or a few grid points and this is shown here to produce severely dissipative dynamics. In order to validate such simulations with observations, enhancements of the plume characteristics to offset the dissipation is necessary. Using a numerically converged simulation, sensitivity tests show that to approximate the observed stratospheric lifetime, a reduction in BC fraction by 50% is necessary for external mixtures. The vorticity dynamics of the plume is also analyzed with a Lagrangian budget to understand the mechanisms responsible for the evolution of a collocated anticyclonic vortex. The results can be distilled down into a simple conceptual model. As the plume rises, the air diverges at the top of the updraft where the largest concentrations of smoke are found. This divergence induces a dilution of the background cyclonic absolute vorticity producing an anticyclonic vortex. Vortex decay occurs from opposite arguments.

54 ENVIRONMENTAL SCIENCES↗

Simulation of the modern arctic climate by the NCAR CCM1

The National Center of Atmospheric Research (NCAR) Community Climate Model Version 1 (CCM1's) simulation of the modern arctic climate is evaluated by comparing a five-year seasonal cycle simulation with the European Center for Medium-Range Weather Forecasts (ECMWF) global analyses. The sea level pressure (SLP), storm tracks, vertical cross section of height, 500-hPa height, total energy budget, and moisture budget are analyzed to investigate the biases in the simulated arctic climate. The results show that the model simulates anomalously low SLP, too much storm activity, and anomalously strong baroclinicity to the west of Greenland and vice versa to the east of Greenland. This bias is mainly attributed to the model's topographic representation of Greenland. First, the broadened Greenland topography in the model distorts the path of cyclone waves over the North Atlantic Ocean. Second, the model oversimulates the ridge over Greenland, which intensifies its blocking effect and steers the cyclone waves clockwise around it and hence produces an artificial circum-Greenland trough. These biases are significantly alleviated when the horizontal resolution increases to T42. Over the Arctic basin, the model simulates large amounts of low-level (stratus) clouds in winter and almost no stratus in summer, which is opposite to the observations. This bias is mainly due to the location of the simulated SLP features and the negative anomaly of storm activity, which prevent the transport of moisture into this region during summer but favor this transport in winter. The moisture budget analysis shows that the model's net annual precipitation (P-E) between 70 deg N and the North Pole is 6.6 times larger than the observations and the model transports six times more moisture into this region. The bias in the advection term is attributed to the positive moisture fixer scheme and the distorted flow pattern. However, the excessive moisture transport into the Arctic basin does not solely result from the advection term. The contribution by the moisture fixer is as large as from advection. By contrast, the semi-Lagrangian transport scheme used in the CCM2 significantly improves the moisture simulation for this region; however, globally the error is as serious as for the positive moisture fixer scheme. Finally, because the model has such serious problems in simulating the present arctic climate, its simulations of past and future climate change for this region are questionable.

Bromwich, David H.↗

Nonlinear Shell Modeling of Thin Membranes with Emphasis on Structural Wrinkling

Thin solar sail membranes of very large span are being envisioned for near-term space missions. One major design issue that is inherent to these very flexible structures is the formation of wrinkling patterns. Structural wrinkles may deteriorate a solar sail's performance and, in certain cases, structural integrity. In this paper, a geometrically nonlinear, updated Lagrangian shell formulation is employed using the ABAQUS finite element code to simulate the formation of wrinkled deformations in thin-film membranes. The restrictive assumptions of true membranes, i.e. Tension Field theory (TF), are not invoked. Two effective modeling strategies are introduced to facilitate convergent solutions of wrinkled equilibrium states. Several numerical studies are carried out, and the results are compared with recent experimental data. Good agreement is observed between the numerical simulations and experimental data.

Tessler, Alexander↗

Advances in reduction techniques for tire contact problems

Some recent developments in reduction techniques, as applied to predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities, are reviewed. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of reduction techniques, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface. Also, the research topics which have high potential for enhancing the effectiveness of reduction techniques are outlined.

Noor, Ahmed K.↗

Analysis of liquid rocket combustion chamber turbulence levels from diffusion data

The intensity of turbulence and the Lagrangian correlation coefficient in a liquid-rocket combustion chamber have been analytically determined from an analysis of experimental diffusion data obtained in a small rocket engine which operated at 300-psia chamber pressure and produced approximately 250 pounds thrust. Results of gas-sample measurements obtained by Orsat and gas-chromatograph techniques to determine helium-concentration profiles were analyzed on the basis of Taylor's (1921) turbulent diffusion theory to obtain turbulence flow-field parameters. The results of the analysis indicate that turbulent diffusion in a combustion chamber can be adequately modeled by the one-dimensional Taylor theory, which assumes that the intensity of turbulence is a function only of axial distance in the chamber and that the Lagrangian correlation coefficient is expressed by a power law. The results indicate a higher intensity of turbulence and lower correlation than previously expected.

Smith, L. O., Jr.↗

Reduction technique for tire contact problems

A reduction technique and a computational procedure are presented for predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of the reduction technique, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface.

Noor, Ahmed K.↗

Meshfree simulation and experimental validation of extreme thermomechanical conditions in friction stir extrusion

Abstract Friction stir extrusion (FSE) is a novel solid-phase processing technique that consolidates and extrudes metal powders, flakes, chips, or billets into high-performance parts by plastic deformation, which has the potential to save substantial processing time and energy. Currently, most studies on FSE are experimental and only a few numerical models have been developed to explain and predict the complex physics of the process. In this work, a meshfree simulation framework based on smoothed particle hydrodynamics (SPH) was developed for FSE. Unlike traditional grid-based methods, SPH is a Lagrangian particle-based method that can handle severe material deformations, capture moving interfaces and surfaces, and monitor the field variable histories explicitly without complicated tracking schemes. These aspects of SPH make it attractive for the FSE process, where in situ evolution of field variables is difficult to observe experimentally. To this end, a 3-D, fully thermomechanically coupled SPH model was developed to simulate the FSE of aluminum wires. The developed model was thoroughly validated by comparing the numerically predicted material flow, strain, temperature history, and extrusion force with experimental results for a certain set of process parameters. The validated SPH model can serve as an effective tool to predict and better understand the extreme thermomechanical conditions during the FSE process.

99 GENERAL AND MISCELLANEOUS↗

Neutrino masses from low scale partial compositeness

We consider a class of models in which the neutrinos acquire Majorana masses through mixing with singlet neutrinos that emerge as composite states of a strongly coupled hidden sector. In this framework, the light neutrinos are partially composite particles that obtain their masses through the inverse seesaw mechanism. We focus on the scenario in which the strong dynamics is approximately conformal in the ultraviolet, and the compositeness scale lies at or below the weak scale. The small parameters in the Lagrangian necessary to realize the observed neutrino masses can naturally arise as a consequence of the scaling dimensions of operators in the conformal field theory. We show that this class of models has interesting implications for a wide variety of experiments, including colliders and beam dumps, searches for lepton flavor violation and neutrinoless double beta decay, and cosmological observations. At colliders and beam dumps, this scenario can give rise to striking signals involving multiple displaced vertices. The exchange of hidden sector states can lead to observable rates for flavor violating processes such as μ → eγ and μ → e conversion. If the compositeness scale lies at or below a hundred MeV, the rate for neutrinoless double beta decay is suppressed by form factors and may be reduced by an order of magnitude or more. The late decays of relic singlet neutrinos can give rise to spectral distortions in the cosmic microwave background that are large enough to be observed in future experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A LaGrangian Perspective of the Surface Energy Budget Contributions to Sea Ice Evolution in the ‘New Arctic’

Recent Arctic warming outpaces global warming and is characterized by dramatic changes in sea ice (record sea ice loss, thinner and more seasonal ice cover), warming SSTs, and a lengthening melt season. These changes impact the surface energy budget and the sensitivity of sea ice to radiative forcing, namely an increased sensitivity of thinner ice to local radiative forcings and a weaker buffering of surface-atmosphere exchanges of heat and moisture. However, these processes are a large source of uncertainty in global climate models with prior studies finding that poor representations of downwelling longwave radiation and sea ice albedo lead to model spread in sea ice conditions. A limitation of prior work is examining sea ice-atmosphere interactions through an Eulerian framework. This study utilizes a sea ice parcel dataset in a LaGrangian framework, allowing for tracking of surface energy budget-sea ice interactions over the lifespan of individual sea ice parcels; this framework also allows for the quantification of the sea ice response to episodic weather events along with seasonal drivers, including the drivers of rapid growth and melt events. First, lag-lead relationships between relevant surface energy budget terms and sea ice thickness changes are explored to understand the atmospheric conditions preceding rapid growth and melt. The impact of episodic events- moisture intrusion, warm air outbreaks, cyclones- on sea ice growth are also investigated by stratifying first-year parcels by sea ice thickness at the end of the growth season and analyzing the effect of atmospheric variability (number and intensity of episodic events) on parcel thickness change. Lastly, parcel growth and melt and the associated radiative anomalies are compared to snowfall events/snow thickness to determine how snowfall modifies the interaction between the local surface energy budget and sea ice mass balance.

Patrick Taylor↗

Evaluation of Stratocumulus Evolution Under Contrasting Temperature Advections in CESM2 Through a Lagrangian Framework

Abstract This study leveraged a Lagrangian framework to examine the evolution of stratocumulus clouds under cold and warm advections (CADV and WADV) in the Community Earth System Model 2 (CESM2) against observations. We found that CESM2 simulates a too rapid decline in low‐cloud fraction (LCF) and cloud liquid water path (CLWP) under CADV conditions, while it better aligns closely with observed LCF under WADV conditions but overestimates the increase in CLWP. Employing an explainable machine learning approach, we found that too rapid decreases in LCF and CLWP under CADV conditions are related to overestimated drying effects induced by sea surface temperature, whereas the substantial increase in CLWP under WADV conditions is associated with the overestimated moistening effects due to free‐tropospheric moisture and surface winds. Our findings suggest that overestimated drying effects of sea surface temperature on cloud properties might be one of crucial causes for the high equilibrium climate sensitivity in CESM2.

54 ENVIRONMENTAL SCIENCES↗

Orbital stability analysis and photometric characterization of the second Earth Trojan asteroid 2020 XL5

Trojan asteroids are small bodies orbiting around the L 4 or L 5 Lagrangian points of a Sun-planet system. Due to their peculiar orbits, they provide key constraints to the Solar System evolution models. Despite numerous dedicated observational efforts in the last decade, asteroid 2010 TK 7 has been the only known Earth Trojan thus far. Here we confirm that the recently discovered 2020 XL 5 is the second transient Earth Trojan known. To study its orbit, we used archival data from 2012 to 2019 and observed the object in 2021 from three ground-based observatories. Our study of its orbital stability shows that 2020 XL 5 will remain in L 4 for at least 4 000 years. With a photometric analysis we estimate its absolute magnitude to be ${H}_{r}=18.5{8}_{-0.15}^{+0.16}$, and color indices suggestive of a C-complex taxonomy. Assuming an albedo of 0.06 ± 0.03, we obtain a diameter of 1.18 ± 0.08 km, larger than the first known Earth Trojan asteroid.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum disorder, duality, and fractional statistics in 2 + 1 dimensions

A low-energy equivalence between two apparently unrelated Lagrangians with fractional statistics is reported. Exploiting this equivalence, it is possible to study the quantum disordered phase of the nonlinear sigma model with Hopf term. It is found that the quasi-particles in the disordered phase also have fractional statistics. There appears to be a dual relationship between the ordered and disordered phases.

Wen, X. G.↗