Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D temperature 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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Magneto-electroluminescence response in 2D and 3D hybrid organic–inorganic perovskite light emitting diodes

We have studied the magneto-electroluminescence (MEL) response in light emitting diodes based on 2D-(PEA) 2 PbI 4 and 3D-MAPbI 3 hybrid organic–inorganic perovskites at cryogenic temperatures. We found that the MEL is negative, i.e., the EL decreases with the applied field strength, B. In addition, the MEL(B) response has a Lorentzian line shape whose width depends on the perovskite used. We interpret the MEL(B) response using the “Δg mechanism” in which the spin of the injected electron–hole (e–h) pairs oscillates between singlet and triplet configurations due to different precession frequencies of the electron and hole constituents that originate from the difference, Δg, in the electron and hole gyromagnetic constants, g. In this model, the MEL(B) linewidth is inversely proportional to the spin lifetime and Δg. The model used is validated by directly measuring the spin lifetime of photogenerated e–h pairs using the circularly polarized pump–probe transmission technique with 100 fs resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Methodology for Thermodynamic Analysis Coupled with Computational Fluid Dynamics Modeling for Casting a Novel Aluminum–Cerium Alloy

In this paper, we present the results of computational fluid dynamics (CFD) analysis to assess castability, porosity, molten metal fluidity and other technological properties of a new Al-Ce alloy. A thermodynamic analysis of a new Al-Ce alloy is done to obtain the physical properties as a function of temperature across both solid and liquid phases. These properties are then used to build the CFD of model the full casting process from initial pouring through to final solidification of the part, here an example of a heavy duty mor mount is used. While several excellent commercial CFD codes exists this process shows that CFD can be used to assess the capability of the alloy to properly fill the mold, as well as give predictions where scattered porosity or large-scale defects may occur in the casting. Further, unlike the commercially available software (e.g., ProCAST, SOLIDCast , MAGMASOFT®) the complex 3D-analysis of the stress / strain fields in cast parts is not performed at this time. However, the availability of free software for assessing the required thermodynamic / thermo-physical properties of new alloys (OpenCALPHAD) and CFD codes such as OpenFOAM® makes the developed option attractive and economical, especially for the analysis of new Al-Ce alloys, for which the available data does not exist.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pressure tuning of Kitaev spin liquid candidate Na 3 Co 2 SbO 6

The search for Kitaev’s quantum spin liquid in real materials has recently expanded with the prediction that honeycomb lattices of divalent, high-spin cobalt ions could host the dominant bond-dependent exchange interactions required to stabilize the elusive entangled quantum state. The layered honeycomb Na 3 Co 2 SbO 6 has been singled out as a leading candidate provided that the trigonal crystal field acting on Co 3d orbitals, which enhances non-Kitaev exchange interactions between $J$ eff = $\frac{1}{2}$ spin-orbital pseudospins, is reduced. Here we show that applied pressure leads to anisotropic compression of the layered structure, significantly reducing the trigonal distortion of CoO 6 octahedra. Ferromagnetic correlations between pseudospins are enhanced in the spin-polarized (3 Tesla) phase up to about 60 GPa. Higher pressures drive a high-spin to low-spin transition destroying the $J$ eff = $\frac{1}{2}$ moments required to map the spin Hamiltonian into Kitaev’s model. The spin transition strongly suppresses the low-temperature magnetic susceptibility and appears to stabilize a paramagnetic phase driven by frustration. The possible emergence of frustrated magnetism of localized $S$ = $\frac{1}{2}$ moments opens the door for exploration of novel magnetic quantum states in compressed honeycomb lattices of divalent cobaltates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2016

The Inner Edge of the Habitable Zone for Synchronously Rotating Planets Around Low-Mass Stars Using General Circulation Models

Terrestrial planets at the inner edge of the habitable zone (HZ) of late-K and M-dwarf stars are expected to be in synchronous rotation, as a consequence of strong tidal interactions with their host stars. Previous global climate model (GCM) studies have shown that, for slowly rotating planets, strong convection at the substellar point can create optically thick water clouds, increasing the planetary albedo, and thus stabilizing the climate against a thermal runaway. However these studies did not use self-consistent orbital/rotational periods for synchronously rotating planets placed at different distances from the host star. Here we provide new estimates of the inner edge of the HZ for synchronously rotating terrestrial planets around late-K and M-dwarf stars using a 3D Earth-analog GCM with self-consistent relationships between stellar metallicity, stellar effective temperature, and the planetary orbital/rotational period. We find that both atmospheric dynamics and the efficacy of the substellar cloud deck are sensitive to the precise rotation rate of the planet. Around mid-to-late M-dwarf stars with low metallicity, planetary rotation rates at the inner edge of the HZ become faster, and the inner edge of the HZ is farther away from the host stars than in previous GCM studies. For an Earth-sized planet, the dynamical regime of the substellar clouds begins to transition as the rotation rate approaches ∼10 days. These faster rotation rates produce stronger zonal winds that encircle the planet and smear the substellar clouds around it, lowering the planetary albedo, and causing the onset of the water-vapor greenhouse climatic instability to occur at up to ∼25% lower incident stellar fluxes than found in previous GCM studies. For mid-to-late M-dwarf stars with high metallicity and for mid-K to early-M stars, we agree with previous studies.

