Search NASA⌕ Search

SEARCH · Search NASA

Results for “single phase”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

High temperature stabilization of ultrafine grain tungsten alloys through synergistic compositional complexities

Thermally-stabilized fine-grained microstructures in tungsten provide a pathway to harnessing enhanced properties such as a reduced ductile-to-brittle transition temperature and improved strength while mitigating the adverse effects of grain growth and recrystallization. Here, in this study, we employ a material design strategy that relies on grain boundary segregation in the nanocrystalline state driven by alloy thermodynamics balanced with impurity scavenging through in situ formation of kinetically-stabilizing metal carbides. A W-Ti-Cr alloy is designed through lattice Monte Carlo methods and subsequently synthesized in a single-phase nanocrystalline state, which upon annealing, evolves into an ultrafine grained microstructure containing chromium grain boundary segregation collectively with a dispersed TiX (X=C,O) phase through the reaction of titanium with carbon and oxygen impurities. In situ synchrotron X-ray diffraction experiments demonstrate that increased Ti promotes stabilization across a larger temperature range but with diminishing returns above 10 at.% Ti. Long-term stability was confirmed through the retention of the ultrafine grained microstructure for a total of 8 days (192 h) at 1300 °C without grain/carbide growth and/or recrystallization. Our results demonstrate that, through strategic tailoring of composition and microstructure, one can harness the benefits of both thermodynamic and kinetic stabilization mechanisms, opening pathways for future alloy formulations that expand the window of stability.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Probabilistic Model-Based Diagnostic Framework for Nuclear Engineering Systems

A fault diagnostic framework was investigated in this study for applications in thermal–hydraulic systems of nuclear power plants. The proposed framework consists of quantitative model-based diagnosis, statistical change detection and probabilistic reasoning. The use of physics-based diagnostic models provides high detection sensitivity and allows noise and measurement uncertainty to be incorporated robustly. Performance-related parametric models for each component are constructed based on first principles. Numerical model residuals are generated using the concept of analytical redundancy. Statistical change detection methods are employed to detect non-zero residuals in the presence of uncertainty. The diagnosis task is performed using Bayesian inference to detect and localize possible faults. Application to a single-phase heat exchanger for demonstration showed that the proposed probabilistic framework can provide improved results in comparison with traditional approaches while remaining less sensitive to false alarms in the presence of measurement and modeling uncertainty.

Bayesian network↗
DOE OSTI2022

Effect of PVD-coated chromium on the subcooled flow boiling performance of nuclear reactor cladding materials

Here we elucidate the separate effect of a thin Cr coating deposited by physical vapor deposition (PVD) on the subcooled flow boiling performance of zircaloy-4. First, we run flow boiling experiments on prototypical zircaloy-4 surfaces mimicking the scratch pattern and surface roughness of nuclear reactor claddings. Then, we PVD-coat a 0.3 µm thick chromium layer on the same exact surface and repeat the same flow boiling investigations. All experiments are run using deionized water at atmospheric pressure, flowing on a 1 × 3 cm 2 rectangular cross section channel at a rate of 1000 kg/m 2 /s and a subcooling of 10 K. We measure the average temperature of the boiling surface at increasing surface heat fluxes, covering a wide range of heat transfer regimes, from single-phase forced convection to the boiling crisis. We also record high-speed videos of the boiling process, which we postprocess to measure bubble nucleation site density, growth time, departure diameter and frequency. The surface analysis reveals that, while the chromium coating does not seem change the surface roughness and morphology, it improves surface wettability. However, it decreases the critical heat flux. The chromium coating causes an increase of nucleation temperature, bubble departure diameter and growth time, and a reduction of the nucleation site density. The concurrence of these observations indicates that a size reduction of the nucleation sites, conformally covered by the chromium coating, may be the cause of the boiling performance deterioration. We confirm this hypothesis repeating the same analysis on a FeCrAl sample prepared and tested using the same protocol as the zircaloy-4 sample, but with a different initial surface texture.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Reconstructing porous media using generative flow networks

