Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D Icing”

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 91 records · Page 5

Comparison of ice-phase microphysical parameterization schemes using numerical simulations of tropical convection

The performance of several ice parameterizations has been evaluated through a numerical cloud model. Ice effects using different schemes are contrasted with each other and with an ice-free control by incorporating them into the cloud model and by applying them to simulations of tropical squall systems. The latter are simulated in 2D so that a large domain can be used to incorporate a complete anvil. Nonsquall-type convective lines are simulated in 3D owing to their smaller horizontal scale. It is concluded that inclusion of ice microphysics in the cloud model enhanced the agreement of the simulated convection with some features of observed convection, including the proportion of surface rainfall in the anvil region and the intensity and structure of the radar brightband near the melting level in the anvil. In the experiments with bulk microphysics, three ice categories produced much better results than two ice categories, which in turn was better than no ice. For the tropical squall-type and nonsquall-type systems the optimal mix was ice, snow, and graupel.

Mccumber, Michael↗

The Effects of Aircraft Wake Dynamics on Contrail Development

Results of large-eddy simulations of the development of young persistent ice contrails are presented, concentrating on the interactions between the aircraft wake dynamics and the ice cloud evolution over ages front a few seconds to approx. 30 min. The 3D unsteady evolution of the dispersing engine exhausts, trailing vortex pair interaction and breakup, and subsequent Brunt-Vaisala oscillations of the older wake plume are modeled in detail in high-resolution simulations, coupled with it bulk microphysics model for the contrail ice development. The simulations confirm that the early wake dynamics can have a strong influence on the properties of persistent contrails even at late times. The vortex dynamics are the primary determinant of the vertical extent of the contrail (until precipitate ton becomes significant): and this together with the local wind shear largely determines the horizontal extent. The ice density, ice crystal number density, and a conserved exhaust tracer all develop and disperse in different fashions from each other. The total ice crystal number can be significantly reduced due to adiabatic compression resulting from the downward motion of the vortex system, even for ambient conditions that are substantially supersaturated with respect to ice. The fraction of the initial ice crystals surviving, their spatial distribution and the ice mass distribution are all sensitive to the aircraft type, ambient humidity, assumed initial ice crystal number, and ambient turbulence conditions. There is a significant range of conditions for which a smaller transport such as a B737 produces as significant a persistent contrail as a larger transport such as a B747, even though the latter consumes almost five times as much fuel. The difficulties involved in trying to minimize persistent contrail production are discussed.

Lewellen, D. C.↗

Heatshield for Extreme Entry Environment Technology (HEEET) TPS for Ice Giants Probe Missions

This poster provides an overview of the requirements, design, development and testing of the 3D Woven TPS being developed under NASAs Heatshield for Extreme Entry Environment Technology (HEEET) project. Under this current program, NASA is working to develop a Thermal Protection System (TPS) capable of surviving entry into Saturn. A primary goal of the project is to build and test an Engineering Test Unit (ETU) to establish a Technical Readiness Level (TRL) of 6 for this technology by 2018. Poster also discusses use of HEEET TPS for probe missions to the Ice Giants, Uranus and Neptune.

Extreme Entry Environment↗

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

Case study of a low-reflectivity pulsating microburst - Numerical simulation of the Denver, 8 July 1989, storm

An intense, low-reflectivity, and pulsating microburst event, detected on July 8, 1989 by the Low-Level Windshear Alert System north of Denver Stapleton Airport is described using a Terminal Area Simulation System (TASS) model (Proctor, 1987). This 3d nonhydrostatic cloud model includes parametrizations for both liquid- and ice-phase microphysics. Simulation results show that the dissipating cloud induces a type-1 pulsating microburst event, with at least three distinct microburst pulses over a period of 20 min. Virga from the dissipating cloud continued to maintain hazardous windshear for a period much longer than the typical lifetime of a single microburst.

Proctor, Fred H.↗

Simulated Phase-dependent Spectra of Terrestrial Aquaplanets in M Dwarf Systems

Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location ofwater vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.

