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 613 records · Page 34

Facile Electrodeposition and Aging to Generate 3-Dimensional α-MnO 2 Battery Cathodes

Conventional tape casting forms 2-dimensional (2D) electrodes containing active material, conductive additive, and binder with restricted ion access as electrodes increase in thickness. To improve the transport properties, 3D architectures were developed using electrodeposition to ensure contact between the active material with the substrate, and provide enhanced electrolyte access. This paper investigates electrodeposition of cryptomelane ( α -MnO 2 ) as a model cathode material to efficiently accommodate (de)lithation and increase areal capacity vs conventional 2D coatings. Electodeposited samples on titantium (Ti) foil substrates were characterized using X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), and scanning electron microscopy (SEM) and show a linear increase of the average oxidation of Mn (3.5–3.8) and active mass loading (1.27–9.9 mg) with deposition and aging times (0–120 min). The initial deposition is amorphous and forms the crystalline material during the elevated temperature aging step. The active material, α -MnO 2 , was also deposited on C-cloth and these cathodes at deposition times of 3, 6, and 9 min deliver 9, 36, and 69% higher areal capacities, respectively, at 0.2 mA cm −2 compared to conventional 2D electrodes with a mass loading equal to the 3 min sample. These results demonstrate the benefit of α -MnO 2 within a porous architecture providing enhanced transport properties.

Electrochemistry↗

Sodium-Cooled Fast Reactor Reference Plant Model

This report details the progress of Idaho National Laboratory (INL) in creating a reference plant multiphysics model for the Advanced Burner Test Reactor (ABTR). This model was developed under Task 13 of the U.S. Nuclear Regulatory Commission project “Development and Modeling Support for Advanced Non-Light Water Reactors,” and is an extension of the reference plant model developed in Task 4b, which was improved upon in the following ways. (1) The discrete ordinates method was used in lieu of the super-homogenization (SPH)-corrected diffusion approximation in order to better capture the anisotropic scattering contribution and the neutron leakage change due to thermal expansion. (2) The novel neutronic spatial discretization approach, termed the ring-heterogeneous (RH) approximation, was conceptualized and introduced to capture the differential expansion of the materials in the core. This new technique proved capable of preserving fission rates and maintaining the eigenvalue within 2.5% and 266 pcm with 9 neutron energy groups, respectively. Separating the different materials in the core enables the differential expansion of materials to be explicitly accounted for, eliminating the need for problem-specific cross-section functionalization techniques. (3) The SAM model for the core and system thermal-hydraulics analysis was updated to include 61 channels instead of just four representative ones. This enables users to obtain improved spatial resolution for sodium temperature and density scalar fields. (4) All the mesh files were created via the Multiphysics Object-Oriented Simulation Environment (MOOSE) Reactor module, eliminating all reliance on external tools for mesh creation. (5) Finally, the fuel axial expansion now leverages the HT9 and UPuZr material properties that have been validated against experimental data. The reference plant model was used to perform a full-core unprotected loss of flow (ULOF) transient calculation, including neutronics, thermal and mechanical feedback mechanisms. Future work will be devoted to further enhancements of the model. Potential improvements to the model include the addition of the control rod driveline expansion feedback and the upgrading of the support plate model so as to explicitly include 3D effects. Additionally, a Nuclear Energy Advanced Modeling and Simulation funded parallel effort has completely automated the creation of the ring-heterogeneous (RH) mesh from the fully heterogeneous (FH) geometry, thus maximizing user friendliness for the sodium fast reactor sodium-cooled fast reactor (SFR) workflow and will be incorporated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Influence of Heat Flow from the Planetary Interior on Surface Habitability