One area of intense scientific interest for the study of sandstones, carbonates, and shale at the pore scale is the use of limited image and petrophysical data to generate multiple realizations of a rock’s pore structure. Such images aid efforts to quantify uncertainty in petrophysical properties, including porosity–permeability transforms. We develop and evaluate a deep learning-based method to synthesize porous media volumes using a so-called generative flow model trained on x-ray microscope images of rock texture and pore structure. These models are optimized on a log-likelihood objective and they synthesize large and realistic three-dimensional images. Further, we demonstrate the rapid generation of sandstone image volumes that display realism as gauged by quantitative comparison of topological features using Minkowski functionals of porosity, specific surface area, and the Euler–Poincaré characteristic (i.e., pore connectivity). We also evaluate the single-phase permeability using Navier–Stokes and lattice Boltzmann methods and show that transport properties of the generated samples match measured trends.

58 GEOSCIENCES↗
DOE OSTI2026

Entropy-defect synergy for dual luminescence mechanism in spinel: Time-resolved anti-counterfeiting and fingerprint visualization

Multimodal luminescent materials, while promising for anti-counterfeiting, often lack dynamic time-dependent responses and controllable spatial distribution, limiting their encryption capabilities in the spatiotemporal dimension. Here, this work presents a coordinated control strategy based on entropy and defect engineering, and uses a backpropagation (BP) neural network for material screening to successfully prepare spinel Mg 0.8 (Fe 0.04 Co 0.04 Ni 0.04 Cu 0.04 Zn 0.04 )Cr 2 O 4 (MgA 5 CO) phosphors with time-dependent dynamic luminescence behavior. This phosphor simultaneously activated the d-d transition luminescence (∼618 nm) derived from Co 2+ /Cr 3+ and the defect luminescence (∼398 nm) related to zinc vacancies (V Zn ) in a single-phase solid solution. The phosphor exhibits a time-dependent color evolution from pink to purple under fixed-wavelength excitation, due to the different excited-state dynamics and decay lifetimes associated with the d-d transition and defect luminescence. Structural characterization and spectral analysis confirmed the existence of V Zn and its significant role in defect luminescence process. The fluorescent and dynamic luminescent properties of entropy-based spinel oxide enable its use in advanced anti-counterfeiting applications like fingerprint recognition and color-changing dedicated anti-counterfeiting mark, showing promise in high-end and time-dynamic anti-counterfeiting fields. This research not only developed a new type of fluorescent dynamic anti-counterfeiting material, but also provided a new idea for constructing advanced optical functional materials with multiple luminescence mechanisms.

Defect project↗
DOE OSTI2020

A perspective on challenges and opportunities in developing high entropy-ultra high temperature ceramics

The emergence of a new class of high entropy ultra-high temperature ceramics (HE-UHTCs) is expected to possess exceptionally superior mechanical, oxidation, and erosion resistance as compared to conventional UHTCs. The improvement in the properties of HE-UHTCs is attributed to the high configurational entropy of the multi-component system stabilized into a single-phase, but with little discussion on entropy contributions. This perspective article presents a viewpoint on the critical scientific question regarding the actual role entropy plays in stabilizing and enhancing the thermo-mechanical properties of a multi-component system. The high throughput first-principle calculations, which can provide the entropy formation ability (EFA) of thermodynamically stable HE-UHTCs, can prove to be valuable toolset to this puzzle. HE-UHTCs show enhancement in the mechanical properties much higher than the predicted rule of the mixture, which cannot traditionally be explained by a mere solid solution formation effect. Additionally, isothermal oxidation studies of HE-UHTCs limited to furnace testing state remarkable oxidation resistance which is attributed to the passivity of complex, layered structure, and composition of the oxide scale. However, whether all the constituents in the HE-UHTC system react simultaneously or preferential oxidation takes place needs to be addressed. To address all these issues, an integrated computational and experimental approach has been postulated towards designing a thermal protection system (TPS) for re-entry applications. Succeeding the approach will provide new opportunities in the compositional space of HE-UHTCs with tailored properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗
DOE OSTI2020

Localized Corrosion Behavior of Non-Equiatomic NiFeCrMnCo Multi-Principal Element Alloys

A series of single phase Ni 38 Fe 20 Cr x Mn 21–0.5x Co 21–0.5x multi-principal element alloys (MPEAs) were fabricated and their localized corrosion behavior was tested in 0.6 M NaCl solution. These MPEAs were found to be susceptible to localized corrosion, exhibiting pitting at room temperature, as indicated by positive hysteresis loops and low repassivation potentials observed in cyclic potentiodynamic polarization experiments. The pitting resistance increased with increasing Cr content. Interestingly, even the MPEAs with low chromium content, Ni 38 Fe 20 Cr 6 Mn 18 Co 18 and Ni 38 Fe 20 Cr 10 Mn 16 Co 16 , exhibited passivity and localized corrosion. Here, the passivity observed in these low chromium MPEAs is likely derived from their relatively high nickel content. Crystallographic pit morphology was observed in most of these MPEAs, suggesting the absence of a salt film on the dissolving pit surface. Such crystallographic pit morphology indicates that the pit growth in most of these MPEAs was under charge-transfer/ohmic control.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stochastic generation of electrolyzer anode catalyst layers