Planet and satellite atmospheres↗
NASA NTRS1997

Structure and Dynamics of Coronal Plasmas

During the past year this grant has funded research in the interaction between magnetic fields and the hot plasma in the solar outer atmosphere. The following is a brief summary of the published papers, abstracts and talks which have been supported. The paper 'Coronal Structures Observed in X-rays and H-alpha Structures' was published in the Kofu Symposium proceedings. The study analyzes cool and hot behavior of two x-ray events, a small flare and a surge. We find that a large H-alpha surge appears in x-rays as a very weak event, while a weak H-alpha feature corresponds to the brightest x-ray emission on the disk at the time of the observation. Calculations of the heating necessary to produce these signatures, and implications for the driving and heating mechanisms of flares vs. surges are presented. A copy of the paper is appended to this report. The paper 'Differential Magnetic Field Shear in an Active Region' has been published in The Astrophysical Journal. We have compared the 3D extrapolation of magnetic fields with the observed coronal structure in an active region. Based on the fit between observed coronal structure throughout the volume of the region and the calculated magnetic field configurations, we propose a differential magnetic field shear model for this active region. The decreasing field shear in the outer portions of the AR may indicate a continual relaxation of the magnetic field with time, corresponding to a net transport of helicity outward. The paper 'Difficulties in Observing Coronal Structure' has been accepted for publication in the journal Solar Physics. In this paper we discuss the evidence that the temperature and density structure of the corona are far more complicated than had previously been thought. The discussion is based on five studies carried out by our group on coronal plasma properties, showing that any one x-ray instrument does see all of the plasma present in the corona, that hot and cool material may appear to be co-spatial at a given location in the corona, and that simple magnetic field extrapolations provide only a poor fit to the observed structure. A copy of the paper is appended to this report.

Golub, Leon↗
DOE OSTI2021

Three-dimensional diagnostics and measurements of inertial confinement fusion plasmas

Recent inertial confinement fusion measurements have highlighted the importance of 3D asymmetry effects on implosion performance. One prominent example is the bulk drift velocity of the deuterium–tritium plasma undergoing fusion (“hotspot”), v HS . Upgrades to the National Ignition Facility neutron time-of-flight diagnostics now provide v HS to better than 1 part in 10 4 and enable cross correlations with other measurements. Here we present the impact of v HS on the neutron yield, downscatter ratio, apparent ion temperature, electron temperature, and 2D x-ray emission. The necessary improvements to diagnostic suites to take these measurements are also detailed. The benefits of using cross-diagnostic analysis to test hotspot models and theory are discussed, and cross-shot trends are shown.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Numerical Simulations of TRMM LBA, TOGA, COARE, GATE, ARM and PRESTORM Convective Systems: Sensitivity tests on Microphysical Processes

The 3D Goddard Cumulus Ensemble (GCE) model was utilized to examine the behavior and response of simulated deep tropical cloud systems that occurred over the west Pacific warm pool region, the Atlantic ocean and the central United States. The periods chosen for simulation were convectively active periods during TOGA-COARE (February 22 1993, December 11-17, 1992; December 19-28, February 9-13, 1993), GATE (September 4, 1974), LBA (January 26 and February 23, 1998), ARM (1997 IOP) and PRESTORM (June 11, 1985). We will examine differences in the microphysics for both warm rain and ice processes (evaporation /sublimation and condensation/ deposition), Q1 (Temperature), Q2 (Water vapor) and Q3 (momentum both U and V) budgets for these three convective events from different large-scale environments. The contribution of stratiform precipitation and its relationship to the vertical shear of the large-scale horizontal wind will also be examined. New improvements to the GCE model (i.e., microphysics: 4ICE two moments and 3ICE one moment; advection schemes) as well as their sensitivity to the model results will be discussed. Preliminary results indicated that various microphysical schemes could have a major impact on stratiform formation as well as the size of convective systems. However, they do not change the major characteristics of the convective systems, such as: arc shape, strong rotational circulation on both ends of system, heavy precipitation along the leading edge of systems.

Tao, W.-K.↗
DOE OSTI2025

Fuel performance analysis of fully-resolved TRISO compact