an primordial and/or radiogenic heat flow from a rocky planet’s interior help maintain a surface environment amenable to life?On rocky Earth-like planets with shallow oceans (≤5 km depth) that are on the outer edge of the habitable zone, the additional energy provided might make the difference between a world supporting regional habitat space and a world uninhabitable by any known life-form. To test this hypothesis, we introduce heat flux from the planetary interior as a forcing to the NASA/GISSROCKE-3D exoplanet GCM, and explore whether it can enhance the habitability of marginal terrestrial rocky planets (e.g.,paleo Earth in a “snowball” climate state).GCM simulations of modern Earth have not included geothermal heat flux as a climate forcing, as it is more than an order of magnitude less than the 2.9 W/m heating caused by anthropogenic greenhouse gases [Flanner, GRL 2009]. Modern Earth measurements in the oceans [Downes et al., GRL 2019] and results from ocean GCMs [Barnes et al., Ocean Model. 2017]suggest that localized higher heat flux can have measurable thermodynamic and dynamic impacts, such as weakening deep ocean stratification and increased poleward ocean heat transports. A similar result with ROCKE-3D would be important, since our snowball Earth simulations without heat flux produce features like ocean stratification that should have hampered the survival of life. We have created a set of spatially variable geothermal heat fluxes at varying horizontal resolutions based on modern Earth heat flow patterns [Davies, Geochem. Geophys. Geosys. 2013] to test the sensitivity of the ROCKE-3D GCM initially in a modern Earth context. Preliminary results show that even modest heat fluxes have a measurable impact on deep ocean temperatures (upto +0.9°C locally) and that this heat can be transported at depth along ocean circulation paths. Estimates of heat flux impacts on vertical ocean mixing and surface ocean conditions require long simulations at the higher spatial resolutions needed to express heat flow at the scale of most ocean ridges. Ongoing simulations are focused on testing the impacts of horizontal resolution and flux scaling on surface environments, especially for “snowball” climate scenarios.

Heat flow↗

A study of the excitation and radiative decay of the 3s-prime 3D-0 and 3d 3D-0 levels of atomic oxygen

The absolute cross-sections for the excitation of the 989 A, 1027 A, 7990 A, 8446 A, 1.1287 micron and 1.3164 micron multiplets of atomic oxygen by electron impact dissociation of O2 are reported. The radiative branching ratios for these transitions are calculated from these results and compared with the NBS compilation of Wiese et al. (1966) and the recent theoretical calculations of Pradhan and Saraph (1977). The cascade models of O(+) radiative recombination and of electron-impact excitation of the O I(3S) state in the terrestrial airglow are discussed in the light of the laboratory measurements, and the effects of the resonant absorption of components of the lambda 989 A and lambda 1027 A multiplets by the Birge-Hopfield band system of N2 are investigated. This process is shown to depend sensitively on the N2 vibrational temperature and to cause characteristic changes in the OI EUV emission spectrum in auroras and in the sunlit F-region at high exospheric temperatures.

Zipf, E. C.↗

Hydrostatic pressure effect on the Co-based honeycomb magnet BaCo 2 (AsO 4 ) 2

Here the honeycomb antiferromagnet BaCo 2 (AsO 4 ) 2 , in which small in-plane magnetic fields (H1 ≈ 0.26T and H2 ≈ 0.52T at T = 1.8 K < T N ≈ 5.4K) induce two magnetic phase transitions, has attracted attention as a possible candidate material for the realization of Kitaev physics based on the 3d element Co 2+ . Here, we report on the change in the transition temperature T N and the critical fields H1 and H2 of BaCo 2 (AsO 4 ) 2 with hydrostatic pressure up to ~20 kbar, as determined from magnetization and specific heat measurements. Within this pressure range, a marginal increase in the magnetic ordering temperature is observed. At the same time, the critical fields are changed significantly (up to ~25 –35%). Specifically, we find that H1 is increased with hydrostatic pressure, i.e., the antiferromagnetic state is stabilized, whereas H2, which was previously associated with a transition into a proposed Kitaev spin-liquid state, decreases with increasing pressure. We discuss to what extent these results are compatible with suggested models with sizable third-nearest-neighbor exchange. Overall, the results put constraints on the magnetic models that are used to describe the low-temperature magnetic properties of BaCo 2 (AsO 4 ) 2 .

2-dimensional systems↗

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗

Oxidation of 3D-printed SiC in air and steam environments