Wolf, E. T.↗

Case study of a low-reflectivity pulsating microburst - Numerical simulation of the Denver, 8 July 1989, storm

Results from a 3D simulation of a very strong microburst, detected by the Low-Level Windshear Alert System, are presented. The simulation was performed using the Terminal Area Simulation System (TASS) model (Proctor, 1987) which is a 3D nonhydrostatic cloud model that includes parameterizations for both liquid and ice-phase microphysics. Emphasis is placed om the pulsating characteristic and the very-low radar reflectivity of the event. The simulation revealed that the dissipating cloud induces a type-1 pulsating microburst event, with at least three distinct microburst pulses over a period of 20 min.

Proctor, Fred H.↗

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗

3D Carbon Coating Enabled High‐capacity and Stable Micro‐sized Silicon Suboxide‐graphite Blended Anodes for Practical Lithium‐ion Batteries

Abstract Silicon oxide (SiO x ) is a promising anode candidate of lithium‐ion batteries (LIBs) owing to its extremely high specific capacity. However, the low initial Coulombic efficiency (ICE) and rapid capacity degradation of SiO x , triggered by the enormous volume variation upon repeated (de)lithiation, gravely hinder its practical use. Herein, two mass‐produced micro‐sized SiO x @C composites with obviously different morphologies for commercial LIBs are reported. Particularly, the SiO x ‐graphite blended anode (SiO x @3D‐G‐Gr) based on SiO x wrapped by three‐dimensional (3D) carbon layers (SiO x @3D‐G) exhibits a capacity of 519 mAh g −1 , an ICE of 90.0 % and a capacity retention of 83.4 % at 0.2 C over 100 cycles. which is far exceeding its counterpart SiO x @C‐H‐Gr (65.7 %). The obtained impressive properties of SiO x @3D‐G originate from the critical contribution of 3D carbon layers, which serves as the effective stress buffer and protective layer as well as the strong networks for electron/Li + transport. Accordingly, the full‐cell based on SiO x @3D‐G‐Gr anode and commercial LiCoO 2 cathode delivers a capacity of 803 mAh and an excellent capacity retention of 95.6 % (616 mAh, 96.6 % for graphite, respectively) at 1 C over 100 cycles with a stabilized CE of nearly 100 %. The micro‐sized SiO x @3D‐G showing a promising prospect in the commercial‐grade anodes in LIBs.

Electrochemistry↗

Laboratory-Scale Coal-Derived Graphene Process (Final Report)