Here, we introduce a stochastic methodology to reproduce the complex pore structure observed in commercial iridium catalyst layers. This method preserves the α pore (pores smaller than 250 nm) and β pore (pores greater than or equal to 250 nm) regions of the catalyst layer. The morphology of the generated materials was validated by comparing the pore size distributions of generated materials against those obtained from commercial materials imaged using x-ray nano computed tomography. We further demonstrate that the pore size distributions of the generated materials are statistically indistinguishable from the imaged catalyst layers, indicating that the stochastic methodology is capable of accurately reproducing catalyst layer morphology. Pore network modelling was conducted on the generated catalyst materials to simulate single-phase permeability, electrical conductivity, and ionic conductivity, and these properties were found to be within experimentally measured ranges for electrolyzer catalyst layers. Additionally, simulations were performed on the generated materials with varying ionomer and iridium catalyst loadings. As the ionomer loading is added, proton conductivity increases exponentially, which demonstrates the importance of optimizing ionomer loading, considering that these effects will be exacerbated in the hydration and temperature conditions of operating electrolyzers. The stochastic material generation method presented in this work is a powerful tool for the development of novel low loading catalyst layers, where the effect of various structural parameters on electrolyzer performance characteristics can be explored.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Flow regime and Reynolds number variation effects on the mixing behavior of parallel flows

The hydraulic single-phase mixing of three parallel rectangular channels is experimentally investigated at various Reynolds numbers (Re) and flow regime combinations. Particle Image Velocimetry results for seven mixing cases are presented and discussed with varying Re combinations ranging from 1,824 to 20,844. While all cases result in the same Re ratio of ~0.69 between the inner and outer flows, two cases represent multi-regime mixing with the inner-outer regime pair of laminar-transitional and transitional-turbulent, while the other 5 cases are all characteristic of turbulent mixing with varying levels of turbulence. The outer channels initially share characteristics with a backward facing step. The center channel is found to initially behave like a slot jet, but then sees a significant increase in velocity decay. This inner flow velocity decay increased dramatically in the laminar-transitional mixing case, whose centerline velocity decay was ~6 times larger than the decay in the turbulent mixing cases. Second order statistics revealed a consistent mixing layer thickness of ~0.1 hydraulic diameters for all the cases but showed more intense shearing in the multi-regime mixing cases. The combined point and thereby the mixing layer length is determined using centerline velocity decay profiles, which show a much more aggressive mixing in multi-regime flows. Multi-regime mixing demonstrated superior characteristics relative to turbulent mixing due to a more dramatic velocity decay in the inner flow and a shorter mixing length. The contributions of this work include communicating the benefits of multi-regime mixing and providing detailed characterization efforts that can serve future efforts for validating computational models. Here this research also lays the groundwork for future studies aimed at achieving high levels of mixing without a severe penalty in pressure drop.

42 ENGINEERING↗
DOE OSTI2024

Experimental investigation of Multi-Mode heat transfer to a Free-Falling dilute particle cloud in a heated vertical tube

The development of dilute particle heat exchangers and reactors for advanced energy systems requires an understanding of the multi-mode heat transfer from a heated wall to falling particles. This study presents experimental results of the overall heat transfer coefficient for a free-falling, dilute flow of particles with solid volume fraction from 0.0005 to 0.006 corresponding to feed rates from 3.7 kg s -1 m -2 to 44 kg s -1 m -2 in a vertical, heated tube containing quiescent air at atmospheric pressure. Tube wall temperatures are varied between 300°C to 900°C while keeping the particle inlet temperature constant. The experimental results show that the overall heat transfer coefficient is a strong function of particle feed rate and surface temperature. Good agreement was obtained with prior studies conducted at comparable temperatures but lower particle feed rates (< 4 kg m -2 s -1 ). The established correlations for particle-to-wall radiation and particle-to-gas convection were used to estimate the wall-to-gas convective contribution from the measured overall heat transfer coefficient. The experimental results indicated a 4 to 6 times improvement in the wall convection in the solid-gas mixture compared to that expected from natural convection in a single-phase gas. Furthermore, the data presented here are applicable to characterize heat transfer in dilute particle heat exchangers, furnaces, and solar receivers.

