Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface tracking”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)

This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Quantifying disorder one atom at a time using an interpretable graph neural network paradigm

Abstract Quantifying the level of atomic disorder within materials is critical to understanding how evolving local structural environments dictate performance and durability. Here, we leverage graph neural networks to define a physically interpretable metric for local disorder, called SODAS. This metric encodes the diversity of the local atomic configurations as a continuous spectrum between the solid and liquid phases, quantified against a distribution of thermal perturbations. We apply this methodology to four prototypical examples with varying levels of disorder: (1) grain boundaries, (2) solid-liquid interfaces, (3) polycrystalline microstructures, and (4) tensile failure/fracture. We also compare SODAS to several commonly used methods. Using elemental aluminum as a case study, we show how our paradigm can track the spatio-temporal evolution of interfaces, incorporating a mathematically defined description of the spatial boundary between order and disorder. We further show how to extract physics-preserved gradients from our continuous disorder fields, which may be used to understand and predict materials performance and failure. Overall, our framework provides a simple and generalizable pathway to quantify the relationship between complex local atomic structure and coarse-grained materials phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lifting and removal device for additive manufacturing system

An attachment device for lifting and removing build plates and additive manufactured parts built on the build plates from an additive manufacturing machine. The attachment device may fit onto a lift trolley to cooperatively engage with and lift the build plate upward off of a surface of the additive manufacturing machine. The attachment device may include a support bar, two lifting arms, and a mounting bracket. The support bar may have a slider track formed therethrough, so that the lifting arms can be adjusted to correspond with a width of the build plate. The lifting arms may each having a channel formed along a length thereof and an interfacing portion slidably located within the slider track. The channels may be sized and shaped to engage with a bottom surface and/or opposing edges of the build plates. The mounting bracket may attach to a lifting device of the lift trolley.

Barr, Cary B.↗
DOE OSTI2025

Evaluation of color in digital nuclear power plant control room displays

Human system interface design in industrial process control is guided by industry standards, human factors best practices, and domain-specific conventions, and often there is a conflict between one or more of the sources of design input for specific design elements. In the nuclear domain, one design element for which conflict arises is the use of color to represent equipment state. Here, this study evaluates the tradeoffs associated with using color in a process control display versus using white and shades of gray. The performance metrics were response time, accuracy, and eye movement metrics using a simplified experimental task and professional operators. Results revealed that adhering to color conventions in nuclear power yielded small advantages in simple tasks, but did not exist for more complex tasks. The results did not provide strong evidence for or against using a particular color scheme and revealed the need for further research on the use of color for commercial nuclear power plants and other process control industries.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2018

Modulating factors of hydrologic exchanges in a large-scale river reach: Insights from three-dimensional computational fluid dynamics simulations

Hyporheic exchange is a critical mechanism shaping hydrological and biogeochemical processes along a river corridor. Recent studies on quantifying the hyporheic exchange were mostly limited to local scales due to field inaccessibility, computational demand, and complexity of geomorphology and subsurface geology. Surface flow conditions and subsurface physical properties are well known factors on modulating the hyporheic exchange, but quantitative understanding of their impacts on the strength and direction of hyporheic exchanges at reach scales is absent. In this study, a high resolution computational fluid dynamics (CFD) model that couples surface and subsurface flow and transport is employed to simulate hyporheic exchanges in a 7-km long reach along the main-stem of the Columbia River. Assuming that the hyporheic exchange does not affect surface water flow conditions due to its negligible magnitude compared to the volume and velocity of river water, we developed a one-way coupled surface and subsurface water flow model using the commercial CFD software STAR-CCM+. The model integrates the Reynolds-averaged Navier-Stokes (RANS) equation solver with a realizable k-e two-layer turbulence model, a two-layer all y^+ wall treatment, and the volume of fluid (VOF) method, and is used to simulate hyporheic exchanges by tracking the free water-air interface as well as flow in the river and the subsurface porous media. The model is validated against measurements from acoustic Doppler current profiler (ADCP) in the stream water and hyporheic fluxes derived from a set of temperature profilers installed across the riverbed. The validated model is then employed to systematically investigate how hyporheic exchanges are influenced by surface water fluid dynamics strongly regulated by upstream dam operations, as well as subsurface structures (e.g. thickness of riverbed and subsurface formation layers) and hydrogeological properties (e.g. permeability). The results suggest that the thickness of riverbed alluvium layer is the dominant factor for reach-scale hyporheic exchanges, followed by the alluvium permeability, the depth of the underlying impermeable layer, and the assumption of hydrostatic pressure.

