Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

CFD Model of Water Droplet Transport for ISS Hygiene Activity

The goal of the study is to assess the impacts of free water propagation in the Waste and Hygiene Compartment (WHC). Free water can be generated inside the WHC in small quantities due to crew hygiene activity. To mitigate potential impact of free water in Node 3 cabin the WHC doorway is enclosed by a waterproof bump-out, Kabin, with openings at the top and bottom. At the overhead side of the rack, there is a screen that prevents large drops of water from exiting. However, as the avionics fan in the WHC causes airflow toward the deck side of the rack, small quantities of free water may exit at the bottom of the Kabin. A Computational Fluid Dynamics (CFD) analysis of Node 3 cabin airflow made possible to identify the paths of water transport. The Node 3 airflow was computed for several ventilation scenarios. To simulate the droplet transport the Lagrangian discrete phase approach was used. Various initial droplet distributions were considered in the study. The droplet diameter was varied in the range of 2-20 mm. The results of the computations showed that most of the drops fall to the rack surface not far from the WHC curtain. The probability of the droplet transport to the adjacent rack surface with electronic equipment was predicted.

Son, Chang H.↗

Km‐Scale Simulations of Mesoscale Convective Systems Over South America—A Feature Tracker Intercomparison

Mesoscale convective systems (MCSs) are clusters of thunderstorms that are important in Earth's water and energy cycle. Additionally, they are responsible for extreme events such as large hail, strong winds, and extreme precipitation. Automated object-based analyses that track MCSs have become popular since they allow us to identify and follow MCSs over their entire life cycle in a Lagrangian framework. This rise in popularity was accompanied by an increasing number of MCS tracking algorithms, however, little is known about how sensitive analyses are concerning the MCS tracker formulation. Here, we assess differences between six MCS tracking algorithms on South American MCS characteristics and evaluate MCSs in kilometer-scale simulations with observational-based MCSs over 3 years. All trackers are run with a common set of MCS classification criteria to isolate tracker formulation differences. The tracker formulation substantially impacts MCS characteristics such as frequency, size, duration, and contribution to total precipitation. The evaluation of simulated MCS characteristics is less sensitive to the tracker formulation and all trackers agree that the model can capture MCS characteristics well across different South American climate zones. Dominant sources of uncertainty are the segmentation of cloud systems in space and time and the treatment of how MCSs are linked in time. Our results highlight that comparing MCS analyses that use different tracking algorithms is challenging. We provide general guidelines on how MCS characteristics compare between trackers to facilitate a more robust assessment of MCS statistics in future studies.

54 ENVIRONMENTAL SCIENCES↗

Color symmetry and confinement as an underlying superconformal structure in holographic QCD

Dedicated to the memory of our colleague, Harald Fritzsch, who, together with Murray Gell-Mann, introduced the color quantum number as the exact symmetry responsible for the strong interaction, thus establishing quantum chromodynamics (QCD) as a fundamental non-Abelian gauge theory. A basic understanding of hadron properties, however, such as confinement and the emergence of a mass scale, from first principles QCD has remained elusive: Hadronic characteristics are not explicit properties of the QCD Lagrangian and perturbative QCD, so successful in the large transverse momentum domain, is not applicable at large distances. In this article, we shall examine how this daunting obstacle is overcome in holographic QCD with the introduction of a superconformal symmetry in anti de Sitter (AdS) space which is responsible for confinement and the introduction of a mass scale within the superconformal group. When mapped to light-front coordinates in physical spacetime, this approach incorporates supersymmetric relations between the Regge trajectories of meson, baryon and tetraquark states which can be visualized in terms of specific SU(3) C color representations of quarks. Finally, we will also briefly discuss here the implications of holographic models for QCD color transparency in view of the present experimental interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Synoptic Mapping of Chemical Composition, Thermal Structure, and Air Motion from UARS Observations