The Energy & Environmental Research Center (EERC) conducted a laboratory-scale coal-derived graphene (CDG) project focused on developing a technological process for making graphene from four U.S. domestic coal or coal wastes, including lignite from North Dakota, subbituminous coal from Wyoming, bituminous coal from Utah, and anthracite from Pennsylvania. The project was divided into two performance or budget periods (BPs), with BP1 comprising the up-front laboratory experiments to make graphene materials from coal beginning on May 1, 2020, to April 30, 2022. BP2 was conducted from May 1, 2022, to April 30, 2023, and was focused on analyzing the CDG process economic feasibility and the technical gaps for technological scale-up and commercialization. During this project, a few different coal-derived high-value products have been demonstrated, including graphite, graphene oxide (GO), reduced graphene oxide (rGO), and graphene quantum dots (GQDs). A new graphite microstructure was discovered and named “croissant graphite” because of the exterior morphological and textural resemblance to croissant food items sold in commercial groceries stores. The new graphite structure and the associated preparation from coal or coal waste feedstocks has been the subject of a U.S. patent application. The systematic experimental processes involving coal cleaning, upgrading, and conversion to high-value carbon products culminated into a developed upgraded coal-to-products (UCP) technology that is being pursued for potential fast-track commercialization, if funding is available. It is envisioned that commercialization of the UCP technology would increase consumption of U.S. domestic coals or coal wastes to make environmentally sustainable high-value products for the electronics industry, high-energy-storage applications, and clean energy technologies such as electric vehicle (EV) lithium-ion batteries (LIBs), for which graphite has become a critical mineral commodity. Croissant graphite microstructures, when observed by field emission scanning electron microscopy (FESEM), display wavy surface morphology and often grow from a base that is made of graphitized particles with honeycomb-like layers, which are believed to be graphene layers. While more studies are needed to fully ascertain the mechanisms of the croissant graphite microstructure formation, it is postulated that their growth may begin from curling of the graphene sheets into ribbon-like structures, and continuous growth and densification of the ribbon-like structures forms croissant microstructures. Additional studies are ongoing to evaluate the electrochemical performance of croissant graphite for LIB applications and to determine the experimental conditions necessary to tune on/off croissant formation so that it can be either optimized or suppressed depending on performance evaluation results. In addition to the discovery of croissant graphite, the graphitization process from the four coal ranks in general was successful. X-ray diffraction (XRD) analysis showed that the degree of graphitization (DoG) ranged from 12% to 80% in an early sample set, and further optimization on lignite coal produces a DoG of about 92%, which was spectacular to see as lignite is the lowest-rank coal. Thus, it is expected that the graphitization performance for higher-rank coals will be similar or better when optimized as well. The coal-derived graphite was used to make GO and rGO. Analytical characterization, e.g., by methods such as Raman spectroscopy, XRD, Fourier transform infrared (FTIR) spectroscopy and FESEM, showed that the sequence of converting the coal to graphite, exfoliating it to GO, and then chemically reducing the GO to rGO was successful. Although coal naturally contains aromatic compounds and some relatively small-sized condensed aromatic units, it does not contain graphene sheets. In the UCP process, the aromatic domains in the coals, particularly low-rank coals, are concentrated and condensed further into graphene sheets, which are ordered into a 3D stack during graphitization. The synthesized graphite is then unpacked by methods such as exfoliation to various graphene products. GQDs were synthesized from all four coal types, and their optical properties were demonstrated to be tunable by the coal precursor preprocessing treatments. In all four coal types, enhanced optical properties were observed for the produced GDQs with incremental improvements made to the coal precursors. GQDs produced from raw coal samples displayed lower ultraviolet–visible (UV–Vis) spectroscopy absorbance intensity compared to those obtained from cleaned and upgraded coal residues. The photoluminescence (PL) intensities also varied with pretreatment conditions and with the concentration of GQDs in aqueous solutions. GQDs obtained from anthracite show longer emission wavelengths and can be excited by visible light as opposed to GQDs derived from the other coal ranks. UV fluorescence 3D maps and spectra revealed that the emission wavelength at which the GQDs solutions display the highest intensity was slightly redshifted based on the coal precursor pretreatments. In low-rank coal (lignite and subbituminous) samples, two clusters were observed in the maps for GQDs, which may suggest that there are potentially two types of fluorophores in solution or two main size populations. The ability to tune the properties of GQDs based on processing methods can be exploited to make GQDs for various optical display or optoelectronics applications. The results also highlight the importance of removing coal-borne impurities to improve the quality of the coal precursor for preparation of graphene products. Coal and/or coal wastes preprocessing methods were developed and applied to clean and upgrade the coal precursors prior to graphitization and subsequent conversion to graphene products. The preprocessing methods involve high specific-gravity separations, mineral acid cleaning (no hydrofluoric acid), and subsequent upgrading by reducing the coal-borne heteroatom (nitrogen, sulfur, and oxygen) content using proprietary chemical agents. Analytical characterization revealed that the preprocessing steps were successful, with ash reductions that range from 38% to 80% and residual ash content that was below the 5 wt% initial target. Based on proximate and ultimate analysis, the heteroatom reduction reactions produced upgraded coal residues with the oxygen content reduced by 8% to 24%, with additional reductions in the nitrogen and sulfur contents. An initial assessment of the waste streams from the UCP process shows very small to negligible environmental impact due to CO 2 , NO x , and SO x because most process steps are performed under inert atmosphere with argon. Consequently, reactive oxygen environments that tend to create these species are avoided. The inorganic and potentially hazardous species are released into aqueous waste streams that are easy to handle for proper disposal. The liquid waste streams were found to contain low-level concentrations of rare-earth elements (REEs), which could be concentrated and recovered as value-added by-products. Additionally, the volatile and gaseous fractions from carbonization and heat treatment contain useful organic compounds that can also be recovered as potential valuable by-products. Thus, the UCP technology is considered an environmentally sustainable and promising emerging technology for making high-value products from coal and coal wastes, with potential additional value-added by-products. Analysis of potential markets for the coal-derived carbon products shows a strong demand in both niche market sectors and across a wide variety of other industrial sectors. Graphite is currently considered a critical mineral commodity that has a large and growing demand in the LIB industry for EV applications. Based on data from Fortune Business Insights (2022) and Marketwatch (2023) reports, the average global graphite market is projected to reach about 33 billion by 2028, growing at a compound annual growth rate (CAGR) of about 7%, with much of this growth expected to be in the LIB industry. GO and rGO have strong market potentials in various application areas, such as coatings for anticorrosion, anti-icing, and antimicrobial protection, thermal barriers, wear resistance, sensors, additive manufacturing such as 3D inks, and others. GQDs are the emerging key player in the bioimaging, photovoltaics, and light-emitting diodes (LEDs) applications, with the potential to replace traditional semiconductor quantum dots (SQDs), which are based on metallic systems that are more toxic and more expensive. Biomedical applications of GQDs are becoming more attractive because of low to no toxicity and extremely low cost compared to SQDs. The major challenges for scale-up and commercialization of coal-derived carbon products such as graphene vary from the inherent attributes of graphene itself to reluctance to accept graphene in new manufacturing processes because of the uncertainty of the unknown. The 2D nature of graphene materials with a thickness of one atom presents significant challenges to proper handling/processing, and process scale-up becomes difficult because it requires high-end, expensive equipment, even for routine handling and analysis for quality assurance and control. Pristine graphene can also be extremely difficult to work into other matrices, thus hindering downstream processibility, especially at large scale. Currently, the cost of graphene and graphene products is still high and presents an economic risk that tends to slow down investment in scaling up emerging technologies. The lack of a standard for graphene materials for quality assurance and quality control poses a great challenge not only for the markets but also for commercialization efforts. A first-look economic feasibility analysis of the UCP technology provided valuable information that suggests the UCP process would be feasible, especially when it is scaled to a pilot scale and could be more competitive at the full scale. Graphitization was found to be the most energy-consuming and most capital-intensive step in the overall process. In small laboratory- and bench-scale experiments, labor is a significant contributor to the total process costs. Although these energy, capital, and labor constraints contribute to a higher selling price for the product, a preliminary economic model suggests that the process would be feasible at large scale when the process is fully integrated, optimized, and automated.

