Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian integration”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Adaptation of multidimensional group particle tracking and particle wall-boundary condition model to the FDNS code

A particulate two-phase flow CFD model was developed based on the FDNS code which is a pressure based predictor plus multi-corrector Navier-Stokes flow solver. Turbulence models with compressibility correction and the wall function models were employed as submodels. A finite-rate chemistry model was used for reacting flow simulation. For particulate two-phase flow simulations, a Eulerian-Lagrangian solution method using an efficient implicit particle trajectory integration scheme was developed in this study. Effects of particle-gas reaction and particle size change to agglomeration or fragmentation were not considered in this investigation. At the onset of the present study, a two-dimensional version of FDNS which had been modified to treat Lagrangian tracking of particles (FDNS-2DEL) had already been written and was operational. The FDNS-2DEL code was too slow for practical use, mainly because it had not been written in a form amenable to vectorization on the Cray, nor was the full three-dimensional form of FDNS utilized. The specific objective of this study was to reorder to calculations into long single arrays for automatic vectorization on the Cray and to implement the full three-dimensional version of FDNS to produce the FDNS-3DEL code. Since the FDNS-2DEL code was slow, a very limited number of test cases had been run with it. This study was also intended to increase the number of cases simulated to verify and improve, as necessary, the particle tracking methodology coded in FDNS.

Chen, Y. S.↗
NASA NTRS2020

Reaching for 20 Years with the IMERG Multi-Satellite Products

The latest releases of Global Precipitation Measurement (GPM) mission products cap five years of vigorous development cycle since the launch of the GPM Core Observatory, and these now provide datasets that are relatively homogeneous across the joint Tropical Rainfall Measuring Mission (TRMM) and GPM eras. Version 06 of the U.S. GPM team's Integrated Multi-satellitE Retrievals for GPM (IMERG) merged precipitation product enforces a consistent intercalibration for all precipitation products computed from individual satellites with the TRMM and GPM Core Observatory sensors as the TRMM- and GPM era calibrators, respectively, and incorporates monthly surface gauge data. The basic IMERG algorithm now features precipitation motion vectors (used to drive the Lagrangian interpolation, or "morphing") that arecomputed by tracking vertically integrated vapor fields analyzed in MERRA2 and GEOS5. This innovation provides globally complete coverage, expanding IMERG's coverage beyond the 60°N-S latitude band provided by IR-based vectors, although we continue to mask out precipitation over snowy/icy surfaces as unreliable. A second innovation is the Quality Index (QI) data field. The half-hourly QI is taken as the approximate Kalman Filter correlation computed in the morphing calculation.

Huffman, George J.↗
NASA NTRS2019

Status and Examples for the Version 06 IMERG Multi-Satellite Products

After five years of development following the launch of the Global Precipitation Measurement (GPM) missionCore Observatory, the GPM data products are now being extended across the joint Tropical Rainfall MeasuringMission (TRMM) and GPM eras. Version 06 of the U.S. GPM team's Integrated Multi-satellitE Retrievals forGPM (IMERG) merged precipitation product provides a consistent intercalibration for all precipitation productscomputed from individual satellites with the TRMM and GPM Core Observatory sensors as the TRMM- andGPM-era calibrators, respectively, and incorporates monthly surface gauge data. One major change in the basicIMERG algorithm for V06 is that precipitation motion vectors (used to drive the quasi-Lagrangian interpolation,or "morphing") are computed by tracking vertically integrated vapor (TQV) fields analyzed in MERRA2 andGEOS5. This innovation provides globally complete coverage, expanding IMERG's coverage beyond the 60N-Slatitude band previously provided by IR-based vectors, although precipitation over snowy/icy surfaces is stillmasked out as unreliable. A second innovation is that the Quality Index (QI) data field computed for the half-hourlydatasets has been refined to include estimates of correlation at microwave overpass times.We will summarize the processing status for V06 IMERG, for which the retrospective processing shouldbe actively advancing at meeting time. We will show early examples of performance. For example, the TQVmotion vectors are typically slightly better than the IR-based vectors at all latitudes. The transition across theTRMM/GPM data boundary will be discussed, including the necessity of filling in the TRMM-based calibrationsover the latitude band 35-65 in each hemisphere. The notional schedule for the eventual retirement of thepredecessor TRMM Multi-satellite Precipitation Analysis (TMPA) multi-satellite dataset will be updated as well.