An operational mapping algorithm was developed to process measurements of individual species observed by different satellite instruments on board UARS. Based on Fast Fourier Synoptic Mapping (FFSM), the algorithm accounts for the precessing orbit of UARS, the error of individual instruments, and gaps associated with instrument duty cycle and the satellite yaw maneuver. It provides synoptic structure and evolution on periods as short as 1 day, derived collectively from all observations of an individual species. The algorithm was applied to synoptically map temperature, thickness, and several chemical constituents observed by the instruments: MLS, CLAES, ISAMS, and HALOE. Each field variable observed by these instruments was mapped twice-daily in continuous global time series several months long. Mapped behavior produced via FFSM was compared against standard archived products generated via Kalman filtering. The standard map products reveal similar behavior, but are limited to features of larger scale. Tracer structure that develops through flow deformation and attending transport is therefore represented in those products more coarsely, eventually being sheared down to scales that are no longer represented properly. The synoptic time series also reveal a diurnal cycle for several of the constituents, one that emerges clearly in their space-time spectra. The zonal and meridional structure of diurnal variability was mapped by filtering the space-time spectrum, an intermediate product of FFSM, to those scales resolved by the UARS sampling. Geographical variations of the diurnal cycle have also been evaluated in its seasonal-mean structure by compositing distributions at individual local times. Both were compared to diurnal variations in chemical models and in dynamical models that account for transport by the diurnal tide. Diurnal variations were found to introduce spurious behavior into the archived products generated via Kalman filtering, behavior that is misrepresented in those synoptic time series. Aliasing by under-sampled diurnal variations introduces a bias into the time-mean structure of archived products. That systematic error is manifested as anomalous wave structure in the time-mean fields. By contrast, synoptic time series generated via FFSM faithfully represent such behavior as varying diurnally, so it does not corrupt time-mean structure. Tracer structure mapped via the above algorithm was being applied to determine the Brewer-Dobson circulation of the stratosphere through Lagrangian considerations. Preliminary results are encouraging: They reveal diabatic ascent in the tropics, balanced by diabatic descent over the winter hemisphere. This Lagrangian description of transport derived from UARS tracer measurements would be potentially valuable for diagnosing interannual changes of transport, specifically, in relation to interannual changes of total ozone. Regrettably, support was reduced to 50% of the approved funding level, leaving this component incomplete.

Salby, Murry L.↗

A Simple Derivation of the Acoustic Boundary Condition in the Presence of Flow

