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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Dissipative MHD solutions for resonant Alfven waves in 1-dimensional magnetic flux tubes

The present paper extends the analysis by Sakurai, Goossens, and Hollweg (1991) on resonant Alfven waves in nonuniform magnetic flux tubes. It proves that the fundamental conservation law for resonant Alfven waves found in ideal MHD by Sakurai, Goossens, and Hollweg remains valid in dissipative MHD. This guarantees that the jump conditions of Sakurai, Goossens, and Hollweg, that connect the ideal MHD solutions for xi(sub r), and P' across the dissipative layer, are correct. In addition, the present paper replaces the complicated dissipative MHD solutions obtained by Sakurai, Goossens, and Hollweg for xi(sub r), and P' in terms of double integrals of Hankel functions of complex argument of order 1/3 with compact analytical solutions that allow a straight- forward mathematical and physical interpretation. Finally, it presents an analytical dissipative MHD solution for the component of the Lagrangian displacement in the magnetic surfaces perpen- dicular to the magnetic field lines xi(sub perpendicular) which enables us to determine the dominant dynamics of resonant Alfven waves in dissipative MHD.

Goossens, Marcel↗
DOE OSTI2024

A High-Order Eulerian–Lagrangian Runge–Kutta Finite Volume (EL–RK–FV) Method for Scalar Nonlinear Conservation Laws

Abstract We present a class of high-order Eulerian–Lagrangian Runge–Kutta finite volume methods that can numerically solve Burgers’ equation with shock formations, which could be extended to general scalar conservation laws. Eulerian–Lagrangian (EL) and semi-Lagrangian (SL) methods have recently seen increased development and have become a staple for allowing large time-stepping sizes. Yet, maintaining relatively large time-stepping sizes post shock formation remains quite challenging. Our proposed scheme integrates the partial differential equation on a space-time region partitioned by linear approximations to the characteristics determined by the Rankine–Hugoniot jump condition. We trace the characteristics forward in time and present a merging procedure for the mesh cells to handle intersecting characteristics due to shocks. Following this partitioning, we write the equation in a time-differential form and evolve with Runge–Kutta methods in a method-of-lines fashion. High-resolution methods such as ENO and WENO-AO schemes are used for spatial reconstruction. Extension to higher dimensions is done via dimensional splitting. Numerical experiments demonstrate our scheme’s high-order accuracy and ability to sharply capture post-shock solutions with large time-stepping sizes.

Chen, Jiajie↗
DOE OSTI2022

CoLoRe: fast cosmological realisations over large volumes with multiple tracers

We present CoLoRe, a public software package to efficiently generate synthetic realisations of multiple cosmological surveys.CoLoRe can simulate the growth of structure with different degrees of accuracy, with the current implementation supporting lognormal fields, first, and second order Lagrangian perturbation theory.CoLoRe simulates the density field on an all-sky light-cone up to a desired maximum redshift, and uses it to generate multiple 2D and 3D maps: galaxy positions and velocities, lensing (shear, magnification, convergence), integrated Sachs-Wolfe effect, line intensity mapping, and line of sight skewers for simulations of the Lyman-α forest. We test the accuracy of the simulated maps against analytical theoretical predictions, and showcase its performance with a multi-survey simulation including DESI galaxies and quasars, LSST galaxies and lensing, and SKA intensity mapping and radio galaxies. We expect CoLoRe to be particularly useful in studies aiming to characterise the impact of systematics in multi-experiment analyses, quantify the covariance between different datasets, and test cross-correlation pipelines for near-future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Demonstration of a new unstructured mesh IMC x-ray transport capability in LAP codes

In this document, the Advanced Simulation and Computing (ASC) Transport Project’s Jayenne team presents evidence that the Los Alamos National Laboratory (LANL) Level 2 milestone statement, Demonstrate a new unstructured mesh IMC x-ray transport capability in LAP codes, due Q4 FY20, has been satisfied with the release and integration of the Jayenne team’s Implicit Monte Carlo (IMC) solver libraries into the lumos multiphysics solver. The Jayenne project’s software includes the algorithms and features prescribed by the milestone description. The lumos software is maintained by the Lagrangian Applications Project (LAP) along with the hydrodynamics code flag and operates under the same ASC program as Jayenne. A set of demonstration problems has been specified, executed, and analyzed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

