SEARCH · Search NASA

Results for “Image coding”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,387 records · Page 77

All sources NASA NTRS DOE OSTI

DOE OSTI2015

SU-E-J-249: Characterization of Gynecological Tumor Heterogeneity Using Texture Analysis in the Context of An 18F-FDG PET Adaptive Protocol

Purpose: We propose a method to examine gynecological tumor heterogeneity using texture analysis in the context of an adaptive PET protocol in order to establish if texture metrics from baseline PET-CT predict tumor response better than SUV metrics alone as well as determine texture features correlating with tumor response during radiation therapy. Methods: This IRB approved protocol included 29 women with node positive gynecological cancers visible on FDG-PET treated with EBRT to the PET positive nodes. A baseline and intra-treatment PET-CT was obtained. Tumor outcome was determined based on RECIST on posttreatment PET-CT. Primary GTVs were segmented using 40% threshold and a semi-automatic gradient-based contouring tool, PET Edge (MIM Software Inc., Cleveland, OH). SUV histogram features, Metabolic Volume (MV), and Total Lesion Glycolysis (TLG) were calculated. Four 3D texture matrices describing local and regional relationships between voxel intensities in the GTV were generated: co-occurrence, run length, size zone, and neighborhood difference. From these, 39 texture features were calculated. Prognostic power of baseline features derived from gradientbased and threshold GTVs were determined using the Wilcoxon rank-sum test. Receiver Operating Characteristics and logistic regression was performed using JMP (SAS Institute Inc., Cary, NC) to find probabilities of predicting response. Changes in features during treatment were determined using the Wilcoxon signed-rank test. Results: Of the 29 patients, there were 16 complete responders, 7 partial responders, and 6 non-responders. Comparing CR/PR vs. NR for gradient-based GTVs, 7 texture values, TLG, and SUV kurtosis had a p < 0.05. Threshold GTVs yielded 4 texture features and TLG with p < 0.05. From baseline to intra-treatment, 14 texture features, SUVmean, SUVmax, MV, and TLG changed with p < 0.05. Conclusion: Texture analysis of PET imaged gynecological tumors is an effective method for early prognosis and should be used complimentary to SUV metrics, especially when using gradient based segmentation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

MO-AB-206-02: Testing Gamma Cameras Based On TG177 WG Report

This education session will cover the physics and operation principles of gamma cameras and PET scanners. The first talk will focus on PET imaging. An overview of the principles of PET imaging will be provided, including positron decay physics, and the transition from 2D to 3D imaging. More recent advances in hardware and software will be discussed, such as time-of-flight imaging, and improvements in reconstruction algorithms that provide for options such as depth-of-interaction corrections. Quantitative applications of PET will be discussed, as well as the requirements for doing accurate quantitation. Relevant performance tests will also be described. Learning Objectives: Be able to describe basic physics principles of PET and operation of PET scanners. Learn about recent advances in PET scanner hardware technology. Be able to describe advances in reconstruction techniques and improvements Be able to list relevant performance tests. The second talk will focus on gamma cameras. The Nuclear Medicine subcommittee has charged a task group (TG177) to develop a report on the current state of physics testing of gamma cameras, SPECT, and SPECT/CT systems. The report makes recommendations for performance tests to be done for routine quality assurance, annual physics testing, and acceptance tests, and identifies those needed satisfy the ACR accreditation program and The Joint Commission imaging standards. The report is also intended to be used as a manual with detailed instructions on how to perform tests under widely varying conditions. Learning Objectives: At the end of the presentation members of the audience will: Be familiar with the tests recommended for routine quality assurance, annual physics testing, and acceptance tests of gamma cameras for planar imaging. Be familiar with the tests recommended for routine quality assurance, annual physics testing, and acceptance tests of SPECT systems. Be familiar with the tests of a SPECT/CT system that include the CT images for SPECT reconstructions. Become knowledgeable of items to be included in annual acceptance testing reports including CT dosimetry and PACS monitor measurements. T. Turkington, GE Healthcare.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2022

