SEARCH · Engineering Papers

Results for “dissipative 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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Hawaii Ocean Mixing Experiment: Program Summary

It is becoming apparent that insufficient mixing occurs in the pelagic ocean to maintain the large scale thermohaline circulation. Observed mixing rates fall a factor of ten short of classical indices such as Munk's "Abyssal Recipe." The growing suspicion is that most of the mixing in the sea occurs near topography. Exciting recent observations by Polzin et al., among others, fuel this speculation. If topographic mixing is indeed important, it must be acknowledged that its geographic distribution, both laterally and vertically, is presently unknown. The vertical distribution of mixing plays a critical role in the Stommel Arons model of the ocean interior circulation. In recent numerical studies, Samelson demonstrates the extreme sensitivity of flow in the abyssal ocean to the spatial distribution of mixing. We propose to study the topographic mixing problem through an integrated program of modeling and observation. We focus on tidally forced mixing as the global energetics of this process have received (and are receiving) considerable study. Also, the well defined frequency of the forcing and the unique geometry of tidal scattering serve to focus the experiment design. The Hawaiian Ridge is selected as a study site. Strong interaction between the barotropic tide and the Ridge is known to take place. The goals of the Hawaiian Ocean Mixing Experiment (HOME) are to quantify the rate of tidal energy loss to mixing at the Ridge and to identify the mechanisms by which energy is lost and mixing generated. We are challenged to develop a sufficiently comprehensive picture that results can be generalized from Hawaii to the global ocean. To achieve these goals, investigators from five institutions have designed HOME, a program of historic data analysis, modeling and field observation. The Analysis and Modeling efforts support the design of the field experiments. As the program progresses, a global model of the barotropic (depth independent) tide, and two models of the baroclinic (depth varying) tide, all validated with near-Ridge data, will be applied, to reveal the mechanisms of tidal energy conversion along the Ridge, and allow spatial and temporal integration of the rate of conversion. Field experiments include a survey to identify "hot spots" of enhanced mixing and barotropic to baroclinic conversion, a Nearfield study identifying the dominant mechanisms responsible for topographic mixing, and a Farfield program which quantifies the barotropic energy flux convergence at the Ridge and the flux divergence associated with low mode baroclinic waves radiation. The difference is a measure of the tidal power available for mixing at the Ridge. Field work is planned from years 2000 through 2002, with analysis and modeling efforts extending through early 2006. If successful, HOME will yield an understanding of the dominant topographic mixing processes applicable throughout the global ocean. It will advance understanding of two central problems in ocean science, the maintenance of the abyssal stratification, and the dissipation of the tides. HOME data will be used to improve the parameterization of dissipation in models which presently assimilate TOPEX-POSEIDON observations. The improved understanding of the dynamics and spatial distribution of mixing processes will benefit future long-term programs such as CLIVAR.

Ray, Richard D.↗

DOE OSTI2021

Rotational Dynamics of an Amphidynamic Zirconium Metal–Organic Framework Determined by Dielectric Spectroscopy

A zirconium metal–organic framework with a difluorophenylene rotator bearing a permanent electric dipole of ~3.2 D was synthesized, and its rotational motion was analyzed by temperature- and frequency-dependent broadband dielectric spectroscopy. While solid-state NMR confirms fast rotation qualitatively, the dissipation factors measured between 113 and 153 K suggested an activation energy E a = 2.6 kcal/mol, but deviations from a single Debye relaxation suggested a dynamic process that cannot be accounted for by a well-defined potential with a single activation barrier. The dynamic heterogeneity of the dipolar rotor was confirmed by analysis in terms of a Cole–Cole relaxation, which suggested a mean barrier of ~1.9 kcal/mol, with a heterogeneity that decreases as temperature increases. Based on the single-crystal structure, we propose that the kinetic heterogeneity results from a temperature-dependent potential where rotation motion is mediated by the escape of the rotator from an energy well created by a double Ph–H···F–Ph hydrogen bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Flow dynamics and heat transfer in simplified battery energy storage systems with heated battery modules

