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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Advancing Geophysical Techniques to Image a Stratigraphic Hydrothermal Resource

Sedimentary-hosted geothermal energy systems are permeable structural, structural-stratigraphic, and/or stratigraphic horizons with sufficient temperature for direct use and/or electricity generation. Sedimentary-hosted (i.e., stratigraphic) geothermal reservoirs may be present in multiple locations across the central and eastern Great Basin of the USA, thereby constituting a potentially large base of untapped, economically accessible energy resources. Sandia National Laboratories has partnered with a multi disciplinary group of collaborators to evaluate a stratigraphic system in Steptoe Valley, Nevada using both established and novel geophysical imaging techniques. The goal of this study is to inform an optimized strategy for subsequent exploration and development of this resource and analogous ones. Building from prior Nevada Play Fairway Analysis (PFA), this team is primarily 1) collecting additional geophysical data, 2) employing novel joint geophysical inversion/modeling techniques to update existing 3D geologic models, and 3) integrating the geophysical results to produce a working, geologically constrained thermo-hydrological reservoir model. Prior PFA work highlights Steptoe Valley as a favorable resource basin that likely has both sedimentary and hydrothermal characteristics. However, there remains significant uncertainty on the nature and architecture of the resource(s) at depth, which increases the risk in exploratory drilling. Newly acquired gravity, magnetic, magnetotelluric, and controlled-source electromagnetic data products, in conjunction with new and preexisting geoscientific measurements and observations, are being integrated and evaluated for efficacy in understanding stratigraphic geothermal resources and mitigating exploration risk. Furthermore, the influence of hydrothermal activity on sedimentary-hosted reservoirs in favorable structural settings, and whether fault-controlled systems may locally enhance temperature and permeability in some deep stratigraphic reservoirs, will also be evaluated.

Geothermal, Sedimentary Heat, Geophysics, Seismic,↗

DOE OSTI2026

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 for pin lattice models 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. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. 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 a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. 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 OSTI2025

Interface‐Induced Stability of Nontrivial Topological Spin Textures: Unveiling Room‐Temperature Hopfions and Skyrmions

Topological spin configurations, such as soliton-like spin texture and Dirac electron assemblies, have recently emerged in fundamental science and technology. Achieving stable topological spin textures at room temperature is crucial for their use as long-range information carriers. However, their creation and manipulation are hindered by multi-step field training and competing interactions. Thus, a spontaneous ground state for multidimensional topological spin textures is desirable, with skyrmions forming swirling, hedgehog-like spin structures in two dimensions and hopfions as their twisted 3D counterparts. Here, the first observation of robust and reproducible topological spin textures of hopfions and skyrmions observed at room temperature and in zero magnetic field is reported, which are stabilized by geometric confinement and protected by interfacial magnetism in a ferromagnet/topological insulator/ferromagnet trilayer heterostructure. These skyrmion-hopfion configurations are directly observed at room temperature with Lorenz transmission electron microscopy. Using micromagnetic modeling, the experimental observations of hopfion-skyrmion assemblies are reproduced. This model reveals a complete picture of how spontaneously organized skyrmion lattices encircled by hopfion rings are controlled by surface electrons, uniaxial anisotropy, and Dzyaloshinskii-Moriya interaction. This study provides evidence that topological chiral spin textures can facilitate the development of magnetic topological carriers, paving the way for ultralow-power and high-density information processing.

hopfions↗

DOE OSTI2025

Harnessing on-machine metrology data for prints with a surrogate model for laser powder directed energy deposition

