GEOS Visualization And Lagrangian dynamics Immersive eXtended Reality Tool (VALIXR) for Scientific Discovery Lightning Talk
Explore the source record for details and available documents.
Thomas Grubb↗
SEARCH · Engineering Papers
Results for “Extended Lagrangian dynamics”
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.
GEOS Visualization and Lagrangian Dynamics Immersive Extended Reality Tool (VALIXR) for Scientific Discovery
Explore the source record for details and available documents.
Thomas G. Grubb↗
GEOS Visualization And Lagrangian dynamics Immersive eXtended Reality Tool (VALIXR) for Scientific Discovery
Explore the source record for details and available documents.
Thomas G. Grubb↗
Using XR for Improving Scientific Discovery With Numerical Weather Models
Our work explores the use of extended reality (XR) to im- prove scientific discovery with numerical weather/climate models that inform Earth science digital twins, specifically the NASA Goddard Earth Observing System (GEOS) global atmospheric model. The overall project is named the Vi- sualization And Lagrangian dynamics Immersive eXtended Reality Toolkit (VALIXR), which has two main areas of focus: (1) enhancing the understanding of and interaction with model output data through advanced visualizations in the XR environment, and (2) the integration of Lagrangian dynamics into the GEOS model, which allows a natural, feature-specific analysis of Earth science phenomena as op- posed to traditional, fixed-point Eulerian dynamics. Here, we report initial work on these focus areas.
Thomas Grubb↗
Gigatraj: An Atmospheric Trajectory Model
Atmospheric trajectory models have a long history of success in tracking air motions in the lower stratosphere and upper troposphere over periods of up to a few days. Parcels have been traced backwards from observations to identify whatever phenomena (strong convection, volcanic eruptions, rocket launches, etc.) put their signature on them. Parcels have also been initialized at a known event and traced forward to examine their subsequent physical and chemical evolution. We describe a new trajectory model, "gigatraj," that aims to increase exibility by (a) making it straightforward to use new meteorological data sources, including those not based on regular latitude-longitude grids; (b) enabling a run-time choice of vertical coordinate system for kinematic and/or quasi-isentropic calculations; (c) allowing for the output of arbitrary meteorological products, selectable by the user and interpolated to the parcels' locations and times. The model can be run in a serial or parallel processing environment, so that large numbers of parcels can be traced in a reasonable time. Information is presented on model accuracy and performance. The former is demonstrated by runs using both test-pattern winds (comparing expected paths with actual output) and real-world winds (comparing forward and backward runs to characterize how well parcels retrace their paths). Sample cases are also shown, including a reverse domain lling (RDF) calculation illustrating a tropopause fold event. Model output can be displayed using the new Visualization And Lagrangian dynamics Immersive eXtended Reality (VALIXR) system, and an example will be shown. In addition, we describe work to incorporate a version of gigatraj into the Goddard Earth Observing System (GEOS) of NASA's Global Modeling and Assimilation O ce (GMAO) at Goddard Space Flight Center. This enables trajectory calculations to be performed within the running GEOS model at the latter's native time resolution, instead of using the winds from every few hours. It also provides access to all of GEOS's internal variables as they are calculated. This module may be useful, for example, for tracking rapid chemical changes in a Lagrangian framework.
dynamics↗
An extended Lagrangian method for subsonic flows
It is well known that fluid motion can be specified by either the Eulerian of Lagrangian description. Most of Computational Fluid Dynamics (CFD) developments over the last three decades have been based on the Eulerian description and considerable progress has been made. In particular, the upwind methods, inspired and guided by the work of Gudonov, have met with many successes in dealing with complex flows, especially where discontinuities exist. However, this shock capturing property has proven to be accurate only when the discontinuity is aligned with one of the grid lines since most upwind methods are strictly formulated in 1-D framework and only formally extended to multi-dimensions. Consequently, the attractive property of crisp resolution of these discontinuities is lost and research on genuine multi-dimensional approach has just been undertaken by several leading researchers. Nevertheless they are still based on the Eulerian description.
Liou, Meng-Sing↗
The dynamics of the Venus ionosphere. II - The effects of the time scale of the solar wind dynamic pressure variations
The effects on the upper dayside Venus ionosphere of a slow increase in solar wind dynamic pressure are simulated numerically with a one-dimensional (spherically symmetric) Lagrangian hydrodynamical code. The simulation is started with an extended ionosphere in pressure equilibrium with the solar wind at the ionopause. The pressure at the ionopause is gradually increased to five times the initial pressure with rise times of 5, 15, and 30 min. It is found that, for rise times greater than about 10 min, the compression of the ionopause is nearly adiabatic, with the ionopause moving downward at velocities of approximately 1-2 km/sec until it reaches a maximally compressed state, at which time the motion reverses. For short rise times the compression produces a shock wave similar to that occurring in the case of a sudden increase in pressure. The global implications of these processes are discussed within the context of Pioneer Venus observations and future theoretical work on this problem is outlined.
Stein, R. F.↗
The influence of dynamic inflow and torsional flexibility on rotor damping in forward flight from symbolically generated equations
The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The Helicopter Equations for Stability and Loads (HESL) program is extended to derive the governing equations of motion for the blade, and a Lagrangian formulation is used to obtain the equations in generalized coordinates. The program generates the steady-state and linearized perturbation equations in symbolic form and then codes them into FORTRAN subroutines. The coefficients for each equation and for each mode are identified through a numerical program; the latter can also be used to obtain the harmonic balance equations. The governing multiblade equations are derived explicitly using HESL. These equations can accommodate any number of elastic blade modes. Stability results are presented for several hingeless rotor blade structural models, and the influence of dynamic inflow in forward flight with an elastic hingeless rotor is investigated.
Reddy, T. S. R.↗
Correlated dynamical changes in the near-earth and distant magnetotail regions - ISEE 3
It is pointed out that the redeployment of ISEE 3 from its continuous monitoring of the solar wind in a large orbit about the upstream Lagrangian point to an extended magnetotail orbit has afforded an opportunity for deep-tail passage of October 1982, and in the radial range from 200 to 220 R(E) during the near-apogee part of the second tail passage in January-February 1983. Attention is given to instrumentation and data sets, spacecraft positions, and observational data.
Baker, D. N.↗
Dynamic Analysis and Control of Lightweight Manipulators with Flexible Parallel Link Mechanisms
The objective is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. The resulting equation of motion have a structure which is useful to reduce the number of terms calculated, to check correctness, or to extend the model to higher order. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. Elastic motion is expressed by the assumed mode method. Mode shape functions of each link are chosen using the load interfaced component mode synthesis. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model.
Lee, Jeh Won↗
Application of a novel finite difference method to dynamic crack problems
A versatile finite difference method (HEMP and HEMP 3D computer programs) was developed originally for solving dynamic problems in continuum mechanics. It was extended to analyze the stress field around cracks in a solid with finite geometry subjected to dynamic loads and to simulate numerically the dynamic fracture phenomena with success. This method is an explicit finite difference method applied to the Lagrangian formulation of the equations of continuum mechanics in two and three space dimensions and time. The calculational grid moves with the material and in this way it gives a more detailed description of the physics of the problem than the Eulerian formulation.
Chen, Y. M.↗
Multivariable control of the Space Shuttle Remote Manipulator System using linearization by state feedback
This paper demonstrates an approach to nonlinear control system design that uses linearization by state feedback to allow faster maneuvering of payloads by the Shuttle Remote Manipulator System (SRMS). A nonlinear feedback law is defined to cancel the nonlinear plant dynamics so that a linear controller can be designed for the SRMS. First a nonlinear design model was generated via SIMULINK. This design model included nonlinear arm dynamics derived from the Lagrangian approach, linearized servo model, and linearized gearbox model. The current SRMS position hold controller was implemented on this system. Next, a trajectory was defined using a rigid body kinematics SRMS tool, KRMS. The maneuver was simulated. Finally, higher bandwidth controllers were developed. Results of the new controllers were compared with the existing SRMS automatic control modes for the Space Station Freedom Mission Build 4 Payload extended on the SRMS.
Gettman, Chang-Ching L.↗
Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design
Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.
Fahrenthold, Eric P.↗
Distributed Constrained Optimization with Semicoordinate Transformations
Recent work has shown how information theory extends conventional full-rationality game theory to allow bounded rational agents. The associated mathematical framework can be used to solve constrained optimization problems. This is done by translating the problem into an iterated game, where each agent controls a different variable of the problem, so that the joint probability distribution across the agents moves gives an expected value of the objective function. The dynamics of the agents is designed to minimize a Lagrangian function of that joint distribution. Here we illustrate how the updating of the Lagrange parameters in the Lagrangian is a form of automated annealing, which focuses the joint distribution more and more tightly about the joint moves that optimize the objective function. We then investigate the use of "semicoordinate" variable transformations. These separate the joint state of the agents from the variables of the optimization problem, with the two connected by an onto mapping. We present experiments illustrating the ability of such transformations to facilitate optimization. We focus on the special kind of transformation in which the statistically independent states of the agents induces a mixture distribution over the optimization variables. Computer experiment illustrate this for &sat constraint satisfaction problems and for unconstrained minimization of NK functions.
Macready, William↗
Mature thunderstorm cloud top structure - Three-dimensional numerical simulation versus satellite observations
Preliminary results of four runs with a three-dimensional model of the effects of vertical wind shear on cloud top height/temperature structure and the internal properties of isolate midlatitude thunderstorms are reported. The model is being developed as an aid to analyses of GEO remote sensing satellite data. The grid is a 27 x 27 x 20 mesh with 2 km horizontal resolution and 0.9 vertical resolution. The total grid is 54 km on a side and 18 km deep. A second-order Crowley scheme for advecting momentum is extended with a third-order correction for spatial truncation error, and the earth-relative horizontal surface wind components are decreased to 50 percent of their values at 0.45 km. A temperature increase with height is included, together with an initial impulse consisting of a nonrotating cylindrical weak buoyant updraft 10 km in radius. The results of the runs are discussed in terms of the time variation of the vertical velocity extrema, the effects of strong and weak shear on a storm, the cloud top height, the Lagrangian dynamics of a thermal couplet, and data from a real storm.
Schlesinger, R. E.↗
LSPRAY: Lagrangian Spray Solver for Applications With Parallel Computing and Unstructured Gas-Phase Flow Solvers
Sprays occur in a wide variety of industrial and power applications and in the processing of materials. A liquid spray is a phase flow with a gas as the continuous phase and a liquid as the dispersed phase (in the form of droplets or ligaments). Interactions between the two phases, which are coupled through exchanges of mass, momentum, and energy, can occur in different ways at different times and locations involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the rate-controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as other phenomena. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on sprays to unstructured grids and parallel computing. LSPRAY, which was developed by M.S. Raju of Nyma, Inc., is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo probability density function (PDF) solver. The LSPRAY solver accommodates the use of an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements in the gas-phase solvers. It is used specifically for fuel sprays within gas turbine combustors, but it has many other uses. The spray model used in LSPRAY provided favorable results when applied to stratified-charge rotary combustion (Wankel) engines and several other confined and unconfined spray flames. The source code will be available with the National Combustion Code (NCC) as a complete package.
Raju, Manthena S.↗
Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows
The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.
Source record↗
Development of an Extended Reality (XR) Tool for Earth Science Visualization
In this presentation, we will discuss our work in adapting the NASA open source XR software, the Mixed Reality Exploration Toolkit (MRET), to an earth science domain. MRET is a NASA open source XR software for rapidly building extended reality (XR) environments for NASA domain problems, e.g., pulling in CAD models of thermal vac chamber and Roman Space Telescope to do fit checks. Primarily used for hardware integration & test, we have been adapting and extending MRET for science problems. Traditionally, scientists view and analyze the result of calculated or measured observables with static 1-D, 2-D or 3-D plots. It can be difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. Additionally, numerical models, such as the NASA GEOS climate model, are almost exclusively formulated and analyzed on Eulerian grids with points fixed in space and time. However, atmospheric phenomena such as convective clouds, hurricanes and wildfire smoke plumes move with the 3-D flow field, and it is often difficult and unnatural to understand these phenomena in an Eulerian reference frame as opposed to the Lagrangian reference frame in which nature operates. As part of an Earth Science Technology Office (ESTO) proposal, we have been adapting MRET to be a scientific exploration and analysis XR tool with integrated Lagrangian Dynamics (LD) for the Goddard Earth Observing System (GEOS) numerical weather prediction model. We believe this will help scientists identify, track, and understand the evolution of Earth Science phenomena. This presentation will discuss current results in our work in developing this XR tool for Earth Science.
Thomas G Grubb↗