Large-scale energy storage systems (ESSs) composed of batteries show promise in addressing current energy challenges, but dissipation of generated heat is important. Here, this paper focuses on buoyant convective flows in simplified ESS battery racks. Natural convection is not generally the primary cooling strategy but can be important in abnormal scenarios where there is module overheat or potentially thermal runaway. We use computational fluid dynamics to investigate the flow dynamics and heat transfer mechanisms in a simplified parameterized rack design. Despite its simplicity, this configuration produces many of the relevant features expected in real ESSs without details of module geometry or hardware, allowing broad conclusions independent of manufacture-specific designs. We start by providing visualizations of the flowfield and measurements of entrainment, heat flux, and pressure. To characterize the dependence on the system parameters, we develop an integral-scale analysis of the average temperature equation to highlight the dominant source terms. We use results from this analysis to derive a steady network model composed of simple algebraic expressions to provide first-order predictions of entrainment through the rack. The network model leads to a linear scaling of the Reynolds number based on convective mass flux with respect to the Grashof number based on the heat source. We deduce empirical relationships that relate the heat exchanged between modules using a surface-averaged Nusselt number as a function of the local Reynolds and Rayleigh numbers. Lastly, we investigate how space between the modules and rack in the spanwise direction creates flow bypass, resulting in different flow pathways.

Battery thermal management↗

DOE OSTI2020

$${\mathcal{P}}{\mathcal{T}}$$-Symmetric Effective Model for Nonequilibrium Phase Transitions in a Dissipative Fermionic Mott Insulator Chain

Nonequilibrium phase transitions in open dissipative systems can be described as instabilities in the spectra and wavefunctions of effective non-Hermitian Hamiltonians invariant under simultaneous parity (P) and time-reversal (T) transformations. The degree of non-Hermiticity reflects the strength of the external drive and dissipation, and the transition is described as a loss of the PT symmetry of the solutions corresponding to stationary low-drive dynamics. This approach has been successfully applied to spin, superconducting, and Mott insulator systems. However, the microscopic foundations for the employed phenomenological models are currently lacking. Here we propose a microscopic mechanism leading to the PT-symmetric effective model in the context of the nonequilibrium Mott transition in a dissipative Hubbard chain. Our model comprises a half-filled fermionic Hubbard chain subject to a constant electric field. The dissipation is introduced via the electron-phonon coupling. We obtain the explicit expressions for the non-Hermitian parameter in terms of the electron-phonon coupling strength and driving field. Analyzing the implications of microscopic model, we find a re-entrant Mott insulator with the increasing electric field for phonon density of states that increases slower than the square of the energy (such as in one or two dimensions), or varies non-monotonously with energy.

36 MATERIALS SCIENCE↗

NASA NTRS1987

Statistics and dynamics of persistent anomalies

Persistent anomalies with recurrent spatial patterns play an important role in the atmosphere's low-frequency variability. A connection between statistical and dynamical methods of description and prediction of persistent anomalies is established by computing and analyzing the empirical orthogonal functions (EOFs) in a simple deterministic model, on the one hand, and in Southern Hemisphere geopotential heights, on the other. The dynamical model is governed by the fully nonlinear, equivalent-barotropic vorticity equation on the sphere, with simplified forcing, dissipation and topography. The Southern Hemisphere data consist in gridded daily maps of 500 mb heights from June 1972 to July 1983. Two types of persistent anomalies appear in this time series, both having a strong wavenumber-three component; they differ by the value of the constant phase of this wave and by the strength of the wavenumber-one component. The first two EOFs bear a striking resemblance to these two patterns. It is concluded that the dynamical interpretation of EOFs is their pointing from the time mean to the most populated regions of the system's phase space. Pursuing this interpretation, a Markov-chain formulation of transitions from one persistent anomaly regime to another is introduced, and the implications for long-range forecasting are discussed.

Mo, Kingtse C.↗

NASA NTRS2019

Dynamic Smagorinsky Modeled Large-Eddy Simulations of Turbulence Using Tetrahedral Meshes

Eddy-resolving numerical computations of turbulent flows are emerging as viable alternatives to Reynolds Averaged Navier-Stokes (RANS) calculations for flows with an intrinsically steady mean state due to the advances in large-scale parallel computing. In these computations, medium to large turbulent eddies are resolved by the numerics while the smaller or subgrid scales are either modeled or taken care of by the inherent numerical dissipation. To advance the state of the art of unstructured-mesh turbulence simulation capabilities, large eddy simulations (LES) using the dynamic Smagorinsky model (DSM) on tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method. In contrast to what has been reported in the literature, the present implementation of dynamic models allows for active backscattering without any ad-hoc limiting of the eddy viscosity calculated from the subgrid-scale model. For the benchmark problems involving compressible isotropic turbulence decay as well as the shock/turbulent boundary layer interaction benchmark problems, no numerical instability associated with kinetic energy growth is observed and the volume percentage of the backscattering portion accounts for about 38-40% of the simulation domain. A slip-wall model in conjunction with the implemented DSM is used to simulate a relatively high Reynolds number Mach 2.85 turbulent boundary layer over a 30° ramp with several tetrahedral meshes and a wall-normal spacing of either Δγ& = 10 or Δγ& = 20. The computed mean wall pressure distribution, separation region size, mean velocity profiles, and Reynolds stress agree reasonably well with experimental data.