Microscale THMC Modeling of Pressure Solution in Salt Rock: Impacts of Geometry and Temperature

Abstract Pressure solution, a mechanism that involves tight coupling between the geometry and thermal-hydro-mechanical-chemical (THMC) processes, plays an important role in diagenesis. In this study, we make the first attempt to conduct microscale THMC modeling to understand and quantify the impacts of geometry and temperature on pressure solution, taking natural salt rock as an example. This modeling capability is achieved by expanding a novel MC code that we developed previously (Hu et al. J Geophys Res: Solid Earth 126:e2021JB023112, 2021) to include temperature effects. We first conduct a simulation of an example that involves a single brine inclusion within a single halite grain and find that the temperature impact is limited for that case. We then extract geometry from an image of a natural salt rock and conduct simulations with different cases: (A) only temperature and no stress, (B) only stress and no temperature, and (C) with both stress and temperature. These different cases result in quite different phenomena. In case A, dissolution and precipitation occur across the entire system due to isolated pore space reaching a localized mass balance between dissolution, precipitation, and diffusion. In case B, intense geometric features (e.g., major asperities, inclusions) in one area undergo stress concentration, thus dominating pressure solution in that area. In case C, pressure solution is spread out at contacting highly stressed geometric features close to the hotter side. We conclude that geometric features dominate stress distribution, thus dominating pressure solution in a natural salt rock that may be affected by the temperature if a sufficient temperature gradient is applied.

58 GEOSCIENCES↗

DOE OSTI2012

Creation and diagnosis of a solid-density plasma with an X-ray free-electron laser

Matter with a high energy density (>10 5 joules per cm 3 ) is prevalent throughout the Universe, being present in all types of stars and towards the centre of the giant planets; it is also relevant for inertial confinement fusion. Its thermodynamic and transport properties are challenging to measure, requiring the creation of sufficiently long-lived samples at homogeneous temperatures and densities. With the advent of the Linac Coherent Light Source (LCLS) X-ray laser, high-intensity radiation (>10 17 watts per cm 2 , previously the domain of optical lasers) can be produced at X-ray wavelengths. The interaction of single atoms with such intense X-rays has recently been investigated. An understanding of the contrasting case of intense X-ray interaction with dense systems is important from a fundamental viewpoint and for applications. Here we report the experimental creation of a solid-density plasma at temperatures in excess of 10 6 kelvin on inertial-confinement timescales using an X-ray free-electron laser. We discuss the pertinent physics of the intense X-ray–matter interactions, and illustrate the importance of electron–ion collisions. Detailed simulations of the interaction process conducted with a radiative-collisional code show good qualitative agreement with the experimental results. We obtain insights into the evolution of the charge state distribution of the system, the electron density and temperature, and the timescales of collisional processes. Furthermore, our results should inform future high-intensity X-ray experiments involving dense samples, such as X-ray diffractive imaging of biological systems, material science investigations, and the study of matter in extreme conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

SU-E-T-58: Calculation of Dose Distribution of Accuboost Brachytherapy in Deformable Polyvinil Alcohol Breast Phantom Using Biomechanical Modeling and Monte Carlo Simulation

Purpose: the accuboost is the most modern method of breast brachytherapy that is a boost method in compressed tissue by a mammography unit. the dose distribution in uncompressed tissue, as compressed tissue is important that should be characterized. Methods: In this study, the mechanical behavior of breast in mammography loading, the displacement of breast tissue and the dose distribution in compressed and uncompressed tissue, are investigated. Dosimetry was performed by two dosimeter methods of Monte Carlo simulations using MCNP5 code and thermoluminescence dosimeters. For Monte Carlo simulations, the dose values in cubical lattice were calculated using tally F6. The displacement of the breast elements was simulated by Finite element model and calculated using ABAQUS software, from which the 3D dose distribution in uncompressed tissue was determined. The geometry of the model is constructed from MR images of 6 volunteers. Experimental dosimetery was performed by placing the thermoluminescence dosimeters into the polyvinyl alcohol breast equivalent phantom and on the proximal edge of compression plates to the chest. Results: The results indicate that using the cone applicators would deliver more than 95% of dose to the depth of 5 to 17mm, while round applicator will increase the skin dose. Nodal displacement, in presence of gravity and 60N forces, i.e. in mammography compression, was determined with 43% contraction in the loading direction and 37% expansion in orthogonal orientation. Finally, in comparison of the acquired from thermoluminescence dosimeters with MCNP5, they are consistent with each other in breast phantom and in chest's skin with average different percentage of 13.7±5.7 and 7.7±2.3, respectively. Conclusion: The major advantage of this kind of dosimetry is the ability of 3D dose calculation by FE Modeling. Finally, polyvinyl alcohol is a reliable material as a breast tissue equivalent dosimetric phantom that provides the ability of TLD dosimetry for validation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Plasma heating in solar flares and their soft and hard X-ray emissions

