Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian approach”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Laser-driven flyer plate impact: Computational studies guided by experiments

We present a computational approach using a multimaterial, arbitrary Lagrangian–Eulerian code termed ALE3D to model the nanosecond/micrometer dynamics of the launch of 0.5–4.5 km/s laser-driven metal flyer plates and the impact with stationary targets of Pyrex and fused silica glasses, and Lexan and Plexiglas polymers, producing pressures in the target in the range of 5–20 GPa. The simulations are compared to experimental results where the flyer velocity profile and the velocity profile imparted to the target material were measured with high-speed velocimetry. The experimental flyer launch by a high-intensity pulsed laser is modeled by depositing heat into a thin vaporizable layer under the flyer plate. This model produces a flyer plate that has not been exposed to the laser pulse, allowing us to compare the properties of the real flyer to a simulated ideal flyer. The simulations of target impact are generally in good agreement with the experiment except at the highest impact velocities where the shock release process in the experiment is slower than that in the simulation. The cause of this disagreement is attributed to an inadequate description of the shock viscosity during the nanosecond unloading process.

Stekovic, Svjetlana (ORCID:0000000163878070)↗
DOE OSTI2023

Rotational symmetry protected edge and corner states in Abelian topological phases

Spatial symmetries can enrich the topological classification of interacting quantum matter and endow systems with non-trivial strong topological invariants (protected by internal symmetries) with additional "weak" topological indices. In this paper, we study the edge physics of systems with a non-trivial shift invariant, which is protected by either a continuous $\text{U}(1)_r$ or discrete $\text{C}_n$ rotation symmetry, along with internal $\text{U}(1)_c$ charge conservation. Specifically, we construct an interface between two systems which have the same Chern number but are distinguished by their Wen-Zee shift and, through analytic arguments supported by numerics, show that the interface hosts counter-propagating gapless edge modes which cannot be gapped by arbitrary local symmetry-preserving perturbations. Using the Chern-Simons field theory description of two-dimensional Abelian topological orders, we then prove sufficient conditions for continuous rotation symmetry protected gapless edge states using two complementary approaches. One relies on the algebraic Lagrangian sub-algebra framework for gapped boundaries while the other uses a more physical flux insertion argument. For the case of discrete rotation symmetries, we extend the field theory approach to show the presence of fractional corner charges for Abelian topological orders with gappable edges, and compute them in the case where the Abelian topological order is placed on the two-dimensional surface of a Platonic solid. Finally, our work paves the way for studying the edge physics associated with spatial symmetries in symmetry enriched topological phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1981

Atmospheric studies related to aerospace activities and remote sensing technology

Parallel investigations were conducted relating to: the sensitivity of 1-D photochemical model simulated column ozone perturbations due to a projected fleet of 1000 aircraft cruising 7 hours per day at altitudes of 15-16 and 18-19 km to uncertainties in kinetic rate constant data determining modeled OH concentrations and eddy diffusivity profile parameterization and a comparison of the inherent strengths and weaknesses of Eulerian and Langrangian averaging processes in the development of multidimensional models and investigation of approaches to applying the Generalized Lagrangian Mean (GLM) formalism to zonal-mean models. The role of multiple scattering and Earth curvature in the evaluation of diurnally dependent photodissociation rates and trace species variations was examined.

Sze, N. D.↗
NASA NTRS1989

Blade mistuning coupled with shaft flexibility effects in rotor aeroelasticity

The effect of bladed-disk polar dissymmetry, resulting from variations in mass from one blade to another, on aeroelastic stability boundaries for a fan stage is presented. In addition to both in-plane and out-of-plane deformations of the bladed-disk, bending of the supporting shaft in two planes is considered, and the resulting Coriolis forces and gyroscopic moments are included in the analysis. A quasi-steady aerodynamics approach is combined with the Lagrangian method to develop the governing equations of motion for the flexible bladed-disk-shaft assembly. Calculations are performed for an actual fan stage.

Khader, Naim↗
NASA NTRS1992

A new Lagrangian method for real gases at supersonic speed