Here, the high-temperature oxidation of additively manufactured and chemically vapor infiltrated (3D-printed SiC) has been compared to chemical vapor deposited (CVD) SiC. 100-h isothermal exposures were conducted at 1425° and 1300°C at 1 atm under both dry air and steam environments. A SiC reaction tube was utilized to reduce silica volatility. After steam oxidation at 1425° and 1300°C, on the 3D-printed SiC surface, which was intrinsically rougher than the CVD surface, scales were 70%–90% thicker at the convex regions compared to concave/flat regions. In the convex regions, large cracks perpendicular to the oxidizing interface were observed. After dry air oxidation, scale thicknesses were comparable between 3D-printed SiC and CVD SiC, regardless of geometry. Finite element modeling, conducted to elucidate the relationship between SiC geometry and ß- to α-cristobalite transformation stress, determined cristobalite transformation tensile stresses to be on the order of 103 MPa during cool down, assuming a 6 vol% reduction. Compared to flat SiC substrates, tensile transformation stresses were elevated at concave regions and relaxed at convex regions. Combined with specimen mass gain (accounting for the rougher surface) of 3D-printed SiC being 15%–32% higher for 3D-printed SiC after 1300°C and 1425°C steam oxidation, the work presented concludes that the increased oxidation of 3D-printed SiC is primarily caused by tensile hoop stresses driven by oxidation volume expansion. Lastly, the efficacy of the 3D-printing method is demonstrated through the production of tristructural isotropic imbedded 3D-printed SiC fuel forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The climatic effects of nuclear war

The effects of various US-USSR nuclear-exchange scenarios on global climate are investigated by means of computer simulations, summarizing the results of Turco et al. (1983) and follow-up studies using 3D global-circulation models. A nuclear-scenario model is used to determine the amounts of dust, smoke, radioactivity, and pyrotoxins generated by a particular type of nuclear exchange (such as a general 5,000-Mt exchange, a 1,000-Mt limited exchange, a 5,000-Mt hard-target counterforce attack, and a 100-Mt attack on cities only): a particle-microphysics model predicts the evolution of the dust and smoke particles; and a radiative-convective climate model estimates the effects of the dust and smoke clouds on the global radiation budget. The findings are presented in graphs, diagrams, and a table. Thick clouds blocking most sunlight over the Northern Hemisphere midlatitudes for weeks or months and producing ground-temperature reductions of 20-40 C, disruption of global circulation patterns, and rapid spread of clouds to the Southern Hemisphere are among the 'nuclear-winter' effects predicted for the 5,000-Mt baseline case. The catastrophic consequences for plant, animal, and human populations are considered, and the revision of superpower nuclear strategies is urged.

Turco, R. P.↗

The QBO as Potential Amplifier and Conduit to Lower Altitudes of Solar Cycle Influence

In several papers, the solar cycle (SC) effect in the lower atmosphere has been linked to the Quasi-biennial Oscillation (QBO), which is generated primarily by small-scale gravity waves. Salby and Callaghan (2000) analyzed the observed zonal winds of the QBO over more than 40 years and found that it contains a relatively large SC signature at 20 km. Following up on an earlier 2D study with our global-scale Numerical Spectral Model (NSM), we discuss here a 3D study with the QBO under the influence of the SC. For a SC period of 10 years, the amplitude of the relative variations of radiative forcing is taken to vary for simplicity from 0.2% at the surface to 2% at 50 km to 20% at 100 km and above. Covering a limited time span of 40 years, this model produces in the lower stratosphere a relatively large modulation of the QBO, which appears to be related to the SC and is in qualitative agreement with the observations. Some of the energy in the QBO, confined to low latitudes primarily, is redistributed globally by the meridional circulation and planetary waves presumably, so that a measurable SC modulation is generated in the tropospheric temperatures of the polar regions. Further studies are needed, (1) to determine whether the effect is real and prevails in more extensive simulations and whether the results are robust when shorter integration steps are employed, and (2) to explore the mechanism(s) that may ample the apparent SC influence of the UV radiation extending into the lower atmosphere. Quasi-decadal oscillations, generated internally by the QBO interacting with the seasonal cycles, may interfere with or aid the SC effect.

Mayr, Hans G.↗

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗

Footprints of Kitaev spin liquid in the Fano lineshape of Raman-active optical phonons

