Search NASA⌕ Search

SEARCH · Search NASA

Results for “test particle simulations”

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

Selective capture and recovery of uranium oxide colloids from aqueous soil suspensions using high gradient magnetic filtration

High Gradient Magnetic Filtration (HGMF) is a promising method for the selective capture and recovery of uranium oxide from surface soils. To date, however, magnetic filtration of uranium oxide has only been demonstrated at a proof-of-principle scale using relatively small filters (<5 cm 3 ) at low flowrates (<60 mL/min). Here, to explore the efficacy of magnetic filtration of uranium oxide at a larger scale, a newly designed HGMF apparatus that is more than an order of magnitude larger than our earlier filters (106 cm 3 ) was designed, fabricated, and tested at relatively high flowrates. Filtration experiments were performed using aqueous uranium oxide particle suspensions with Arizona Road Dust (ARD) as a soil simulant. At a flowrate of 125 mL/min, the apparatus’ uranium capture rate was exceptionally high (96 %), but selectivity was poor due to the high rate of capture for diamagnetic soil constituents (e.g., 77 % for silicon). All particles were captured at a lower rate when the flowrate was increased to 250 mL/min, but uranium selectivity was significantly increased due to the more substantial reduction in diamagnetic particle capture (i.e., capture rate of 77 % and 15 % for uranium and silicon, respectively). When backwashing the apparatus at the same flowrates used during filtration experiments, the rate of uranium recovery tended to be fairly low. Nevertheless, higher flowrates (1 L/min) and sonication were both shown to be highly effective methods of increasing uranium recovery. Magnetic field simulations were also performed to investigate potential optimizations to the design of the apparatus. These simulations showed that the intensity of the applied magnetic field could be increased by increasing the thickness of the steel magnetic housing. Additionally, stochastic trajectory simulations were performed to investigate the potential mechanisms of particle capture.

HGMF↗

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↗

Chimera Events for Performance Studies of the MicroBooNE Deep Learning-based Low Energy Excess Search

MicroBooNE is a short baseline neutrino oscillation experiment based at Fermilab that employs Liquid Argon Time Projection Chamber (LArTPC) technology. One of its target measurements is to investigate the nature of the excess of low energy electron-like events observed by MiniBooNE. This measurement will require an excellent understanding of systematic uncertainties, obtained through testing the performance of reconstruction algorithms on samples with known properties. However, using exclusively Monte Carlo events for this task is limited by how well the discrepancies between simulation and data are understood. An alternative is to test against samples of “chimera” events, which are made up of separate single-particle components from data that are combined to create neutrino-like events. These chimera events can be used to help quantify systematic uncertainties. This note covers the performance and status of creating and using chimera events that match a target neutrino topology in MicroBooNE.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Energy Characterization of the UMD Electron Accelerator

Multiple different groups at Los Alamos National Laboratory have a cited interest in understanding the effects of ionizing radiation in materials. Of especial interest to the ISR division is understanding the effects in materials that are used in the space environment, where radiation is omnipresent and impossible to fully shield from. Development is under way to simulate the effects of radiation damage in detectors and electronics, primarily with particle transport codes such as Geant4 and RAMPART, a code that provides a user friendly interface to Geant4. However, these simulations need benchmarking to experimental tests, some of which have been performed at a variety of small electron accelerator facilities. One of these facilities is the Linear Electron Accelerator at the University of Maryland Radiation Facilities. This facility can provide electron beams on the order of 100 mA delivered in tight macropulse bunches, which allows for experiments where the behavior of electronics is monitored as a function of pulse number. While the beam rates per macropulse of this accelerator is monitored during an experiment via a current pick-off connected to a large graphite beamstop, the energy distribution of this electron beam is not measured during the experiment and has not been well characterized. The energy distribution of the beam is pertinent for these types of experiments, since dose is directly related to energy deposited in a medium. Further, these experiments are performed in open air and sample placement in the room can dictate what dose rate they receive per macropulse. We have created simulated dosemaps for the UMD vault room, but these simulations assume a given energy of the beam, both for calculating the dose and to simulate how much the beam diverges in the air. There have been two attempts now to characterize this accelerator beam energy with a plate type spectrometer, including the most recent experiment earlier this year to quantify the beam energy distribution. These results will also be compared to a dosimetry experiment performed by collaborators at the UMD facility, the details of which are described in Reference. Understanding the electron beam energy is important for future experiments, because the energy directly impacts the dose and radiation damage given to a sample and is important for properly simulating the electron beam when comparing simulation to experiment.

