Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Inter-species molecular attraction effect in the development of a two-species mixing layer

Simulations of fluid mixing occurring at high-pressures mandate the use of real-fluid equations of state to capture the effect of boh intermolecular repulsive forces and attractive forces which are neglected in perfect gases. Dedicated mixing rules account for these forces by providing a generic mapping of the mixture state-space onto a pure fluid domain. However, the specific inter-species molecular forces depend on the particular characteristics of the molecular pairs under consideration, and typically require the knowledge of binary interaction coefficients k_{αβ}. The impact of k_{αβ} on a flow field is as of now unclear; further, k_{αβ} is unknown for many molecular pairs and when unknown, set to be null. This study addresses the impact of k_{αβ} on the temporal evolution of a real-fluid mixing layer and uses Direct Numerical Simulation (DNS) to explore this impact. The results show differences between augmented attraction (k_{αβ}<0) and diminished attraction (k_{αβ}>0), affecting both fluid dynamics and diffusion processes. Specifically, in a binary mixing layer, it is observed that augmented attraction leads to delayed transition and mixing layer growth, manifested in the momentum layer thickness. The fractal dimension, which is a manifestation of interface corrugation, is calculated. These results provide an hitherto undocumented mechanism affecting the development of mixing layers.

Banuti, Daniel↗
DOE OSTI2021

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Learning in Modal Space: Solving Time-Dependent Stochastic PDEs Using Physics-Informed Neural Networks

One of the open problems in scientific computing is the long-time integration of nonlinear stochastic partial differential equations (SPDEs), especially with arbitrary initial data. We address this problem by taking advantage of recent advances in scientific machine learning and the spectral dynamically orthogonal (DO) and borthogonal (BO) methods for representing stochastic processes. The recently introduced DO/BO methods reduce the SPDE to solving a system of deterministic PDEs and a system of stochastic ordinary differential equations. Specifically, we propose two new physics-informed neural networks (PINNs) for solving time-dependent SPDEs, namely the neural network (NN)-DO/BO methods. The proposed methods incorporate the DO/BO constraints into the loss function (along with the modal decomposition of the SPDE) with an implicit form instead of generating explicit expressions for the temporal derivatives of the DO/BO modes. Hence, the NN-DO/BO methods can overcome some of the drawbacks of the original DO/BO methods. For example, we do not need the assumption that the covariance matrix of the random coefficients is invertible as in the original DO method, and we can remove the assumption of no eigenvalue crossing as in the original BO method. Moreover, the NN-DO/BO methods can be used to solve time-dependent stochastic inverse problems with the same formulation and same computational complexity as for forward problems. Furthermore, we demonstrate the capability of the proposed methods via several numerical examples, namely: (1) A linear stochastic advection equation with deterministic initial condition: we obtain good results with the proposed methods, while the original DO/BO methods cannot be applied directly in this case. (2) Long-time integration of the stochastic Burgers' equation: we show the good performance of NN-DO/BO methods, especially the effectiveness of the NN-BO approach for such problems with many eigenvalue crossings during the whole time evolution, while the original BO method fails. (3) Nonlinear reaction diffusion equation: we consider both the forward problem and the inverse problems, including very noisy initial point values, to investigate the flexibility of the NN-DO/BO methods in handling inverse and mixed type problems. Taken together, these simulation results demonstrate that the NN-DO/BO methods can be employed to effectively quantify uncertainty propagation in a wide range of physical problems, but future work should address the efficiency issue of PINNs for forward problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Generation of localized reactor point kinetics parameters using coupled neutronic and thermal fluid models for pebble-bed reactor transient analysis

The systems analysis of anticipated operating occurrences and design basis accidents for pebble-bed reactor systems requires knowledge of neutron point kinetics equations (PKE) parameters. Typically, the generation of PKE parameters is performed in a global manner using standalone neutronics calculations, without the inclusion of thermal fluid distributions. We utilize Griffin and Pronghorn for generating global and local PKE parameters which includes the use of thermal fluid distributions to account for localized effects. This work establishes a methodology for calculating PKE parameters for a pebble bed reactor with a coupled neutronics/thermal fluids analysis. PKE parameters generated on a global and local basis are compared against a diffusion solve for a typical load-following transient. Locally-generated neutron kinetic parameters are able to reduce the maximum error in the transient power level from 5% to below 1.5%; along with this, bulk temperature errors were reduced from 12 K to 4 K.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

The advancing wave front on a sloping channel covered by a rod canopy following an instantaneous dam break