Huffman, George J.↗
NASA NTRS1974

Numerical calculations of velocity and pressure distribution around oscillating airfoils

An analytical procedure based on the Navier-Stokes equations was developed for analyzing and representing properties of unsteady viscous flow around oscillating obstacles. A variational formulation of the vorticity transport equation was discretized in finite element form and integrated numerically. At each time step of the numerical integration, the velocity field around the obstacle was determined for the instantaneous vorticity distribution from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints, using the Lagrangian Multiplier Technique, at each time step of the numerical integration. The procedure was then applied for determining pressures around obstacles oscillating in unsteady flow. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗
NASA NTRS1983

Ignition of confined gaseous mixtures by hot surfaces and hot wires

Ignition times and spatial and temporal variations of temperature and concentration in gaseous mixtures confined between two infinite parallel walls or two infinite cylinders have been obtained by numerical integration of the appropriate conservation equations written in Lagrangian coordinates. Ignition times and ignition energies are presented for the case of an isothermal wall in terms of the initial mixture pressure and equivalence ratio for both one and two-step chemical reaction mechanisms. The numerical results indicate that there is a critical mixture pressure for which the ignition time is minimum. The values of this critical pressure are larger (smaller) than 1 atm for the one- (two-) step reaction mechanism. The critical pressure for the ignition time is not equal to the critical pressure for the ignition energy. The ignition time and energy decrease with the equivalence ratio within a certain range and then remain constant.

Ramos, J. I.↗
NASA NTRS1992

Finite element Euler calculations of unsteady transonic cascade flows

A Galerkin finite element procedure incorporating an explicit Runge-Kutta time-stepping scheme has been developed in this work to solve unsteady transonic flow in cascades. The computational domain is discretized by a globally unstructured but locally structured blade-fitted deformable mesh. The Galerkin approximation is applied to the unsteady Euler equations based on a mixed Eulerian-Lagrangian description. The semi-discretized equations are integrated forward in time using a multistage Runge-Kutta scheme. An artificial dissipation operator of the type proposed by Jameson is adapted in the current scheme to capture shocks and suppress nonphysical oscillations. Phase-shifted boundary conditions are used to reduce the computational domain to a single reference passage. Results for both steady and unsteady transonic flows through cascades are presented and compared to existing finite volume solutions.

Hsiao, Chingteng↗
NASA NTRS1994

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗
NASA NTRS2003

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↗
NASA NTRS1990

Behaviour of plate elements based on the first-order shear deformation theory

A new analytical technique to assess the performance of shear deformable elements is presented, which makes it possible to determine a priori whether a given element will lock when used to model thin structures. The role that shear constraints play in determining the behavior of thin elements was established by comparing the results of key numerical tests with the predictions of element behavior made by studying the form of the shear constraints. Conclusions regarding locking behavior and the effects of reduced integration in thin shear deformable elements are presented, including the findings (1) that singularity of the shear stiffness matrix is not sufficient to avoid locking; (2) that the effect of mesh refinement on an element that contains spurious constraints is two-fold; and (3) that reduced integration does not remove spurious constraints but rather relaxes them. The results of the study are in agreement with previous studies of Mindlin plate elements in regarding Lagrangian elements as superior to serendipity elements when either full or reduced integration is employed.

Averill, R. C.↗
NASA NTRS2018

Extratropical Cyclone Precipitation Life Cycles: A Satellite-Based Analysis