With the renewed interest in high speed flights, the real gas effect is of theoretical as well as practical importance. In the past decade, upwind splittings or Godunov-type Riemann solutions have received tremendous attention and as a result significant progress has been made both in the ideal and non-ideal gas. In this paper, we propose a new approach that is formulated using the Lagrangian description, for the calculation of supersonic/hypersonic real gas inviscid flows. This new formulation avoids the grid generation step which is automatically obtained as the solution procedure marches in the 'time-like' direction. As a result, no remapping is required and the accuracy is faithfully maintained in the Lagrangian level. In this paper, we give numerical results for a variety of real gas problems consisting of essential elements in high speed flows, such as shock waves, expansion waves, slip surfaces and their interactions. Finally, calculations for flows in a generic inlet and nozzle are presented.

Loh, C. Y.↗
NASA NTRS1994

Eulerian-Lagrangian Simulations of Transonic Flutter Instabilities

This paper presents an overview of recent applications of Eulerian-Lagrangian computational schemes in simulating transonic flutter instabilities. This approach, the fluid-structure system is treated as a single continuum dynamics problem, by switching from an Eulerian to a Lagrangian formulation at the fluid-structure boundary. This computational approach effectively eliminates the phase integration errors associated with previous methods, where the fluid and structure are integrated sequentially using different schemes. The formulation is based on Hamilton's Principle in mixed coordinates, and both finite volume and finite element discretization schemes are considered. Results from numerical simulations of transonic flutter instabilities are presented for isolated wings, thin panels, and turbomachinery blades. The results suggest that the method is capable of reproducing the energy exchange between the fluid and the structure with significantly less error than existing methods. Localized flutter modes and panel flutter modes involving traveling waves can also be simulated effectively with no a priori knowledge of the type of instability involved.

Bendiksen, Oddvar O.↗
NASA NTRS1980

Lagrangian methods in nonlinear plasma wave interaction

Analysis of nonlinear plasma wave interactions is usually very complicated, and simplifying mathematical approaches are highly desirable. The application of averaged-Lagrangian methods offers a considerable reduction in effort, with improved insight into synchronism and conservation (Manley-Rowe) relations. This chapter indicates how suitable Lagrangian densities have been defined, expanded, and manipulated to describe nonlinear wave-wave and wave-particle interactions in the microscopic, macroscopic and cold plasma models. Recently, further simplifications have been introduced by the use of techniques derived from Lie algebra. These and likely future developments are reviewed briefly.

Crawford, F. W.↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
DOE OSTI2024

Weak boundary conditions for Lagrangian shock hydrodynamics: A high-order finite element implementation on curved boundaries

Here, we propose a new Nitsche-type approach for weak enforcement of normal velocity boundary conditions for a Lagrangian discretization of the compressible shock-hydrodynamics equations using high-order finite elements on curved boundaries. Specifically, the variational formulation is appropriately modified to enforce free-slip wall boundary conditions, without perturbing the structure of the function spaces used to represent the solution, with a considerable simplification with respect to traditional approaches. Total energy is conserved and the resulting mass matrices are constant in time. The robustness and accuracy of the proposed method are validated with an extensive set of tests involving nontrivial curved boundaries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Quantitative proton radiography and shadowgraphy for arbitrary intensities

Charged-particle radiography and shadowgraphy data can be directly inverted to obtain a line-integrated transverse Lorentz force or a line-integrated transverse refractive index gradient if intensity modulations due to scattering and absorption are negligible, and angular deflections are small. We develop a new direct-inversion algorithm based on plasma physics and compare it to a new Monge–Ampère code and an existing power diagram code. The measured or source intensity is represented by electrons subject to drag, and the other intensity by fixed ions. The decrease in kinetic plus electrostatic energy determines convergence. The displacement of the electrons from their initial to their equilibrium positions determines the line-integrated force or refractive index gradient. We have implemented two approaches: PIC (particle in cell) and Lagrangian fluid, in 1-D and 2-D. The PIC code works for arbitrary intensities, can work efficiently in parallel, and can make use of existing codes. The Lagrangian code requires less memory and is faster than the PIC code without massively parallel processing, but fails in 2-D for large intensity modulations. The Monge–Ampère code is by far the fastest in 2-D, without massively parallel processing, but fails for intensities with large voids, high contrast ratios and large deflections across the boundaries, and could not obtain the degree of convergence possible with the PIC code. As a result, the power diagram code was by far the slowest and most memory intensive, and failed for large peaks in the measured intensity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Dynamic Subgrid Turbulence Modeling for Shallow Cumulus Convection Simulations beyond LES Resolutions