43 PARTICLE ACCELERATORS↗

Evaluation of High Level Waste Sludge Processing Behavior

The U.S. Department of Energy’s (DOE) Hanford Site has 177 underground storage tanks that contain wastes from past nuclear fuel reprocessing and waste-management operations. Over 20% of this waste is in the form of an insoluble sludge that will require slurry modification before its transfer to the Waste Treatment and Immobilization Plant (WTP). Specific WTP acceptance criteria for waste feed delivery describe the physical and chemical characteristics of the waste that must be met before the waste is transferred to the WTP. One challenging requirement relates to the undissolved solids (UDS) composition in a waste feed because the waste contains solid particles that settle, and their concentration and relative proportion can change during the transfer of the waste in individual batches. A key uncertainty is the ability to transfer and mix wastes with large variations in UDS concentrations and resulting settling rates. To address this uncertainty, a number of small scale mixing and settling tests have been conducted to determine the mobilization performance of variable chemistry simulants. Comparison of the size and density of the particulate for each simulant to that of southeast area Hanford sludge was made using metrics for particle mobilization, suspension, settling, and pipeline transfer where dependance on particle size and density may be different, including: 1. Settling velocity, 2. Critical shear stress for erosion, 3. Just-suspended impeller speed, and 4. Pipeline critical transport velocity. Existing high-level waste sludge data has shown the effect that increasing Al concentration has on resulting settled solids. This differential settling of particles in the sludge has the possibility of resulting in solids segregation during feed preparation and uneven particle distribution during pipeline transportation or mixer jet pump operations. Understanding the predictive capabilities of HLW solids settling and transport as well as potential remedies for addressing disparate sludge behaviors can help provide technical guidance during HLW flowsheet planning.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Spent Crystalline Silicotitanate Storage Study

Washington River Protection Solutions is working to support initial production of immobilized low-activity waste (LAW) by feeding Hanford tank supernate from tank farms to the Hanford Waste Treatment and Immobilization Plant (WTP) LAW Facility. This goal incorporates the design of a Tank-Side Cesium Removal (TSCR) system, which filters tank waste supernate to remove suspended solids and then removes Cs by processing it through crystalline silicotitanate (CST) ion exchange media manufactured by Honeywell UOP, LLC. The 137 Cs-depleted product is intended to be sent to the WTP for vitrification. The Cs-loaded CST columns will be stored indefinitely, with a goal of eventual CST removal and treatment. Thus, the spent CST needs to be recoverable. The testing described herein looks to potential upset process conditions where CST storage may be required before various rinse steps are completed. This study evaluated upset conditions at three sequential processing steps envisioned for TSCR (feed, 0.1 M NaOH rinse, water rinse), and in-column drying with compressed air (normal end step) to assess impacts on the nature of the CST bed. Testing was conducted at the small scale (12-mL bed volume); simulated AP-105 tank waste was used as the feed. Following process disruption, the CST bed was dried in place at 70 °C. Post-dried CST bed physical properties (angle of repose, penetration depth, particle morphological changes) were measured to evaluate how CST moved and flowed. The testing is intended to provide a preliminary assessment of issues that may arise from desiccation of CST with the indicated salt solutions in place. Since these were small-scale tests, the processing conditions will not match full scale conditions exactly; however, the tests do provide insight into the impact of stopping processing at an earlier step than normal. Except for the feed that was dried in-place, all other process stop-conditions showed the CST bed flowed well after drying. At this small scale, CST beds would not present an issue for retrievability. The feed that was dried in place had solidified into a rock-hard monolith with no movement possible. Samples had to be chipped from the surface.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Nonlinear dielectric relaxation of polar liquids