Bao, Jie↗
DOE OSTI2022

Synchrotron X-ray based particle image velocimetry to measure multiphase streamflow and densitometry

Synchrotron-based radiation techniques are increasingly used in the field of flow condition measurements. As a complicated multiphase flow, cavitation has been investigated in a wide variety of industrial fields, to prevent damages on spillways at large dams or intensify the bacteria eradication rate during water treatments. To control the detrimental effects of cavitating flows, the dynamics of and mechanisms affecting cavitation development shall be visualized and identified in real-time. Due to the large reflection and scattering effects of multiphase flow, standard visible-light imaging technologies cannot penetrate vapor bubbles and provide extremely limited information. In this study, the synchrotron source emitted high-flux X-ray pulses are used to capture flow motion and visualize internal structures with sufficiently high spatial and temporal resolutions. The image processing procedures on raw x-ray images are developed using 2D Fourier transform (FT) and wavelet transform (WT) to implement the contrast enhancement and de-noising. Through the edge detection algorithm, tracking either seeded particles or phase interfaces inside the opaque multiphase flow can be achieved to perform the particle image velocimetry. Further, with decomposing particle-only and particle-eliminated subplots from the X-ray phase-contrast image, the time-resolved velocity and void fraction fields are obtained simultaneously, paving a way for further flow condition analysis such as densitometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Development of Modified Perturbation Solutions to the One-Phase Stefan Problems With a Convective Boundary

The classical Stefan problem is used to track the moving solid-liquid interface during the freezing process. Perturbation theory has often been applied to find an approximate analytical solution due to the nonlinearity of the moving interface. However, the Stefan number (i.e., the sensible over latent heat) must be small and usually less than 0.01 to assume the perturbation expansion, which in turn limits the thermal engineering applications. In this study, a modified perturbation solution is developed by adding a correction term after the leading-order solution to be valid for a much wider range of Stefan numbers (i.e., 0.01 less than or equal to Ste less than or equal to 1). Specifically, a one-phase Stefan problem is first formulated subjected to a convective boundary in the Cartesian, cylindrical, and spherical coordinate systems for inward solidification. The leading-order solution is calculated based on the regular perturbation theory, while the correction term is obtained using the Monte-Carlo method and a multi-variant regression. Results show that the correction term has a linear relationship with the Stefan number and is not significantly influenced by the Biot number. The proposed modified perturbation solution can accurately and rapidly predict the nonlinear moving interface motion for the freezing process.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A Numerical Study of Two-Phase Flow and Interfacial Mass Transfer in a Wetted Wall Column for Counter-Current Gas Absorption

We present a numerical study to investigate the effects of surface waves of falling films on the mass transfer from flowing gas phase into the liquid film. The model employs the volume of fluid method to explicitly track the gas-liquid phase interface and a one-fluid formulation to model the multiphase flow coupled with mass transfer between the two phases. The surface waves are generated by the time-varying injection rate of the falling film. The effects of surface wave frequencies and amplitudes on the overall mass transfer also are investigated for a counter-current gas-liquid flow of oxygen and water. Because of the solubility difference in the two phases, the model can capture the concentration discontinuity across the liquid-gas interface. Good agreement is obtained between the numerical results and experimental data reported in the literature. By systematically varying the frequencies and amplitudes of flow rate for the falling film, the surface waves generated on the film surface enhance mass transfer in general. Here, the model prediction shows a proportional increase in mass transfer for small frequencies of surface waves, but saturation at higher frequencies.

