Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 883 records · Page 49

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Demonstration of CFD to support the coupled analysis of a reactor pressure vessel subjected to pressurized thermal shock

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended duration or under accident condition. As one of the most significant components of a reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Aging degradation mechanisms (e.g. thermo-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The prediction of the temperature field is achieved by using computational fluid dynamics (CFD) simulation. In this report, a demonstration CFD standalone simulation is performed to support coupled analysis for Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS). The study use a simplified computational domain to represents a real RPV. The purpose of the study is to demonstrate the transient temperature response of RPV to ECC injection. The CFD model is built in a robust and efficient way for further coupled calculation. The next steps of this work, including the coupled thermal and tensor mechanics capabilities using Cardinal are expected to be complete by the end of FY21 for the demo problem. After this, into FY22, the capability will be demonstrated for a realistic RPV.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational Imaging for Intelligence in Highly Scattering Aerosols (Final Report)

Natural and man-made degraded visual environments pose major threats to national security. The random scattering and absorption of light by tiny particles suspended in the air reduces situational awareness and causes unacceptable down-time for critical systems and operations. To improve the situation, we have developed several approaches to interpret the information contained within scattered light to enhance sensing and imaging in scattering media. These approaches were tested at the Sandia National Laboratory Fog Chamber facility and with tabletop fog chambers. Computationally efficient light transport models were developed and leveraged for computational sensing. The models are based on a weak angular dependence approximation to the Boltzmann or radiative transfer equation that appears to be applicable in both the moderate and highly scattering regimes. After the new model was experimentally validated, statistical approaches for detection, localization, and imaging of objects hidden in fog were developed and demonstrated. A binary hypothesis test and the Neyman-Pearson lemma provided the highest theoretically possible probability of detection for a specified false alarm rate and signal-to-noise ratio. Maximum likelihood estimation allowed estimation of the fog optical properties as well as the position, size, and reflection coefficient of an object in fog. A computational dehazing approach was implemented to reduce the effects of scatter on images, making object features more readily discernible. We have developed, characterized, and deployed a new Tabletop Fog Chamber capable of repeatably generating multiple unique fog-analogues for optical testing in degraded visual environments. We characterized this chamber using both optical and microphysical techniques. In doing so we have explored the ability of droplet nucleation theory to describe the aerosols generated within the chamber, as well as Mie scattering theory to describe the attenuation of light by said aerosols, and correlated the aerosol microphysics to optical properties such as transmission and meteorological optical range (MOR). This chamber has proved highly valuable and has supported multiple efforts inclusive to and exclusive of this LDRD project to test optics in degraded visual environments. Circularly polarized light has been found to maintain its polarization state better than linearly polarized light when propagating through fog. This was demonstrated experimentally in both the visible and short-wave infrared (SWIR) by imaging targets made of different commercially available retroreflective films. It was found that active circularly polarized imaging can increase contrast and range compared to linearly polarized imaging. We have completed an initial investigation of the capability for machine learning methods to reduce the effects of light scattering when imaging through fog. Previously acquired experimental long-wave images were used to train an autoencoder denoising architecture. Overfitting was found to be a problem because of lack of variability in the object type in this data set. The lessons learned were used to collect a well labeled dataset with much more variability using the Tabletop Fog Chamber that will be available for future studies. We have developed several new sensing methods using speckle intensity correlations. First, the ability to image moving objects in fog was shown, establishing that our unique speckle imaging method can be implemented in dynamic scattering media. Second, the speckle decorrelation over time was found to be sensitive to fog composition, implying extensions to fog characterization. Third, the ability to distinguish macroscopically identical objects on a far-subwavelength scale was demonstrated, suggesting numerous applications ranging from nanoscale defect detection to security. Fourth, we have shown the capability to simultaneously image and localize hidden objects, allowing the speckle imaging method to be effective without prior object positional information. Finally, an interferometric effect was presented that illustrates a new approach for analyzing speckle intensity correlations that may lead to more effective ways to localize and image moving objects. All of these results represent significant developments that challenge the limits of the application of speckle imaging and open important application spaces. A theory was developed and simulations were performed to assess the potential transverse resolution benefit of relative motion in structured illumination for radar systems. Results for a simplified radar system model indicate that significant resolution benefits are possible using data from scanning a structured beam over the target, with the use of appropriate signal processing.

58 GEOSCIENCES↗

