Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion behavior”

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 181 records · Page 10

Impact of hydrophobic pendant phenyl groups on transport and co-transport of methanol and acetate in PEGDA-SPMAK cation exchange membranes

CO 2 reduction cells are innovative devices that reduce CO 2 into valuable chemicals (i.e. methanol (MeOH) and acetate (OAc – )) on the cathode. A major challenge with such devices is to develop ion exchange membranes that allow ion-selective transport (i.e., protons for cation exchange membranes, CEM) and suppress the crossover of CO 2 reduction products. To design such membranes, it is important to understand multi-solute transport behavior of these solutes. Previously, our group reported acetate diffusivities in co-transport with MeOH increase in sulfonated CEMs, where we speculated charge screening of the electrostatic interactions by co-diffusing MeOH has an impact. Here, crosslinked membranes fabricated by photopolymerization of poly(ethylene glycol) diacrylate (PEGDA), 3-sulfopropyl methacrylate potassium (SPMAK, SO 3 – -containing ionomer), and a phenyl-containing comonomer either phenoxyethyl acrylate (PEA) or poly(ethylene glycol) phenyl ether acrylate (PEGPEA)) are investigated. We observe OAc – diffusivities to both (1) PEGDA-SPMAK and (2) PEGDA-PEGPEA increase in co-diffusion with MeOH and those to (3) PEGDA-SPMAK/PEGPEA decrease. To rationalize this emergent co-transport behavior, we speculate (1) electrostatic interactions are interfered with by co-diffusing MeOH, (2) chain mobility (segmental dynamics) increases in the presence of MeOH and (3) for films with both hydrophobic (PEGPEA) and hydrophilic (SPMAK) comonomers chain mobility (segmental dynamics) decrease due to interactions between the comonomer sidechains depressing overall solute diffusivities. Finally, while further investigations are needed, this work contributes to improving our fundamental understanding of the relationships between polymer film chemistry, solute chemistry, and emergent cotransport behavior observed and described in this work (and others) towards enabling the design of improved membrane materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of controlled magnetic island bifurcation on electron diffusion

Magnetic islands strongly influence cross-field electron transport in magnetized plasmas. In particular, bifurcations of the island topology modify the number and location of O-points, X-points, and separatrix boundaries, thereby altering diffusion pathways. In recent DIII-D experiments, external magnetic perturbations were used to rotate and periodically bifurcate the island on the q = 2 surface, causing a switchback between a q = 2/1-dominated structure and a narrower q = 4/2-dominated structure. To investigate how this topological change affects electron transport, we employ the field line tracing code TRIP3D with an implemented collisional operator. Thermal, tracer electrons launched from O-points, X-points, and outside separatrix boundaries reveal distinct diffusion regimes, including classical, subdiffusive, and superdiffusive behavior, depending on both the dominant island mode and launch location. These results suggest that island bifurcation can alter electron diffusion across rational surfaces, with direct implications for particle confinement. While the present work emphasizes diffusion as a general framework, the findings provide insight into the conditions under which electron trapping into an island or stochastization of the island's separatrix can enable additional mechanisms, such as the generation of energetic electrons.

DIII-D↗

Atomic Diffusion, Segregation, and Grain Boundary Migration in Nickel-Based Alloys from Molecular Dynamics Simulations

Grain boundary diffusion and metal mobility in alloys control material performance in many applications and yet remain poorly understood at a mechanistic level. With advances in accessible time and length scales for computational molecular simulations, and recent force field developments, we now possess tools to help unravel those mechanisms. Using large-scale molecular dynamics simulations, here we examined vacancy-mediated diffusion processes in Ni-5Cr alloy with low and high-energy grain boundaries. We show that atomic diffusion inside the grain boundary plane is about four times higher than bulk diffusion, at any temperature, and exhibits a typical Arrhenius behavior with a very small energy barrier (0~.8 eV for Cr and 0.7 eV for Ni within 1300-1600 K). Additionally, the fastest diffusing species inverts; Cr diffusion was faster than Ni in the bulk but slower in the grain boundaries. This is attributed to the creation of high cohesive energy clusters of Cr at the grain boundary. Grain boundary migration was also observed to be temperature dependent and appears to be two times higher in the 5% Cr alloy than in pure Ni, highlighting the important role of the alloying element on grain boundary motion.