The drag coefficient Cd for a rigid and uniformly distributed rod canopy covering a sloping channel following the instantaneous collapse of a dam was examined using flume experiments. The measurements included space x and time t high resolution images of the water surface h(x, t) for multiple channel bed slopes So and water depths behind the dam Ho along with drag estimates provided by sequential load cells. Using these data, an analysis of the Saint-Venant equation (SVE) for the front speed was conducted using the diffusive wave approximation. An inferred Cd=0.4 from the h(x, t) data near the advancing front region, also confirmed by load cell measurements, is much reduced relative to its independently measured steady-uniform flow case. This finding suggests that drag reduction mechanisms associated with transients and flow disturbances are more likely to play a dominant role when compared to conventional sheltering or blocking effects on Cd examined in uniform flow. The increased air volume entrained into the advancing wave front region as determined from an inflow–outflow volume balance partly explains the Cd reduction from unity.

Mechanics↗
NASA NTRS1972

Effects of vaporization and condensation on Apollo 11 glass spherules - Implications for cooling rates.

Fourteen of 40 glass spherules present in a section of the Apollo 11 microbreccia 10019,22 were found to exhibit steep concentration gradients at their rims; the oxides of the relatively volatile elements Na, K, and P increase by, respectively, factors of up to 67, 16, and greater than 54 at the spherule rims in comparison to the homogeneous centers. These gradients usually extend over 25-30 microns. Furthermore, the natural surfaces of 11 unpolished glass spherules separated from the fines 10084,97, in 60-80% of all cases, show higher concentrations in Na2O, K2O, and P2O5 than the averages for the central portions of 45 independently analyzed glass spherules. It is suggested that the concentration gradients observed are diffusion gradients and are the result of the impact event which produced the glasses; in the impact melting process, splash drops of melt formed that lost volatiles by vaporization. A procedure is described to determine the cooling rate of a spherule from measured concentration gradients, provided the diffusion coefficients are known for the particular glass compositions in question.

Kurat, G.↗
NASA NTRS1989

An embedded grid formulation applied to delta wings

An embedded grid algorithm for the Euler and/or Navier-Stokes equations is developed and applied to delta wings at high angles of attack in low speed flow. The Navier-Stokes code is an implicit, finite volume algorithm, using flux difference splitting for the convective and pressure terms and central differencing for the viscous and heat transfer terms. Calculations are compared with detailed experimental results over an angle of attack range up to and beyond the maximum lift coefficient, corresponding to vortex breakdown at the trailing edge, for a delta wing nominally of unit aspect ratio. The results indicate that the overall flowfield, including surface pressures, surface streamlines, and vortex trajectories, can be simulated accurately with the global grid version of the present algorithm. However, comparison of computed velocities and vorticity with experimentally measured off-body values at an angle of attack of 20.5 deg indicates the core region is substantially more diffuse in the computations than that measured with either a five-hole probe or a laser velocimeter. Embedded grids, used to improve the numerical discretization in the core region, are formulated within the framework of the implicit, upwind-biased multi-grid algorithm. Structured levels of local nested refinements are made. Three-dimensional results for both Euler and Navier-Stokes calculations are shown, with up to 3 levels of embedded refinement. The embedding procedure was effective in eliminating a crossflow secondary separation produced in the Euler solutions on coarse grids.

Thomas, James L.↗
NASA NTRS1996

Radiative Effects of Aerosols

During the Atlantic Stratocumulus Transition Experiment (ASTEX) in June 1992, two descents in cloud-free regions allowed comparison of the change in aerosol optical depth as determined by an onboard total-direct-diffuse radiometer (TDDR) to the change calculated from measured size-resolved aerosol microphysics and chemistry. Both profiles included a pollution haze from Europe but the second also included the effect of a Saharan dust layer above the haze. The separate contributions of supermicrometer (coarse) and submicrometer (fine) aerosol were determined and thermal analysis of the pollution haze indicated that the fine aerosol was composed primarily of a sulfate/water mixture with a refractory soot-like core. The soot core increased the calculated extinction by about 10% in the most polluted drier layer relative to a pure sulfate aerosol but had significantly less effect at higher humidities. A 3 km descent through a boundary layer air mass dominated by pollutant aerosol with relative humidities (RH) 10-77% yielded a close agreement between the measured and calculated aerosol optical depths (550 nm) of 0.160 (+/- 0.07) and 0. 157 (+/- 0.034) respectively. During descent the aerosol mass scattering coefficient per unit sulfate mass varied from about 5 to 16 m(exp 2)/g and primarily dependent upon ambient RH. However, the total scattering coefficient per total fine mass was far less variable at about 4+/- 0.7 m(exp 2)/g. A subsequent descent through a Saharan dust layer located above the pollution aerosol layer revealed that both layers contributed similarly to aerosol optical depth. The scattering per unit mass of the coarse aged dust was estimated at 1.1 +/- 0.2 m(exp 2)/g. The large difference (50%) in measured and calculated optical depth for the dust layer exceeded measurements.