Mechanisms of Pacific decadal variability in ESMs: the roles of stochastic forcing, feedbacks & external forcing

This project explored how and why climate patterns in the Pacific Ocean—specifically Pacific Decadal Variability (PDV)—change over time, and how well leading Earth System Models (ESMs) capture these changes under different climate scenarios. By analyzing large climate model datasets and running new simulations, the team investigated how ocean and atmosphere processes interact to shape PDV. The study tested a leading scientific hypothesis using both simplified and complex models, helping to uncover the key mechanisms that drive long-term climate swings in the Pacific. This work improves our ability to detect, predict, and understand how human-driven climate change may influence major ocean-climate patterns in the future.

54 ENVIRONMENTAL SCIENCES↗

The 15 meter hoop-column antenna dynamics: Test and analysis

A 15 meter model of the hoop-column antenna concept has been vibration tested for model characterization and analytical model verification. Linear finite element analysis predicted the global vibration frequencies accurately. Good agreement between analysis and test data was obtained only after the analytical model was refined using static test data. As structures become more flexible, structural properties determined from static data become more accurate and should be used to update analytical models. Global vibration modes are not significantly affected by the surface mesh which permits simplified analytical models to be used for prediction of global behavior. These reduced models are believed sufficient for preliminary design and controls simulations where only global behavior is desired. The mesh modes were highly damped due to the knit mesh used for the reflector surface. These modes were also highly coupled and very difficult to measure in the laboratory. The inability to fully characterize the antenna mesh modes in the laboratory indicates robust methods for active surface vibration suppression will be needed. Fortunately, the surface mesh exhibits high passive damping which should be beneficial to active control systems.

Belvin, W. Keith↗

Earth observing system instrument pointing control modeling for polar orbiting platforms

An approach to instrument pointing control performance assessment for large multi-instrument platforms is described. First, instrument pointing requirements and reference platform control systems for the Eos Polar Platforms are reviewed. Performance modeling tools including NASTRAN models of two large platforms, a modal selection procedure utilizing a balanced realization method, and reduced order platform models with core and instrument pointing control loops added are then described. Time history simulations of instrument pointing and stability performance in response to commanded slewing of adjacent instruments demonstrates the limits of tolerable slew activity. Simplified models of rigid body responses are also developed for comparison. Instrument pointing control methods required in addition to the core platform control system to meet instrument pointing requirements are considered.

Briggs, H. C.↗

Theoretical studies of solar lasers and converters

The previously constructed one dimensional model for the simulated operation of an iodine laser assumed that the perfluoroalkyl iodide gas n-C3F7I was incompressible. The present study removes this simplifying assumption and considers n-C3F7I as a compressible fluid.

Heinbockel, John H.↗

Development of comprehensive numerical schemes for predicting evaporating gas-droplets flow processes of a liquid-fueled combustor

An existing Computational Fluid Dynamics code for simulating complex turbulent flows inside a liquid rocket combustion chamber was validated and further developed. The Advanced Rocket Injector/Combustor Code (ARICC) is simplified and validated against benchmark flow situations for laminar and turbulent flows. The numerical method used in ARICC Code is re-examined for incompressible flow calculations. For turbulent flows, both the subgrid and the two equation k-epsilon turbulence models are studied. Cases tested include idealized Burger's equation in complex geometries and boundaries, a laminar pipe flow, a high Reynolds number turbulent flow, and a confined coaxial jet with recirculations. The accuracy of the algorithm is examined by comparing the numerical results with the analytical solutions as well as experimented data with different grid sizes.

Chen, C. P.↗

Ground effects on the low-speed aerodynamics of a powered, generic hypersonic configuration

A study was undertaken in the NASA Langley 14- by 22-foot subsonic tunnel to determine the low-speed aerodynamic characteristics of a powered, generic, hypersonic configuration in ground effect. The model was a simplified configuration consisting of a triangular wedge forebody, a rectangular mid-section which housed the flow through, an ejector type propulsion simulation system, and a rectangular wedge afterbody. Additional model components included a delta wing, a rectangular wedge forebody, inlet fences, exhaust flow deflectors, and afterbody fences. Aerodynamic force and moment data were obtaind over an angle of attack range from -4 to 18 degrees while model height above the tunnel floor was varied from 1/4 inch to 6 feet. Variations in freestream dynamic pressure, from 10 psf to 80 psf, and engine ejector pressure yielded a range of thrust coefficients from 0 to 0.8.