The TRi-structural ISOtropic (TRISO) fuel multilayered coating structure offers multiple barriers to fission product release, enhancing safety and performance. The heterogeneous nature of TRISO fuel compacts, comprising thousands of randomly distributed coated fuel particles embedded in a graphite matrix, creates intricate stress fields and thermal gradients that cannot be accurately modeled using simplified one-dimensional or homogenized approaches. Consequently, three-dimensional modeling enables the prediction of fuel compact dimensional changes, internal pressure buildup, and fission product transport pathways under diverse irradiation and thermal conditions. This capability facilitates detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating failures, which directly impact fuel performance and safety margins. This capability is particularly critical for advanced reactors, such as high-temperature gas-cooled reactors and other Generation IV reactor designs where TRISO fuel operates at elevated temperatures and burn-up levels. This work introduces a novel method to generate an optimized packing of TRISO compacts and a complete 3D mesh with random distribution of TRISO particles, which are discretized into each coating component layer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Resolving atomic transport through iron oxide under irradiation using isotopic tracers

Structural materials in nuclear environments are exposed to complex and harsh conditions, including high temperatures, corrosive solutions, and irradiation. Material stability and likewise degradation are fundamentally controlled by atomic transport mechanisms. Here, we investigate the influence of radiation-induced, non-equilibrium point defect populations on mass transport within model oxide systems using isotopic tracers and atom probe tomography. An isotopically labelled hematite (a-Fe2O3) film was synthesized with 16O- and 18O-rich layers (>30nm and 10 nm, respectively) sequentially deposited via molecular beam epitaxy. One region of the film was proton irradiated to ~0.1 dpa at 450°C, and regions outside the beam spot provided a non-irradiated control of thermal annealing only. 3D APT was used to directly visualize the elemental and O-isotopic redistribution following both thermal and combined thermal and irradiation conditions. Complementary scanning transmission electron microscopy (STEM) was also employed to study the microstructural evolution. Error function fits of the 18O-concentration profiles enabled quantification of diffusion coefficients.

isotopic tracers, diffusion, proton irradiation, h↗
DOE OSTI2023

Long‐Range Electron Transport Rates Depend on Wire Dimensions in Cytochrome Nanowires

Abstract The ability to redirect electron transport to new reactions in living systems opens possibilities to store energy, generate new products, or probe physiological processes. Recent work by Huang et al. showed that 3D crystals of small tetraheme cytochromes (STC) can transport electrons over nanoscopic to mesoscopic distances by an electron hopping mechanism, making them promising materials for nanowires. However, fluctuations at room temperature may distort the nanostructure, hindering efficient electron transport. Classical molecular dynamics simulations of these fluctuations at the nano‐ and mesoscopic scales allowed us to develop a graph network representation to estimate maximum electron flow that can be driven through STC wires. In longer nanowires, transient structural fluctuations at protein‐protein interfaces tended to obstruct efficient electron transfer, but these blockages are ameliorated in thicker crystals where alternative electron transfer pathways become more efficient. The model implies that more flexible proteinprotein interfaces limit the required minimum diameter to carry currents commensurate with conventional electronics.

Chemistry↗
DOE OSTI2025

Boosting H I -Galaxy Cross-Clustering Signal through Higher-Order Cross-Correlations

After reionization, neutral hydrogen (${\rm H\, \small {I}}$) traces the large-scale structure (LSS) of the Universe, enabling ${\rm H\, \small {I}}$ intensity mapping (IM) to capture the LSS in 3D and constrain key cosmological parameters. We present a new framework utilizing higher-order cross-correlations to study ${\rm H\, \small {I}}$ clustering around galaxies, tested using real-space data from the IllustrisTNG300 simulation. This approach computes the joint distributions of k-nearest neighbor (kNN) optical galaxies and the ${\rm H\, \small {I}}$ brightness temperature field smoothed at relevant scales (the kNN-field framework), providing sensitivity to all higher-order cross-correlations, unlike two-point statistics. To simulate ${\rm H\, \small {I}}$ data from actual surveys, we add random thermal noise and apply a simple foreground cleaning model, filtering out Fourier modes of the brightness temperature field with k ∥ < k min,∥ . Under current levels of thermal noise and foreground cleaning, typical of a Canadian Hydrogen Intensity Mapping Experiment (CHIME)-like survey, the ${\rm H\, \small {I}}$-galaxy cross-correlation signal in our simulations, using the kNN-field framework, is detectable at >30σ across r = [3, 12] h –1 Mpc. In contrast, the detectability of the standard two-point correlation function (2PCF) over the same scales depends strongly on the foreground filter: a sharp k ∥ filter can spuriously boost detection to 8σ due to position-space ringing, whereas a less sharp filter yields no detection. Nonetheless, we conclude that kNN-field cross-correlations are robustly detectable across a broad range of foreground filtering and thermal noise conditions, suggesting their potential for enhanced constraining power over 2PCFs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

4D structural changes and pore network model of biomass during pyrolysis