42 ENGINEERING↗
DOE OSTI2021

Tribocatalytically-activated formation of protective friction and wear reducing carbon coatings from alkane environment

Minimizing the wear of the surfaces exposed to mechanical shear stresses is a critical challenge for maximizing the lifespan of rotary mechanical parts. In this study, we have discovered the anti-wear capability of a series of metal nitride-copper nanocomposite coatings tested in a liquid hydrocarbon environment. The results indicate substantial reduction of the wear in comparison to the uncoated steel substrate. Analysis of the wear tracks indicates the formation of carbon-based protective films directly at the sliding interface during the tribological tests. Raman spectroscopy mapping of the wear track suggests the amorphous carbon (a-C) nature of the formed tribofilm. Further analysis of the tribocatalytic activity of the best coating candidate, MoN-Cu, as a function of load (0.25–1 N) and temperature (25 °C and 50 °C) was performed in three alkane solutions, decane, dodecane, and hexadecane. Results indicated that elevated temperature and high contact pressure lead to different tribological characteristics of the coating tested in different environments. The elemental energy dispersive x-ray spectroscopy analysis and Raman analysis revealed formation of the amorphous carbon film that facilitates easy shearing at the contact interface thus enabling more stable friction behavior and lower wear of the tribocatalytic coating. These findings provide new insights into the tribocatalysis mechanism that enables the formation of zero-wear coatings.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Dynamics and unsteady morphologies at ice interfaces driven by D 2 O–H 2 O exchange

The growth dynamics of D 2 O ice in liquid H 2 O in a microfluidic device were investigated between the melting points of D 2 O ice (3.8 °C) and H 2 O ice (0 °C). As the temperature was decreased at rates between 0.002 °C/s and 0.1 °C/s, the ice front advanced but retreated immediately upon cessation of cooling, regardless of the temperature. This is a consequence of the competition between diffusion of H 2 O into the D 2 O ice, which favors melting of the interface, and the driving force for growth supplied by cooling. Raman microscopy tracked H/D exchange across the solid H 2 O–solid D 2 O interface, with diffusion coefficients consistent with transport of intact H 2 O molecules at the D 2 O ice interface. At fixed temperatures below 3 °C, the D 2 O ice front melted continuously, but at temperatures near 0 °C a scalloped interface morphology appeared with convex and concave sections that cycled between growth and retreat. Finally, this behavior, not observed for D 2 O ice in contact with D 2 O liquid or H 2 O ice in contact with H 2 O liquid, reflects a complex set of cooperative phenomena, including H/D exchange across the solid–liquid interface, latent heat exchange, local thermal gradients, and the Gibbs–Thomson effect on the melting points of the convex and concave features.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Chapter 18: Learning and Tracking Ad Hoc Fiducial Markers in Spatial Augmented Reality

We describe a spatial augmented reality system with a tangible user interface used to control computer simulations of complex systems. In spatial augmented reality, the user’s physical space is augmented with projected imagery, blending real objects with projected information, and a tangible user interface enables users to manipulate physical objects as controllers for interactive visualizations. Our system learns ad hoc objects in the user’s environment as fiducial markers (i.e., objects that are visually recognized and tracked). When combined with simulation and visualization tools, these interfaces allow the user to control simulations or ensembles of simulations via physical objects using apt metaphors. While other research has leveraged the use of depth cameras, our system enables the use of standard cameras in readily available smartphones and webcams and has an implementation that runs completely in JavaScript in the web browser. We discuss the prerequisite object-recognition requirements for such tangible user interfaces and describe computer-vision and machine-learning algorithms meeting those requirements. We conclude by presenting example applications, which are also available online.

energy distribution↗
DOE OSTI2025

Ultrafast low-temperature metal–insulator interface phonon dynamics and heat transport in a Pt/Gd 3 Fe 5 O 12 heterostructure