Gatlin, Gregory M.↗

The performance of trellis-coded MDPSK with multiple symbol detection

The idea of using a multiple (more than two) symbol observation interval to improve error probability performance is applied to differential detection of trellis-coded multiple phase-shift keying (MPSK) over an additive white Gaussian noise (AWGN) channel. An equivalent Euclidean distance measure per trellis branch is determined for this detection scheme. This is used to define an augmented (larger multiplicity) trellis code whose distance measure is the conventional squared Euclidean distance typical of conventional trellis-coded modulation on the AWGN. Such an augmented multiple trellis code is a convenient mathematical tool for simplifying the analysis. Results are obtained by a combination of analysis (upper Chernoff bounds and asymptotic large-SNR approximations) and computer simulation. It is shown that only a slight increase (e.g., one symbol) in the length of the observation interval will provide a significant improvement in bit error probability performance.

Divsalar, Dariush↗

A simplified wind vector algorithm for satellite scatterometers

A simplified algorithm is derived for retrieving wind vectors from microwave scatterometer observations. The azimuthal dependence of the sea-surface radar cross section is modeled by a double-cosine function, rather than the traditional second-order cosine expansion. The algorithm is tested using the aircraft scatterometer observations obtained during the AAFE-RADSCAT Experiment in 1975 and 1976. There is little difference between the performance of the simplified algorithm and the more involved least-squares searching algorithm. The AAFE-RADSCAT aircraft observations are sampled so as to simulate the three-look satellite scatterometer NSCAT that will be launched on the Japanese Advanced Earth Observing Satellite. The results of the retrievals indicate that it is probably not possible to uniquely determine the wind direction using just NSCAT observations because of a 180 degree ambiguity. However, if forecast models can predict the wind direction to an accuracy of +/- 90 deg, thereby eliminating the ambiguity, then the results indicate that NSCAT can determine the wind direction to an accuracy of about +/- 20 deg. Better performance is obtained when the three observations are the same polarization (either v-pol or h-pol), as compared to using a mix of v-pol and h-pol observations.

Wentz, Frank J.↗

Reactions during conversion of simplified low-activity waste glass feeds

The mechanisms that affect the incorporation of 99Tc, a volatile radioactive component of concern, into glass melt during vitrification of low-activity waste (LAW) are being investigated to develop the method to increase the retention of 99Tc in glass waste form. Previous studies with simulated LAW glass feeds (slurry mixture of liquid waste and chemical/mineral additives) demonstrated that the early stage feed-to-glass conversion reactions below 800°C are critical for the Re (used as a nonradioactive surrogate of 99Tc) retention in glass. To examine the effect of feed composition on the feed-to-glass conversion reactions, simplified systems containing major LAW components (NaNO3 and NaOH) and representative additive components (SiO2 and H3BO3) were designed and tested. The ratio of H3BO3 to NaNO3 was varied in three-component system without NaOH and that of NaOH to NaNO3 was varied in the four-component system at a fixed H3BO3 to NaNO3 ratio. As a first step of testing with simplified feeds, this study applied thermal analyses and phase characterization of the reacting feeds, which were performed without the addition of Re, to investigate the evolution of salt phases during slurry drying process and upon heating of dried feeds.

Jin, Tongan↗

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗

High-Temperature Gas-Cooled Reactors Multiphysics Simulation Demonstration and Code Validation

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor applications. The first part of this effort focuses on the running-in verification of Griffin’s multiphysics capabilities, specifically for simulating the evolution of pebble-bed reactor cores from startup to equilibrium. Since Fiscal Year 2024, improvements and enhancements have been implemented in Griffin, including simplifying the process to specify streamlines and developing the online cross-section generation capability. In the absence of validation data, code-to-code comparisons are conducted with kugelpy, showing good agreement for integral quantities like k-eff predictions and predictions for maximum power density. However, accuracy issues are noted for more detailed quantities like the spatial distribution of fission rate densities which will require further work to address. The second part of this report presents an improved System Analysis Module (SAM) core channel model where the effects of cross flow are considered during the pressurized loss of forced cooling transient, resulting in an improved agreement of the predicted pebble temperature with respect to the predictions from the SAM 2D porous media model. Additionally, the wall channeling effect due to variable porosity at the near wall region of the core is also investigated. Furthermore, to demonstrate Griffin’s online cross-section generation capability, a Multiphysics simulation is performed by coupling Griffin to the SAM core channel model. In the third part of the report, as a part of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA) thermal-hydraulic code validation benchmark activity for a high-temperature gas-cooled reactor, the High Temperature Test Facility (HTTF) is investigated first using the NekRS computational fluid dynamics (CFD) code to study the flow mixing phenomenon in the lower plenum of the facility. Then, code-to-code and code-to-data comparisons are performed for Test PG27, which is a pressurized conduction cooldown (PCC) test, using five different codes by six organizations from five countries. The different simulations show good agreements in terms of the general trend but there are differences in some results such as the peak temperatures of different regions and heat removal rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Simplified Model of VIPER Thermal Management System. Part II: Integrated Vehicle

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) thermal management system (TMS) relies on four loop heat pipes (LHPs) to transport electronic waste heat to the vehicle cooling radiative surface and avoid overheating. The TMS has also ten constant conductance heat pipes (CCHPs) dedicated to balance the thermal load within the internal environment where the avionics boxes are mounted, also called warm electronic box (WEB), and to transport the heat from two of the science payload instruments. The TMS also uses two thermal straps to thermally link the batteries to the WEB. These thermal components, in addition to heaters, thermostat, multi-layer insulation (MLIs), and isolators forms the core of the VIPER TMS. The complex heat transport balance managed by the TMS is challenging to characterize and model. The more fidelity and granularity of a model, the more costly the computational resources needed and the longer the simulation and modeling time. When the priority is to provide quick but reliable assessments of the thermal performance or real time thermal feedback for training of console operators, simplified modeling tools are needed. To satisfy that need, this paper describes the effort to develop and correlate a model of VIPER TMS based on control volume approach. The correlation effort in particular focuses on hibernation, cold thermal balance, and hot thermal balance data from the integrated vehicle thermal vacuum (TVAC) test. Thus, the correlated model captures the heat leaks during hibernations and the performance at two extremes, bounding, operating scenarios.

Loop Heat Pipe↗

A Simplified Model of VIPER Thermal Management System. Part II: Integrated Vehicle

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) thermal management system (TMS) relies on four loop heat pipes (LHPs) to transport electronic waste heat to the vehicle cooling radiative surface and avoid overheating. The TMS has also ten constant conductance heat pipes (CCHPs) dedicated to balance the thermal load within the internal environment where the avionics boxes are mounted, also called warm electronic box (WEB), and to transport the heat from two of the science payload instruments. The TMS also uses two thermal straps to thermally link the batteries to the WEB. These thermal components, in addition to heaters, thermostat, multi-layer insulation (MLIs), and isolators forms the core of the VIPER TMS. The complex heat transport balance managed by the TMS is challenging to characterize and model. The more fidelity and granularity of a model, the more costly the computational resources needed and the longer the simulation and modeling time. When the priority is to provide quick but reliable assessments of the thermal performance or real time thermal feedback for training of console operators, simplified modeling tools are needed. To satisfy that need, this paper describes the effort to develop and correlate a model of VIPER TMS based on control volume approach. The correlation effort in particular focuses on hibernation, cold thermal balance, and hot thermal balance data from the integrated vehicle thermal vacuum (TVAC) test. Thus, the correlated model captures the heat leaks during hibernations and the performance at two extremes, bounding, operating scenarios.

Thermal Modeling↗

Diagnostic analysis of two-dimensional monthly average ozone balance with Chapman chemistry

Chapman chemistry has been used in a two-dimensional model to simulate ozone balance phenomenology. The similarity between regions of ozone production and loss calculated using Chapman chemistry and those computed using LIMS and SAMS data with a photochemical equilibrium model indicate that such simplified chemistry is useful in studying gross features in stratospheric ozone balance. Net ozone production or loss rates are brought about by departures from the photochemical equilibrium (PCE) condition. If transport drives ozone above its PCE condition, then photochemical loss dominates production. If transport drives ozone below its PCE condition, then photochemical production dominates loss. Gross features of ozone loss/production (L/P) inferred for the real atmosphere from data are also simulated using only eddy diffusion. This indicates that one must be careful in assigning a transport scheme for a two-dimensional model that mimics only behavior of the observed ozone L/P.

Stolarski, Richard S.↗