Biochar is an engineered carbon-rich substance used for soil improvement, environmental management, and other diverse applications. To date, the understanding of how biomass affects biochar microstructure has been limited due to the complexity of analysis involved in tracing the changes in the physical structure of biomass as it undergoes thermochemical conversion. In this study, we used synchrotron x-ray micro-tomography to visualize changes in the internal structure of biochar from diverse feedstock (miscanthus straw pellets, wheat straw pellets, oilseed rape straw pellets, and rice husk) during pyrolysis by collecting a sequence of 3D scans at 50 °C intervals during progressive heating from 50 °C to 800 °C. The results show a strong dependence of biochar porosity on feedstock as well as pyrolysis temperature, with observed porosity in the range of 7.41–60.56%. Our results show that the porosity, total surface area, pore volume, and equivalent diameter of the largest pore increases with increasing pyrolysis temperature up to about 550 °C. The most dramatic development of pore structure occurred in the temperature range of 350–450 °C. This understanding is pivotal for optimizing biochar’s properties for specific applications in soil improvement, environmental management, and beyond. By elucidating the nuanced variations in biochar’s physical characteristics across different production temperatures and feedstocks, this research advances the practical application of biochar, offering significant benefits in agricultural, environmental, and engineering contexts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Electronic structure, dimer physics, orbital-selective behavior, and magnetic tendencies in the bilayer nickelate superconductor $\mathrm{La_3Ni_2O_7}$under pressure

Motivated by the recently reported high-temperature superconductivity in the bilayer La 3 Ni 2 O 7 (LNO) under pressure, we comprehensively study this system using ab initio techniques. The Ni 3d orbitals have a large bandwidth at ambient pressure, increasing by ~22% at 29.5 GPa. Without electronic interactions, the Ni d 3z 2 -r 2 orbitals form a bonding-antibonding molecular orbital state via the O p z inducing a “dimer” lattice in the LNO bilayers. The Fermi surface consists of two-electron sheets with mixed e g orbitals and a hole pocket defined by the d 3z 2 -r 2 orbital, suggesting a Ni two-orbital minimum model. Different from the infinite-layer nickelate, we obtained a large interorbital hopping between d 3z 2 -r 2 and d x 2 -y 2 states in LNO, caused by the ligand “bridge” of in-plane O p x or p y orbitals connecting those two orbitals, inducing d-p σ-bonding characteristics. The competition between the intraorbital and interorbital hoppings leads to an interesting dominant spin stripe (π,0) order because of bond ferromagnetic tendencies via the recently discussed “half-empty” mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗
DOE OSTI2022

Resolution Study of Thermonuclear Initiation in White Dwarf Tidal Disruption Events

Abstract We study the initiation of thermonuclear detonations in tidally disrupted white dwarf stars by intermediate-mass (10 3 M ⊙ ) black holes. The length scales required to resolve the initiation mechanism are not easily reached in 3D, so instead we have devised 2D proxy models, which, together with a logarithmic gridding strategy, can adequately capture detonation wave fronts as material undergoes simultaneous compression and stretching from tidal forces. We consider 0.15 and 0.6 M ⊙ white dwarf stars parameterized by tidal strengths in the range β = 4–23. High spatial resolution elucidates the manner and conditions leading to thermonuclear detonation, linking the initiation sequence to stellar composition and tidal strength. All of our models suffer sustained detonations triggered by a combination of adiabatic compression, mild thermonuclear preconditioning, and collisional heating, in degrees depending primarily on tidal strength. We find that many diagnostics, such as temperature, total released energy, and iron-group products, are fairly well converged (better than 10%) at resolutions below 10 km along the scale height of the orbital plane. The exceptions are intermediate-mass transients like calcium, which remain uncertain up to factors of 2, even at 1 km resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Radiative effects of polar stratospheric clouds during the Airborne Antarctic Ozone Experiment and the Airborne Arctic Stratospheric Expedition

Results are presented of a study of the radiative effects of polar stratospheric clouds during the Airborne Antarctic Ozone Experiment (AAOE) and the Airborne Arctic Stratospheric Expedition (AASE) in which daily 3D Type I nitric acid trihydrate (NAT) and Type II water ice polar stratospheric clouds (PSCs) were generated in the polar regions during AAOE and the AASE aircraft missions. Mission data on particular composition and size, together with NMC-analyzed temperatures, are used. For AAOE, both Type I and Type II clouds were formed for the time period August 23 to September 17, after which only Type I clouds formed. During AASE, while Type I clouds were formed for each day between January 3 and February 10, Type II clouds formed on only two days, January 24 and 31. Mie theory and a radiative transfer model are used to compute the radiative heating rates during the mission periods, for clear and cloudy lower sky cases. Only the Type II water ice clouds have a significant radiative effect, with the Type I NATO PSCs generating a net heating or cooling of 0.1 K/d or less.

Rosenfield, Joan E.↗