Simonnin, Pauline GN↗

Electron-electron interactions in low-dimensional Si:P delta layers

Scientists have long studied the physics of highly disordered conducting systems, seeking to understand the multitude of quantum phenomena that govern how electrons move through material systems. Recently, research into silicon-based quantum computing has made disordered conducting systems, such as Si:P monolayers embedded in isotopically pure Si, technically relevant. Consequently, applying and advancing the theoretical frameworks developed to describe electron behavior in disordered systems is a necessary objective in this field of research. This study investigates key components of dopant-based Si quantum computing devices: embedded regions of highly doped delta layers (δ layers). We examine the transport behavior and the electron-electron interaction (EEI) physics in embedded Si:P δ layers by means of self-consistent magnetotransport measurements. Parameters associated with the electronic transport offer a meaningful quantitative characterization of δ-layer quality and dopant diffusion. In addition, by examining EEI behaviors in a set of samples with embedded Si:P δ layers produced with different PH3 exposure procedures prior to Si encapsulation, we show how details of material synthesis affect the dimensionality of charge carrier interactions in embedded Si:P δ layers. Furthermore, the relationship between δ-layer confinement and EEI screening lengths is established here. This understanding will help validate important models used for device simulation and design and lead to improvements in the control of electrostatic gating of and tunneling transport through Si:P single atom transistors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mesoporous Thin Film Architectures: Addressing Material Demands through Molecular Self-Assembly

Mesoporous thin films spark interest across a wide range of disciplines due to their tunable nanostructures, large internal surface areas, and strong compatibility with planar optical, electronic, and microfluidic devices. While attention in the porous materials community has shifted toward macroporous or disordered nanoporous systems, a resurgence in mesoporous thin film research is underway, driven by new molecular self-assembly methods, advanced materials chemistry, and improved characterization techniques. The integration of high-χN block copolymer design, kinetically persistent micelle templating, and postdeposition processing protocols now allows control over structural parameters such as pore size, wall thickness, porosity, and connectivity. These advances have overcome many of the thermodynamic and processing constraints that previously limited widespread adoption. Rather than serving only as high-surface-area supports, mesoporous thin films are engineered as active interfaces where responsive chemistries and nanoscale confinement act in tandem. Embedding switchable ligands, thermoresponsive polymers, redox mediators, or ion-selective groups directly within the pore walls enables real-time control over transport, optical, and electrochemical properties. These capabilities open up new directions in adaptive coatings, gated membranes, and fast-response biosensors. To further expand their functional scope, mesoporous films are integrated into hierarchical and multicomponent architectures. Techniques such as triblock terpolymer templating, crack-directed assembly, and nanoimprint lithography allow for control over spatial organization on the micron and submicron scale and pore system orientation. This enables programmable anisotropy, enhanced molecular diffusion, and wavelength-selective photonic behavior, essential for next-generation sensing, catalysis, and energy applications. Such structural and functional complexity requires equally sophisticated characterization. Multimodal and in situ techniques can track material dynamics under operational conditions. Recent progress includes extended-range ellipsometric porosimetry (EP) for hierarchical architectures, vacuum EP for interface energetics, time-resolved EP for diffusion kinetics, and correlative AFM-SAXS mapping. The introduction of advanced neutron-based spectroscopies, particularly quasielastic neutron scattering (QENS), promises to provide real-time access to ion transport dynamics and segmental motion under nanoscale confinement, offering a path toward deeper mechanistic understanding of structure-performance correlations in mesoporous systems. This Account reflects the technical advances made and the interdisciplinary collaborations that have shaped our collective vision. The particular dimensions of mesopores enable us to subtly tune interactions at the molecular, interfacial, and mesoscopic levels that permit us to harness nanoconfinement. What emerges is a versatile, modular platform capable of chemical gating, energy transduction, and sensing with a level of tunability unmatched by other porous materials. We highlight critical challenges including the need for more robust large-area processing, a deeper understanding of dynamic behavior under cycling, and better integration with device-level architectures. Our strategies support the transition of mesoporous thin films into active high-performance components in next-generation energy, environmental, and biomedical systems.