Valero, Francisco P. J.↗
NASA NTRS2000

Passive, Direct-Read Monitoring System for Selective Detection and Quantification of Hydrogen Chloride

Monitoring the exposure of an employee to hydrogen chloride or hydrochloric acid in the presence of other acids has been a challenge to the industrial hygiene community. The capability of a device to differentiate the levels of acid vapors would allow for more accurate determinations of exposure and therefore improved occupational health. In this work, a selective direct-read colorimetric badge system was validated for Short Term Exposure Limit (STEL) monitoring of hydrogen chloride. The passive colorimetric badge system consists of a direct reading badge and a color scale. The badge has a coated indicator layer with a diffusive resistance in the shape of an exclamation mark. An exclamation mark will appear if hydrogen chloride is present in the atmosphere at concentrations at or above 2.0 ppm. By using the color scale, the intensity of the color formed on the badge can be further quantified up to 25 ppm. The system was validated according to a protocol based on the NIOSH Protocol for the Evaluation of Passive Monitors. The badge was exposed to relative humidities ranging from 11% to 92%, temperatures ranging from 7 C to 400 C and air velocities ranging from 5 cm/sec to 170 cm/sec. All experiments were conducted in a laboratory vapor generation system. Hydrofluoric acid, nitric acid, sulfuric acid, chlorine, hydrogen sulfide and organic acids showed no effect on the performance of the hydrogen chloride monitoring system. The passive badge and color scale system exceeded the accuracy requirements as defined by NIOSH. At ambient conditions, the mean coefficient of variation was 10.86 and the mean bias was 1.3%. This data was presented previously at the American Industrial Hygiene Conference and Exposition in Toronto, Canada in June 1999.

Chapman, K. B.↗
DOE OSTI2022

Oxygen and Proton Transport in Flooded Graphene Pores with N-Dopants and Defects

Reactant transport is an important consideration in the design of ideal electrode structures. For the oxygen reduction reaction catalyzed by Pt/C in proton exchange membrane fuel cell cathodes, O 2 and H + delivery to Pt surfaces and H 2 O transport away are required. Some Pt nanoparticles may only be accessible via micropores that are too small for ionomer molecules to enter, possibly requiring flooding for H + transport. Here, to test if these “buried” Pt particles can play a role in activity through this proposed transport mechanism, we have performed atomic-scale simulations based on reactive force field molecular dynamics with an emphasis on confinement below 20 Å. Diffusion coefficients as a function of the molar concentration and local environment have been evaluated in water domains confined in two-dimensional graphene nanochannels of various channel heights representing a morphological model for micropores in proton exchange fuel cell cathodes. Our study shows that local atomic-scale structures can strongly modify H + , O 2 , and H 2 O transport rates in flooded micropores less than 20 Å in size. We find that there is a critical crossover in diffusion behavior around the 20 Å spacing with larger pores having bulk-like diffusion properties and confinement below 20 Å monotonically decreases diffusion rates. As pore size decreases, we observe locally dispersed water regions that ultimately strand reactants from long-distance transport. These findings suggest that flooded micropores may in fact be viable transport pathways for relevant reactants and products if the pore walls on opposite sides remain separated by ≥10 Å separation. Furthermore, the confinement effect is so strong that N-doping and C-vacancy defects in the C pore wall have only a minimal impact on diffusion rates and their effects are counterintuitively more apparent at larger spacings. These findings provide valuable insight regarding cathode performance and the role “stranded” catalyst particles may play in fuel cell cathodes.

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

The Effect of Acetylation on Iron Uptake and Diffusion in Water Saturated Wood Cell Walls and Implications for Decay