CESE method↗

NASA NTRSDate not supplied

Boundary-Layer Cloud Modeling Challenges on the North Slope of Alaska

The accurate modeling and prediction of cloud base heights is critical for energy balance calculations and aviation operations, alike. Low-level (i.e., boundary-layer) Arctic clouds can be difficult to model, making prediction of formation and dissipation challenging. Primarily mixed-phase, these clouds typically contain low quantities of supercooled liquid water and often slowly precipitate relatively small amounts of moderately and heavily rimed snow particles. While this appears to be the predominant cloudy state on the North Slope of Alaska (NSA), the delicate balance of microphysical, dynamical, radiative, surface coupling, and advective processes can rapidly shift to heavy snow (with various degrees of riming) or to a complete dissipation of the cloud layer without any precipitation, depending on the dominant processes. Here we strive to disentangle these various processes. First, we compare the predictive performances of four different numerical weather models in forecasting the presence and base-heights of low-level clouds: the High-Resolution Rapid Refresh - Alaska (HRRR-AK) model, the Polar Weather Research and Forecasting (Polar WRF) model, the Unified Model (UM), and the European Centre for Medium-range Weather Forecasting (ECMWF) model. Initial results comparing model output at two U.S. Department of Energy Atmospheric Radiation Measurement (AMT) NSA sites, during the fall season in 2019 and 2022, show that the UM slightly outperforms the HRRR-AK in terms of accurately forecasting the presence of a low-level cloud layer (89% of the time). All models have a significant bias of 300 to 800 meters in forecasting cloud base height (lower than is observed); however, the UM and ECMWF models have the lowest biases. Finally, a case study for a particularly challenging April 2017 thin-cloud event is presented, wherein we compare the performance of four different bulk microphysical parameterization schemes using a higher-resolution large eddy simulation (LES) model, the WRF-LES. Initial results show that the Thompson scheme was the only one able to reproduce and sustain a substantial supercooled liquid layer, but it was unable to reproduce the transition from a deep, liquid-rich cloud to a thin layer with moderately and heavily rimed precipitation. This is the first step in linking simulated LES-scale riming processes with those parameterized at a coarser mesoscale model scale. This has important implications for forecasting low-level clouds in an operational environment, given the efficiency of the riming process.

cloud base heights↗

DOE OSTI2024

Monodomain Liquid‐Crystal Elastomer Lattices for Broad Strain‐Rate Mechanical Damping

Designing structures that effectively dissipate energy across a broad range of mechanical loading rates, including those from compression, shock, and vibration, poses a significant engineering challenge. In this study, liquid‐crystal elastomers (LCEs), which possess anisotropic properties due to the alignment of their polymer backbone, are explored. As a result, LCEs exhibit a soft elastic response under mechanical loading, making them ideal for energy dissipation. Advances in additive manufacturing (AM) enable simple fabrication of foamlike dissipative structures with complex lattice geometries. Herein, direct ink write 3D printing, an extrusion AM method, is used to fabricate aligned, monodomain LCE lattice structures for broad strain‐rate mechanical damping. In this work, it is shown that these structures can dissipate strain energy in quasi‐static environments, comparable to traditional elastomeric lattices, and provide improved damping under high strain‐rate drop testing due to LCE soft elasticity. Additionally, under dynamic mechanical vibration, monodomain LCE lattices enhance damping at structural natural frequencies and provide high‐frequency attenuation. Finally, in these findings, a promising method is presented for fabricating mechanical damping devices that effectively dissipate energy across a broad range of loading rates.

3D printings↗

NASA NTRS2016

Improving Fidelity of Launch Vehicle Liftoff Acoustic Simulations