Although extratropical cyclones are the most common midlatitude storms, the relationship between their precipitation life cycle and dynamical strength life cycle has not been thoroughly analyzed. Given that thermodynamic heat exchanges associated with precipitation impact cyclone circulation, there is a need to understand the precipitation/dynamics relationship. Based on Integrated MultisatellitE Retrievals for Global Precipitation MeasurementGPM precipitation and Lagrangian cyclone tracks, the precipitation maximum occurs prior to the dynamical strength maximum 70% of the time. The lag in timing is consistent with the difference in cyclone precipitable water vapor at the two peaks. Conditional subsetting of the cyclone composites shows that if the precipitable water vapor distribution is constrained to be equal throughout the composite life cycle, the precipitation peak occurs very near the time of the peak in cyclone dynamical strength. Thus, the boost in dynamical strength caused by latent heat associated with precipitation manifest itself with little to no time lag.

Extratropical cyclone↗
NASA NTRSDate not supplied

IMERG Multi-Satellite Products Across Two Decades

The Version 06 Global Precipitation Measurement (GPM) mission products were completed over the last year, capping five years of development since the launch of the GPM Core Observatory, and covering the joint Tropical Rainfall Measuring Mission (TRMM) and GPM eras with consistently processed algorithms. The U.S. GPM team’s Integrated Multi-satellitE Retrievals for GPM (IMERG) merged precipitation product enforces a consistent intercalibration for all precipitation products computed from individual satellites with the TRMM and GPM Core Observatory sensors as the TRMM- and GPM-era calibrators, respectively, and incorporates monthly surface gauge data in the Final (research) product. Mid-latitude calibrations during the TRMM era necessarily are more approximate because TRMM only covered the latitude band 35°N-S, while GPM covers 65°N-S. Starting in V06, IMERG employs precipitation motion vectors (used to drive the quasi-Lagrangian interpolation, or “morphing”) that are computed by tracking the vertically integrated vapor as analyzed in MERRA2 and GEOS FP. This approach covers the entire globe, expanding coverage beyond the 60°N-S latitude band provided by IR-based vectors in previous versions, although we choose to mask out microwave-based precipitation over snowy/icy surfaces as unreliable. We will provide examples of performance for the V06 IMERG products, including comparison with the long-term record of GPCP and TMPA, showing higher values by about 8% in the latitude band 50°N-S over oceans; diurnal cycle, demonstrating improvement over previous versions; and daily precipitation PDFs for the entire record, showing a shift at the TRMM/GPM boundary, as well as interannual variations. These analyses have important implications for the utility of V06 IMERG data for long-record calculations. Finally, we will review the retirement of the predecessor TMPA multi-satellite dataset.

George J. Huffman↗
NASA NTRS2013

Large-Deformation Displacement Transfer Functions for Shape Predictions of Highly Flexible Slender Aerospace Structures

Large deformation displacement transfer functions were formulated for deformed shape predictions of highly flexible slender structures like aircraft wings. In the formulation, the embedded beam (depth wise cross section of structure along the surface strain sensing line) was first evenly discretized into multiple small domains, with surface strain sensing stations located at the domain junctures. Thus, the surface strain (bending strains) variation within each domain could be expressed with linear of nonlinear function. Such piecewise approach enabled piecewise integrations of the embedded beam curvature equations [classical (Eulerian), physical (Lagrangian), and shifted curvature equations] to yield closed form slope and deflection equations in recursive forms.

Deflection equations↗
NASA NTRS1975

Monitoring estuarine circulation and ocean waste dispersion using an integrated satellite-aircraft-drogue approach

The author has identified the following significant results. An inexpensive, integrated drogue-aircraft-satellite approach was developed which is based on the Lagrangian technique and employs remotely tracked drogues and dyes together with satellite observation of natural tracers, such as suspended sediment. Results include current circulation studies in Delaware Bay in support of an oil slick movement model; investigations of the dispersion and movement of acid wastes dumped 40 miles off the Delaware coast; and coastal current circulation. In each case, the integrated drogue-aircraft-satellite approach compares favorably with other techniques on the basis of accuracy, cost effectiveness, and performance under severe weather conditions.

Klemas, V.↗
NASA NTRS1979

Vortex simulation of three-dimensional, spotlike disturbances in a laminar boundary layer