DYCAST: A finite element program for the crash analysis of structures

DYCAST is a nonlinear structural dynamic finite element computer code developed for crash simulation. The element library contains stringers, beams, membrane skin triangles, plate bending triangles and spring elements. Changing stiffnesses in the structure are accounted for by plasticity and very large deflections. Material nonlinearities are accommodated by one of three options: elastic-perfectly plastic, elastic-linear hardening plastic, or elastic-nonlinear hardening plastic of the Ramberg-Osgood type. Geometric nonlinearities are handled in an updated Lagrangian formulation by reforming the structure into its deformed shape after small time increments while accumulating deformations, strains, and forces. The nonlinearities due to combined loadings are maintained, and stiffness variation due to structural failures are computed. Numerical time integrators available are fixed-step central difference, modified Adams, Newmark-beta, and Wilson-theta. The last three have a variable time step capability, which is controlled internally by a solution convergence error measure. Other features include: multiple time-load history tables to subject the structure to time dependent loading; gravity loading; initial pitch, roll, yaw, and translation of the structural model with respect to the global system; a bandwidth optimizer as a pre-processor; and deformed plots and graphics as post-processors.

Pifko, A. B.↗
NASA NTRS2015

Resurrected DSCOVR Propulsion System - Challenges and Lessons Learned

The Deep Space Climate Observatory (DSCOVR), formerly known as Triana, is a unique mission, not because of its objectives but because of how long it was in storage before launch. The Triana spacecraft was built in the late 90s and later renamed as DSCOVR, but the project was canceled before the spacecraft was launched. The nearly-complete spacecraft was put in controlled storage for 10 years, until the National Oceanic and Atmospheric Administration (NOAA) provided funding for the National Aeronautics and Space Administration (NASA) to refurbish the spacecraft. On February 11, 2015, DSCOVR was launched on a Falcon 9 v1.1 from launch complex 40 at Cape Canaveral Air Force Station. This paper describes the DSCOVR propulsion system, which utilizes ten 4.5 N thrusters in blowdown mode to perform Midcourse Correction (MCC) maneuvers, Lissajous Orbit Insertion (LOI) at Lagrangian point L1, momentum unloading maneuvers, and station keeping delta-v maneuvers at L1. This paper also describes the testing that was performed, including susbsystem-level and spacecraft-level tests, to verify the propulsion system's integrity for flight. Finally, this paper concludes with a discussion of the challenges and lessons learned during this unique mission, including replacement of a bent thruster and installation of an auxiliary heater over existing propellant line heaters.

System↗
NASA NTRS2013

NASA Tech Briefs, January 2013

Topics include: Single-Photon-Sensitive HgCdTe Avalanche Photodiode Detector; Surface-Enhanced Raman Scattering Using Silica Whispering-Gallery Mode Resonators; 3D Hail Size Distribution Interpolation/Extrapolation Algorithm; Color-Changing Sensors for Detecting the Presence of Hypergolic Fuels; Artificial Intelligence Software for Assessing Postural Stability; Transformers: Shape-Changing Space Systems Built with Robotic Textiles; Fibrillar Adhesive for Climbing Robots; Using Pre-Melted Phase Change Material to Keep Payloads in Space Warm for Hours without Power; Development of a Centrifugal Technique for the Microbial Bioburden Analysis of Freon (CFC-11); Microwave Sinterator Freeform Additive Construction System (MS-FACS); DSP/FPGA Design for a High-Speed Programmable S-Band Space Transceiver; On-Chip Power-Combining for High-Power Schottky Diode-Based Frequency Multipliers; FPGA Vision Data Architecture; Memory Circuit Fault Simulator; Ultra-Compact Transputer-Based Controller for High-Level, Multi-Axis Coordination; Regolith Advanced Surface Systems Operations Robot Excavator; Magnetically Actuated Seal; Hybrid Electrostatic/Flextensional Mirror for Lightweight, Large-Aperture, and Cryogenic Space Telescopes; System for Contributing and Discovering Derived Mission and Science Data; Remote Viewer for Maritime Robotics Software; Stackfile Database; Reachability Maps for In Situ Operations; JPL Space Telecommunications Radio System Operating Environment; RFI-SIM: RFI Simulation Package; ION Configuration Editor; Dtest Testing Software; IMPaCT - Integration of Missions, Programs, and Core Technologies; Integrated Systems Health Management (ISHM) Toolkit; Wind-Driven Wireless Networked System of Mobile Sensors for Mars Exploration; In Situ Solid Particle Generator; Analysis of the Effects of Streamwise Lift Distribution on Sonic Boom Signature; Rad-Tolerant, Thermally Stable, High-Speed Fiber-Optic Network for Harsh Environments; Towed Subsurface Optical Communications Buoy; High-Collection-Efficiency Fluorescence Detection Cell; Ultra-Compact, Superconducting Spectrometer-on-a-Chip at Submillimeter Wavelengths; UV Resonant Raman Spectrometer with Multi-Line Laser Excitation; Medicine Delivery Device with Integrated Sterilization and Detection; Ionospheric Simulation System for Satellite Observations and Global Assimilative Model Experiments - ISOGAME; Airborne Tomographic Swath Ice Sounding Processing System; flexplan: Mission Planning System for the Lunar Reconnaissance Orbiter; Estimating Torque Imparted on Spacecraft Using Telemetry; PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils; Multiple-Frame Detection of Subpixel Targets in Thermal Image Sequences; Metric Learning to Enhance Hyperspectral Image Segmentation; Basic Operational Robotics Instructional System; Sheet Membrane Spacesuit Water Membrane Evaporator; Advanced Materials and Manufacturing for Low-Cost, High-Performance Liquid Rocket Combustion Chambers; Motor Qualification for Long-Duration Mars Missions.