Acetylation is widely used as a wood modification process that protects wood from fungal decay. The mechanisms by which acetylation protects wood are not fully understood. With these experiments, we expand upon the literature and test whether previously observed differences in iron uptake by wood were a result of decreased iron binding capacity or slower diffusion. We measured the concentration of iron in 2 mm thick wood sections at 0, 10, and 20% acetylation as a function of time after exposure to iron solutions. The iron was introduced either strongly chelated with oxalate or weakly chelated with acetate. The concentrations of iron and oxalate in solution were chosen to be similar to those found during brown rot decay, while the concentration of iron and acetate matched previous work. The iron content of oxalate-exposed wood increased only slightly and was complete within an hour, suggesting little absorption and fast diffusion, or only slight surface adsorption. The increase in iron concentration from acetate solutions with time was consistent with Fickian diffusion, with a diffusion coefficient on the order of 10 -16 m 2 s -1 . The rather slow diffusion rate was likely due to significant binding of iron within the wood cell wall. The diffusion coefficient did not depend on the acetylation level; however, the capacity for iron absorption from acetate solution was greatly reduced in the acetylated wood, likely due to the loss of OH groups. We explored several hypotheses that might explain why the diffusion rate appears to be independent of the acetylation level and found none of them convincing. Implications for brown rot decay mechanisms and future research are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Propagation of elastic waves in a flexomagnetic solid

Flexomagnetism is the coupling between magnetism and strain gradients and is a technologically relevant phenomenon. We present a theory of elastic wave propagation in a linear elastic flexomagnetic material with microstructure and strain gradient elastic interactions. The expressions of frequency, phase velocity, and group velocity of longitudinal and transverse waves are derived and are shown to depend on the flexomagnetic coefficient and microstructure. We also show that the effect of flexomagnetism and microstructure can lead to some interesting phenomena in wave propagation, which are not observed in classical linear elasticity theory of waves. Specifically, in contrast to classical linear elastic materials, where wave propagation is non-dispersive, flexomagnetic materials with microstructure can exhibit both normal and abnormal dispersion. It is also noteworthy that, in flexomagnetic materials with gradient elasticity, the phase velocities of transverse waves can exceed those of longitudinal waves, which is atypical in classical elasticity. Furthermore, waves can also attenuate for a certain range of wavenumbers that depend on the flexomagnetic coefficient and microstructural parameters. Finally, we explore the possibility of waves exhibiting zero group velocity modes, where waves are non-propagating but have strong local energy confinement, negative group velocity modes, where the wave packet moves in the opposite direction to that of wave propagation, and the phenomenon of wave freezing, where a propagating wave stops in space without diffusing or spreading.

Ghosh, Swarnava [ORNL] (ORCID:0000000338005264)↗
NASA NTRS2001

Modelling of Trapped Radiation Near Jupiter

Energetic (62 to approx. 130 MeV) proton fluxes measured with the Galileo Probe inside Jupiter's main ring (radii 122,500 to 128,940 km) show a modest increase as Jupiter is approached. Solutions of a reduced (1-d) equatorial diffusion equation with-constant losses described by a lifetime tau match those data for tau greater than or equal to 10(exp 9) s for diffusion coefficient D(sub LL) = 10(exp -9) L(exp 4)/s, 10(exp -9) L(exp -3)/s, or 10(exp -10) L(exp 4)/s (this particle population may be undergoing nearly loss-free inward radial diffusion). Tau would increase as Jupiter is approached if its value were determined by pitch angle scattering due to EM wave-particle interactions. if the absolute amplitudes of the waves' magnetic fluctuations did not vary with radial distance. Exploration of numerical solutions of the diffusive transport equation for Jupiter's magnetospheric ions at die planetary ring and closer to Jupiter has been done. Explicit range-energy relationships for energy loss in ring matter (SiO2), modeled as a continuous disk, are incorporated. A spatial resolution of approx. 1000 km in the radial direction is used in the vicinity of the main ring; off-equatorial ion fluxes have been included (assuming the atmosphere is a perfect absorber). For protons the maximum effect of energy loss in the microscopic size range of ring particles is expected at energies approx. 0.1 Mev. Presumably there could be sufficiently large, unobserved bodies within the ring that would have to be modeled as discrete objects, as planetary satellites would be. Within limitations of the model and computational resources it's planned to characterize solutions for transport of magnetospheric ions past the planetary ring; the procedures that are used are fairly well known and may be applied generally to the simpler magnetospheric models.

Mihalov, John D.↗
DOE OSTI2025

Preparation and electrochemical properties of high-entropy oxide Li x (CrMnCoNiZn) 3-x O 4