14 SOLAR ENERGY↗
DOE OSTI2023

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Microstructure, compression properties, and oxidation behavior of Hf-25Ta-5Me alloys (Me is Mo, Nb, W, 0.5Mo + 0.5 W, Cr, or Zr)

Hf–Ta based alloys have recently been investigated as potential candidates for high-temperature structural applications. While most attentions have been given to the properties above ~1200 °C where the alloys are mainly single-phase BCC structures and have excellent oxidation performance due to formation of super-oxides, structural properties at lower temperatures are equally important for applications in which an alloy may experience a range of temperatures. In the present work, microstructure, phase composition, mechanical properties and oxidation behavior of six Hf-25Ta-5Me alloys (Me is Mo (HTM alloy), W (HTW), 0.5Mo + 0.5 W (HTMW), Cr (HTC) or Zr (HTZ), the compositions are in at.%) are reported at temperatures below the eutectoid transformation. The alloys were prepared by arc melting followed by hot isostatic pressing for 3 h at 1400 °C and 207 MPa. All the alloys display coarse grains of partially or fully transformed (by a eutectoid reaction) high-temperature BCC phase. The eutectoid regions consist of fine lamellae of Hf-rich HCP and Ta-rich BCC phases. The ternary alloys containing W or Cr also contained small amounts of a cubic Laves phase. At 25 °C, the HTW alloy was the strongest (yield stress σ y = 1966 MPa) but brittle and HTZ was the weakest (σ y = 1120 MPa) but ductile among the studied alloys. Other alloys showed intermediate behaviors. In general, the room temperature ductility of the alloys increased with decreasing σy. All alloys maintained high strength up to 800 °C, but displayed a noticeable strength decrease at 1000 °C. At 1000 °C, HTMW was the strongest alloy (σ y = 468 MPa) and HTC was the weakest alloy (σ y = 368 MPa). All alloys had excellent deformability at 1000 °C. Oxidation behavior of the alloys was studied at 800 °C and 1000 °C and compared with that of Hf–27Ta binary alloy. Although the ternary additions improved oxidation resistance, overall oxidation performance of the studied alloys at 800 °C and 1000 °C was poor.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effects of HCP/BCC element ratios on the room-temperature tensile properties of Ti-Zr-Hf-Nb-Ta refractory high-entropy alloys

Equiatomic and non-equiatomic Ti-Zr-Hf-Nb-Ta refractory high-entropy alloys (RHEAs) were arc melted, homogenized, cold rolled, and recrystallized to produce single-phase, body-centered cubic (BCC), microstructures with weak texture and equiaxed grains 76–199 μm in size. Here, the non-equiatomic alloys had either a 60:40 or 80:20 atomic ratio of hexagonal close-packed (HCP) elements (Ti + Zr + Hf) to BCC elements (Nb + Ta). Alloy compositions were measured after thermomechanical processing to determine the concentrations of the major (substitutional) and minor (interstitial) elements. We investigated how elastic constants and uniaxial tensile properties were affected by changes in the relative concentrations of the constituent elements at fixed HCP:BCC ratios. Yield strengths ranged from 801 to 922 MPa and ultimate tensile strengths from 815 to 933 MPa. Good agreement is obtained between the experimental yield strengths and those predicted by a strength theory based on edge dislocations indicating that the observed compositional effects are due to their effects on shear modulus and volume misfit. Fracture occurred by dimpled rupture with fracture strains of 19.4%–25.7%, but uniform strains were an order of magnitude lower at 1.1%–3.2%, calling into question the useable ductility (prior to necking) of RHEAs considered to be ductile based on their fracture strain. Contrary to predictions in the literature that HCP elements promote ductility, our present results show that increasing the HCP:BCC ratio decreases both the total strain and the uniform strain. Similar trends were not evident in the yield or ultimate strengths; rather, strengths were affected mainly by shear modulus and volume misfit.

BCC high-entropy alloys↗
DOE OSTI2021

Thermal Processing Conditions for the Synthesis of Near Theoretical Density Li 5 La 3 Ta 2 O 12 Ceramics for Ceramic Dual-Mode Detectors

