Search NASA⌕ Search

SEARCH · Search NASA

Results for “SIMULATION METHODS”

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

Evaluating the Use of High-Fidelity Simulator Research Methods to Study Airline Flight Crew Resilience

As it evolves, aviation will continue to require integration of a wide range of safety systems and practices, some of which are already in place and others that are yet to be developed. New concepts in system safety thinking have emerged to consider not only what may go wrong, but also what can be learned when things go right during commercial flight operations. Taken together, these complementary perspectives form a more comprehensive approach to systemsafety thinking that can help to recognize and preserve the resilient performance capabilities currently provided by humans. A need exists, however, for research methods to enable better understanding of the human contributions to aviation safety. NASA’s System-Wide Safety Project supports research on using flight simulation methods to study operator resilience and safety-producing behaviors. Building on prior NASA efforts investigating procedural non-adherences during area navigation standard terminal route arrivals, a high-fidelity commercial aviation line operational simulation (LOS) experiment has been designed to study how flight crews anticipate, monitor for, respond to, and learn from expected and unexpected disturbances during these operations. A diverse set of LOS scenarios were developed to simulate highly realistic, complex, but routinely encountered operational situations. Each scenario provided multiple opportunities to collect data on how flight crews manage threats and errors, as well as novel opportunities to observe resilient and safety-producing behaviors. The experimental design, implications for the study of safety-producing behaviors using simulation, and considerations for airline pilot training will be discussed.

Chad L Stephens↗

SU-D-BRF-06: A Brachytherapy Simulator with Realistic Haptic Force Feedback and Real-Time Ultrasounds Image Simulation for Training and Teaching

Purpose: Surgical procedures require dexterity, expertise and repetition to reach optimal patient outcomes. However, efficient training opportunities are usually limited. This work presents a simulator system with realistic haptic force-feedback and full, real-time ultrasounds image simulation. Methods: The simulator is composed of a custom-made Linear-DELTA force-feedback robotic platform. The needle tip is mounted on a force gauge at the end effector of the robot, which responds to needle insertion by providing reaction forces. 3D geometry of the tissue is using a tetrahedral finite element mesh (FEM) mimicking tissue properties. As the needle is inserted/retracted, tissue deformation is computed using a mass-tensor nonlinear visco-elastic FEM. The real-time deformation is fed to the L-DELTA to take into account the force imparted to the needle, providing feedback to the end-user when crossing tissue boundaries or needle bending. Real-time 2D US image is also generated synchronously showing anatomy, needle insertion and tissue deformation. The simulator is running on an Intel I7 6- core CPU at 3.26 MHz. 3D tissue rendering and ultrasound display are performed on a Windows 7 computer; the FEM computation and L-DELTA control are executed on a similar PC using the Neutrino real-time OS. Both machines communicate through an Ethernet link. Results: The system runs at 500 Hz for a 8333-tetrahedron tissue mesh and a 100-node angular spring needle model. This frame rate ensures a relatively smooth displacement of the needle when pushed or retracted (±20 N in all directions at speeds of up to 2 m/s). Unlike commercially-available haptic platforms, the oblong workspace of the L-DELTA robot complies with that required for brachytherapy needle displacements of 0.1m by 0.1m by 0.25m. Conclusion: We have demonstrated a real-life, realistic brachytherapy simulator developed for prostate implants (LDR/HDR). The platform could be adapted to other sites or training for other types of needle-based procedures.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-T-439: An Improved Formula of Scatter-To-Primary Ratio for Photon Dose Calculation