Molecular dynamics of two water models, SPC/E and TIP3P, at a number of temperatures are used to test the Kivelson-Madden equation connecting single-particle and collective dielectric relaxation times through the Kirkwood factor. The relation is confirmed by simulations and used to estimate the nonlinear effect of the electric field on the dielectric relaxation time. We show that the main effect of the field comes through slowing down of the single-particle rotational dynamics and the relative contribution of the field-induced alteration of the Kirkwood factor is insignificant for water. Theories of nonlinear dielectric relaxation need to mostly account for the effect of the field on rotations of a single dipole in a polar liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Verification and Validation of a Magnetic Plasma Fluid Model Utilizing the MOOSE (Multiphysics Object Oriented Simulation Environment) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification, benchmarking, and validation studies have been conducted. Verification studies involved utilizing the method of manufactured solutions and comparing the convergence slope of a known solution to the theoretical slope. Benchmarking consists of comparisons to existing edge codes, namely BOUT++ and SOLEDGE3X. Validation efforts focused on comparisons against open-source data from the TCV tokamak.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SR19021 Tritium Aging of Regenerated LANA.75 (Final Project Report)

The Savannah River Tritium Enterprise (SRTE) has used the metal hydride LaNi 4.25 Al 0.75 (LANA.75) in the Tritium Facilities for over two decades. LANA.75 beds store significant quantities of tritium but have a limited service life due to the radiolytic decay of tritium to He-3 within the metal matrix. It has been shown that heating tritium-aged LANA.75 under vacuum can reverse tritium aging effects, eliminating the heel of trapped hydrogen, and restoring the reversible capacity. Additional investigation is needed to ensure there are no unexpected changes to the hydride before this restoration technique is employed in full scale beds in the Tritium Facilities. This project was to be comprised of three distinct scopes: obtain tritium aging data on the sample regenerated in 2018, regenerate a tritium-aged sample at 600 °C, and perform thermal stability testing on a non-tritiated sample. Isotherms were collected on the previously regenerated sample after approximately two years of tritium aging. Isotherms were collected at 80, 100, and 120 °C. As expected, there was a decrease in the plateau pressure, an increase in plateau slope, and a portion of the “heel” of tritium trapped in the metal had been reestablished. Unexpectedly, it appeared that the plateau had shortened at the higher tritium to metal ratios as well. This is typically seen in older samples. The second scope, to regenerate a second tritium-aged LANA.75 sample, was not completed. A Task Technical and Quality Assurance Plan was written and approved, a high temperature test cell was fabricated, several pre-job briefs were held, and the hydride sample was passivated with air. Despite these successes, the hydride sample was not recovered from the legacy test cell. The third scope was to perform thermal stability testing on a “cold” LANA.75 sample. A non-tritiated sample of LANA.75 was held at 750 °C under vacuum for 200 hours to simulate exposure to multiple regeneration evolutions. Hydride isotherm performance, chemical composition, crystallinity, particle size, and morphology are compared between the pre- and post-regeneration samples. No significant changes were observed in composition, crystallinity, or particle size. Comparison of before and after isotherms showed that performance improved rather than deteriorated during the evolution. Scanning Electron Microscopy (SEM) analysis showed small growths on the particle surface after exposure to regeneration conditions. Additional testing will be required to determine the cause of these growths.

07 ISOTOPE AND RADIATION SOURCES↗

Geant4 Simulations and Experimental Results for the Pacific-Ocean Neutrino Experiment (P-ONE)

The Pacific-Ocean Neutrino Experiment (P-ONE) is an upcoming large-volume neutrino telescope planned to be built off the coast of Vancouver Island, Canada. It aims to detect high-energy neutrinos via the Cherenkov radiation produced in water by secondary particles born from high-energy neutrino interactions. Its first set of detectors is planned to be commissioned this year. However, since the telescope is still in its early stages of development, many tests and simulations must still be run to characterize the efficacy of the detectors. One necessary simulation must be devoted towards examining whether the polarization of the incoming Cherenkov light impacts the strength of the signals by the PMTs, and it must be supported by experimental data. In this report, I detail what I have completed during my SULI internship at the Los Alamos National Lab (LANL); I present the preliminary results of a Geant4 polarization simulation as well as preliminary laboratory measurements that will later be used to inform the future P-ONE design.