In this paper, the energy budgets of two single-loop-like flares observed in X-ray are analyzed under the assumption that nonthermal electrons (NTEs) are the only source of plasma heating during all phases of both events. The flares were observed by RHESSI and GOES on 2002 February 20 and June 2, respectively. Using a one-dimensional (1D) hydrodynamic code for both flares, the energy deposited in the chromosphere was derived applying RHESSI observational data. The use of the Fokker-Planck formalism permits the calculation of distributions of the NTEs in flaring loops and thus spatial distributions of the X-ray nonthermal emissions and integral fluxes for the selected energy ranges that were compared with the observed ones. Additionally, a comparative analysis of the spatial distributions of the signals in the RHESSI images was conducted for the footpoints and for all the flare loops in selected energy ranges with these quantities' fluxes obtained from the models. The best compatibility of the model and observations was obtained for the 2002 June 2 event in the 0.5-4 Å GOES range and total fluxes in the 6-12 keV, 12-25 keV, 20-25 keV, and 50-100 keV energy bands. Results of photometry of the individual flaring structures in a high energy range show that the best compliance occurred for the 2002 June 2 flare, where the synthesized emissions were at least 30% higher than the observed emissions. For the 2002 February 20 flare, synthesized emission is about four times lower than the observed one. However, in the low energy range the best conformity was obtained for the 2002 February 20 flare, where emission from the model is about 11% lower than the observed one. The larger inconsistency occurs for the 2002 June 2 solar flare, where synthesized emission is about 12 times greater or even more than the observed emission. Some part of these differences may be caused by inevitable flaws of the applied methodology, like by an assumption that the model of the flare is symmetric and there are no differences in the emissions originating from the feet of the flares loop and by relative simplicity of the applied numerical 1D code and procedures. No doubt a significant refinement of the applied numerical models and more sophisticated implementation of the various physical mechanisms involved are required to achieve a better agreement. Despite these problems, a collation of modeled results with observations shows that soft and hard X-ray emissions observed for analyzed single-loop-like events may be fully explained by electron-beam-driven evaporation only.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