The growth of a turbulent spot in a laminar boundary layer as the spot evolves from a localized disturbance in the layer, is simulated numerically using a three-dimensional vortex filament description of the vorticity field. The filaments are marked with a sequence of node points which are tracked in a Lagrangian reference frame. Velocity computation is done by Biot-Savart integration. Although some discrepancies with experiment appear to exist in the near wall region, the gross properties of the spot, including the velocities of the leading and trailing edges and the velocity perturbations away from the wall, are in good agreement with experiment.

Leonard, A.↗
NASA NTRS1979

Vortex simulation of three-dimensional, spotlike disturbances in a laminar boundary layer

The growth of a turbulent spot in a laminar boundary layer, as the spot evolves from a localized disturbance in the layer, is simulated numerically using a three-dimensional vortex filament description of the vorticity field. The filaments are marked with a sequence of mode points which are tracked in a Lagrangian reference frame. Velocity computation is done by Biot-Savart integration. Although some discrepancies with experiment appear to exist in the near wall region, the gross properties of the spot, including the velocities of the leading and trailing edges and the velocity perturbations away from the wall, are in good agreement with experiment.

Leonard, A.↗
NASA NTRS2003

A Vertically Lagrangian Finite-Volume Dynamical Core for Global Models

A finite-volume dynamical core with a terrain-following Lagrangian control-volume discretization is described. The vertically Lagrangian discretization reduces the dimensionality of the physical problem from three to two with the resulting dynamical system closely resembling that of the shallow water dynamical system. The 2D horizontal-to-Lagrangian-surface transport and dynamical processes are then discretized using the genuinely conservative flux-form semi-Lagrangian algorithm. Time marching is split- explicit, with large-time-step for scalar transport, and small fractional time step for the Lagrangian dynamics, which permits the accurate propagation of fast waves. A mass, momentum, and total energy conserving algorithm is developed for mapping the state variables periodically from the floating Lagrangian control-volume to an Eulerian terrain-following coordinate for dealing with physical parameterizations and to prevent severe distortion of the Lagrangian surfaces. Deterministic baroclinic wave growth tests and long-term integrations using the Held-Suarez forcing are presented. Impact of the monotonicity constraint is discussed.

Lin, Shian-Jiann↗
NASA NTRSDate not supplied

Effect of Space Radiation on TES Detectors Performance

The Athena mission and its X-IFU instrument (X-ray Integral Field Unit) will be positioned at the Sun-Earth Lagrangian point L1, where it will be subject to solar wind (low flux plasma of 95% protons and 5% alpha particles) with energy below 0.1 MeV, and to galactic cosmic rays and solar flares (energetic protons, alpha particles, and electrons) with energies up to hundreds of MeV for protons and the GeV for heavier ions. Some of these particles will go through the satellite and hit the focal plane assembly and hence the detectors. These detectors will be TES (Transition-edge sensor) microcalorimeters, flown for the first time in such an environment. In order to ensure the performance of this type of detectors throughout the duration of such mission, it is critical to study the impact of the radiation on their behavior. Indeed, although a lot of reference material exist for semiconductor detectors such as CCDs, little is currently known about the impact of radiation on TES detectors. These energetic events could cause local heating or damage to the detectors and affect their performance. In this work, we describe how we designed a test campaign to assess the impact of L1 radiation on TES detectors for Athena/X-IFU-like missions and present the results of the tests. Analyses includes assessing changes in the pulse shapes and energy resolution of the detectors measured at 55 mK after several radiation dose steps performed at 4 K.

Sophie Beaumont↗
NASA NTRS1977

Variational energy principle for compressible, baroclinic flow. 1: First and second variations of total kinetic action

The case of a cold gas in the absence of external force fields is considered. Since the only energy involved is kinetic energy, the total kinetic action (i.e., the space-time integral of the kinetic energy density) should serve as the total free-energy functional in this case, and as such should be a local minimum for all possible fluctuations about stable flow. This conjecture is tested by calculating explicit, manifestly covariant expressions for the first and second variations of the total kinetic action in the context of Lagrangian kinematics. The general question of the correlation between physical stability and the convexity of any action integral that can be interpreted as the total free-energy functional of the flow is discussed and illustrated for the cases of rectillinear and rotating shearing flows.

Schmid, L. A.↗