Interfacial thermal and acoustic phenomena have an important role in quantum science and technology, including in spintronic and spincaloritronic materials and devices. Simultaneous measurements of the low-temperature thermal and acoustic properties of a metal/insulator heterostructure reveal distinct dynamics in the characteristic phonon frequency ranges of acoustic and thermal transport. The measurements probed a heterostructure consisting of a thin film of Pt on the ferrimagnetic insulator gadolinium iron garnet (Gd 3 Fe 5 O 12 , GdIG) grown epitaxially on a gadolinium gallium garnet substrate. Ultrafast structural dynamics within the Pt layer were tracked using time-resolved ultrafast x-ray diffraction and analyzed to probe interfacial acoustic and thermal properties. The rapid heating of the Pt layer by a 400 nm wavelength femtosecond-duration optical pulse produced transient structural changes that provided the stimulus for these measurements. Rapid heating produced a broadband acoustic pulse that was partially reflected by the Pt/GdIG interface. Temporal frequencies up to 740 GHz, corresponding to angular frequencies of several THz, were detected in a wavelet analysis of the acoustic oscillations of the strain in the Pt layer. The structural results were analyzed to determine (i) the acoustic damping coefficient and phonon mean free path in Pt at frequencies of hundreds of GHz and (ii) the Grüneisen anharmonicity parameter. The thermal conductance of the Pt/GdIG interface was tracked using the slower, tens-of-picosecond-scale, dynamics of the initial cooling of the heated Pt layer. Analysis using a model based on the Boltzmann transport equation shows that the phonon transmission is lower at the phonon frequencies relevant to thermal transport than for subterahertz regime acoustics.

Acoustic phenomena↗
DOE OSTI2019

Monovalent and Divalent Cations at the a-Al2O3(0001)/Water Interface: How Cation Identity Affects Interfacial Ordering and Vibrational Dynamics

Mineral oxide-water interfaces are important for a wide range of industrial, geochemical, and biological processes. The reactivity of these interfaces is strongly impacted by the presence of ions. Thus, it is critical to understand how ions alter the interfacial environment. This can be achieved by measuring the changes in the structure and vibrational dynamics of interfacial water induced by the presence of ions in close vicinity to the mineral surface. The a-Al2O3(0001) surface represents a flexible platform to study the effect of ions on interfacial aqueous environments at positive, neutral and negative surface charge. Using vibrational sum frequency generation (vSFG) in the frequency and time domain, we investigate how monovalent and divalent cations affect the hydrogen bonding environment of the first few layers of interfacial water next to a-Al2O3(0001). Our results indicate that monovalent cations, such as Li+, Na+, K+, and Cs+, appear to have lower binding affinities for the interface compared to Ca2+, Sr2+, and Ba2+. This leads to an interfacial region that is structured in a cation valence dependent manner. Addition of divalent cations at the negatively charged interface (pH 10) increases the spectral intensity in the 3400 cm-1 region compared to neat pH 10 H2O, in contrast to monovalent cations that only attenuate the vSFG signal. Time resolved vSFG measurements reveal that the O-H vibrational lifetime (T1) of interfacial species at pH 10 in the presence of NaCl and BaCl2 remains similar, but restructuring of the surface seen in steady state vSFG is manifested in the degree to which strongly hydrogen bonded species recover to their original populations post excitation. By tracking the accumulation of ions at the interface via the vSFG response, we can characterize the unique surface arrangements of interfacial water molecules induced by a range of monovalent and divalent cations at the a-Al2O3(0001)/water interface.

Piontek, Stefan↗
DOE OSTI2019

Production Operations Management System

The Production Operations Management System (POMS) is a set of software tools which allows production teams and analysis groups across multiple Fermilab experiments to launch, modify and monitor large scale campaigns of related Monte Carlo or data processing jobs. POMS provides a web service interface that enables automated jobs submission on distributed resources according to customers’ requests and subsequent monitoring and recovery of failed submissions, debugging and record keeping. POMS interfaces with existing HEP data access, processing, movement and monitoring tools at Fermilab, including Jobsub, dCache, SAM, and FIFEmon, and in combination with them handles the creation and tracking of unique filenames for multistage campaigns. A flexible, interactive GUI interface exists as part of the system to visualize campaign stages, parallelization tracks and data flow. An important feature of POMS is a one-to-one connection between the GUI campaign visualizer and a text based representation of the complete campaign configuration. An extensible library of template campaign configurations is available. The templates are user modifiable and map cleanly to and from the GUI description of the campaign. Future releases of POMS will interface with Rucio for file movement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Three-Dimensional Tracking of Interfacial Hopping Diffusion