SIMULATION OF ASTRONOMICAL IMAGES FROM OPTICAL SURVEY TELESCOPES USING A COMPREHENSIVE PHOTON MONTE CARLO APPROACH

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons s{sup −1}, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to (1) construct kilometer scale phase screens necessary for wide-field telescopes, (2) reproduce atmospheric point-spread function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, (3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and (4) recover system effective area predicted from analytic photometry integrals. This new code, the Photon Simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

SU-E-CAMPUS-I-03: Dosimetric Comparison of the Hypoxia Agent Iodoazomycin Arabinoside (IAZA) Labeled with the Radioisotopes I-123, I-131 and I-124

Purpose: To compare the radiation dose to normal organs from the radio-iodinated, hypoxia-binding radiosensitizer iodoazomycin arabinoside (IAZA) for three different isotopes of iodine. Methods: Dosimety studies with normal volunteers had been carried out with [{sup 123}I]IAZA, a drug binding selectively to hypoxic sites. Two other isotopes of iodine, {sup 131}I and {sup 124}I, offer the opportunity to use IAZA as an agent for radioisotope therapy and as an imaging tracer for Positron Emission Tomography. Radioisotope dosimetry for {sup 131}I and {sup 124}I was performed by first deriving from the [{sup 123}I]IAZA studies biological uptake and excretion data. The cumulated activities for {sup 131}I or {sup 124}I where obtained by including their half-lives when integrating the biological data and then extrapolating to infinite time points considering a) physical decay only or b) physical and biological excretion. Doses were calculated using the Medical Internal Radiation Dose (MIRD) schema (OLINDA1.1 code, Vanderbilt 2007). Results: Compared to {sup 123}I, organ doses were elevated on average by a factor 6 and 9 for {sup 131}I and {sup 124}I, respectively, if both physical decay and biological excretion were modeled. If only physical decay is considered, doses increase by a factor 18 ({sup 131}I) and 19 ({sup 124}I). Highest organ doses were observed in intestinal walls, urinary bladder and thyroid. Effective doses increased by a factor 11 and 14 for {sup 131}I and {sup 124}I, respectively, if biological and physical decay are present. Purely physical decay yields a 23-fold increase over {sup 123}I for both, {sup 131}I and {sup 124}I. Conclusion: Owing to the significant dose increase, caused by their longer half life and the approximately 10 times larger electronic dose deposited in tissue per nuclear decay, normal tissue doses of IAZA labeled with {sup 131}I and {sup 124}I need to be carefully considered when designing imaging and therapy protocols for clinical trials. Effective blocking of iodine uptake in the thyroid is essential. Alberta Innovates - Health Solutions (AIHS) and Canadian Institutes of Health Research (CIHR)

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Modeling the effect of laser heating on the strength and failure of 7075-T6 aluminum

The effect of rapid laser heating on the response of 7075-T6 aluminum has been characterized using 3-D digital image correlation and a series of thermocouples. The experimental results indicate that as the samples are held under a constant load, the heating from the laser profile causes non-uniform temperature and strain fields, and the strain-rate increases dramatically as the sample nears failure. Simulations have been conducted using the LLNL multi-physics code ALE3D, and compared to the experiments. The strength and failure of the material was modeled using the Johnson–Cook strength and damage models. Here, in order to capture the response, a dual-condition criterion was utilized which calibrated one set of parameters to low temperature quasi-static strain rate data, while the other parameter set is calibrated to high temperature high strain rate data. The thermal effects were captured using temperature dependent thermal constants and invoking thermal transport with conduction, convection, and thermal radiation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Generative diffusion model surrogates for mechanistic agent-based biological models

Mechanistic, multicellular, agent-based models are commonly used to investigate tissue, organ, and organism-scale biology at single-cell resolution. The Cellular-Potts Model (CPM) is a powerful and popular framework for developing and interrogating these models. CPMs become computationally expensive at large space- and time- scales making application and investigation of developed models difficult. Surrogate models may allow for the accelerated evaluation of CPMs of complex biological systems. However, the stochastic nature of these models means each set of parameters may give rise to different model configurations, complicating surrogate model development. In this work, we leverage denoising diffusion probabilistic models (DDPMs) to train a generative AI surrogate of a CPM used to investigate in vitro vasculogenesis. We describe the use of an image classifier to learn the characteristics that define unique areas of a 2-dimensional parameter space. We then apply this classifier to aid in surrogate model selection and verification. Our CPM model surrogate generates model configurations 20,000 timesteps ahead of a reference configuration and demonstrates approximately a 22x reduction in computational time as compared to native code execution. Our work represents a step towards the implementation of DDPMs to develop digital twins of stochastic biological systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Power System Waveform Classification Using Time-Frequency and CNN