We develop a theoretical description of the Raman spectroscopy in the spin-phonon-coupled Kitaev system and show that it can provide observable signatures of fractionalized excitations characteristic of the underlying spin-liquid phase. In particular, we obtain the explicit form of the phonon modes and construct the coupling Hamiltonian based on the D 3d symmetry. We then systematically compute the Raman intensity and show that the spin-phonon coupling renormalizes phonon propagators and generates the salient Fano lineshape. We find that the temperature evolution of the Fano lineshape displays two crossovers, and the low-temperature crossover shows pronounced magnetic-field dependence. We thus identify the observable effect of the Majorana fermions and the Z2 gauge fluxes encoded in the Fano lineshape. Our results are consistent with the phonon Raman scattering experiments in the candidate material α-RuCl 3 .

36 MATERIALS SCIENCE↗

Coupled Thermo–Hydro–Mechanical Modeling of Hydro-Shearing Stimulation in an Enhanced Geothermal System in the Raft River Geothermal Field, USA

The Raft River Enhanced Geothermal System (EGS) demonstration project aims to improve the geothermal production by enhancing the reservoir permeability via shear stimulation in well RRG-9. In this paper, we performed a series of 3D thermo– hydro–mechanical (THM) simulations to investigate the influence of hydraulic and thermal effects on the development of this EGS. The model includes synthetic fracture populations based on borehole televiewer images and in-situ stress measurements from well RRG-9. Fracture permeability evolution is determined using an empirical permeability law developed from laboratory experiments. The model was calibrated by comparing the hydraulic response to field observations, including wellhead pressure, injection rate, and well injectivity. Particularly, we analyzed the enhancement of reservoir permeability and the spatial extent of the stimulation zone for the given injection schedule. Our results indicate that the permeability enhancement of fractured geothermal reservoir is caused by the combined effects of injection-induced cooling and fluid pressure increase. The decrease of temperature plays a dominant role in reactivation of natural fractures under the hydro-shearing mechanism, while the higher injection pressure promotes shear failure and enlarges the stimulation zone. Finally, for the specified extensional stress state, the model favors greatest permeability enhancement along the maximum horizontal principal stress, a moderate vertical enhancement, and a smaller gain along the minimum horizontal principal stress.

58 GEOSCIENCES↗

Lattice Boltzmann simulation of the dissolution of slag in alkaline solution using real-shape particles

Highlights: • A dissolution numerical model was proposed to capture the real dissolution kinetics of slag in alkaline solution. • The log forward dissolution rate of Si was described as a function of NBO/T and solution pH. • A threshold solid volume fraction of 0.688 was found for a voxel in 3D, 63.8% larger than that for a pixel in 2D. • The proposed dissolution numerical model provides a reliable alternative to study the dissolution kinetics of slag. A dissolution numerical model was proposed in this study to capture the real dissolution kinetics of slag in alkaline solution. It consists of three modules, i.e. (i) simulation of the initial particle parking structure of slag in alkaline solution using real-shape particles of slag, (ii) simulation of the chemical reactions between slag and solution based on the transition state theory, and (iii) simulation of the physical transport of aqueous ions using the lattice Boltzmann method. This dissolution numerical model was verified using experimental results, showing reasonable accuracy. After verification, the dissolution numerical model was implemented to study the influences of temperature and particle shape using a proper recipe of slag in alkaline solution. This recipe was designed to avoid solid phase precipitation or gel formation via thermodynamic analysis. The simulation results showed faster dissolution kinetics of slag when using higher temperatures and more irregular particle shapes.

36 MATERIALS SCIENCE↗

Wavelength modulation laser-induced fluorescence for plasma characterization

Laser-Induced Fluorescence (LIF) spectroscopy is an essential tool for probing ion and atom velocity distribution functions (VDFs) in complex plasmas. VDFs carry information about the kinetic properties of species that is critical for plasma characterization. Accurate interpretation of these functions is challenging due to factors such as multicomponent distributions, broadening effects, and background emissions. Our research investigates the use of Wavelength Modulation (WM) LIF to enhance the sensitivity of VDF measurements. Unlike standard Amplitude Modulation (AM) methods, WM–LIF measures the derivative of the LIF signal. This approach makes variations in VDF shape more pronounced. VDF measurements with WM–LIF were investigated with both numerical modeling and experimental measurements. The developed model enables the generation of both WM and AM signals, facilitating comparative analysis of fitting outcomes. Experiments were conducted in a weakly collisional argon plasma with magnetized electrons and non-magnetized ions. Measurements of the argon ion VDFs employed a narrow-band tunable diode laser, which scanned the 4p 4 D 7/2 –3d 4 F 9/2 transition centered at 664.553 nm in vacuum. A lock-in amplifier detected the second harmonic WM signal, which was generated by modulating the laser wavelength with an externally controlled piezo-driven mirror of the diode laser. Finally, our findings indicate that the WM–LIF signal is more sensitive to fitting parameters, allowing for better identification of VDF parameters such as the number of distribution components, their temperatures, and velocities. In addition, WM–LIF can serve as an independent method to verify AM measurements and is particularly beneficial in environments with substantial light noise or background emissions, such as those involving thermionic cathodes and reflective surfaces.