In this study, we leverage the massive amount of multi-modal on-machine metrology data generated from Laser Powder Directed Energy Deposition (LP-DED) to construct a comprehensive surrogate model of the 3D printing process. By employing Dynamic Mode Decomposition with Control (DMDc), a data-driven technique, we capture the complex physics inherent in this extensive dataset. This physics-based surrogate model emphasizes thermodynamically significant quantities, enabling us to accurately predict key process outcomes. The model ingests 21 process parameters, including laser power, scan rate, and position, while providing outputs such as melt pool temperature, melt pool size, and other essential observables. Furthermore, it incorporates uncertainty quantification to provide bounds on these predictions, enhancing reliability and confidence in the results. We then deploy the surrogate model on a new, unseen part and monitor the printing process as validation of the method. Our experimental results demonstrate that the predictions align with actual measurements with high accuracy, confirming the effectiveness of our approach. Furthermore, this methodology not only facilitates real-time predictions but also operates at process-relevant speeds, establishing a basis for implementing feedback control in LP-DED.

Digital twins↗

DOE OSTI2023

Chemical structure and complex growth modes of magnesium silicate hydrate: Nanoparticle orientation, aggregation, and fusion

The extent of utilization of magnesium silicate hydrate (MSH) in construction is limited partly because of insufficient data ascertaining the kinetics of its precipitation. In this report MSH grown homogeneously or heterogeneously in the range of magnesium to silicon solution molar concentration ratios, [Mg]/[Si] = 0.5-1.5, and temperature = 25-80 °C was analyzed for its chemical structure and morphology. Infrared spectroscopy and geochemical modeling show that increasing the [Mg]/[Si] ratio results in less silicate polymerization in MSH. Atomic force microscopy reveals the presence of nanoparticles that aggregate and fuse to form films. Oriented attachment of nanoparticles and features indicative of enhanced crystallinity were observed at higher temperatures and longer reaction times. These data provide direct evidence for the persistence of hierarchical structures (i.e., nanoparticles forming 3D microparticles and 2D film layers) at the nanoscale to mesoscale. These findings offer insights into the precise chemical synthesis of MSH and its widespread use as a binder for construction purposes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

RF sheath induced sputtering on Proto-MPEX. I. Sheath equivalent dielectric layer for modeling the RF sheath

The pulsed linear plasma device Prototype Material Plasma Exposure eXperiment (Proto-MPEX) uses a radio frequency (RF) helicon antenna with an aluminum nitride ceramic window for plasma production. The RF sheath created under the helicon antenna is sufficient to cause ion impact energies to be greater than the sputtering threshold of the AlN helicon window material and for impurities to be created. In this study, we investigate the RF sheath on the inner diameter of the helicon window and its impact on the impurity production rates in Proto-MPEX. Three models—a 3D COMSOL finite element RF model of the Proto-MPEX helicon region, a rectified DC sheath potential model, and the 3D Global Impurity TRansport code—are coupled together to study impurity production and transportation. This novel method of impurity generation and transport modeling spans length scales ranging from the sheath (millimeters or less) up to the full device (meters) and can be applied to other radio frequency sources and antennas in a wide range of plasma physics studies, including basic plasmas, low-temperature processing plasmas, plasma thrusters, and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Factors regulating ozone over the United States and its export to the global atmosphere

Attention is given to the factors regulating summertime O3 over the U.S. and its export to the global atmosphere, which are examined via a 3-mo simulation involving a continental-scale 3D photochemical model. It is found that reducing NO(x) emissions by 50 percent from 1985 levels would decrease rural O3 concentrations over the eastern U.S. by about 15 percent under almost all meteorological conditions, while reducing anthropogenic hydrocarbon emissions by 50 percent would have less than a 4 percent effect except in the largest urban plumes. The correlation between O3 concentrations and temperature observed at eastern U.S. sites is attributed in part to the association of high temperatures with regional stagnation, and in part to an actual dependence of O3 production on temperature driven primarily by conversion of NO(x) to PAN.

Jacob, Daniel J.↗

DOE OSTI2023

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Multiband effects on the upper critical field angular dependence of 122-family iron pnictide superconductors