In modeling the noise propagation in a duct of a ducted fan engine, the authors needed the boundary condition satisfied by the acoustic pressure on a lined duct wall in the presence of flow. The results, in the form used by the authors, was derived by M.K. Myers, who presented a formal derivation of a related results by K. Taylor. Here, a brief and simple derivation of the acoustic boundary condition of Myers is given. The main difference in the derivation from that by Myers is that the Gaussian coordinate system (q(exp 1, q(exp 2), q(exp 3)(q(exp 1, q(exp 2), t)) is used to specify the instantaneous position, i.e. Lagrangian variable, of the point on the mean position of the wall with curvilinear (Gaussian) coordinates (q(exp 1), q(exp 2)). Myers uses a locally orthogonal coordinate system which is somewhat less specific than what is used here.

Farassat, F.↗

Application of the Marsupial Paradigm to Tropical Cyclone Formation from Northwestward-Propagating Disturbances

A wave-tracking algorithm is developed for northwestward-propagating waves that, on occasion, play a role in tropical cyclogenesis over the western oceans. To obtain the Lagrangian flow structure, the frame of reference is translated obliquely at the same propagation speed with the precursor disturbance. Trajectory analysis suggests that streamlines in the obliquely translated frame of reference can be used to approximate flow trajectories. The algorithm was applied to Super Typhoon Nakri (2008), Tropical Cyclone Erika (2009), and a few other examples. Diagnoses of meteorological analyses and satellite-derived moisture and precipitation fields show that the marsupial framework for tropical cyclogenesis in tropical easterly waves is relevant also for northwestward-propagating disturbances as are commonly observed in the tropical western Atlantic, the Gulf of Mexico, and the western North Pacific. Finally, it is suggested that analysis of the global model data and satellite observations in the marsupial framework can provide useful guidance on early tropical cyclone advisories.

Wang, Zhuo↗

A Generalized Eulerian-Lagrangian Analysis, with Application to Liquid Flows with Vapor Bubbles

Under a NASA MSFC SBIR Phase 2 effort an analysis has been developed for liquid flows with vapor bubbles such as those in liquid rocket engine components. The analysis is based on a combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it combines features of the so-called particle distribution approach with those of the so-called particle trajectory approach and can, in fact, be considered as a generalization of both of those traditional methods. The result of this generalization is a reduction in CPU time and memory requirements. Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle equations of motion (and, for certain applications, the particle temperature equation), although practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the simulation of bubble formation and growth have been included, as well as models for bubble drag and heat transfer. The results indicate that bubble formation is more or less 'explosive'. for a given flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and the surface tension. The bubble motion, on the other hand, is much less sensitive to the properties, but is affected strongly by the local pressure gradients in the flow field. In situations where either the material properties or the flow field are not known with sufficient accuracy, parametric studies can be carried out rapidly to assess the effect of the important variables. Future work will include application of the analysis to cavitation in inducer flow fields.

Dejong, Frederik J.↗

Shock-Layer Radiation Insights Available Through Flowfield Property Binning

This paper introduces a method for assessing similarity between available shock-layer radiation measurements and the radiative heating to a vehicle of interest. This method will be referred to as the flowfield property binning (FPB) approach. The challenge in developing this approach is that radiometer measurements from a surface normal line of sight, observed radiation measurements of the entire shock-layer from a distance, or measured shock-tube nonequilibrium profiles are notably different geometrically than the radiative flux to a vehicle surface. Therefore, the relative contributions of various regions of the shock-layer to the measurement, and how these relative contributions compare with the radiative heating to a vehicle of interest, are unclear. By tracking the radiation contribution from each flowfield location, and then grouping this contribution based on two local flowfield properties, these geometric differences are removed. This provide a two dimensional group, or bin, distribution that represents the radiation environment in the same format for various types of radiation simulations, such as observed radiation or radiative flux to a vehicle surface. In the limits of optically-thin radiation transport and thermochemical equilibrium flow, the choice of pressure and temperature as the binning flowfield properties is shown to provide a bin distribution that uniquely defines the radiation simulation. This concept is found to be sufficient for establishing similarity between two cases even when the optically-thin radiation transport and thermochemical equilibrium assumptions are not satisfied. A simple metric of comparison between bin distributions from different cases is provided by the histogram intersection algorithm. This algorithm results in a histogram intersection factor (H) between 0 and 1, which represents the fraction of the radiation that is similar between the two cases. Therefore, if H evaluated between bin distributions for two completely different radiation problems is near 1, then they represent essentially the same radiation problem. This means that they depend on the same weighting of various radiation modeling parameters evaluated at the same pressure and temperature. If one of these cases agrees well with a radiation measurement, then the matching bin distributions provide justification for applying this experimental validation to the other case. Even for H values as low as 0.2, the measurement case could be considered relevant if there are no other useful measurements. This approach allows observed radiation measurements to directly inform radiative heating simulations, which allows NASA’s past and future investments in observed radiation measurements to directly inform the radiative heating uncertainty or margins for future vehicle design. This FPB approach is also used to determine shock-tube conditions for informing backshell radiative heating simulations, where the normal shock velocity and Lagrangian time from the shock is used for the binning variables instead of pressure and temperature.

Shock Layer Radiation↗

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗

Analytical Dynamics and Nonrigid Spacecraft Simulation

Application to the simulation of idealized spacecraft are considered both for multiple-rigid-body models and for models consisting of combination of rigid bodies and elastic bodies, with the elastic bodies being defined either as continua, as finite-element systems, or as a collection of given modal data. Several specific examples are developed in detail by alternative methods of analytical mechanics, and results are compared to a Newton-Euler formulation. The following methods are developed from d'Alembert's principle in vector form: (1) Lagrange's form of d'Alembert's principle for independent generalized coordinates; (2) Lagrange's form of d'Alembert's principle for simply constrained systems; (3) Kane's quasi-coordinate formulation of D'Alembert's principle; (4) Lagrange's equations for independent generalized coordinates; (5) Lagrange's equations for simply constrained systems; (6) Lagrangian quasi-coordinate equations (or the Boltzmann-Hamel equations); (7) Hamilton's equations for simply constrained systems; and (8) Hamilton's equations for independent generalized coordinates.

Likins, P. W.↗

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗

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

A Simple Derivation of the Acoustic Boundary Condition in the Presence of Flow

In modelling the noise propagation in a duct of a ducted fan engine, the authors needed the boundary condition satisfied by the acoustic pressure on a lined duct wall in the presence of flow. The result in the form used by the authors was derived by M. K. Myers who presented a formal derivation of a related result by K. Taylor. A good review of the subject can be found in reference. Here, a brief and simple derivation of the acoustic boundary condition of Myers is given. The main difference in the derivation from that in reference is that the Gaussian co-ordinate system (q(sup 1), q(sup 2), q(sup 3)(q(sup 1), q(sup 2), t)) is used to specify the instantaneous position, i.e., Lagrangian variable, of the point on the mean position of the wall with curvilinear (Gaussian) co-ordinates (q(sup 1), q(sup 2)). Myers uses a locally orthogonal coordinate system which is somewhat less specific than what is used here.

Farassat, F.↗

Energy-Efficient Driving in Connected Corridors via Minimum Principle Control: Vehicle-in-the-Loop Experimental Verification in Mixed Fleets

Connected and automated vehicles (CAVs) can plan and actuate control that explicitly considers performance, system safety, and actuation constraints in a manner more efficient than their human-driven counterparts. In particular, eco-driving is enabled through connected exchange of information from signalized corridors that share their upcoming signal phase and timing (SPaT). This is accomplished in the proposed control approach, which follows first principles to plan a free-flow acceleration-optimal trajectory through green traffic light intervals by Pontryagin's Minimum Principle in a feedback manner. Urban conditions are then imposed from exogeneous traffic comprised of a mixture of human-driven vehicles (HVs) - as well as other CAVs. As such, safe disturbance compensation is achieved by implementing a model predictive controller (MPC) to anticipate and avoid collisions by issuing braking commands as necessary. The control strategy is experimentally vetted through vehicle-in-the-loop (VIL) of a prototype CAV that is embedded into a virtual traffic corridor realized through microsimulation. Up to 36% fuel savings are measured with the proposed control approach over a human-modelled driver, and it was found connectivity in the automation approach improved fuel economy by up to 26% over automation without. Additionally, the passive energy benefits realizable for human drivers when driving behind downstream CAVs are measured, showing up to 22% fuel savings in a HV when driving behind a small penetration of connectivity-enabled automated vehicles.

33 ADVANCED PROPULSION SYSTEMS↗

Traffic Control via Connected and Automated Vehicles (CAVs): An Open-Road Field Experiment with 100 CAVs

The CIRCLES project aims to reduce instabilities in traffic flow, which are naturally occurring phenomena due to human driving behavior. Also called “phantom jams” or “stop-and-go waves,” these instabilities are a significant source of wasted energy. Toward this goal, the CIRCLES project designed a control system, referred to as the MegaController by the CIRCLES team, that could be deployed in real traffic. Our field experiment, the MegaVanderTest (MVT), leveraged a heterogeneous fleet of 100 longitudinally controlled vehicles as Lagrangian traffic actuators, each of which ran a controller with the architecture described in this article. The MegaController is a hierarchical control architecture that consists of two main layers. The upper layer is called the Speed Planner and is a centralized optimal control algorithm. It assigns speed targets to the vehicles, conveyed through the LTE cellular network. The lower layer is a control layer, running on each vehicle. It performs local actuation by overriding the stock adaptive cruise controller, using the stock onboard sensors. The Speed Planner ingests live data feeds provided by third parties as well as data from our own control vehicles and uses both to perform the speed assignment. The architecture of the Speed Planner allows for the modular use of standard control techniques, such as optimal control, model predictive control (MPC), kernel methods, and others. The architecture of the local controller allows for the flexible implementation of local controllers. Corresponding techniques include deep reinforcement learning (RL), MPC, and explicit controllers. Depending on the vehicle architecture, all onboard sensing data can be accessed by the local controllers or only some. Likewise, control inputs vary across different automakers, with inputs ranging from torque or acceleration requests for some cars to electronic selection of adaptive cruise control (ACC) setpoints in others. The proposed architecture technically allows for the combination of all possible settings proposed previously, that is {Speed Planner algorithms} × {local Vehicle Controller algorithms} × {full or partial sensing} × {torque or speed control}. As a result, most configurations were tested throughout the ramp up to the MegaVandertest (MVT).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

"Traffic Control via Connected and Automated Vehicles: An Open-Road Field Experiment with 100 CAVs"

The CIRCLES project aims to reduce instabilities in traffic flow, which are naturally occurring phenomena due to human driving behavior. These "phantom jams" or "stop-and-go waves,"are a significant source of wasted energy. Toward this goal, the CIRCLES project designed a control system referred to as the MegaController by the CIRCLES team, that could be deployed in real traffic. Our field experiment leveraged a heterogeneous fleet of 100 longitudinally-controlled vehicles as Lagrangian traffic actuators, each of which ran a controller with the architecture described in this paper. The MegaController is a hierarchical control architecture, which consists of two main layers. The upper layer is called Speed Planner, and is a centralized optimal control algorithm. It assigns speed targets to the vehicles, conveyed through the LTE cellular network. The lower layer is a control layer, running on each vehicle. It performs local actuation by overriding the stock adaptive cruise controller, using the stock on-board sensors. The Speed Planner ingests live data feeds provided by third parties, as well as data from our own control vehicles, and uses both to perform the speed assignment. The architecture of the speed planner allows for modular use of standard control techniques, such as optimal control, model predictive control, kernel methods and others, including Deep RL, model predictive control and explicit controllers. Depending on the vehicle architecture, all onboard sensing data can be accessed by the local controllers, or only some. Control inputs vary across different automakers, with inputs ranging from torque or acceleration requests for some cars, and electronic selection of ACC set points in others. The proposed architecture allows for the combination of all possible settings proposed above. Most configurations were tested throughout the ramp up to the MegaVandertest.

Lee, Jonathan↗

The hidden structure of hydrodynamic transport in random fracture networks

We study the large-scale dynamics and prediction of hydrodynamic transport in random fracture networks. The flow and transport behaviour is characterized by first passage times and displacement statistics, which show heavy tails and anomalous dispersion with a strong dependence on the injection condition. The origin of these behaviours is investigated in terms of Lagrangian velocities sampled equidistantly along particle trajectories, unlike classical sampling strategies at a constant rate. The velocity series are analysed by their copula density, the joint distribution of the velocity unit scores, which reveals a simple, albeit hidden, correlation structure that can be described by a Gaussian copula. Based on this insight, we derive a Langevin equation for the evolution of equidistant particle speeds. In this framework, particle motion is quantified by a stochastic time-domain random walk, the joint density of particle position, and speed satisfies a Klein–Kramers equation. The upscaled theory quantifies particle motion in terms of the characteristic fracture length scale and the distribution of Eulerian flow velocities. That is, it is predictive in the sense that it does not require the a priori knowledge of transport attributes. The upscaled model captures non-Fickian transport features, and their dependence on the injection conditions in terms of the velocity point statistics and average fracture length. It shows that the first passage times and displacement moments are dominated by extremes occurring at the first step. The presented approach integrates the interaction of flow and structure into a predictive model for large-scale transport in random fracture networks.

42 ENGINEERING↗

"Designing, simulating, and performing the 100-AV field test for the CIRCLES consortium: Methodology and Implementation of the Largest mobile traffic control experiment to date"

Previous controlled experiments on single-lane ring roads have shown that a single partially autonomous vehicle (AV) can effectively mitigate traffic waves. This naturally prompts the question of how these findings can be generalized to field operational, high-density traffic conditions. To address this question, the Congestion Impacts Reduction via CAV-in-the-loop Lagrangian Energy Smoothing (CIRCLES) Consortium conducted MegaVanderTest (MVT), a live traffic control experiment involving 100 vehicles near Nashville, TN, USA. This article is a tutorial for developing analytical and simulation-based tools essential for designing and executing a live traffic control experiment like the MVT. It presents an overview of the proposed roadmap and various procedures used in designing, monitoring, and conducting the MVT, which is the largest mobile traffic control experiment at the time. The design process is aimed at evaluating the impact of the CIRCLES AVs on surrounding traffic. The article discusses the agent-based traffic simulation framework created for this evaluation. A novel methodological framework is introduced to calibrate this microsimulation, aiming to accurately capture traffic dynamics and assess the impact of adding 100 vehicles to existing traffic. The calibration model's effectiveness is verified using data from a six-mile section of Nashville's I-24 highway. The results indicate that the proposed model establishes an effective feedback loop between the optimizer and the simulator, thereby calibrating flow and speed with different spatiotemporal characteristics to minimize the error between simulated and real-world data. Finally, We simulate AVs in multiple scenarios to assess their effect on traffic congestion. This evaluation validates the AV routes, thereby contributing to the execution of a safe and successful live traffic control experiment via AVs.

Ameli, Mostafa↗