01 COAL, LIGNITE, AND PEAT↗

Microphysical properties of the Shuttle exhaust cloud

A data base describing the properties of the exhaust cloud produced by the launch of the STS has been developed based on data from a series of ground and aircraft based measurements made during the launches of STS 2, 3, and 4. Aircraft observations were performed during the STS-3 launch with a NOAA WP-3D Orion hurricane research aircraft which contained instrumentation for cloud condensation nucleus and ice nucleus counting, Aitken particle counting, and pH determination. Ground observations were conducted at 50 different sites, as well as in the direct exhaust from the solid rocket booster flame trench at all three launches. The data is analyzed in order to determine any possible adverse impacts of the exhaust products on human health and/or the environment. Analyses of the exhaust cloud measurements indicate that in the case of the ground cloud where plenty of large water drops are present and considerable scavenging and fallout of aerosol takes place, possible adverse impacts of the remaining aerosols (CCN and IN) on natural precipitation processes which may occur in the launch area hours after the launch are remote. However, it is determined that under certain atmospheric conditions there could be short term adverse effects on visibility.

Keller, V. W.↗

A Theory for the Balance between Warm Rain and Ice Crystal Processes of Precipitation in Mixed-Phase Clouds

Abstract Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow (“ice crystal process”), forming “cold” precipitation. Equally, cloud droplets, when large enough, coalesce to form “warm” precipitation by the “warm rain process.” Warm rain could be supercooled and freeze as “warm” graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, “cold ice” (crystals, snow), and “warm rain” (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system’s precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and “wins” the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space.

Meteorology & Atmospheric Sciences↗

Three-dimensional calculation of the mixing of radial jets from slanted slots with a reactive cylindrical crossflow

A numerical model that calculates the 3D chemically reacting flowfield in an experimental low emission combustor is described. The ICEd-ALE finite volume computational methodology is employed in this study. Radial jets issuing from slanted slots interact with a cylindrical axially flowing mainstream to produce a bulk swirl velocity downstream of the slot region. The swirl pattern at a given axial station is composed of a clockwise rotating region near the wall and a counterclockwise rotating region extending from the combustor centerline. The jet radial penetration and downstream swirl velocity axial development are shown to depend principally on the jet-to-mainstream momentum flux ratio.

Winowich, N. S.↗

LDV flowfield measurements on a straight and swept wing with a simulated ice accretion

Simulated glaze ice accretion effects on a flowfield are presently studied for the case of a 3D semispan wing with NACA 0012 airfoil section on a rectangular untwisted planform. The model was tested at both zero-sweep and 30-deg sweep setting with four-beam/two-color LDV flowfield mapping. The comparison of LDV-measure velocity profiles with 3D Navier-Stokes predictions revealed correct trends, but with several differences that are attributable to the turbulence models and grid resolution used in the computations.

Bragg, M. B.↗

Effect of Time-Step Duration and Other Parameters on Ice Accretion

This work addresses the ice prediction cases studies for the 1st AIAA Ice Prediction Workshop generated using a mix of commercial and in-house tools. All predictions are compared with those by LEWICE. Results from both 2D NACA 23012 airfoil cases and 3D NACA 0012 swept wing cases have been included. It is observed that the computational tools have a much better correlation with experimental data for rime icing conditions. For glaze icing conditions, the ice horns are not accurately predicted. The impact of parameters such as surface temperature, heat transfer coefficient, and time-step for flow and droplet field re-computation is explored.

Icing, Rotorcraft↗

Target fabrication fields first NIF shot with new 3D-printed microfluidic target design

Capable of producing ultrafine features, two-photon polymerization (2PP) 3D printing made it possible to field a first-of-its-kind leak-proof liquid-filled diffraction target. These targets are used for shots aimed at using laser pressure instead of low temperatures turn water into ice. The shot is the product of months of design, assembly, and installation work by LSEO’s target fabrication team, in collaboration with 3D printing experts Widi Moestopo (MED) and Xiaoxing Xia (MED) and PLS researcher Neal Bhandarkar.

42 ENGINEERING↗

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of 3D Simulation Capabilities within xRAGE using a Suite of Test Problems

This presentation describes a 3D verification and validation suite for xRAGE. Three verification test problems for hydrodynamics are presented: the Kidder ball problem, the Verney shell problem, and a 5-material compression problem. There is excellent agreement between 2D and 3D XRAGE simulation results for these problems, and between the xRAGE results and the benchmark solutions. In addition, two 3D ICF test problems are presented, based on an OMEGA direct drive capsule experiment, and on a NIF indirect drive capsule experiment. These two ICF test problems are used to evaluate two hydrodynamic methods within xRAGE. In addition to hydrodynamics the simulations also include three temperature plasma physics, thermal conduction and radiation diffusion. When simulating either capsule in 3D the newer unsplit hydrodynamic method in xRAGE produces more vorticity relative to the older default method. For the indirect drive capsule the 3D simulations are in reasonable agreement with the experimental values of ion temperature and neutron production. Unstable 3D vortex rings are present near the DT ice/carbon interface which enhance the confinement of the DT ice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