Launch vehicles experience high acoustic loads during ignition and liftoff affected by the interaction of rocket plume generated acoustic waves with launch pad structures. Application of highly parallelized Computational Fluid Dynamics (CFD) analysis tools optimized for application on the NAS computer systems such as the Loci/CHEM program now enable simulation of time-accurate, turbulent, multi-species plume formation and interaction with launch pad geometry and capture the generation of acoustic noise at the source regions in the plume shear layers and impingement regions. These CFD solvers are robust in capturing the acoustic fluctuations, but they are too dissipative to accurately resolve the propagation of the acoustic waves throughout the launch environment domain along the vehicle. A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed to improve such liftoff acoustic environment predictions. The framework combines the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin (DG) solver, Loci/THRUST, developed in the same computational framework. Loci/THRUST employs a low dissipation, high-order, unstructured DG method to accurately propagate acoustic waves away from the source regions across large distances. The DG solver is currently capable of solving up to 4th order solutions for non-linear, conservative acoustic field propagation. Higher order boundary conditions are implemented to accurately model the reflection and refraction of acoustic waves on launch pad components. The DG solver accepts generalized unstructured meshes, enabling efficient application of common mesh generation tools for CHEM and THRUST simulations. The DG solution is coupled with the CFD solution at interface boundaries placed near the CFD acoustic source regions. Both simulations are executed simultaneously with coordinated boundary condition data exchange.

Liever, Peter↗

DOE OSTI2026

Muon-induced fission of actinide nuclei

A negative muon captured by an actinide cascades down through the muonic atomic levels; deeply bound transitions can proceed via inverse internal conversion, depositing the muonic transition energy directly into the nucleus and, when the deposited energy exceeds the fission barrier, inducing prompt fission. Because the muon mean lifetime exceeds the saddle-to-scission timescale by orders of magnitude, the muon can survive the entire fission event as a 1⁢𝑠 spectator and ultimately attach to one or both of the emerging fragments. Its postscission attachment probability to the light fragment, 𝑃 𝐿 , can be used as a direct electromagnetic probe of fission dynamics on a timescale of 10 −21 s. In previous work, we introduced a three-dimensional lattice solution of the time-dependent Dirac equation coupled to the electromagnetic field generated by a fissioning nucleus and reported 𝑃 𝐿 for several actinides at a single dissipation strength. In this work, we extend that framework to a systematic survey of 232 Th , 238 U , and 240 Pu and implement a more realistic fission model which incorporates dynamic pairing correlations. We find that 𝑃 𝐿 falls steeply with the fragment charge asymmetry, a robust structural fingerprint of the fissioning system, while its dependence on nuclear dissipation is secondary and sensitive to the phenomenological friction prescription. These results establish 𝑃 𝐿 as a clean electromagnetic probe of fragment charge asymmetry and motivate a self-consistent, coordinate- and time-dependent treatment of nuclear dissipation as the natural next step.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1996

Computational Analysis of Static and Dynamic Behaviour of Magnetic Suspensions and Magnetic Bearings

Static modelling of magnetic bearings is often carried out using magnetic circuit theory. This theory cannot easily include nonlinear effects such as magnetic saturation or the fringing of flux in air-gaps. Modern computational tools are able to accurately model complex magnetic bearing geometries, provided some care is exercised. In magnetic suspension applications, the magnetic fields are highly three-dimensional and require computational tools for the solution of most problems of interest. The dynamics of a magnetic bearing or magnetic suspension system can be strongly affected by eddy currents. Eddy currents are present whenever a time-varying magnetic flux penetrates a conducting medium. The direction of flow of the eddy current is such as to reduce the rate-of-change of flux. Analytic solutions for eddy currents are available for some simplified geometries, but complex geometries must be solved by computation. It is only in recent years that such computations have been considered truly practical. At NASA Langley Research Center, state-of-the-art finite-element computer codes, 'OPERA', 'TOSCA' and 'ELEKTRA' have recently been installed and applied to the magnetostatic and eddy current problems. This paper reviews results of theoretical analyses which suggest general forms of mathematical models for eddy currents, together with computational results. A simplified circuit-based eddy current model proposed appears to predict the observed trends in the case of large eddy current circuits in conducting non-magnetic material. A much more difficult case is seen to be that of eddy currents in magnetic material, or in non-magnetic material at higher frequencies, due to the lower skin depths. Even here, the dissipative behavior has been shown to yield at least somewhat to linear modelling. Magnetostatic and eddy current computations have been carried out relating to the Annular Suspension and Pointing System, a prototype for a space payload pointing and vibration isolation system, where the magnetic actuator geometry resembles a conventional magnetic bearing. Magnetostatic computations provide estimates of flux density within airgaps and the iron core material, fringing at the pole faces and the net force generated. Eddy current computations provide coil inductance, power dissipation and the phase lag in the magnetic field, all as functions of excitation frequency. Here, the dynamics of the magnetic bearings, notably the rise time of forces with changing currents, are found to be very strongly affected by eddy currents, even at quite low frequencies. Results are also compared to experimental measurements of the performance of a large-gap magnetic suspension system, the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Eddy current effects are again shown to significantly affect the dynamics of the system. Some consideration is given to the ease and accuracy of computation, specifically relating to OPERA/TOSCA/ELEKTRA.