A family of Li-containing garnet compounds, extensively investigated for Li-ion solid electrolyte applications, show great potential to be developed into transparent ceramic dual-mode scintillators. In particular, Li 5 La 3 Ta 2 O 12 (LLTaO) shows promise as an efficient dual-mode scintillator due to its high Li content, high density, non-hygroscopicity, and high effective-Z. Often the synthesis of LLTaO ceramics involves solid state processing techniques starting with oxide and carbonate powders, calcining, sintering for the removal of residual porosity, and post-processing. This study focuses on two objectives: (1) establishing thermal processing boundary conditions for a flexible alkoxide-based sol-gel method for synthesizing fine-grained powders; and (2), consolidating powders via hot-pressing to form high-density ceramics. In this study, near single-phase ceramics were achieved by hot-pressing sol-gel powders with varying amounts of excess Li. Rietveld refinements were used to analyze X-ray diffraction data and provide a quantitative evaluation of the phases present in the hot pressed ceramics. Some of the diffraction data showed large low angle peak asymmetry, suggesting the presence of two cubic Li 5 La 3 Ta 2 O 12 garnet phases with different lattice parameters. Depending on the synthesis and processing conditions, the refinement revealed ~ 90 wt% concentration of a phase with a ≈ 12.8 Å and ~ 10 wt% concentration of a phase with a ≈ 12.9 Å, indicating a remaining excess of Li in the hot pressed ceramics in the large lattice parameter phase. These synthesis experiments allowed for the determination of a baseline for Li loss during hot-pressing. X-ray excited radioluminescence and photoluminescence excitation and emission spectra were obtained from the hot pressed ceramics. In conclusion, the measured emission spectra featured broad tantalate emission and line emissions likely from a rare-earth impurity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Synthesis and magnetic properties of LaFe 13- x - y M x Si y N 3 nitrides

Here, single-phase LaFe 13-x-y M x Si y N 3 compounds (M = Mn) were synthesized by subjecting intermetallic precursors to the ammonia flow at 623 K for 4 h. The nitrides derived from Si-poor (y ≤ 3) cubic intermetallics retain the NaZn 13 -type crystal structure and show considerable (~3 %) lattice expansion compared to the nitrogen-free materials. The temperatures of the magnetic ordering transitions, all of which are second order, increase by ~ 200 K as a result of nitrogenation, and all nitrides exhibit rather weak and broad magnetocaloric effects. The nitrides are stable below 750 K, but they begin to lose nitrogen above this temperature, fully decomposing above 1050 K. The nitrogenation of the tetragonal LaFe 9 Si 4 results in the decomposition of the intermetallic phase and formation of the Fe 3 N compound together with unidentified amorphous phase(s).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tetragonal structure and uniaxial magnetic anisotropy in the arc-melted (Ce,Zr) 2 (Fe,M) 17 (M = Mo, W, Co) alloys

The tetragonal ThMn 12 -type structure is stabilized in (Ce 1-x Zr x ) 2 T 16 M (x = 0.2 – 0.3; T = Fe or Fe/Co, M = Mo or W) arc-melted alloys. Approximately 5 at% of Mo or W admixture is sufficient to transform the hexagonal Th 2 Ni 17 -type structure of (Ce 1-x Zr x ) 2 Fe 17 into the tetragonal ThMn 12 -type structure of nearly single-phase (Ce 1-x Zr x ) 2 Fe 16 M and/or (Ce 1-x Zr x ) 2 Fe 15 CoM bulk alloys. X-ray Rietveld refinements reveal that Zr and Mo (W) substitute different sites in the tetragonal crystal structure. Zirconium preferentially replaces Ce, whereas Mo (W) substitutes Fe. At room temperature, the tetragonal phases exhibit strong ferromagnetism and a uniaxial magneto-crystalline anisotropy with anisotropy fields of 15 – 17 kOe. These materials possess room-temperature saturation magnetizations of 95 – 105 emu/g and Curie temperature of 420 – 505 K. At some concentrations, concurrent Zr and Mo (W) site occupancy facilitates 1:12 structure formation in bulk alloys with minimal presence of non-magnetic element in the Fe sublattice, thus securing the highest known magnetic moment per Fe atom (∼ 1.37 µ B ) in this type of materials. The intrinsic magnetic characteristics, as well as the absence of critical rare earths (RE), make these compounds interesting for development as low-cost permanent magnets. In conclusion, small Co additions improve the Curie temperature, especially if combined with W.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