Purpose: Scatter-to-primary ratio (SPR) is an important dosimetric quantity that describes the contribution from the scatter photons in an external photon beam. The purpose of this study is to develop an improved analytical formula to describe SPR as a function of circular field size (r) and depth (d) using Monte Carlo (MC) simulation. Methods: MC simulation was performed for Mohan photon spectra (Co-60, 4, 6, 10, 15, 23 MV) using EGSNRC code. Point-spread scatter dose kernels in water are generated. The scatter-to-primary ratio (SPR) is also calculated using MC simulation as a function of field size for circular field size with radius r and depth d. The doses from forward scatter and backscatter photons are calculated using a convolution of the point-spread scatter dose kernel and by accounting for scatter photons contributing to dose before (z' d) reaching the depth of interest, d, where z' is the location of scatter photons, respectively. The depth dependence of the ratio of the forward scatter and backscatter doses is determined as a function of depth and field size. Results: We are able to improve the existing 3-parameter (a, w, d0) empirical formula for SPR by introducing depth dependence for one of the parameter d0, which becomes 0 for deeper depths. The depth dependence of d0 can be directly calculated as a ratio of backscatter-to-forward scatter doses for otherwise the same field and depth. With the improved empirical formula, we can fit SPR for all megavoltage photon beams to within 2%. Existing 3-parameter formula cannot fit SPR data for Co-60 to better than 3.1%. Conclusion: An improved empirical formula is developed to fit SPR for all megavoltage photon energies to within 2%.

60 APPLIED LIFE SCIENCES↗

Method for simulating dose reduction in digital mammography using the Anscombe transformation

Purpose: This work proposes an accurate method for simulating dose reduction in digital mammography starting from a clinical image acquired with a standard dose. Methods: The method developed in this work consists of scaling a mammogram acquired at the standard radiation dose and adding signal-dependent noise. The algorithm accounts for specific issues relevant in digital mammography images, such as anisotropic noise, spatial variations in pixel gain, and the effect of dose reduction on the detective quantum efficiency. The scaling process takes into account the linearity of the system and the offset of the detector elements. The inserted noise is obtained by acquiring images of a flat-field phantom at the standard radiation dose and at the simulated dose. Using the Anscombe transformation, a relationship is created between the calculated noise mask and the scaled image, resulting in a clinical mammogram with the same noise and gray level characteristics as an image acquired at the lower-radiation dose. Results: The performance of the proposed algorithm was validated using real images acquired with an anthropomorphic breast phantom at four different doses, with five exposures for each dose and 256 nonoverlapping ROIs extracted from each image and with uniform images. The authors simulated lower-dose images and compared these with the real images. The authors evaluated the similarity between the normalized noise power spectrum (NNPS) and power spectrum (PS) of simulated images and real images acquired with the same dose. The maximum relative error was less than 2.5% for every ROI. The added noise was also evaluated by measuring the local variance in the real and simulated images. The relative average error for the local variance was smaller than 1%. Conclusions: A new method is proposed for simulating dose reduction in clinical mammograms. In this method, the dependency between image noise and image signal is addressed using a novel application of the Anscombe transformation. NNPS, PS, and local noise metrics confirm that this method is capable of precisely simulating various dose reductions.

60 APPLIED LIFE SCIENCES↗

Direct simulation of electron transfer using ring polymer molecular dynamics: Comparison with semiclassical instanton theory and exact quantum methods

The use of ring polymer molecular dynamics (RPMD) for the direct simulation of electron transfer (ET) reaction dynamics is analyzed in the context of Marcus theory, semiclassical instanton theory, and exact quantum dynamics approaches. For both fully atomistic and system-bath representations of condensed-phase ET, we demonstrate that RPMD accurately predicts both ET reaction rates and mechanisms throughout the normal and activationless regimes of the thermodynamic driving force. Analysis of the ensemble of reactive RPMD trajectories reveals the solvent reorganization mechanism for ET that is anticipated in the Marcus rate theory, and the accuracy of the RPMD rate calculation is understood in terms of its exact description of statistical fluctuations and its formal connection to semiclassical instanton theory for deep-tunneling processes. In the inverted regime of the thermodynamic driving force, neither RPMD nor a related formulation of semiclassical instanton theory capture the characteristic turnover in the reaction rate; comparison with exact quantum dynamics simulations reveals that these methods provide inadequate quantization of the real-time electronic-state dynamics in the inverted regime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unbiased simulation of near-Clifford quantum circuits

Modeling and simulation are essential for predicting and verifying the behavior of fabricated quantum circuits, but existing simulation methods are either impractically costly or require an unrealistic simplification of error processes. In this paper, we present a method of simulating noisy Clifford circuits that is both accurate and practical in experimentally relevant regimes. In particular, the cost is weakly exponential in the size and the degree of non-Cliffordness of the circuit. Our approach is based on the construction of exact representations of quantum channels as quasiprobability distributions over stabilizer operations, which are then sampled, simulated, and weighted to yield unbiased statistical estimates of circuit outputs and other observables. As a demonstration of these techniques, we simulate a Steane [[7,1,3]]-encoded logical operation with non-Clifford errors and compute its fault tolerance error threshold. Finally, we expect that the method presented here will enable studies of much larger and more realistic quantum circuits than was previously possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-D-19A-06: The Effect of Beam Parameters On Very High-Energy Electron Radiotherapy: A Planning Study

Purpose: We evaluated the effect of very high-energy electron (VHEE) beam parameters on the planning of a lung cancer case by means of Monte Carlo simulations. Methods: We simulated VHEE radiotherapy plans using the EGSnrc/BEAMnrc-DOSXYZnrc code. We selected a lung cancer case that was treated with 6MV photon VMAT to be planned with VHEE. We studied the effect of beam energy (80 MeV, 100 MeV, and 120 MeV), number of equidistant beams (16 or 32), and beamlets sizes (3 mm, 5 mm or 7 mm) on PTV coverage, sparing of organs at risk (OARs) and dose conformity. Inverse-planning optimization was performed in a research version of RayStation (RaySearch Laboratories AB) using identical objective functions and constraints for all VHEE plans. Results: Similar PTV coverage and dose conformity was achieved by all the VHEE plans. The 100 MeV and 120 MeV VHEE plans were equivalent amongst them and were superior to the 80 MeV plan in terms of OARs sparing. The effect of using 16 or 32 equidistant beams was a mean difference in average dose of 2.4% (0%–7.7%) between the two plans. The use of 3 mm beamlet size systematically reduced the dose to all the OARs. Based on these results we selected the 100MeV-16beams-3mm-beamlet-size plan to compare it against VMAT. The selected VHEE plan was more conformal than VMAT and improved OAR sparing (heart and trachea received 125% and 177% lower dose, respectively) especially in the low-dose region. Conclusion: We determined the VHEE beam parameters that maximized the OAR dose sparing and dose conformity of the actually delivered VMAT plan of a lung cancer case. The selected parameters could be used for the planning of other treatment sites with similar size, shape, and location. For larger targets, a larger beamlet size might be used without significantly increasing the dose. B Palma: None. M Bazalova: None. B Hardemark: Employee, RaySearch Americas. E Hynning: Employee, RaySearch Americas. B Qu: None. B Loo Jr.: Research support, RaySearch, Varian. P Maxim: Research support, RaySearch, Varian.

60 APPLIED LIFE SCIENCES↗

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗

Microstructure of vapor deposited coatings on curved substrates

Thermal barrier coating systems consisting of a metallic bond coat and ceramic over layer are widely used to extend the life of gas turbine engine components. They are applied using either high-vacuum physical vapor deposition techniques in which vapor atoms rarely experience scattering collisions during propagation to a substrate, or by gas jet assisted (low-vacuum) vapor deposition techniques that utilize scattering from streamlines to enable non-line-of-sight deposition. Both approaches require substrate motion to coat a substrate of complex shape. Here, direct simulation Monte Carlo and kinetic Monte Carlo simulation methods are combined to simulate the deposition of a nickel coating over the concave and convex surfaces of a model airfoil, and the simulation results are compared with those from experimental depositions. The simulation method successfully predicted variations in coating thickness, columnar growth angle, and porosity during both stationary and substrate rotated deposition. It was then used to investigate a wide range of vapor deposition conditions spanning high-vacuum physical vapor deposition to low-vacuum gas jet assisted vapor deposition. The average coating thickness was found to increase initially with gas pressure reaching a maximum at a chamber pressure of 8–10 Pa, but the best coating thickness uniformity was achieved under high vacuum deposition conditions. However, high vacuum conditions increased the variation in the coatings pore volume fraction over the surface of the airfoil. The simulation approach was combined with an optimization algorithm and used to investigate novel deposition concepts to tailor the local coating thickness.

42 ENGINEERING↗

SU-F-T-62: Three-Dimensional Dosimetric Gamma Analysis for Impacts of Tissue Inhomogeneity Using Monte Carlo Simulation in Intracavitary Brachytheray for Cervix Carcinoma

Purpose: The aim of this study was to investigate the impacts of tissue inhomogeneity on dose distributions using a three-dimensional (3D) gamma analysis in cervical intracavitary brachytherapy using Monte Carlo (MC) simulations. Methods: MC simulations for comparison of dose calculations were performed in a water phantom and a series of CT images of a cervical cancer patient (stage: Ib; age: 27) by employing a MC code, Particle and Heavy Ion Transport Code System (PHIT) version 2.73. The {sup 192}Ir source was set at fifteen dwell positions, according to clinical practice, in an applicator consisting of a tandem and two ovoids. Dosimetric comparisons were performed for the dose distributions in the water phantom and CT images by using gamma index image and gamma pass rate (%). The gamma index is the minimum Euclidean distance between two 3D spatial dose distributions of the water phantom and CT images in a same space. The gamma pass rates (%) indicate the percentage of agreement points, which mean that two dose distributions are similar, within an acceptance criteria (3 mm/3%). The volumes of physical and clinical interests for the gamma analysis were a whole calculated volume and a region larger than t% of a dose (close to a target), respectively. Results: The gamma pass rates were 77.1% for a whole calculated volume and 92.1% for a region within 1% dose region. The differences of 7.7% to 22.9 % between two dose distributions in the water phantom and CT images were found around the applicator region and near the target. Conclusion: This work revealed the large difference on the dose distributions near the target in the presence of the tissue inhomogeneity. Therefore, the tissue inhomogeneity should be corrected in the dose calculation for clinical treatment.

60 APPLIED LIFE SCIENCES↗

The DESI N -body simulation project – I. Testing the robustness of simulations for the DESI dark time survey

ABSTRACT Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogues to validate data analysis pipelines and identify potential systematics. We compare three N-body simulation codes, abacus, gadget-2, and swift, to investigate the regimes in which their results agree. We run N-body simulations at three different mass resolutions, 6.25 × 108, 2.11 × 109, and 5.00 × 109 h−1 M⊙, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the Dark Energy Spectroscopic Instrument (DESI) statistical error for $s\ \gt\ 20\ h^{-1}$ Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of 2.1 × 109 h−1 M⊙ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than $20\ h^{-1}$ Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1 per cent for k ≤ 10 h Mpc−1. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi-N-body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy–halo relation will be presented in a follow-up study.

79 ASTRONOMY AND ASTROPHYSICS↗

Pore-Scale and Multiscale Numerical Simulation of Flow and Transport in a Laboratory-Scale Column

Pore-scale models are useful for studying relationships between fundamental processes and phenomena at larger (i.e., Darcy) scales. However, the size of domains that can be simulated with explicit pore-scale resolution is limited by computational and observational constraints. Direct numerical simulation of pore-scale flow and transport is typically performed on millimeter-scale volumes at which X-ray computed tomography (XCT), often used to characterize pore geometry, can achieve micrometer resolution. In contrast, the scale at which a continuum approximation of a porous medium is valid is usually larger, on the order of centimeters to decimeters. Furthermore, laboratory experiments that measure continuum properties are typically performed on decimeter-scale columns. At this scale, XCT resolution is coarse (tens to hundreds of micrometers) and prohibits characterization of small pores and grains. We performed simulations of pore-scale processes over a decimeter-scale volume of natural porous media with a wide range of grain sizes, and compared to results of column experiments using the same sample. Simulations were conducted using high-performance codes executed on a supercomputer. Two approaches to XCT image segmentation were evaluated, a binary (pores and solids) segmentation and a ternary segmentation that resolved a third category (porous solids with pores smaller than the imaged resolution). We used a mixed Stokes-Darcy simulation method to simulate the combination of Stokes flow in large open pores and Darcy-like flow in porous solid regions. Simulations based on the ternary segmentation provided results that were consistent with experimental observations, demonstrating our ability to successfully model pore-scale flow over a column-scale domain.

CFD↗

Investigating the atomic level influencing factors of glass forming ability in NiAl and CuZr metallic glasses

Bulk metallic glasses are a relatively new class of amorphous metal alloy which possess unique mechanical and magnetic properties. The specific concentrations and combinations of alloy elements needed to prevent crystallization during melt quenching remains poorly understood. A correlation between atomic properties that can explain some of the previously identified glass forming ability (GFA) anomalies of the NiAl and CuZr systems has been identified, with these findings likely extensible to other transition metal–transition metal and transition metal–metalloid (TM–M) alloy classes as a whole. In this work, molecular dynamics simulation methods are utilized to study thermodynamic, kinetic, and structural properties of equiatomic CuZr and NiAl metallic glasses in an attempt to further understand the underlying connections between glass forming ability, nature of atomic level bonding, short and medium range ordering, and the evolution of structure and relaxation properties in the disordered phase. The anomalous breakdown of the fragility parameter as a useful GFA indicator in TM–M alloy systems is addressed through an in-depth investigation of bulk stiffness properties and the evolution of (pseudo)Gruneisen parameters over the quench domain, with the efficacy of other common glass forming ability indicators similarly being analyzed through direct computation in respective CuZr and NiAl systems. Comparison of fractional liquid-crystal density differences in the two systems revealed 2-3 times higher values for the NiAl system, providing further support for its efficacy as a general purpose GFA indicator.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THE COSMIC-RAY ENERGY SPECTRUM OBSERVED WITH THE SURFACE DETECTOR OF THE TELESCOPE ARRAY EXPERIMENT

The Telescope Array (TA) collaboration has measured the energy spectrum of ultra-high energy cosmic rays (UHECRs) with primary energies above 1.6 Multiplication-Sign 10{sup 18} eV. This measurement is based upon four years of observation by the surface detector component of TA. The spectrum shows a dip at an energy of 4.6 Multiplication-Sign 10{sup 18} eV and a steepening at 5.4 Multiplication-Sign 10{sup 19} eV which is consistent with the expectation from the GZK cutoff. We present the results of a technique, new to the analysis of UHECR surface detector data, that involves generating a complete simulation of UHECRs striking the TA surface detector. The procedure starts with shower simulations using the CORSIKA Monte Carlo program where we have solved the problems caused by use of the ''thinning'' approximation. This simulation method allows us to make an accurate calculation of the acceptance of the detector for the energies concerned.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-CAMPUS-I-05: Internal Dosimetric Calculations for Several Imaging Radiopharmaceuticals in Preclinical Studies and Quantitative Assessment of the Mouse Size Impact On Them. Realistic Monte Carlo Simulations Based On the 4D-MOBY Model

Purpose: Commonly used radiopharmaceuticals were tested to define the most important dosimetric factors in preclinical studies. Dosimetric calculations were applied in two different whole-body mouse models, with varying organ size, so as to determine their impact on absorbed doses and S-values. Organ mass influence was evaluated with computational models and Monte Carlo(MC) simulations. Methods: MC simulations were executed on GATE to determine dose distribution in the 4D digital MOBY mouse phantom. Two mouse models, 28 and 34 g respectively, were constructed based on realistic preclinical exams to calculate the absorbed doses and S-values of five commonly used radionuclides in SPECT/PET studies (18F, 68Ga, 177Lu, 111In and 99mTc).Radionuclide biodistributions were obtained from literature. Realistic statistics (uncertainty lower than 4.5%) were acquired using the standard physical model in Geant4. Comparisons of the dosimetric calculations on the two different phantoms for each radiopharmaceutical are presented. Results: Dose per organ in mGy was calculated for all radiopharmaceuticals. The two models introduced a difference of 0.69% in their brain masses, while the largest differences were observed in the marrow 18.98% and in the thyroid 18.65% masses.Furthermore, S-values of the most important target-organs were calculated for each isotope. Source-organ was selected to be the whole mouse body.Differences on the S-factors were observed in the 6.0–30.0% range. Tables with all the calculations as reference dosimetric data were developed. Conclusion: Accurate dose per organ and the most appropriate S-values are derived for specific preclinical studies. The impact of the mouse model size is rather high (up to 30% for a 17.65% difference in the total mass), and thus accurate definition of the organ mass is a crucial parameter for self-absorbed S values calculation.Our goal is to extent the study for accurate estimations in small animal imaging, whereas it is known that there is a large variety in the anatomy of the organs.

60 APPLIED LIFE SCIENCES↗

SU-E-J-217: Multiparametric MR Imaging of Cranial Tumors On a Dedicated 1.0T MR Simulator Prior to Stereotactic Radiosurgery

Purpose: Quantitative magnetic resonance imaging (MRI) of cranial lesions prior to stereotactic radiosurgery (SRS) may improve treatment planning and provide potential prognostic value. The practicality and logistics of acquiring advanced multiparametric MRI sequences to measure vascular and cellular properties of cerebral tumors are explored on a 1.0 Tesla MR Simulator. Methods: MR simulation was performed immediately following routine CT simulation on a 1T MR Simulator. MR sequences used were in the order they were performed: T2-Weighted Turbo Spin Echo (T2W-TSE), T2 FLAIR, Diffusion-weighted (DWI, b = 0, 800 to generate an apparent diffusion coefficient (ADC) map), 3D T1-Weighted Fast Field Echo (T1W-FFE), Dynamic Contrast Enhanced (DCE) and Post Gadolinium Contrast Enhanced 3D T1W-FFE images. T1 pre-contrast values was generated by acquiring six different flip angles. The arterial input function was derived from arterial pixels in the perfusion images selected manually. The extended Tofts model was used to generate the permeability maps. Routine MRI scans took about 30 minutes to complete; the additional scans added 12 minutes. Results: To date, seven patients with cerebral tumors have been imaged and tumor physiology characterized. For example, on a glioblastoma patient, the volume contoured on T1 Gd images, ADC map and the pharmacokinetic map (Ktrans) were 1.9, 1.4, and 1.5 cc respectively with strong spatial correlation. The mean ADC value of the entire volume was 1141 μm2/s while the value in the white matter was 811 μm2/s. The mean value of Ktrans was 0.02 min-1 in the tumor volume and 0.00 in the normal white matter. Conclusion: Our initial results suggest that multiparametric MRI sequences may provide a more quantitative evaluation of vascular and tumor properties. Implementing functional imaging during MR-SIM may be particularly beneficial in assessing tumor extent, differentiating radiation necrosis from tumor recurrence, and establishing reliable bio-markers for treatment response evaluation. The Department of Radiation Oncology at Henry Ford Health System has research agreement with Varian Medical System and Philips Health Care.

60 APPLIED LIFE SCIENCES↗

Pebble Fuel Handling and Reactivity Control for Salt-Cooled High Temperature Reactors

This report documents the work completed on the X-PREX facility under NEUP Project 11- 3172. This project seeks to demonstrate the viability of pebble fuel handling and reactivity control for fluoride salt-cooled high-temperature reactors (FHRs). The research results also improve the understanding of pebble motion in helium-cooled reactors, as well as the general, fundamental understanding of low-velocity granular flows. Successful use of pebble fuels in with salt coolants would bring major benefits for high-temperature reactor technology. Pebble fuels enable on-line refueling and operation with low excess reactivity, and thus simpler reactivity control and improved fuel utilization. If fixed fuel designs are used, the power density of salt- cooled reactors is limited to 10 MW/m 3 to obtain adequate duration between refueling, but pebble fuels allow power densities in the range of 20 to 30 MW/m 3 . This can be compared to the typical modular helium reactor power density of 5 MW/m3. Pebble fuels also permit radial zoning in annular cores and use of thorium or graphite pebble blankets to reduce neutron fluences to outer radial reflectors and increase total power production. Combined with high power conversion efficiency, compact low-pressure primary and containment systems, and unique safety characteristics including very large thermal margins (>500°C) to fuel damage during transients and accidents, salt-cooled pebble fuel cores offer the potential to meet the major goals of the Advanced Reactor Concepts Development program to provide electricity at lower cost than light water reactors with improved safety and system performance.This report presents the facility description, experimental results, and supporting simulation methods of the new X-Ray Pebble Recirculation Experiment (X-PREX), which is now operational and being used to collect data on the behavior of slow dense granular flows relevant to pebble bed reactor core designs. The X-PREX facility uses novel digital x-ray tomography methods to track both the translational and rotational motion of spherical pebbles, which provides unique experimental results that can be used to validate discrete element method (DEM) simulations of pebble motion. The validation effort supported by the X-PREX facility provides a means to build confidence in analysis of pebble bed configuration and residence time distributions that impact the neutronics, thermal hydraulics, and safety analysis of pebble bed reactor cores. Experimental and DEM simulation results are reported for silo drainage, a classical problem in the granular flow literature, at several hopper angles. These studies include conventional converging and novel diverging geometries that provide additional flexibility in the design of pebble bed reactor cores. Excellent agreement is found between the X-PREX experimental and DEM simulation results. This report also includes results for additional studies relevant to the design and analysis of pebble bed reactor cores including the study of forces on shut down blades inserted directly into a packed bed and pebble flow in a cylindrical hopper that is representative of a small test reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Monte-Carlo modeling of polar wind photoelectron distributions with anomalous heat flux

In situ measurements above the polar cap by the DE satellites show asymmetric field-aligned electron velocity distributions in the photoelectron energy range 5-60 eV. The formation of these anisotropic distributions is described by a collisional kinetic equation that includes the effects of Coulomb collisions (via a Fokker-Planck operator), those of a parallel electric field, and the magnetic mirror effects of the geomagnetic field. This equation is solved using a Monte Carlo simulation method. The simulation illustrates the formation of the portion of the observed electron distributions originating in the ionosphere and makes it possible to evaluate the significance of the suprathermal electrons at higher altitude. Both observations and calculations indicate that these electrons possess an energy flux which, by its effect on the polar-wind electric field, may influence the dynamics of the polar-wind outflow.

Yasseen, F.↗