Theoretical predictions have suggested that molecular motion at interfaces—which influences processes including heterogeneous catalysis, (bio)chemical sensing, lubrication and adhesion, and nanomaterial self-assembly—may be dominated by hypothetical “hops” through the adjacent liquid phase, where a diffusing molecule readsorbs after a given hop according to a probabilistic “sticking coefficient.” Here, we use three-dimensional (3D) single-molecule tracking to explicitly visualize this process for human serum albumin at solid-liquid interfaces that exert varying electrostatic interactions on the biomacromolecule. Following desorption from the interface, a molecule experiences multiple unproductive surface encounters before readsorption. An average of approximately seven surface collisions is required for the repulsive surfaces, decreasing to approximately two and a half for surfaces that are more attractive. The hops themselves are also influenced by long-range interactions, with increased electrostatic repulsion causing hops of longer duration and distance. These findings explicitly demonstrate that interfacial diffusion is dominated by biased 3D Brownian motion involving bulk-surface coupling and that it can be controlled by influencing short- and long-range adsorbate-surface interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

TOPAS-nBio: An Extension to the TOPAS Simulation Toolkit for Cellular and Sub-cellular Radiobiology

The TOPAS Monte Carlo (MC) system is used in radiation therapy and medical imaging research, having played a significant role in making Monte Carlo simulations widely available for proton therapy related research. While TOPAS provides detailed simulations of patient scale properties, the fundamental unit of the biological response to radiation is a cell. Thus, our goal was to develop TOPAS-nBio, an extension of TOPAS dedicated to advance understanding of radiobiological effects at the (sub-)cellular, (i.e., the cellular and sub-cellular) scale. TOPAS-nBio was designed as a set of open source classes that extends TOPAS to model radiobiological experiments. TOPAS-nBio is based on and extends Geant4-DNA, which extends the Geant4 toolkit, the basis of TOPAS, to include very low-energy interactions of particles down to vibrational energies, explicitly simulates every particle interaction (i.e., without using condensed histories) and propagates radiolysis products. To further facilitate the use of TOPAS-nBio, a graphical user interface was developed. TOPAS-nBio offers full track-structure Monte Carlo simulations, integration of chemical reactions within the first millisecond, an extensive catalogue of specialized cell geometries as well as sub-cellular structures such as DNA and mitochondria, and interfaces to mechanistic models of DNA repair kinetics. We compared TOPAS-nBio simulations to measured and published data of energy deposition patterns and chemical reaction rates (G values). Our simulations agreed well within the experimental uncertainties. Additionally, we expanded the chemical reactions and species provided in Geant4-DNA and developed a new method based on independent reaction times (IRT), including a total of 72 reactions classified into 6 types between neutral and charged species. Chemical stage simulations using IRT were a factor of 145 faster than with step-by-step tracking. Finally, we applied the geometric/chemical modeling to obtain initial yields of double-strand breaks (DSBs) in DNA fibers for proton irradiations of 3 and 50 MeV and compared the effect of including chemical reactions on the number and complexity of DSB induction. Over half of the DSBs were found to include chemical reactions with approximately 5% of DSBs caused only by chemical reactions. In conclusion, the TOPAS-nBio extension to the TOPAS MC application offers access to accurate and detailed multiscale simulations, from a macroscopic description of the radiation field to microscopic description of biological outcome for selected cells. TOPAS-nBio offers detailed physics and chemistry simulations of radiobiological experiments on cells simulating the initially induced damage and links to models of DNA repair kinetics.

62 RADIOLOGY AND NUCLEAR MEDICINE↗