SEARCH · Engineering Papers
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.
Using our open-source programming environment Morpho, an energy-optimization method is developed to compute equilibrium shapes of an arbitrarily confined swollen hydrogel. The internal strains and contact forces are measured, probing its mechanics.
Process monitoring in laser powder bed fusion additive manufacturing can provide insights into stochastic anomalies, melt pool and plume dynamics, and part quality. Swelling, a build anomaly where overbuilt material protrudes through the powder layer after recoating, is readily detectable in post-recoat visible light images of the powder bed. Here, this work identifies several of the underlying mechanisms driving swelling formation by analyzing the influence of processing parameters, laser scan paths, and build plate locations on the presence of swelling detected in situ. Swelling near the edge of the part and swelling in the internal region of the part are shown to correlate with different process conditions. Edge and internal swelling may be driven by different phenomena, with edge swelling predominately occurring on the edge of a part facing the laser module and correlated to clusters of near-surface voids (detected with X-ray computed tomography). A larger spot size, higher laser power, and lower scan velocity also increased the presence of edge swelling. Laser spot size and scan path influenced internal swelling, which occurred preferentially with a larger spot size and in regions with large melt pools, caused by localized heat accumulation due to non-optimal processing parameters or scan path strategies. For coupons processed with a slicer-defined maximum scan vector length, swelling seldom occurred at internal vector-stripe boundaries. These results provide a mechanistic understanding of how swelling can be linked to material flaws, insight into how some instances of swelling can be avoided, and evidence supporting the use of swelling as an in situ indicator for quality assurance and part qualification.
Polymeric gels crosslinked by DNA sequences can exploit DNA strand-displacement reactions to promote swelling through dynamic polymerization. The degree of swelling and the rate of swelling must be directly tunable to achieve the promise of programmable soft matter. Though the kinetics of the strand-displacement reaction provide insertion rates up to 10 4 /Molar/second as measured in bulk solution, DNA hydrogel swelling can take upwards of 30 h to complete. Computational modeling of the reaction-induced swelling of these gels with our recently-developed reactive electrochemomechanical theory (Zimmerman et al., 2024) suggests that their extraordinarily slow swelling is partly due to a scaling mismatch between the addition of charge and the addition of fluid volume, leading to a large transient increase in the fixed charge density. The significant increase in the gel’s fixed charge density, due to the binding of negatively charged DNA, sharply restricts the concentration of mobile hairpins through the phenomenon of Donnan charge exclusion, an effect commonly exploited in nanofiltration applications using polymeric membranes. The scaling problem is overcome when the mean additional swelling provided to the hydrogel by addition of a crosslink is above a critical value, thus the swelling outpaces the charge accumulation, leading the fixed charge density to drop and significantly accelerating the swelling process. This study shows that Donnan exclusion can explain the kinetics of DNA hydrogel swelling, and studies ways to modulate the reaction speed by either modifying the salt concentration or increasing or decreasing the number of base pairs in each DNA sequence.
Ferritic/martensitic (F/M) steels are being considered as potential structural materials for next generation nuclear reactors, and variants of the alloy HT9 are some of the most promising candidates. In this study, two conventional and two composition engineered HT9 alloys were irradiated using 3.5 MeV Fe 2+ up to 600 peak displacement-per-atom (dpa) at 450 °C. Void swelling and microstructure evolution were characterized for each alloy and compared. The two conventional HT9 alloys (INL and ACO3) showed similar void swelling behavior due to their similar elemental composition and processing conditions. The INL HT9 exhibited a maximum of 2.4% swelling and the ACO3 HT9 showed a maximum of 2.8% swelling at 342 and 393 average local dpa, respectively. On the other hand, the two-composition engineered HT9 alloys with varying N contents (10 ppm for low N and 440 ppm for high N) showed disparate swelling behavior. The low N HT9 exhibited a maximum of 4.6% swelling, while the high N HT9 showed a maximum of 0.7% swelling at 342 average local dpa. Changes in the N content also affected Ni/Si rich G-phase formation. The low N HT9 showed a larger size and lower density of G-phase precipitates compared with the high N HT9 after 600 peak dpa irradiation. This study compares the void swelling behavior of the ion irradiated four current HT9 alloys to extremely high doses, with the void swelling data from neutron irradiated HT9 alloys. Furthermore, the comparison lends critical insights into how well these current alloys can withstand high neutron fluxes in future reactors, especially since the low N and high N HT9 alloys have never been exposed to such high doses before.
Many rate theory models of cavity (void) swelling have been published over the past 50 years, all having the same, or similar, structures. A rigorous validation of the models has not been possible because of the dearth of information concerning the microstructures that correspond with the swelling data. Whereas the lack of microstructure information is still an issue for historical swelling data, in the past 10–20 years data have been published on the evolution of the microstructure (point defect yields from collision cascades, cavity number densities, and dislocation densities/yield strengths) allowing certain gaps in information to be filled when considering historic swelling data. With reasonable estimates of key microstructure parameters, a standard rate theory model can be applied, and the model parameter space explored, in connection with historical swelling data. By using published data on: (i) yield strength as a function of dose and temperature (to establish an empirical expression for dislocation density evolution); (ii) cavity number densities as a function of temperature; and (iii) freely migrating defect (FMD) production as a function of primary knock-on atom (PKA) spectrum, the necessary parameter and microstructure inputs that were previously unknown can be used in model development. This paper describes a rate-theory model for void swelling of 316 stainless steel irradiated in the EBR-2 reactor as a function of irradiation temperature and neutron dose.
This work describes the development of a micromechanical-based constitutive model accounting for the effect of internal expansion on the residual elasticity of a Hashin composite material. This material is made of spherical inclusions that are subjected to gradual swelling within a quasi-brittle matrix. The main focus of this work is to describe and analyze the material mechanical response, with an additional focus on the internal swelling’s effect on the stress–strain response and residual elasticity. Microstructural features and parameters of major importance for the mechanical responses were identified. The innovative characteristics of the proposed approach are summarized as follows: (1) a full determination of the physics of a complete-damage problem throughout the whole process of inclusion swelling with upscaling techniques, which transfers the microcrack-related properties from the lower scale to upper scale; and (2) an evolution of the mechanical fields and the corresponding residual elasticity for various inclusion swelling levels. Concerning the matrix–inclusion composite, it was hypothesized that only the matrix was susceptible to cracking, with varied degrees of damage, whereas the inclusions behave elastically and the elastic modulus of the expanding inclusions remains constant. It is to emphasize that the gradual swelling of inclusions is modeled by an increasing strain in the current micromechanical-based constitutive model, and the microcracks are represented by a set of randomly oriented penny-shaped microcracks with identical radii (namely, closed cracks). The main contribution of the current research is to establish the exact mathematical solutions for the mechanical fields (stress, strain) caused by the swelling of the inclusions (mechanical loading) and derive the effective residual elasticity of a composite (structural response) subject to internal expansion and quasi-brittle damage. Based on the assumption of closed cracks that could be extended to open cracks in the upcoming work, the results proposed by the present paper help to understand the non-linear mechanical behavior of quasi-brittle materials subject to microcracking and provide a theoretical framework to be used as academical benchmark for numerical simulations.
Machine learning methods are applied to predict the volumetric swelling rate of the nuclear fuel uranium nitride (UN) over various temperatures, irradiation conditions, and power densities. Both kernel-based methods and symbolic regression models for UN swelling are developed and compared with multiple experimental datasets. We find that the UN pellet geometry and dimensions must be taken into account to accurately model swelling behavior. Strong agreement is observed between the developed machine learning models and the data. The predictive error generated by the machine learning models improves on empirical models taken from the literature. Sensitivity analysis is performed to determine which properties such as temperature, burnup, and power density, are most important in the swelling process. We find that machine learning can be used to quickly develop accurate swelling models for nuclear materials. In conclusion, the presented results illustrate the potential of machine learning to determine volumetric swelling in UN.
With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.
Abstract A long-standing and critical issue in the field of irradiated structural materials is that void swelling is significantly higher in face-centered cubic-structured (fcc) materials (1% dpa −1 ) as compared to that of body-centered cubic-structured (bcc) materials (0.2% dpa −1 ). Despite extensive research in this area, the underlying mechanism of the difference in swelling resistance between these two types of materials is not yet fully understood. Here, by combining atomistic simulations and STEM imaging, we find stacking fault tetrahedra (SFTs) are the primary cause of the high swelling rate in pure fcc copper. We reveal that SFTs in fcc copper are not neutral sinks, different from the conventional knowledge. On the contrary, they are highly biased compared to other types of sinks because of the SFT-point defect interaction mechanism. SFTs show strong absorption of mobile self-interstitial atoms (SIAs) from the faces and vertices, and weak absorption of mobile vacancies from the edges. We compare the predicted swelling rates with experimental findings under varying conditions, demonstrating the distinct contributions of each type of sink. These findings will contribute to understanding the swelling of irradiated structural materials, which may facilitate the design of materials with high swelling resistance.
Void swelling can develop in materials under persistent irradiation when nonequilibrium vacancy and self-interstitial populations migrate under sufficiently asymmetric interaction biases. In conventional metals, the propensity is determined to the first approximation by comparing point-defect relaxation strains. We thus present DFT-based calculations of structures and formation energies and volumes of point defects in the α and the δ phases of plutonium. Here we discuss the pros and cons of various levels of electronic structure theory: spin polarization, spin-orbit coupling, and orbital polarization. Our results show that lattice defects in δ-Pu, in contrast to most fcc metals, have surprisingly small formation volumes. Equally unexpected are the large defect formation volumes found in the low-symmetry α-Pu phase. Both these unusual properties can be satisfactorily explained from defect-induced spin/orbital moment formation and destruction in the Pu phases. Surprisingly, the point defects in α-Pu are found to induce far larger transformation of the local electronic structure than in δ-Pu. When we use the calculated defect properties to estimate the classic void swelling bias in each of the phases, we find it to be unusually small in δ-Pu but likely much larger in α-Pu. Hence, swelling rates and mechanisms can diverge dramatically between the different phases of Pu. Especially in the transient regime before the formation of large defect clusters, the swelling rate of α-Pu can reliably be expected to be much larger than δ-Pu. However, accurate forecasts over longer times will require the conventional void-swelling theory to be modified to handle the complexities presented by the different Pu phases. As a case in point, we show the possible anomalous temperature dependence of vacancy properties in δ-Pu, caused by entropic contributions from defect-induced spin-lattice fluctuations. Such complications may affect defect-defect interactions and thus alter the void swelling bias.
The VTR core has six reactivity control assemblies and three safety assemblies. The control assemblies or primary control rods are adjusted during the normal operation to balance the core reactivity and to control the reactor power. A typical control assembly radial layout is presented in Figure 1. The figure shows the swelled absorber (B 4 C) rod. Initially, helium gas fills the gap between the pin and the cladding before irradiation swelling takes place. For VTR, HT9 steel was selected as the cladding and duct material. The main neutron absorbing material used in the VTR is B 4 C. When residing in the core, the neutronics, thermophysical, and mechanical properties of the materials used in a control assembly will degrade due to accumulated neutron damage. Material degradation limits how long a control assembly can reside in the core. Many phenomena affect the control assembly lifetime, such as the loss of reactivity worth due to B 4 C depletion, the mechanical interaction of the absorber rod and the cladding due to B 4 C swelling, the helium gas buildup in the pin due to B-10 capture, etc. B 4 C swelling, which causes closure of the gap between the absorber rod and the cladding, is usually considered as the main limiting factor from past experience. An initial study was conducted at PNNL to evaluate the irradiation behavior of a VTR control assembly. The evaluation was performed using the CNRD2 code that was initially developed for the FFTF. The study also included an assessment of the VTR control assembly and focused on a 61-pin control assembly design, which is different from that used (37-pin design) in the core design study. The study conducted by PNNL was reviewed independently by ANL. A Python script referred to as the Control Assembly Evaluation Script (CAES) was developed for the independent review and additional assessment of 37-pin control assembly design. The script has focused on the assessment of the absorber rod swelling for its importance in determining the control assembly lifetime. CAES uses geometry, neutronics, materials data as input to predict the swelling of the absorber rod during its residence in the reactor core. The results from CAES showed some non-negligible differences against the PNNL results. Some of the differences can be attributed to the different interpretation of the control rod assembly dimensions. To resolve this issue, a benchmark exercise was proposed. The benchmark specification was developed by PNNL. The benchmark exercise was performed independently at PNNL and ANL using different codes/scripts (CRND2 and CAES). This memo documents the results calculated using the different codes. However, this report is limited to presenting the results obtained. Further investigation of the cause of the observed difference will be performed as part of future activities, pending continuation of the VTR program.
This report describes the work in NEAMS (Nuclear Energy Advanced Modeling and Simulation) to develop a mechanistic fission gas model for uranium nitride fuels in BISON. The existing Sifgrs (Simple integrated fission gas release and swelling) model tracks the average properties of two bubble populations in the bulk and at the grain boundaries. It was recognized that dislocations play a crucial role in the fission gas swelling of UN, and an irradiation-induced dislocation density model was needed, as well as a model describing how fission gas interacts with dislocations creating a third population of bubbles along dislocations. A mechanistic model that tracks dislocation bubbles was implemented in Sifgrs. This model was used to simulate fission gas swelling and release in UN. Lower-length-scale calculations and experimental observations from carbide fuel were leveraged to help populate the model with essential parameters. As a placeholder, an empirical function was formulated for the evolution of the dislocation network. Because of the difference in evolution of this network at different temperatures, the dislocation bubbles are able to capture behavior that the bulk intragranular bubbles cannot. It was found that the bulk bubbles dominate microscopic swelling at low temperatures, and the dislocation bubbles dominate at higher temperatures. Based on this model, the transition between the two bubble types is the main factor behind the breakaway swelling phenomenon in UN. In order to test the model, a couple of assessments were run. For the lower temperature JOYO pins, the model produces reasonable fission gas release and swelling values. For the higher temperature SP1 pin, the model dramatically underestimates the fission gas release. This issue is attributed to the gas being trapped inside dislocation bubbles, unable to escape to grain boundaries to cause fission gas release, and could be remedied by a mechanistic dislocation model that allows the dislocation density to decrease at very high temperatures. A preliminary mechanistic dislocation model was developed supported by first-principle calculations. These calculations provided valuable insight into how interstitial defects cluster in UN in the {110} orientation, and may eventually form dislocation loops, leading to the conclusion that dislocation loops may nucleate from these clusters. An estimate of the dislocation line energy in UN was obtained and will be improved in future work. The free-energy cluster dynamics code Centipede was used to track interstitial and vacancy absorption at dislocation loops and calculate their growth. In addition, improvements to Centipede were made including new convergence criteria for transient simulations and sink driving force updates.
We report Ferritic-martensitic (FM) steels for in-vessel components in proposed fusion reactors are expected to suffer from high levels of displacement damage and helium (He) generation by neutron transmutation. However, a thorough understanding of the He synergistic effects on the cavity swelling in FM steels is still not well established. To gain fundamental insights into the He effect on cavity swelling, high purity Fe and Fe-10 wt.% Cr ferritic model alloys were irradiated with 8 MeV Ni ions and co-implanted He ions at 400 to 550 °C up to 30 displacements per atom (dpa) with He implantation rates of 0.1, 10 and 50 appm He/dpa. The current study focuses on the 50 appm He/dpa, 500 °C behavior vs. 0.1 and 10 appm He/dpa. Irradiation-induced defects, including cavities, dislocation loops, and dislocation networks were characterized using transmission electron microscopy (TEM). In the grain interior, a bimodal cavity size distribution was observed in the 10 and 50 appm He/dpa samples, but not for 0.1 appm He/dpa. Cavity swelling was maximized at intermediate He implantation rates of ~10 appm He/dpa for both ion-irradiated Fe and Fe-10Cr alloys. The cavity swelling behavior as a function of He implantation rate appears to be controlled by the He/dpa-dependent variation of cavity sink strengths. Treating small bubbles as biased sinks for interstitial absorption can significantly increase the ratio of biased to unbiased sink strengths (Q) and results in maximized cavity swelling for a Q ratio close to one.
This work is the third and final part in an initial series on addressing the behavior of MX precipitate stability in an advanced Fe-9Cr reduced activation ferritic/martensitic (RAFM) alloy under fusion-relevant ion irradiation conditions. Here, the helium trapping properties of MX precipitates are investigated across varying damage levels (15–100 dpa), temperatures (400–600 °C), and helium doses (10–25 appm He/dpa) using sophisticated dual ion beam experiments and electron microscopy. Results indicate that individual MX precipitates efficiently sequester helium in the form of nanoscale bubbles at the precipitate-matrix interfaces near the peak swelling temperature (∼5 bubbles/precipitate at 500 °C). Swelling was primarily due to matrix cavities. The Fe-9Cr alloy reached 2% swelling by 100 dpa, suggesting a shift to steady-state swelling around 50 dpa at 500 °C. Furthermore, MX precipitate dissolution beginning at 15 dpa did not coincide with this onset of steady-state swelling.
Experiments indicated that metallic fuel in sodium-cooled fast reactors (SFRs) rapidly swells radially and axially at low burnup. Despite that, prior studies have been focused on describing high burnup axial fuel elongation. With recent conventional and non-conventional metallic fuel concepts being considered for license applications, understanding multidimensional fuel swelling at a wide range of burnup levels is important to fuel analysis and qualification activities. Here, we report the development and demonstration efforts of a low-burnup SFR metallic fuel swelling model evaluation framework using the BISON advanced fuel performance code. The framework leverages the Integral Fast Reactor (IFR) program X423 experiment data through the ongoing integration project to enable standardized and automated use of legacy metallic fuel irradiation data maintained in the SFR fuel irradiation and physics database (FIPD) for BISON metallic fuel model verification and validation. In conclusion, the performance of the framework was demonstrated using the two representative metallic fuel swelling model sets implemented in BISON, with a series of insights about future advanced swelling model development.
The development of biomolecular stimuli-responsive hydrogels is important for biomimetic structures, soft robots, tissue engineering, and drug delivery. DNA polymerization gels are a new class of soft materials composed of polymer gel backbones with DNA duplex crosslinks that can be swollen by sequential strand displacement using hairpin-shaped DNA strands. The extensive swelling can be tuned using physical parameters such as salt concentration and biomolecule design. Previously, DNA polymerization gels have been used to create shape-changing gel automata with a large design space and high programmability. Here we systematically investigate how the swelling response of DNA polymerization gels can be tuned by adjusting the design and concentration of DNA crosslinks in the hydrogels or DNA hairpin triggers, and the ionic strength of the solution in which swelling takes place. Finally, we also explore the effect hydrogel size and shape have on the swelling response. Tuning these variables can alter the swelling rate and extent across a broad range and provide a quantitative connection between biochemical reactions and macroscopic material behaviour.
Semiconducting conjugated polymers bearing glycol side chains can simultaneously transport both electronic and ionic charges with high charge mobilities, making them ideal electrode materials for a range of bioelectronic devices. However, heavily glycolated conjugated polymer films have been observed to swell irreversibly when subjected to an electrochemical bias in an aqueous electrolyte. The excessive swelling can lead to the degradation of their microstructure, and subsequently reduced device performance. An effective strategy to control polymer film swelling is to copolymerize glycolated repeat units with a fraction of monomers bearing alkyl side chains, although the microscopic mechanism that constrains swelling is unknown. Here we investigate, experimentally and computationally, a series of archetypal mixed transporting copolymers with varying ratios of glycolated and alkylated repeat units. Experimentally we observe that exchanging 10% of the glycol side chains for alkyl leads to significantly reduced film swelling and an increase in electrochemical stability. Through molecular dynamics simulation of the amorphous phase of the materials, we observe the formation of polymer networks mediated by alkyl side-chain interactions. When in the presence of water, the network becomes increasingly connected, counteracting the volumetric expansion of the polymer film.
Bilayer systems comprising lipid mixtures are the most well-studied model of biological membranes. While the plasma membrane of the cell is a single bilayer, many intra- and extra-cellular biomembranes comprise stacks of bilayers. Most bilayer stacks in nature are periodic, maintaining a precise water layer separation between bilayers. That equilibrium water separation is governed by multiple inter-bilayer forces and is highly responsive. Biomembranes re-configure inter-bilayer spacing in response to temperature, composition, or mass transport cues. In synthetic bilayer systems for applications in cosmetics or topical treatments, control of the hydration level is a critical design handle. Herein we investigate a binary lipid system that leverages key inter-bilayer forces leading to unprecedented levels of aqueous swelling while maintaining a coherent multilamellar form. We found that combining cationic lipids with bicontinuous cubic phase-forming lipids (lipids with positive Gaussian modulus), results in the stabilization of multilamellar phases against repulsive steric forces that typically lead to bilayer delamination at high degrees of swelling. Using ultra-small-angle X-ray scattering alongside confocal laser scanning microscopy, we characterized various super-swelled states of 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP) and glycerol monooleate (GMO) lipids, as well as other analogous systems, at varied concentration and molar ratios. Through these experiments we established swelling profiles of various binary lipid systems that were near-linear with decreasing lipid volume fraction, showing maximum swelling with periodicity well above 200 nanometers.