Detailed measurements of the in-plane resistivity were performed in a high-quality Ba(Fe 1-x Co x ) 2 As 2 (x=0.065) single crystal, in magnetic fields up to 9 T and with different orientations θ relative to the crystal c axis. A significant ρ(T) H,θ rounding is observed just above the superconducting critical temperature T c due to Cooper pairs created by superconducting fluctuations. These data are analyzed in terms of a generalization of the Aslamazov-Larkin approach, that extends its applicability to high reduced-temperatures and magnetic fields. This method allows us to carry out a criterion-independent determination of the angular dependence of the upper critical field, H c2 (θ). In spite of the relatively small anisotropy of this compound, it is found that H c2 (θ) presents a significant deviation from the single-band 3D anisotropic Ginzburg-Landau (3D-aGL) approach, particularly for large θ (typically above ~60o). These results are interpreted in terms of the multiband nature of these materials, in contrast with other proposals for similar H c2 (θ) anomalies. Our results are also consistent with an effective anisotropy factor almost temperature independent near T c , a result that differs from the ones obtained by using a single-band model.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Estimating the Electric Fields Driving Lightning Dart Leader Development With BIMAP‐3D Observations

Abstract In this paper, a numerical dart leader model has been implemented to understand the leader's development and the corresponding electric field changes observed by the 3D broadband mapping and polarization (BIMAP‐3D) system. The model assumes the extending leader channel is equipotential and has a linear charge distribution induced by an ambient electric field. The charge distribution induced by the ambient field can be used to model the electric field change at the ground. We then find the ambient electric field which best fits the field change measurements at the two BIMAP stations. The estimated ambient electric field decreases in the direction of dart leader propagation. Our observations and modeling results are consistent with our earlier hypothesis that dart leader speed is proportional to the electric field at the leader tip. The model also supports our earlier analysis that leader speed variations near branch junctions were due to previous charge deposits near the junctions. The modeled tip electric field is generally lower than the breakdown field unless the pre‐dart‐leader channel has a significant temperature of 3,000 K. This is consistent with the fact that dart leaders typically do not form new branches into the virgin air. Furthermore, the tip field is generally close to the negative streamer stability field at ambient temperatures, explaining the nature of the narrow and well‐defined channel structure. In addition to the charge distribution and the ambient and tip electric field, the development of the channel potential and current distribution are also presented.

3D lightning mapping↗

DOE OSTI2024

Three-dimensional reconstruction of implosion stagnation in laser direct drive on OMEGA

Multidimensional effects on hot-spot formation must be considered to better understand the current limits on the performance of direct-drive inertial confinement fusion experiments on OMEGA with cryogenically layered solid deuterium–tritium targets. A comprehensive reconstruction effort has been established at the Laboratory for Laser Energetics to infer hot-spot and shell conditions at stagnation from a large collection of x-ray, neutron, and particle detectors along multiple lines of sight. Several time-gated and time-integrated x-ray imagers are being used to record the shape of the hot-spot plasma. A 3D hot-spot x-ray emission tomography technique has been developed to infer low-mode drive asymmetries from the hot-spot shape. A suite of neutron diagnostics is used to provide measurements of hot-spot flow velocity, ion temperature, and areal density. Here, the information obtained from the x-ray and neutron detectors will be combined into a coherent model of the shape of the hot spot and shell assembly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water

The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Development of a Particle Flow Control Mechanism

The next generation of Concentrated Solar Power (CSP) systems are utilizing solid particles as the heat transfer medium. The state-of-the-art particle flow control mechanism currently utilized is the slide-gate mechanism which functions by linearly actuating a gate across a tear-drop shaped opening. The sensitivity, which is a measure of the change in mass flow rate per unit movement of the slide-gate, is 0.354 g/s-mm[1]. The slide-gate’s main shortcomings include particle lodging and particle leakage between the plates, and difficulty actuating under the pressure from the particles. The goal of a new device is to minimize or eliminate these flaws. The new proposed device functions similar to a chuck mechanism in a drill where a rotation of an outer ring produces a change in outlet area through moving “jaws” diagonally. This mechanism aims to solve the issues of the slide-gate by having the jaws on the outside of the particle flow volume to reduce the probability of particle lodging, and the jaw tips are designed to overlap as they actuate to reduce particle leakage. Additionally, the vertical orientation of the device as well as the inclined movement of the jaws allows for easier actuation of the device under pressure. Experiments involving 3D printed prototypes have been carried out, producing a max mass flowrate of about 50 g/s and a sensitivity curve with the equation 𝑦 (𝑔/𝑠) = 0.0851𝑥^2(𝑔/𝑠−𝑚𝑚2) +3.6993𝑥(𝑔/𝑠−𝑚𝑚). These experiments validated the device’s capability to address the shortcomings of the slide-gate. Currently, an aluminum model of the device is being tested at room temperature, and a stainless-steel model is being developed to be tested at 700oC, which is the expected operating temperature of the CSP system. Furthermore, simulations are on-going to predict the mass flowrate of different device configurations at room and high temperatures.