A scale-dependent dynamic Smagorinsky model is implemented in the Met Office/NERC Cloud (MONC) model using two averaging flavors, along Lagrangian pathlines and local moving averages. The dynamic approaches were compared against the conventional Smagorinsky–Lilly scheme in simulating the diurnal cycle of shallow cumulus convection. The simulations spanned from the LES to the near-gray-zone and gray-zone resolutions and revealed the adaptability of the dynamic model across the scales and different stability regimes. The dynamic model can produce a scale- and stability-dependent profile of the subfilter turbulence length scale across the chosen resolution range. At gray-zone resolutions the adaptive length scales can better represent the early precloud boundary layer leading to temperature and moisture profiles closer to the LES compared to the standard Smagorinsky. As a result, the initialization and general representation of the cloud field in the dynamic model is in good agreement with the LES. In contrast, the standard Smagorinsky produces a less well-mixed boundary layer, which fails to ventilate moisture from the boundary layer, resulting in the delayed spinup of the cloud layer. Moreover, strong downgradient diffusion controls the turbulent transport of scalars in the cloud layer. However, the dynamic approaches rely on the resolved field to account for nonlocal transports, leading to overenergetic structures when the boundary layer is fully developed and the Lagrangian model is used. Introducing the local averaging version of the model or adopting a new Lagrangian time scale provides stronger dissipation without significantly affecting model behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

An Improved Convection Parameterization with Detailed Aerosol–Cloud Microphysics for a Global Model

Abstract A new microphysical treatment that includes aerosol–cloud interactions and secondary ice production (SIP) mechanisms is implemented in the convection scheme of the Community Atmosphere Model, version 6 (CAM6). The approach is to embed a 1D Lagrangian parcel model in the bulk convective plume of the existing deep convection parameterization. Aerosol activation, growth processes including collision/coalescence, and three processes of SIP mechanisms, two of which are normally overlooked in atmospheric models, are represented in this embedded parcel model. These microphysical processes are treated with a hybrid bin/bulk scheme and a high spatial and temporal resolution for the integration of the embedded parcel in 1D, allowing vertical velocity to determine the microphysical evolution following the in-cloud motion during ascent. Simulations of an observed case (Midlatitude Continental Convective Clouds Experiment) of a mesoscale convective system in Oklahoma, United States, with a single-column model (SCAM) version of CAM, are compared with aircraft in situ and ground-based observations of microphysical properties from the convection and precipitation. Results from the validation show the new microphysical scheme has a good representation of the ice initiation in the bulk convective plume, including the known and empirically quantified pathways of primary and secondary initiation, with benefits for the accuracy of properties of its supercooled cloud liquid. The sensitivity simulations and use of tagging tracers for the validated simulation confirm that the newly included SIP mechanisms are of paramount importance for convective microphysics and can be successfully treated in the global model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2017

Osiris-REx Spacecraft Current Status and Forward Plans

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

Messenger, Scott↗
NASA NTRS1994

Floating shock fitting via Lagrangian adaptive meshes

In recent works we have formulated a new approach to compressible flow simulation, combining the advantages of shock-fitting and shock-capturing. Using a cell-centered Roe scheme discretization on unstructured meshes, we warp the mesh while marching to steady state, so that mesh edges align with shocks and other discontinuities. This new algorithm, the Shock-fitting Lagrangian Adaptive Method (SLAM) is, in effect, a reliable shock-capturing algorithm which yields shock-fitted accuracy at convergence. Shock-capturing algorithms like this, which warp the mesh to yield shock-fitted accuracy, are new and relatively untried. However, their potential is clear. In the context of sonic booms, accurate calculation of near-field sonic boom signatures is critical to the design of the High Speed Civil Transport (HSCT). SLAM should allow computation of accurate N-wave pressure signatures on comparatively coarse meshes, significantly enhancing our ability to design low-boom configurations for high-speed aircraft.

Vanrosendale, John↗