74 ATOMIC AND MOLECULAR PHYSICS↗

Software Quality Assurance for the MOOSE-Based Open-Source Multiphysics Code Cardinal - An Expanded CI Testing Suite

Cardinal is a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and the Monte Carlo particle transport code OpenMC within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Cardinal provides high-resolution thermal-hydraulics and/or radiation transport feedback to MOOSE multiphysics simulations. Multiphysics feedback is implemented in a geometry-agnostic manner which eliminates the need for rigid one-to-one mappings. A generic data transfer implementation also allows NekRS and OpenMC to couple to any MOOSE application, enabling a broad set of multiphysics capabilities. Cardinal simulations can also leverage combinations of MPI, OpenMP, and GPU resources. Cardinal continuous development and improvement efforts have led to the software being considered as a high-fidelity design and licensing tool for key areas of nuclear reactor relevant physics, including neutron transport, fluid flow, heat transfer, and mechanical processes. The fast development and expansion of the software from a pure R&D framework towards its application in the nuclear industry and regulation require a focus on developing, enhancing and, maintaining Cardinal’s software quality through strict adherence to a Software Quality Assurance (SQA) framework and SQA program. To facilitate compliance with SQA standards, the Cardinal SQA Program has been initiated during Fiscal Year 2023 (FY23). During the development of the Cardinal SQA Program, multiple gaps have been identified. These gaps are primarily related to model verification and code pedigree as they relate to the use of Cardinal as a safety analysis tool. These gaps have been captured in a report published in 2023. A second report highlighted the progress made during Fiscal Year 2024 (FY24) and described Argonne’s effort to document and integrate software verification within Cardinal’s software development process. This report documents a snapshot of the verification test cases currently available for Cardinal and NekRS in their assimilation into a Continuous Integration (CI) platform. Following the CI practice permits the integrating of source code changes frequently and ensuring that the integrated codebase clears the verification testing for the software. It should be noted that the SQA program itself, including the program plans, procedures, configuration management, and testing strategies, need to be developed in a future step of this task.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental and theoretical study of weakly coherent mode in I-mode edge plasmas in the EAST tokamak

The I-mode is a promising operation mode for fusion in the future, featuring high-temperature and low-density confinement, but the reason why the temperature and density are decoupled remains an important aspect to be explored. The experimental results from the experimental advanced superconducting tokamak (EAST) showed that the weakly coherent mode (WCM) is directly related to sustaining the I-mode and that the peak amplitude of the WCM is proportional to the temperature in the pedestal. Here, simulating the experimental data from EAST with the six-field model of BOUT++, we find a density perturbation close to the frequency of the WCM observed in the experiment. By testing all the physical terms in this model, we find that the density perturbation and particle transport are directly related to the drift Alfvén wave (DAW) mode. Additionally we use the SymPIC program (Xiao et al Plasma Sci. Technol. 20 110501; Phys. Plasmas 22 112504; Plasma Sci. Technol. 23 055102) to simulate the same experimental data and find that the frequency range of the WCM is close to both experimental and BOUT++ results. Therefore, the WCM of the I-mode can be considered to be driven by the DAW, which helps improve the transport of the I-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Point containment algorithms for constructive solid geometry with unbounded primitives

Here, we present several algorithms for evaluating point containment in constructive solid geometry (CSG) trees with unbounded primitives. Three algorithms are presented based on postfix, prefix, and infix notations of the CSG binary expression tree. We show that prefix and infix notations enable short-circuiting logic, which reduces the number of primitives that must be checked during point containment. To evaluate the performance of the algorithms, each algorithm was implemented in the OpenMC Monte Carlo particle transport code, which relies on CSG to represent solid bodies through which subatomic particles travel. Two sets of tests were carried out. First, the execution time to generate a rasterized image of a 2D slice of three CSG models of varying complexity was measured. Use of both prefix and infix notations offered significant speedup over the postfix notation that has traditionally been used in particle transport codes, with infix resulting in a 6 x reduction in execution time relative to postfix for a model of a tokamak fusion device. We then measured the execution time of neutron transport simulations of the same three models using each of the algorithms. The results and performance improvements reveal the same trends as for the rasterization test, with a 5.52 x overall speedup using the infix notation relative to the original postfix notation in OpenMC for the tokamak model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Practical Insights on Applying Simulation-Based Control Methods in Experimental Studies