14 SOLAR ENERGY↗

NASA NTRSDate not supplied

Secondary Volcanically-Induced Lunar Atmosphere and Lunar Volatiles: 3-D Modeling and Analysis

The fact that the Moon could have a transient secondary atmosphere due to volcanic outgassing has been known for some time, though typically such an atmosphere was believed to be extremely thin (~10-8 bar) [1]. But recent research by Needham and Kring (NK) [2] suggests that during the peak of volcanic activity ~3.5 Ga such a volcanically-outgassed atmosphere could reach ~10-2 bar of surface pressure. In similar research Wilson et al. [3] proposed a more conservative estimate, arguing that the thickness of such an atmosphere would depend on the intervals between major eruptions and may not exceed microbar densities. In either case a collisional atmosphere could be present, which would control transport of outgassed volatiles (such as H2O) and their deposition in polar regions, where they could be preserved until modern day frozen in permanently shadowed regions (PSR) or buried beneath the regolith. Here we study such a hypothetical atmosphere to investigate its stability, meteorological properties and the effect on transport of volatiles. We use the ROCKE-3D planetary 3-D General Circulation Model (GCM)[4]. The insolation and orbital parameters were set to conditions 3.5 Ga. The atmospheric composition, based on the list of outgassed species presented by NK in combination with our estimates for atmospheric escape, condensation and the results from our 1-D chemistry model, was chosen to be either CO-dominated or CO2-dominated (depending on atmospheric temperature). In this study we restricted ourselves to relatively "thick" lunar atmospheres of 1-10 mb, though we believe that our results will scale to thinner atmospheres as well. We present the results for ground and atmospheric temperature for modeled atmospheres over a wide parameter space. In particular we consider different atmospheric compositions (CO or CO2 dominated), a set of atmospheric pressures from 1 mb to 10 mb and a set of obliquities from 0o to 40o. We also present an experiment of a single major eruption [5] and show that in just 3 years ~80% of the outgassed water is deposited in polar regions. This demonstrates the efficiency of such an atmosphere in delivering volatiles. We argue that a secondary lunar atmosphere could play a significant role in forming volatile deposits currently observed in the polar regions of the Moon.

Lunar atmospheres↗

NASA NTRS2019

Manufacturing of Thermoset Polyimide Composites by Laser Sintering

Selective Laser Sintering (SLS) is an additive manufacturing technique that builds 3D models layer by layer using a laser to selectively melt cross sections in powdered polymeric materials, following sequential slices of the computer-aided design (CAD) model. SLS generally uses thermoplastic polymeric powders such as polyamides. The resultant 3D-printed objects are often weaker in their strength compared to traditionally processed materials, due to their higher porosity. This paper described the process development of using melt-processable imide oligomers terminated with reactive 4-phenylethynylphthalic anhydride (4-PEPA) to conduct laser sintering (LS). The first successful 3D-printing of high temperature RTM370 thermoset polyimide carbon fiber composites were further post-cured to promote additional crosslinking for achieving higher temperature (T g = 370°C) capability. Another novel imide oligomer, RTM385-SLS, formulated with a complex melt viscosity [ƞ*] of ~104-105 poise is also suitable for LS. RTM385-SLS resin powder was mixed with 20-25% of hexagonal boron nitride (h-BN) and subjected to LS to print out “Green” specimens which could be further post-cured to afford a thermally conductive but electrically insulating composites with high T g of 385 °C. The cured composite specimens were then subjected to mechanical testing, thermal conductivity and porosity measurements as well as SEM characterization.