High-entropy oxides (HEOs) have gained significant attention as anode materials for lithium-ion batteries (LIBs) due to their high theoretical specific capacity, synergistic effects of constituent elements, and enhanced structural stability induced by high entropy. Here, this study investigates the electrochemical properties of a novel series of HEOs, Li x (CrMnCoNiZn) 3-x O 4 , synthesized via the sol-gel method. As lithium content increases, the electrochemical lithium storage performance of the HEOs improves, attributed to the rise in high-valence states and oxygen vacancies. The high valence state of elements facilitates greater electron transfer, enhancing the specific capacity during lithium insertion. Additionally, a higher concentration of oxygen vacancies boosts ionic conductivity and lithium ion diffusion kinetics, leading to superior electrochemical performance. Among the synthesized materials, Li 0.11 (CrMnCoNiZn) 2.89 O 4 exhibited the best performance, delivering a discharge capacity of 365 mAh·g −1 after 100 cycles at 100 mA·g −1 , 257 mAh·g −1 after 1000 cycles at 1 A·g −1 , and a rate capacity of 196 mAh·g −1 at 2 A·g −1 . The Li + diffusion coefficient reached 2.78 × 10 −15 cm 2 ·s −1 . This work not only deepens the understanding of HEOs' electrochemical behavior in LIBs but also provides insights into the impact of lithium ion substitution on their performance.

25 ENERGY STORAGE↗
DOE OSTI2020

Carrier recombination and diffusion in high-purity diamond after electron irradiation and annealing

Carrier-transport mechanisms are studied in high-purity diamond irradiated with 6 MeV electrons in the dose range of 10 12 –10 16 cm -2 and annealed at different temperatures up to 1450°C. Lifetimes and diffusion coefficients are extracted using two pump–probe techniques based on free-carrier absorption and transient-grating principles and then correlated with the corresponding defect evolution from spectroscopic measurements. Furthermore, the neutral monovacancy is revealed as the main carrier recombination center in the as-irradiated diamond, providing bipolar carrier lifetimes of a few nanoseconds at the highest irradiation dose. Carrier-capture cross sections are reduced during annealing as vacancies aggregate into divacancies at ≤1000°C and extended vacancy clusters at 1450°C.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Simulation of CEFR neutronic start-up tests with FENNECS

This paper presents simulation results of selected Neutronic Start-up Tests of the China Experimental Fast Reactor (CEFR) obtained by the neutronics code FENNECS that have been performed within the frame of the IAEA Coordinated Research Program I31032. The Finite Element Neutronics code FENNECS is developed at GRS and solves the few-group steady-state and transient diffusion equation using a Galerkin-based finite element approach. Its main purpose is the safety assessment of Small Modular Reactors and Micro Reactors with complex geometry (e.g., rotating control drums) which gain increased interest internationally. Serpent has been applied to create reference models of the CEFR and cross-section libraries in 10 energy groups for FENNECS. Using these libraries, the following experiments have been simulated with FENNECS: net criticality, control rod integral and differential worth, void reactivity effect, subassembly exchange reactivity effects and reaction rate distribution. The obtained satisfactory agreements with the measurements represent a valuable contribution to the validation of FENNECS. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2022

Effect of Sulfur on siderophile element partitioning between olivine and martian mantle primary melt

Ni and Co variations in primary martian magmas exhibit anomalous incompatible behavior, which has remained an unexplained conundrum. Because martian magmas are S-rich, and some trace metals are reported to have enhanced solubility in S-bearing magmas, we have carried out a series of experiments to evaluate the effect of high-S melts on the olivine/melt partitioning of Ni, Co, Mn, V, and Cr. Near-liquidus experiments on a synthetic primary martian mantle melt (Yamato-980459 [Y98]) were completed in a piston-cylinder apparatus at 0.75 GPa. Previous studies in S-free systems illustrate that the partition coefficients for these elements are dependent chiefly on D(Mg(Ol/melt)) (the partition coefficient defined as wt% Mg in olivine/wt% Mg in melt, a proxy for temperature), and were used to calibrate a predictive expression that includes the effects of temperature [i.e., D(Mg(Ol/melt))], melt composition, and oxygen fugacity. These predictive expressions are then used to isolate any effect in DM olivine/melt due to dissolved sulfur. The results show that S might have a small effect for Co, but not enough to change Co partitioning from compatible to incompatible in our experiments. The addition of a sulfur term to the D(Co) predictive expressions shows that nearly 8000 ppm of sulfur would be required in the melt (at liquidus temperature of Y98) for D(Co) to become <1. These S contents are two times higher than those of a sulfide-saturated melt at the P-T conditions of a martian mantle source region. Therefore, the anomalous incompatible behavior observed in these primary magma suites must be due to another mechanism. High temperature, oxygen fugacity, and diffusion are not viable mechanisms, but magma mixing, assimilation, or kinetic crystallization effects remain possibilities.

shergottite↗