47 OTHER INSTRUMENTATION↗

Quantifying Trapped Magnetic Vortex Losses in Niobium Resonators at mK Temperatures

Trapped magnetic vortices in niobium introduce microwave losses that degrade the performance of superconducting resonators. While such losses have been extensively studied above 1 K, we report here their direct quantification in the millikelvin and low-photon regime relevant to quantum devices. Using a high-quality factor 3D niobium cavity cooled through its superconducting transition in controlled magnetic fields, we isolate vortex-induced losses and find the resistive component of the sensitivity to trapped flux S to be approximately 2 nOhm/mG at 10 mK and 6 GHz. The decay rate is initially dominated by two-level system (TLS) losses from the native niobium pentoxide, with vortex-induced degradation of T1 occurring above Btrap~50 mG. In the absence of the oxide, even 10 mG of trapped flux limits performance, Q0 ~ 1010, or T1 ~350 ms, underscoring the need for stringent magnetic shielding. The resistive sensitivity, S, decreases with temperature and remains largely field-independent, whereas the reactive component, S , exhibits a maximum near 0.8 K. These behaviors are well modeled within the Coffey Clem framework in the zero-creep limit, under the assumption that vortex pinning is enhanced by thermally activated processes. Our results suggest that niobium-based transmon qubits can tolerate vortex-induced dissipation at trapped field levels up to several hundred mG, but achieving long coherence times still requires careful magnetic shielding to suppress lower-field losses from other mechanisms

Bafia, Daniel P. [Fermilab]↗

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m 3 . This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion, Large drops, Wind tunnel↗

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m3. This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion↗

Mesosphere Dynamics with Gravity Wave Forcing: Planetary Waves - 2

We present results from a non-linear, 3D, time dependent numerical spectral model (NSM) which extends from the ground up into the thermosphere and incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Our focal point is the mesosphere where wave interactions are playing a dominant role. We discuss planetary waves in the present paper and diurnal and semi-diurnal tides in the companion paper. Without external time dependent energy or momentum sources, planetary waves (PWs) are generated in the model for zonal wavenumbers 1 to 4, which have amplitudes in the mesosphere above 50 km as large as 30 m/s and periods between 2 and 50 days. The waves are generated primarily during solstice conditions, which indicates that the baroclinic instability (associated with the GW driven reversal in the latitudinal temperature gradient) is playing an important role. Results from a numerical experiment show that GWs are also involved directly in generating the PWs. For the zonal wavenumber m = 1, the predominant wave periods in summer are around 4 days and in winter between 6 and 10 days. For m = 2, the periods are in summer and close to 2.5 and 3.5 days respectively For m = 3, 4 the predominant wave periods are in both seasons close to two days. The latter waves have the characteristics of Rossby gravity waves with meridional winds at equatorial latitudes. A common feature of the PWs (m = 1 to 4) generated in summer and winter is that their vertical wavelengths throughout the mesosphere are large which indicates that the waves are not propagating freely but are generated throughout the region. Another common feature is that the PWs propagate preferentially westward in summer and eastward in winter, being launched from the westward and eastward zonal winds that prevail respectively in summer and winter altitudes below 80 km. During spring and fall, for m = 1 and 2 eastward propagating long period PWs are generated that are launched from the smaller eastward zonal winds that prevail in these seasons. PWs generated in the model produce large amplitude modulations of the diurnal tides at altitudes above 80 km and contribute to their seasonal variations.

Mayr, H. G.↗