Spherical harmonic–based DEM in LAMMPS: Implementation, verification and performance assessment
Not Available
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.
Not Available
Boiling point elevation (BPE) for high salinity solutions is a crucial parameter that significantly affects the overall thermal efficiency of desalination, water treatment, and many other industrial processing systems, especially when minimum or zero liquid discharge is a key requirement. This paper addresses the discrepancies among the empirical correlations reported in the literature that were commonly used to determine BPE for seawater solutions. An evaluation of the experimental data with low solution concentration reveals that the data cannot be extrapolated to determine the BPE at high concentration (X > 140 ppt or g/kg). Besides, disparities are established within different equations used to determine the BPE for high salinity seawater. As a result of the minimal availability of experimental data of BPE for X > 140 ppt and the need to eliminate the discrepancies, an experimental setup was built, which can produce reliable and repeatable BPE data for any solution concentration. Boiling points were obtained experimentally for seawater solutions with a wide range of salinity (X between 20 and 280 ppt) and brine solution temperature (t b between 20 and 80°C). The experimental results of boiling point elevations for seawater are further compared with the experimental data and analytical results reported in the literature. The newly obtained experimental results show that the minimum and maximum BPE values for the brine solutions tested are 0.175°C at t b = 80°C and X = 20 ppt, and 11.6°C at t b = 80°C and X = 280 ppt, respectively. Furthermore, the measured BPE data for seawater concentration over 80 ppt are fitted and presented using an appropriate fitting model. The new correlation can accurately predict the BPE for high salinity seawater brines and thus should be very helpful in the design of high-energy-efficient evaporation systems.
Not Available
Not Available
We verify a recently-developed nonperturbative guiding center formalism to charged particle dynamics in fields with two-parameter continuous symmetry groups. This entails finding exact constants of motion, valid in the nonperturbative regime, that agree with Kruskal’s adiabatic invariant series to all orders in the perturbative regime, when the field scale length is large compared with a typical gyroradius. We demonstrate that the nonperturbative guiding center model makes exact predictions in these cases, even though it eliminates the cyclotron timescale, thereby establishing a theoretical baseline for performance of the nonperturbative formalism.
The mesh-based Rigorous-Two-Step (R2S) method has been widely used in the accurate estimation of the shutdown dose rate (SDR) of fusion systems. Several mesh-based R2S code has been implemented based on Monte Carlo particle transport code MCNP5 or DAG-MCNP5 and inventory calculation code such as FISPACT-II, ALARA or ACAB. Recently, the OpenMC, a community-developed Monte Carlo particle transport code, implemented the photon particle transport capability, which enables OpenMC to be used in R2S. In this paper, PyNE R2S is extended with the support of DAG-OpenMC, i.e., OpenMC with DAGMC. Modifications that allow PyNE to read the flux result in the state point file of OpenMC has been implemented to PyNE R2S workflow. A custom source routine of OpenMC that reads the photon source file from PyNE R2S and sampling photon particle in OpenMC photon transport has been implemented in OpenMC. With these modifications, PyNE R2S is now capable of performing R2S calculations with both DAG-MCNP5 and DAG-OpenMC. The ITER FNG dose rate benchmark problem has been used to validate the code. The FNG neutron source has been modified for neutron transport with DAG-OpenMC. The shutdown dose rate of 19 cooling times has been calculated and compared with the experimental data and the computational results of other R2S codes. The results of PyNE R2S with DAG-OpenMC show satisfactory agreement with the experimental data and other computational results. Therefore, we consider PyNE R2S with DAG-OpenMC is a reliable code to calculate the SDR of fusion systems.
Not Available
This study develops and validates a sequentially coupled thermo-mechanical welding simulation for the three-pass 316L stainless steel NeT TG4 benchmark weldment using the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) and the Nuclear Engineering Material model Library (NEML). A diffused ellipsoidal heat source was calibrated against thermocouple data and weld macrographs to accurately model the fusion zone geometry and transient thermal fields. Material hardening is represented using the Lemaitre-Chaboche mixed isotropic-kinematic hardening model, while four annealing models - no annealing, single-stage at 1050 °C and 1300 °C, and two-stage at 800 °C/1300 °C - were implemented to assess the impact of annealing models on the accuracy of the predicted welding-induced plasticity, distortions, and residual stresses. The predictions were validated against experimental measurements and benchmarked against results from commercial software, demonstrating that thermo-mechanical MOOSE welding simulations achieve comparable accuracy with enhanced computational efficiency. This work highlights the potential of using open-source finite element frameworks like MOOSE for advanced manufacturing simulations.
Numerical simulations become a necessity when experimental approaches cannot cover the required physical and time scale of interest. One of such area is a simulation of long-term host rock behaviors for nuclear waste disposal and simulation tools involved in the assessment must go through rigorous validation tests. The DECOVALEX project (Development of COupled models and their VAlidation against EXperiments) is dedicated to this purpose by international participants. This work is part of the ongoing phase DECOVALEX–2023 (D–2023, Task G) particularly aiming to simulate fracture behaviors under various conditions. Here, in this study, we cross-verified a variety of numerical methods including continuous and discontinuous approaches against four benchmark exercises with emphasis on numerical accuracy and parameterization of the various numerical approaches. The systematic inter-comparisons of test cases highlight advantages and disadvantages of the different numerical models. Numerical details on discretization effects (e.g. mesh density and orientation) and domain size were investigated in detail for practical applications. It became evident that meticulous attention to mesh resolution and domain size is imperative for achieving accurate numerical simulations, even for static cracks. Moreover, when comparing numerical methods to closed-form solutions for static cracks, all models successfully reproduced the maximum crack opening but encountered challenges near the crack tips. Finally, the paper discusses how to convert between and therefore compare parameters of various numerical approaches. Our benchmark studies reveal that each model necessitates a distinct number of parameters, even in simple scenarios like static crack aperture benchmarks. It is generally more practical to employ fewer parameters to mitigate model over-parameterization and enhance experimental feasibility.
Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.
Not Available
Explore the source record for details and available documents.
Standard UO 2 fuel pellets used in light-water reactors fracture during irradiation due to the large thermal gradient in the radial direction. Over the decades, numerous researchers have explored fuel cracking from experimental and modeling points of view. To date, there have been both empirical and mechanistic approaches to predict the number of fragments that form in UO 2 . The empirical models only consider maximum power and burnup as inputs. Existing mechanistic approaches for normal operation have not accounted for irradiation effects. Here, this work employs a mechanistic fuel cracking model using the extended finite element method to explore radial crack formation while including a sensitivity analysis that accounts for the randomization of tensile strength within the fuel, the strength randomization criteria (uniform or volume-weighted Weibull), power ramping rates, computational mesh density, maximum power level, and irradiation (burnup) effects. The results indicate that the uncertainty in this mechanistic modeling approach envelopes the predicted values from three different empirical correlations in almost all cases. This means that, for computationally intensive analyses involving UO 2 fragmentation, the empirical correlations can be used. However, since the mechanistic calculations bound those of the empirical correlations, there is confidence in the applicability of the mechanistic approach developed in this work to generate a correlation for fuel types where limited data exists (e.g., doped-UO 2 , U 3 Si 2 ).
Not Available
A portable In-Field Alpha Spectrometry (IFAS) system was developed to support on-site uranium enrichment determination for samples taken at facilities working with gaseous UF 6 . It consists of a modified commercial alpha spectrometer, a vacuum pump, a laptop computer, and associated connectors and cables. The IFAS uses specially designed Single-Use Destructive Assay (SUDA) samplers to capture gaseous UF 6 and convert and store it as UO 2 F 2 in a thin-film zeolite. To examine the quantitative performance of the IFAS system a large number of replicate SUDA samples was acquired for natural uranium (NU, 0.72 wt%, 23 samples), low enriched uranium 1 (LEU1, 2.33 wt% 235 U, 10 samples), and low enriched uranium 2 (LEU2, 4.62 wt% 235 U, 10 samples). The standard deviation (SD) values for uranium enrichment determination with the IFAS for NU, LEU1, and LEU2, in 20-hour assays, were 0.037 wt% 235 U, 0.063 wt% 235 U, and 0.072 wt% 235 U, respectively. In addition, a separate analysis suggests that comparable SD values can be achieved in less time, approximately 8-10 hours. For context, these SD values are well below the International Atomic Energy Agency’s International Target Values for the determination of uranium enrichment in non-destructive assay measurements.
Not Available
Not Available