Many modern reclosers and circuit breakers have microprocessor relays that record waveforms of system events. In some cases, utilities may record a half-a-dozen event captures for every event. This is thousands of events per year. The industry needs faster, more automated, more conclusive, and easy-to-use systems that can process massive amounts of event recordings without extensive input/support from power system engineers. To address the need for a commercially viable solution that can classify waveform data, energies were directed to develop a universal neural network (NN) structure (deep learning algorithm) that works for a wide variety of system event types. The structure that showed the most promise was one that included the use of spectrograms. The technique has shown positive results in audio engineering, particularly with respect to speech recognition. A waveform signature could be treated as a spoken word like audio waveforms for specific things such as “YES” or “UP”. No two people produce the exact same waveform when speaking each of these words, but audio processing algorithms based on spectrograms and convolutional neural networks (CNN) can still distinguish the word regardless of the speaker. No two circuits produce the exact same waveform for a given event, but the NN can be trained to classify the event type regardless of the circuit or location on the circuit. A Power System Neural Network (PSNN) has been developed to use a CNN to classify events within waveform data for power systems. The waveform is converted to an array of values by way of spectrograms and interpreted as an image. This image is passed into the CNN. The test results on independent simulated test and validation datasets show greater than 99% accuracy. While the results thus far are based on simulated data, the performance of the PSNN is very promising and should work for a wide variety of power system conditions of interest. Ultimately, much of the custom code and tools used today and much of the manual effort expended today may be automated using this PSNN.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2019

Moirai Land Data System (Moirai) v3.0.0

The Moirai Land Data System (Moirai LDS) is designed to produce recent historical land data inputs for the AgLU module of GCAM data system , but the Moirai LDS outputs could also be used by other models/applications. The Moirai are the Greek Fates, and this software is named Moirai to represent the fundamental influence of land data inputs on model outcomes. The primary function of the Moirai LDS is to combine spatially explicit input data (e.g., raster images) with tabular input data (e.g., crop price table) to generate tabular output data for a suite of variables. Some of these outputs replace the data provide by the Global Trade Analysis Project (GTAP), and other data replace and augment the original GCAM GIS processing. The Moirai LDS output data are aggregated by Geographic Land Unit (GLU) within each country. The GLU coverage is an input to the Moirai LDS (as a thematic raster image and an associated CSV file that maps the thematic integers to names), and the GLU boundaries can be determined arbitrarily. Previous versions of GCAM (and Moirai LDS) used only bioclimatic Agro-Ecological Zones (AEZs) and corresponding data that were provided by GTAP, as the GLUs. As a result, some AEZ terminology still exists in the code, but this terminology now refers more generally to GLUs. The Moirai LDS now enables any set of boundaries to be used as GLUs (including AEZs), allowing for more flexible generation of land use region boundaries (defined as the intersection of GLUs with geopolitical regions). The current default set of GLUs is the same set of 235 global watersheds as used by the GCAM water module. The GCAM 5.1 geopolitical regions (32 or 14) are included and used as inputs to Moirai to generate a mapping file between the Moirai outputs, which are at the level of the intersection between the GLUs and the country boundaries, and the geopolitical regions. The diagnostics scripts use this geopolitical region mapping in some cases. Moirai can also recalibrate three of the outputs (crop production, harvested area, and land rent) to a specified year that is the center of a five-year averaging window. No recalibration retains the circa 2000, 7-year average of the source data. The cucrrency-year for land rent can also be specified, and the default is 2001 to match the GTAP data.

Di Vittorio, Alan↗

DOE OSTI2024

Data and scripts associated with the manuscript "A Hydrogeophysical Framework to Assess Infiltration During a Simulated Ecosystem-Scale Flooding Experiment”

This package contains the data and scripts used in "A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment" (Adebayo et al., 2023). The Geophysics folder contains two subfolders containing Ground Penetrating Radar (GPR) and Electrical Resistivity Imaging (ERI) datasets acquired between August 23 and August 26, 2021. The GPR subfolder contains raw data (Tempest_Project.GPZ) and processed results using the Sensors and Software EKKO Project software. The ERI subfolder contains two additional subfolders named Static and Timelapse. The Static subfolder contains background resistivity datasets acquired before infiltration, both raw and inverted results using the ResIpy code. The Timelapse subfolder contains repeated resistivity datasets acquired during infiltration on two profiles, along with raw and inverted results using the ResIpy code. The Flow model folder contains ATS model simulation files. The Geophysics-informed subfolder contains simulation datasets for two profiles generated using information from the geophysics data. The National_database-informed subfolder contains simulation datasets for two profiles generated using information from a national soil database. Each subfolder contains initial simulation results, multiple ensembles, and final simulation files. The Sensitivity_analysis folder contains the observation files exported from the flow model, notebooks used in running the code, and sensitivity analysis results showing the sensitivity of different parameters over a 48-hour period. We used mutual information-based global sensitivity analysis for the sensitivity analysis, following Jiang et al., 2022. This study shows a novel approach that integrates geophysical data with hydrological modeling to understand subsurface flow during coastal flooding.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