oxides↗

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On Relative Importance of Vehicular and Structural Motions in Defining Electrolyte Transport

Molecular motions of ionic and neutral electrolyte species determine transport properties at the continuum scale. These molecular motions can be classified as vehicular (e.g., cations moving with a solvation shell of neighboring solvent molecules) and structural (e.g., cations hopping from one solvation shell to another) motions. While literature studies have described the presence, and relative importance, of each of these motions in various electrolytes, a clear link to macroscopic transport properties has not been made. We herein establish this link by using the fluctuation-dissipation theorem to develop theoretical expressions connecting the molecular displacements to Stefan-Maxwell diffusivities. To illustrate the usefulness of the proposed equations, we study LiPF 6 in propylene carbonate as an exemplar electrolyte. We show that its transport behavior improves at all concentrations when structural diffusion of cations is promoted. On the other hand, boosting the cation vehicular diffusion negatively affects the concentrated compositions. We extend this understanding to a generalized electrolyte of a salt dissolved in a solvent. Our theory suggests that while structural diffusion influences Stefan-Maxwell diffusivities globally, vehicular diffusion is only relevant under certain conditions. Such guidelines are critical for a bottom-up design of electrolyte transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rare Earth and Transition Metal Containing Glasses

Transition metal (TM) and rare earth (RE) ions have been incorporated into many glass systems such as silicate, phosphate, and borosilicate-based oxide glasses, as well as in halide and chalcogenide glasses, that find applications ranging from optical, photonic, and magnetic devices, solid-state battery, to nuclear waste disposal. Understanding the structural role of RE and TM in these glasses can help to develop glass compositions for targeted applications with either high-optical emission efficiency, electrical conductivity, or chemical durability. In this chapter, we first provide a general introduction of the applications and structural features of RE and TM in glasses, then the critical aspects of molecular dynamics (MD) simulations of these glasses such as interatomic potentials, structural analysis tools to study RE and TM ions in glasses and their clustering behaviors, Quantitative Structure–Property Analysis (QSPR), diffusion and dynamic property calculations, and electronic structure calculations to understand electronic defects such as charge trapping and radiation effects are introduced. Three representative case studies are presented: the first one is on MD simulations of erbium- and europium-doped silica and silicate glasses, as well as cerium doped aluminophosphate glasses, that revealed the effect of glass composition on RE ion local structure and clustering behavior. Electronic structure calculations of cerium-doped glass show how the existence of multioxidation states help to mediate radiation-induced damages caused by excited electron–hole pairs was also discussed. The second one focuses on alkali vanadophosphate glasses where the existence of two vanadium oxidation states help to provide electronic conduction in the glasses while alkali ions provide ionic conduction. MD simulations were used to understand vanadium environments and other structural aspects in the phosphate glasses, as well as the ionic transport behaviors of alkali ions. The third case study is on zirconium-containing borosilicate and aluminosilicate glasses that find wide applications in nuclear waste disposal. MD simulations help to provide structural details of zirconium ions that are validated by diffraction and EXAFS spectra. The structural information was used to interpret changes of mechanical properties and chemical durability by using QSPR and other analyses-based MD-generated structure models.

Du, Jincheng↗

Eddy diffusivity operator in homogeneous isotropic turbulence

Here, we use the recently developed macroscopic forcing method to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However, for small scales eddy diffusivity is found to vanish inversely proportional to the wave number. Behavior at all scales is reasonably captured by a nonlocal eddy diffusivity operator modeled as D/$\sqrt{I - l^2∇^2}$, where D is the eddy diffusivity in the Boussinesq limit, and l is a constant on the order of the large-eddy length. These results present a direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Extended families of critical and stationary droplets for nonequilibrium phase transitions in spatially discrete bistable systems