Source record↗
NASA NTRS1985

The combined effects of stress and moisture on the structural integrity of composites

A variational approach is used to model the moisture content and distribution profiles which cause swelling in polymers and composites which are used as aircraft structural components. The moisture intrusion stems from humidity, which makes the problem one of diffusion of mass from the environment into a semiinfinite plate with a given diffusion coefficient and thickness. Lagrangian, ordinary differential equations are defined for the mass diffusion, including driving force, mass potential and volume dissipation expressions. Account is taken of the effects of stress on diffusivity, including the dilatational stress caused by inhomogeneous swelling of the material. It is noted that material near the top and bottom of a part undergoes the most swelling, and the evolution of moisture diffusion in the materials is traced. Finally, experimental data are used to define differential equations for predicting the final moisture distribution within the material.

Hiel, C.↗
DOE OSTI2025

Spatially resolved measurements of plasma ion velocity distributions in a dipole magnetic field

The equilibrium flows of a plasma discharge in a dipole magnetic field are a topic of interest in low temperature plasma physics. Experimentalists typically rely on probe-based and line-integrated diagnostic techniques in these environments to describe plasma behavior. Presented here are measurements of argon ion dynamics with laser induced fluorescence techniques to provide insight into plasma dynamics in dipole magnetic fields with nonperturbative, spatially localized measurements. Simulation results from a Lagrangian approach to track particle orbits are compared to measured density profiles and provide evidence to support the mechanism distinguishing experimental configurations is the initial approach of particles. Applying a negative DC bias to the magnet induces strong E×B flows around the magnet, even exceeding the ion acoustic speed as measured far from the magnet. A strong enough bias also produces two distinct ion populations and provides a method for controlling the density gradient on the equator.

McLaughlin, Jacob W. (ORCID:0000000152661888)↗
NASA NTRS2000

Climatology of the Finite-Volume CCM and Sensitivity Studies

A new model has been developed in the Data Assimilation Office at NASA Goddard Space Flight Center. It is intended for data assimilation and general circulation modeling, including representations of radiatively active trace gases such as stratospheric ozone. The model is based on the flux-form semi-Lagrangian' dynamical core, developed by Lin and Rood. In collaboration with the National Center for Atmospheric Research, the physical parametrization package from the Community Climate model, Version 3, has been included. This paper will discuss the climatology of a multi-year integration of the model, with observed boundary conditions (sea-surface temperature, sea ice, and incoming solar radiation). The standard model version has a resolution of 2.5 degrees longitude, 2 degrees latitude and has 55 levels between the surface and 80km. The mean climate and interannual variability of this model will be discussed, with emphasis on the middle atmosphere. The sensitivity to horizontal resolution will be examined, as will the impact of increasing and decreasing the vertical resolution. Changes which differ significantly from the interannual variations of the control run will be highlighted. Impacts of using long-term mean boundary conditions or values which vary from year to year will be discussed. Changes in the tropospheric circulation induced by lowering the upper boundary of the model will be discussed. The emphasis of the paper is on the climatology of the middle atmosphere and its impacts on the troposphere, with particular attention given to the climatologically important region near the tropopause.