Advanced nuclear reactors are crucial to the future of energy both in the United States and around the globe. In contrast to the current operating fleet, they are characterized as being deployable in remote locations and able to operate in semi-autonomous or autonomous fashion. This leap forward necessitates a new reactor control paradigm. Because advanced nuclear reactors are still under development in the United States, the creation of new control methods to achieve autonomous operations has been based on systems modeling and simulation. However, an important factor in successfully deploying these new control methods is the ability to seamlessly transition from simulation environments to real-world settings. Control methods tested in both simulation and experimental settings need to be investigated in the context of advanced reactor applications. This work developed a series of simple controllers for Idaho National Laboratory (INL)’s Microreactor Applications Research Validation and Evaluation (MARVEL) microreactor operating in load-following scenarios. These controllers were tested in both simulation and experimental settings, and a comparative performance analysis was performed. The simulation tests leveraged the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software developed in a previous stage of the current effort, along with the MARVEL Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The experimental tests leveraged the COMMAND software, MARVEL models, and the U.S. Department of Energy Microreactor Program’s Microreactor Automated Control System (MACS). MACS was developed to serve as a control method testbed. It was customized to mirror the MARVEL microreactor, and COMMAND enabled MACS to emulate the physics of MARVEL. The load-following controller was developed using the simulation platform, with efforts to emulate real systems by introducing actuator saturation and noise. These factors were incrementally accounted for in the controller design. After finalizing the controller design, it was implemented with the experimental setup. The experimental conditions tested included an initial test under conditions similar to the final simulation test, and two additional scenarios. The first scenario introduced additional actuator saturation to account for equipment aging over time, which was unknown to the controller. The second scenario introduced sensor delay, a phenomenon anticipated with the use of remote operations or wireless communication in advanced reactors. These tests revealed several notable differences. While the controller performed well in simulation, it exhibited several limitations when transitioning to hardware. The main challenges involved maintaining the steady-state target power, as evidenced by larger error values between the true reactor power and setpoint power, as well as persistent oscillations in controlled reactor power. These issues could lead to unacceptable transient conditions in real reactor testing. Introducing actuator aging and stochastic delays in the experimental setup significantly impacted controller performance, resulting in increased overshoot and undershoot, and exacerbated error and oscillations previously mentioned. These findings underscore the importance of experimental testbeds for testing and validating control methods, as controllers developed using only theory and/or simulation may perform unexpectedly when applied to actual hardware. This research emphasizes the need for an experimental testbed for achieving such validation.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

Reverse scaling of a bonded-sphere DEM model: Formulation and application to lignocellulosic biomass microstructures

We explore scaling laws for adapting a bonded-sphere discrete element method (BS-DEM) model developed for woody structural mechanics at the millimeter scale to model the mechanics of realistic lignocellulosic biomass microstructures. Two scaling approaches, i.e., the reverse coarse graining (RCG) and equivalent bulk behavior (EBB), are proposed based on the classical mechanics principles and assessed in the single-particle compression and rectangular cuboid block tension tests. The EBB approach is recommended for BS-DEM models with general 3D geometries and is applied to simulate the microindentation test on a realistic microscale pinewood specimen. Simulations are performed to elucidate the impact of specimen thickness and loading position on the specimen’s force–deformation behavior. Furthermore, the range of Young’s modulus obtained from the calibrated BS-DEM simulations can match the experimental measurements. This is the first-of-its-kind study that has explored scale-bridging modeling approaches in biomass micromechanics and has proposed solutions based on BS-DEM models in the microscale.

59 BASIC BIOLOGICAL SCIENCES↗