Bistable nonequilibrium systems are realized in catalytic reaction-diffusion processes, biological transport and regulation, spatial epidemics, etc. Behavior in spatially continuous formulations, described at the mean-field level by reaction-diffusion type equations (RDEs), often mimics that of classic equilibrium van der Waals type systems. When accounting for noise, similarities include a discontinuous phase transition at some value, p eq , of a control parameter, p , with metastability and hysteresis around p eq . For each p , there is a unique critical droplet of the more stable phase embedded in the less stable or metastable phase which is stationary (neither shrinking nor growing), and with size diverging as p → p eq . Spatially discrete analogs of these mean-field formulations, described by lattice differential equations (LDEs), are more appropriate for some applications, but have received less attention. It is recognized that LDEs can exhibit richer behavior than RDEs, specifically propagation failure for planar interphases separating distinct phases. Herein, we show that this feature, together with an orientation dependence of planar interface propagation also deriving from spatial discreteness, results in the occurrence of entire families of stationary droplets. The extent of these families increases approaching the transition and can be infinite if propagation failure is realized. In addition, there can exist a regime of generic two-phase coexistence where arbitrarily large droplets of either phase always shrink. Such rich behavior is qualitatively distinct from that for classic nucleation in equilibrium and spatially continuous nonequilibrium systems.

97 MATHEMATICS AND COMPUTING↗

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗

Vacancy diffusion barrier spectrum and diffusion correlation in multicomponent alloys

Vacancy diffusion serves a crucial role in many important kinetic behaviors and properties of multicomponent alloys. Essential questions, however, persist regarding how chemical complexity affects diffusion and what unique characteristics, if any, set these alloys apart from traditional metals. Using neural network kinetics model, we study vacancy diffusion in NbMoTa alloy across a broad temperature range (2600 to 800 K). Unlike pure metals, the two key diffusion parameters—diffusion correlation factor f and activation energy ΔG m —are not constant in alloys, but instead substantially decrease with decreasing temperature. This temperature dependence arises from a reduced number of active vacancy jump pathways at lower temperatures, leading to more correlated diffusion. Upon examining vacancy diffusion throughout the entire compositional space of the Nb-Mo-Ta system, we discover that the slowest vacancy diffusion surprisingly occurs in the non-equimolar region, rather than the equimolar concentration where the configurational entropy is highest. The diffusion barrier spectrum, characterizing the diffusion energy landscape, is an intrinsic material characteristic, which controls both f and ΔG m and, thereby, the diffusivity. Lastly, we find that the vacancy diffusion rate drops noticeably in the presence of local chemical order in the NbMoTa system, particularly for MoTa alloys with long-range B2 order.

Defects↗

Merging of transport theory with the time-dependent Hartree-Fock approach: Multinucleon transfer in U + U collisions

Multinucleon transfer mechanism in the collision of 238 U + 238 U system is investigated at E c.m. = 833 MeV in the framework of the quantal diffusion description based on the stochastic mean-field approach. Double cross sections σ (N, Z ) as a function of the neutron and proton numbers, the cross sections σ (Z ) and σ (A) as a function of the atomic numbers and the mass numbers are calculated for production of the primary fragments. The calculation indicates the 238 U + 238 U system may be located at an unstable equilibrium state at the potential energy surface with a slightly negative curvature along the β stability line on the (N, Z ) plane. Here, this behavior may lead to rather large diffusion along the β stability direction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Isothermal and Cyclic Oxidation of Haynes 282 Processed by Electron Beam Melting (EBM) and Laser Powder Bed Fusion (LPBF) in Dry Air at 800 and 950°C