Pawson, Steven↗
NASA NTRS2017

Investigating Sources of Ozone over California Using AJAX Airborne Measurements and Models: Assessing the Contribution from Long Range Transport

High ozone (O3) concentrations at low altitudes (1.5e4 km) were detected from airborne Alpha Jet Atmospheric eXperiment (AJAX) measurements on 30 May 2012 off the coast of California (CA). We investigate the causes of those elevated O3 concentrations using airborne measurements and various models. GEOS-Chem simulation shows that the contribution from local sources is likely small. A back trajectory model was used to determine the air mass origins and how much they contributed to the O3 over CA. Low-level potential vorticity (PV) from Modern Era Retrospective analysis for Research and Applications 2 (MERRA-2) reanalysis data appears to be a result of the diabatic heating and mixing of airs in the lower altitudes, rather than be a result of direct transport from stratospheric intrusion. The Q diagnostic, which is a measure of the mixing of the air masses, indicates that there is sufficient mixing along the trajectory to indicate that O3 from the different origins is mixed and transported to the western U.S.The back-trajectory model simulation demonstrates the air masses of interest came mostly from the mid troposphere (MT, 76), but the contribution of the lower troposphere (LT, 19) is also significant compared to those from the upper troposphere/lower stratosphere (UTLS, 5). Air coming from the LT appears to be mostly originating over Asia. The possible surface impact of the high O3 transported aloft on the surface O3 concentration through vertical and horizontal transport within a few days is substantiated by the influence maps determined from the Weather Research and Forecasting Stochastic Time Inverted Lagrangian Transport (WRF-STILT) model and the observed increases in surface ozone mixing ratios. Contrasting this complex case with a stratospheric-dominant event emphasizes the contribution of each source to the high O3 concentration in the lower altitudes over CA. Integrated analyses using models, reanalysis, and diagnostic tools, allows high ozone values detected by in-situ measurements to be attributed to multiple source processes.

Ryoo, Ju-Mee↗
DOE OSTI2019

An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables

We derive a formulation of the nonhydrostatic equations in spherical geometry with a Lorenz staggered vertical discretization. The combination conserves a discrete energy in exact time integration when coupled with a mimetic horizontal discretization. The formulation is a version of Dubos and Tort (2014, https://doi.org/10.1175/MWR-D-14-00069.1) rewritten in terms of primitive variables. It is valid for terrain following mass or height coordinates and for both Eulerian or vertically Lagrangian discretizations. The discretization relies on an extension to Simmons and Burridge (1981, https://doi.org/10.1175/1520-0493(1981)109<0758:AEAAMC>2.0.CO;2) vertical differencing, which we show obeys a discrete derivative product rule. This product rule allows us to simplify the treatment of the vertical transport terms. Energy conservation is obtained via a term-by-term balance in the kinetic, internal, and potential energy budgets, ensuring an energy-consistent discretization up to time truncation error with no spurious sources of energy. We demonstrate convergence with respect to time truncation error in a spectral element code with a horizontal explicit vertically implicit implicit-explicit time stepping algorithm.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Static geometric and material nonlinear analysis.

This paper presents a survey of four aspects of nonlinear analysis - basic formulation, plasticity relations, computational procedures, and methods of solution. This survey is for the most part limited to static geometric and material nonlinear analysis under the assumption of small strains. The two formulations covered are the total Lagrangian and the incremental moving coordinate formulations. A comparison of the two formulations is presented and an attempt is made to clarify the underlying basis of each formulation. The survey of computational procedures reveals that most researchers are now using some form of numerical integration to evaluate nonlinear contributions.

Stricklin, J. A.↗
NASA NTRS1992

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.↗
DOE OSTI2024

Chiral limit of 2d QCD revisited with lightcone conformal truncation

We study the chiral limit of 2d QCD with a single quark flavor at finite Nc using LCT. By modifying the LCT basis according to the quark mass in a manner motivated by ’t Hooft’s analysis, we are able to restore convergence for quark masses much smaller than the QCD strong coupling scale. For such small quark masses, the IR of the theory is expected to be well described by the Sine-Gordon model. We verify that LCT numerics are able to capture in detail the spectrum and correlation functions of the Sine-Gordon model. This opens up the possibility for studying deformations of various integrable CFTs using LCT by considering the chiral limit of QCD like theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Decomposing a renewable energy design and dispatch model

We address a mixed-integer linear programming model which selects a cost-minimizing set of available technologies with which to design a renewable energy system and prescribe their associated dispatch decisions. Realistically sized instances of such models pose computational challenges. To this end, we develop a Lagrangian heuristic based on a decomposition methodology which partitions the model into blocks and optimizes these more manageable, smaller subproblems. It also provides a lower bound to assess solution quality. In conclusion, we apply this methodology to the National Renewable Energy Laboratory's Renewable Energy Integration and Optimization (REopt TM ) model to generate near-optimal solutions to realistic instances containing, on average, approximately 300,000 variables and at least as many constraints, with a mean 30% optimality gap improvement using a five-minute solution time limit, compared to directly solving the original monolith.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

IMERG V06: Changes to the Morphing Algorithm

As the US Science Team’s globally gridded precipitation product from the NASA/JAXA Global Precipitation Measurement (GPM) mission, the Integrated Multi-satellitE Retrievals for GPM (IMERG) estimates the surface precipitation rates at 0.1° every half-hour using spaceborne sensors for various scientific and societal applications. One key component of IMERG is the morphing algorithm, which uses motion vectors to perform quasi-Lagrangian interpolation to fill in gaps in the passive microwave precipitation field using motion vectors. Up to IMERG V05, the motion vectors were derived from the large-scale motions of infrared observations of cloud tops. This study details the changes introduced in IMERG V06 to derive motion vectors from large-scale motions of selected atmospheric variables in numerical models, which allow IMERG estimates to be extended from the 60°N/S latitude band to the entire globe. Evaluation against both instantaneous passive microwave retrievals and ground measurements demonstrates the general improvement in the precipitation field of the new approach. Most of the model variables tested exhibited similar performance, but total precipitable water vapor was chosen as the source of the motion vectors for IMERG V06 due to its competitive performance and global completeness. Continuing assessments will provide further insights into possible refinements of this revised morphing scheme in future versions of IMERG.

Jackson Tan↗
NASA NTRSDate not supplied

Seasonal and Dayurnal Planetary Albedo Variability from Six Years of DSCOVR EPIC Data

Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar radiation from the Lissajous orbital position near the Lagrangian L1 point provide continuous monitoring of the Earth’s sunlit hemisphere. CERES-based angle models were used to convert the near-hourly reflected radiances of the EPIC images into a climate-style planetary albedo data-point over the sunlit hemisphere. Integration over the sunlit hemispheres averages out the meteorological weather noise, but retains the seasonal and planetary-scale variability. As the Earth rotates, this generates variations in the Earth’s planetary albedo that are precisely aligned in longitude, which constitutes the dayurnal cycle. This dayurnal variability in the planetary albedo arises from planetary-scale changes in cloud radiative properties that can be directly compared to similarly sampled climate GCM output data. Six years of EPIC data have been analyzed, showing characteristic patterns in the seasonal and dayurnal variability of the Earth’s planetary albedo. Much of the seasonal change in planetary albedo is associated with the changing DSCOVR viewing geometry and the change in solar declination. But throughout the year, the highest planetary albedos are observed over the Central Asia (Iraq) longitude, while the lowest planetary albedos occur over the Central Pacific longitude. For these longitudes, the relative seasonal changes in the planetary albedo are slowly varying and anti-correlated. Dayurnal amplitude maxima tend to occur during the July-September timeframe, with April-May and December exhibit distinct minima in the dayurnal amplitude. West Africa and the West Pacific longitudes likewise exhibit anti-correlated seasonal variability, while they also undergo anti-correlated short period oscillations. On the other hand, in nearby longitudes, there are short period spikes in planetary albedo of a few-days duration, as well as longer period oscillations that may range from a week to several months, that generally tend to be correlated.

Planetary Albedo↗