Britcher, Colin P.↗

NASA NTRS1992

Simulating the arcjet thruster using a nonlinear active load

Arcjet thrusters are known to have distinctly nonlinear electrical characteristics with negative incremental resistance throughout their normal operating region. The 1-kW arcjet simulator proposed here accurately mimics these characteristics statically, and dynamically up to 40 kHz. Testing with a power-processing unit (PPU) and the simulator, a 1-kW arcjet, and a resistive load bank established that the simulator accurately predicted the transient behavior of the PPU-arcjet combination, while the resistive load bank did not. The simulator uses resistors and insulated-gate bipolar transistors to dissipate the PPU output power; a nonlinear feedback is applied to produce the desired v-i characteristics. This feedback is analyzed dynamically. The arcjet simulator allows development of a PPU to a higher level of confidence before beginning testing at an arcjet facility.

Altenburger, Gene P.↗

DOE OSTI2024

A note on higher-order and nonlinear limiting approaches for continuously bounds-preserving discontinuous Galerkin methods

In Dzanic (2024), a limiting approach for high-order discontinuous Galerkin schemes was introduced which allowed for imposing constraints on the solution continuously (i.e., everywhere within the element). While exact for linear constraint functionals, this approach only imposed a sufficient (but not the minimum necessary) amount of limiting for nonlinear constraint functionals. This short note shows how this limiting approach can be extended to allow exactness for general nonlinear quasiconcave constraint functionals through a nonlinear limiting procedure, reducing unnecessary numerical dissipation. Finally, some examples are shown for nonlinear pressure and entropy constraints in the compressible gas dynamics equations, where both analytic and iterative approaches are used.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗

NASA NTRS1975

An approximation technique for predicting the transient response of a second order nonlinear equation

An approximation technique has been developed for determining the transient response of a nonlinear dynamic system. The nonlinearities in the system which has been considered appear in the system's dissipation function. This function was expressed as a second order polynomial in the system's velocity. The developed approximation is an extension of the classic Kryloff-Bogoliuboff technique. Two examples of the developed approximation are presented for comparative purposes with other approximation methods.

Laurenson, R. M.↗

NASA NTRS1978

Gate-assisted turn-off thyristor

1,000-volt, 200-ampere gate-assisted turn-off thyristor has been developed for power circuits requiring high efficiency, small size, and low weight. Design features include shunted cathode for high dV/dt capability. Cathode in interdigitated with dynamic gate for fast, low-loss switching. Operating frequency exceeds 20 kHz with overall energy dissipation of less that 12 mJ per pulse for typical 20-microsecond half-sine waveform. Device has turn-on time of 2 microseconds and turn-off time as short as 3 microseconds with only 2 amperes of gate drive.

Lowry, L. R.↗

NASA NTRS1980

Stellar convection. II - A multimode numerical solution for convection in spheres

A stable, equilibrium solution for convection in a self-gravitating sphere of Boussinesq fluid was computed by using a modal analysis in which the theta, phi dependence of the fluid is expanded in a set of 168 spherical harmonics. The temperature, velocity, and convective flux of the fluid as well as the kinetic and thermal energy spectra as functions of wavelength are computed. The spectra are found to be in agreement with both experimental observations and analytic scaling laws. The dynamics of the energy cascade is examined by computing the ratio of the amount of energy dissipation at a particular wavelength to the amount of energy produced at that same wavelength. It is found that there is only a slight cascade of kinetic energy to smaller wavelengths but a large cascade of thermal energy.

Marcus, P. S.↗

NASA NTRS1983

A mechanistic model of Eulerian, Lagrangian mean, and Lagrangian ozone transport by steady planetary waves

Ozone transport is calculated for steady, dissipative planetary waves using the Eulerian, Lagrangian mean, and residual circulation. A Lagrangian model of parcel dynamics is used to interpret planetary wave-photochemistry interaction. In chemically active regions the mean field ozone changes are found to be significant only where there are large gradients in chemical sources and sinks along parcel trajectories. The largest changes in the mean field are found in the lower stratosphere and are due to the Lagrangian mean advection. When the Lagrangian mean advection is approximated by the residual circulation, errors in the transport velocities as large as 30 pct may occur.

Rood, R. B.↗