Here, the isothermal and cyclic oxidation behavior of Haynes 282 (H282) processed by laser powder bed fusion and electron beam melting was compared to its cast counterpart during exposures in air at 800 and 950°C. The specific microstructure of the AM alloys compared to coarse-grain cast 282 resulted in differences in oxidation rates, internal oxidation and spallation behavior. At 800°C, faster Ti diffusion to the surface due to the smaller grain size of the AM alloys led to Ti doping of the Cr 2 O 3 scale and faster oxidation rates. The higher density of grain boundaries also resulted in more pronounced internal oxidation. At 950°C, a duplex Cr 2 O 3 scale was observed with a dense inner layer and porous cracked outer layer with embedded TiO 2 particles. The impact of the external and internal oxide compositions and oxidation-induced elemental depletions on the oxidation lifetime, creep properties and spallation behavior is discussed.

36 MATERIALS SCIENCE↗

Fundamental Studies of Tritium Formation and Diffusivity in Pure and Defective Zircaloy-4 Getters

Zirconium (Zr) and its alloys are used as fuel cladding in nuclear reactors. These materials have their low absorption cross-section for thermal neutrons and have good mechanical and thermal properties. In addition, they are used as tritium getters in tritium-producing burnable absorber rods (TPBARs), in the form of nickel-plated Zircaloy-4 tubes. To produce 3 H, the breading blankets such as lithium aluminate are irradiated in pressurized water. 3 H produced in this process is captured by the getter and 3 H reacts chemically with Zr to form metal hydride (ZrTx). The precipitation of metal hydride results in not only the volume increase and hydrogen embrittlement of zircaloy but also changes in getter’s chemical properties. This degrades the getter’s performance. Therefore, understanding the behavior of 3 H (such as diffusion and solubility) in zircaloy is important to understand getter performance.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

U-rich U-Mo Solidus, Liquidus, Enthalpy, and Thermal Diffusivity

Several material properties have been studied as a function of composition as part of the USHPRR effort to use U-10Mo for LEU monolithic fuel foils, with special attention to areas where simulation outcomes can be improved. The liquidus and solidus phase boundaries are important for engineering design of processing parameters and for predicting the degree of microsegregation within initial castings. The numerous phase diagram constructions available in the literature differ significantly in their solid/liquid boundary definitions, therefore, an experimental strategy was employed in order to resolve some of the discord. That strategy focuses not on proving any singular point, but instead to provide data with highly characterized uncertainty. The data presented here is the most consistent with the U-rich liquidus and solidus boundaries provided by H. Okamoto (2012), Berche et al. (2011), and Mardon et al. (1959). The recently suggested change in the peritectic reaction temperature (to 1302 K) suggested by Berchè and adopted by Okamoto is also supported, but with lower confidence. The enthalpy of fusion is a basic characteristic of the solidification reaction property for which is required for both research-oriented predictions of behavior and development of applied design strategies for casting engineering. It has been determined that an admixture model for the heat of fusion between U BCC and Mo BCC is a reasonable assumption within the composition region of interest, though limitations of the measurement accuracy impacts the ability to support any trends more detailed than this baseline behavior. A limited number of thermal diffusivity measurements were also performed into the gamma region to provide additional information for solidification modeling and simulations as well as a check on the component properties of thermal conductivity, density, and heat capacity.

36 MATERIALS SCIENCE↗

Frequency-domain vs time-domain SP{sub N} equations to simulate neutron noise

Inside the reactor core mechanical vibrations of fuel assemblies can produce high fluctuations around a steady-state configuration, known as neutron noise. This effect can cause the triggering of power reduction measures. Classically, diffusion theory has been used to simulate this behavior. However, this equation has some limitations if the materials of the reactor have strong variations. In this work, we use the diffusive time-dependent simplified spherical harmonics equations that improve the previous results without the necessity of using high computational requirements. In particular, two types of analyses with these equations (SP{sub 3}) are made: a frequency-domain and a time-domain. A numerical neutron noise benchmark tests the methodology and compare both formulations. First, numerical results show a good agreement between the amplitudes and phases of the SP3 equations computed with the frequency-domain and time-domain. Therefore, as the frequency-domain computation only requires to solve a linear system, it is a recommendable option for neutron noise computations. Second, one can conclude that for this type of nuclear systems, where the assemblies are not homogenized, the SP{sub 3} approximation results improve considerably the accuracy of the diffusion theory.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