additive Manufacturing↗

NASA NTRSDate not supplied

Study of Martian Polar Caps with GISS ROCKE-3D GCM

Martian polar caps consist of both H 2 O and CO 2 ice. While H 2 O ice is mainly passive on modern Mars, it may have not been the case in recent Martian history, when its obliquity was higher, or when it was changing rapidly. The distribution of ice species in the snowpack affects its physical and thermodynamic properties. In the upper layers, it determines its albedo and thermal emissivity. Thus understanding the mutual effect between these ices and their interaction with the atmosphere is crucial for understanding the evolution of Martian polar regions. In this study, we employ a newly-developed Exotic Ices snow model coupled to the NASA Goddard Institute for Space Studies (GISS) ROCKE-3D planetary General Circulation Model (GCM) [1] to study the behavior of Martian polar caps. ROCKE-3D is a planetary GCM developed at NASA GISS as an extension of its Earth climate model, modelE [2]. It has been extensively used to simulate climate of various planets, including Mars (e.g. [3,4]). The Exotic Ices snow model was specially developed for planetary applications which involve more than one condensable in the atmosphere, in which case snow can contain multiple species of ice (CO2 and H 2 O in the Mars case). For each species of ice, the model uses their proper physical properties and phase diagram, but otherwise it treats all species of ice on an equal footing. The combined effects on albedo, thermal inertia and mutual insulation are treated accordingly. The snowpack interacts with the atmospheric dust cycle, and can accumulate a prognostic amount of dust, though the effect of dust on snow properties is not currently treated explicitly, and is prescribed. In this study, we first validate our model against the modern Martial climate, for which we use mission results from Mars Climate Sounder (atmospheric temperature and dust optical depth), SPICAM on Mars Express (atmospheric water), and Viking 2 (surface pressure). We investigate the effect of snow radiative properties on CO 2 and water cycles and the ability of our model to accurately reproduce those with minimal model tuning. We then perform simulations for several obliquities from a recent Martian past, and investigate the behavior of the Martian polar caps in such conditions.

Mars↗

DOE OSTI2019

An anisotropic damage model based on dislocation-mediated nucleation of cracks under high-rate compression

In this report, we developed a thermodynamically-consistent, rate-dependent micromechanics model for brittle damage nucleated by dislocation plasticity applicable for large deformations. Dislocation substructure evolution was used to inform a nucleation criterion for a microcrack. Under global compression, the sliding of a microcrack induces formation of wing cracks. Effective stress drives dynamic growth of these cracks under a 3D stress state, resulting in an anisotropic material stiffness. The model was also advanced to predict grain size dependence of a polycrystalline solid. Internal variables were constrained based on the laws of thermodynamics. Material constants were calibrated for polycrystalline beryllium to demonstrate the applicability of the model to simulate dynamic failure under compression. We demonstrate the versatility of the model to capture brittle to ductile transition governed by temperature and strain rate. The predictive capability of the model to simulate failure stress and failure strain is compared with dynamic and quasistatic data on beryllium.

36 MATERIALS SCIENCE↗

DOE OSTI2025

3D Multiphysics Model for Large Scale Planar Cell: Inductance Investigation in Impedance Analysis

Developed a comprehensive 3D multiphysics model to analyze impedance behavior in large-scale planar Solid Oxide Fuel Cells (SOFCs). The model investigates how impedance varies across operating conditions and cell components, with a particular focus on inductive loops or negative capacitance at the air electrode. It was found that the Inductive loop disappears in the low-frequency zone when the activation overpotential or faradic current is not temperature-dependent. These findings offer new insights into SOFC impedance behavior.

impedence analysis↗