Abstract The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm , ipole , and igrmonty , paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Modelling the Alfvén eigenmode induced fast-ion flow measured by an imaging neutral particle analyzer

An imaging neutral particle analyzer (INPA) provides energy and radially resolved measurements of the confined fast-ion population ranging from the high-field side to the edge on the midplane of the DIII-D tokamak. In recent experiments, it was used to diagnose fast-ion flow in the INPA-interrogated phase-space driven by multiple, marginally unstable Alfvén eigenmodes (AEs). Here, the key features of this measured fast-ion flow are: (I) a fast-ion flow from q min and the injection energy (81 keV) towards lower energies and plasma periphery.(II) A flow from the same location towards higher energies and the plasma core, (III) a phase-space 'hole' at the injected energy and plasma core and (IV) a pile-up at the plasma core at lower energies (~60 keV). Ad hoc energetic particle diffusivity modelling of TRANSP significantly deviates from the observation. Comparably, a reduced modelling, i.e. a combination of NOVA-K and ASCOT5 code with the measured mode structure and amplitude, generally reproduce some key features of the observed phase-space flow, but largely failed to interpret fast ion depletion near the plasma axis. At last, self-consistent, first-principle multi-phase hybrid simulations that include realistic neutral beam injection and collisions are able to reproduce most features of the time-resolved phase-space flow. During consecutive hybrid phases, an RSAE consistent with the experiment grows and saturates, redistributing the injected fast ions. The resulting synthetic INPA images are in good agreement with the measurement near the injection energy. The simulations track the fast-ion redistribution within the INPA range, confirming that the measured fast-ion flow follows streamlines defined by the intersection of phase-space surfaces of constant magnetic moment μ and constant E' = nE + ωP φ , where n and ω are the instability toroidal mode number and frequency, and E and P φ the ion energy and toroidal canonical momentum. Nonperturbative effects are required to reproduce the depletion of fast ions near the magnetic axis at the injection energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

MO-AB-206-01: PET Physics

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2022

Investigation of the OECD/NEA PWR MOX/UO{sub 2} core transient benchmark using a coupled whole core pin-by-pin route in WIMS

This paper investigates the new CAMELOT coupling route developed in WIMS, which uses prepared tabulated cross-sections combined with the flux solver MERLIN module and the integrated sub-channel thermal hydraulics solver ARTHUR module to solve for the coupled core state. This coupling route is applied to the OECD/NEA pressurized water reactor (PWR) core transient benchmark exercise, involving 3-D neutronics and thermal-hydraulics analysis of a full-sized PWR core featuring MOX fuel at various core states, as well as simulating a control rod ejection transient. Analysis of the results for each part of the benchmark demonstrate that the high-resolution pin-by-pin analysis used in WIMS closely matches the other participants in the benchmark, with differences of between 20-200 pcm for k-effective, 1-2.5% RMS difference for assembly-averaged powers, 2-8% difference in rod worth calculations, < 200 ppm for calculations of critical boron concentration, zero difference in delayed neutron fraction calculation and a very close match for power response to the rod ejection transient. This is the first application of the WIMS-CAMELOT approach for the transient analysis of a full-sized reactor and as such further improvements such as reducing computational cost through parallelization, optimisation and memory reduction, are currently in active development. This paper demonstrates the accuracy and flexibility of the CAMELOT coupling route, providing a straightforward and consistent user image that can be easily used for a variety of modelling problems involving coupled neutronics and thermal-hydraulics. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