Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluids”

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 613 records · Page 34

Spatio-temporal Fourier Transformer for Long-term Dynamics Prediction (StFT) v1.0

We propose a novel machine learning model spatio-temporal Fourier transformer (StFT) to emulate long-term dynamics of multi-scale and multi-physics systems. Our method StFT overcomes the limitations of rapid error accumulation, particularly in long-term forecasting of systems characterized by complex and coupled dynamics. StFT achieves outstanding accuracy and computational efficiency by effectively capturing multi-scale interactions, and quantify the uncertainties inherent in the predictions. Our model leverages a structured hierarchy of StFT blocks, and explicitly captures dynamics across both macro- and micro- spatial scales. Evaluations conducted on three benchmark datasets (plasma, fluid, and atmospheric dynamics) demonstrate the advantages of our approach over state-of-the-art ML methods.

Bai, Zhe [Lawrence Berkeley National Laboratory (L↗

Liquid-induced topological transformations of cellular microstructures

The fundamental topology of cellular structures—the location, number and connectivity of nodes and compartments—can profoundly affect their acoustic, electrical5, chemical, mechanical and optical properties, as well as heat, fluid and particle transport. Approaches that harness swelling, electromagnetic actuation and mechanical instabilities in cellular materials have enabled a variety of interesting wall deformations and compartment shape alterations, but the resulting structures generally preserve the defining connectivity features of the initial topology. Achieving topological transformation presents a distinct challenge for existing strategies: it requires complex reorganization, repacking, and coordinated bending, stretching and folding, particularly around each node, where elastic resistance is highest owing to connectivity. Here we introduce a two-tiered dynamic strategy that achieves systematic reversible transformations of the fundamental topology of cellular microstructures, which can be applied to a wide range of materials and geometries. Our approach requires only exposing the structure to a selected liquid that is able to first infiltrate and plasticize the material at the molecular scale, and then, upon evaporation, form a network of localized capillary forces at the architectural scale that ‘zip’ the edges of the softened lattice into a new topological structure, which subsequently restiffens and remains kinetically trapped. Reversibility is induced by applying a mixture of liquids that act separately at the molecular and architectural scales (thus offering modular temporal control over the softening–evaporation–stiffening sequence) to restore the original topology or provide access to intermediate modes. Guided by a generalized theoretical model that connects cellular geometries, material stiffness and capillary forces, we demonstrate programmed reversible topological transformations of various lattice geometries and responsive materials that undergo fast global or localized deformations. In conclusion, we then harness dynamic topologies to develop active surfaces with information encryption, selective particle trapping and bubble release, as well as tunable mechanical, chemical and acoustic properties.

42 ENGINEERING↗

Continued Validation Studies using the MOOSE Framework for Plasma Simulation with Electromagnetics

Resources and tools for the modeling and simulation of low-temperature plasma (LTP) discharges are increasingly vital to progress in the field in order to properly characterize and study complex source designs and plasma chemistries beyond the scope of traditional diagnostics. Open-source software provides powerful platforms for this work and can enable community-driven LTP R&D. Within the Multiphysics Object-Oriented Simulation Environment (MOOSE) open-source framework [1], capabilities have been demonstrated in the areas of plasma fluids (Zapdos [2]), plasma chemistry (CRANE [3]), and general electromagnetic wave theory (Electromagnetic Library for Kinetics & fluids [ELK] [4]). ELK has since been coupled to Zapdos/CRANE to enable fully coupled electromagnetic plasma simulations. This talk will detail the continued validation efforts and discuss Zapdos-ELK-CRANE code coupling with various-low temperature plasma sources. The impact of fully coupled electromagnetics on process parameters (e.g., temperature and electron/ion energy) versus an electrostatic description will also be discussed. [1] Permann et al., SoftwareX, 11 (2020) 100430. [2] Lindsay et al., J. Phys. D: Appl. Phys. 49 (2016) 235204. [3] Keniley et al., Talk, GEC Session DT3.00002 (2019). [4] Icenhour et al., Talk, GEC Session RR1.00006 (2019). Acknowledgements: INL Graduate Fellowship Program, NSF SI2 Grant 1740300, U.S. Dept. of Energy Office of Science Graduate Student Research Program

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synthesis and Single Crystals of Refractory Oxides of Lanthanides and Thorium

At the completion of this program, we can report that we developed a considerable degree of technical improvements in our ability to perform hydrothermal reactions at high temperatures and pressures. We can now routinely perform reactions at 700-750°C and 200 MPa. Currently we are in the process of exploiting this new technology synthesizing a range of exotic new materials investigating relatively poorly understood materials. Our initial efforts focused on the chemistry of rare earth oxides with tetravalent and pentavalent oxides. We recently published a study of the lanthanides with Nb 5+ and Ta 5+ ions, where we grew oxides such as RENdO 4 and RETaO 4 as high quality single crystals. These compounds were targeted as potential hosts for luminescent and scintillation materials, particularly given that they are among the densest oxide hosts and hence have good potential as absorbers for high energy radiation like X-rays and gamma rays. We also isolated a range of unusual new rare earth tantalates with very complex structures. indicating that the chemistry is very sensitive to conditions. We performed some fairly comprehensive examinations of the solid-state chemistry of rare earth ions with various tetravalent metal ions especially Si 4+ , Ge 4+ , Sn 4+ and Ti 4+ . Given the potential role of rare earth silicates in immobilizing radioactive waste elements in long-term storage, and the similarity of our hydrothermal fluids with known geological conditions, this chemistry continues to be relevant. We prepared an extensive series of new lanthanide germanates (e.g. RE 13 Ge 6 O 31 (OH), BaRE 10 (GeO 4 ) 4 O 8 ). and found that there is there is almost no overlap between the chemistry of the rare earth silicates. Stannic oxide (SnO 2 ) is much more refractory and requires higher temperatures and of mineralizer concentrations. One significant result is the growth of RE 2 Sn 2 O 7 pyrochlore single crystals. These are of interest because the rare earth stannate pyrochlores are known to display a wide range of magnetic frustration such as spin ice behavior. We grew high quality single crystals of rare earth germanate and stannate pyrochlores and this led to a collaboration with Professor Kate Ross at Colorado State. Preliminary measurements, indicate that the crystals contain no detectable defects or site disorder. Initial neutron diffraction on single crystals was performed at Oak Ridge, and more detailed experiments involving the Ross group are underway at both NIST and ORNL. This particular chemistry has turned out to be the most potentially significant work on this project and the collaborative effort with Prof. Ross is the topic of a DoE renewal project on quantum materials. Our initial foray into the hydrothermal chemistry of rare earth titanates has also been very promising and a range of cubic and polar ferroic phases of the light rare earths RE 2 Ti 2 O 7 (RE = La - Pr) in the P2 1 phase. We also discovered an interesting new phase Ce 2 Ti 4 O 11 that can have implications in heavy metal immobilization and storage. along with a series of new rare earth titanates (La 5 Ti 4 O 15 (OH) Sm 3 TiO 5 (OH) 3 and Lu 5 Ti 2 O 11 (OH) with exceptionally complex structures. One interesting sidelight has been high temperature hydrothermal chemistry terbium, including the growth of large crystals of TbO(OH). This is not a new compound but it is the first time it has been grown as large single crystals. The Tb atom density is almost as high as that in Tb 2 O 3 and has a very high Verdet constant (ca. 70), making it a very attractive candidate as a Faraday rotator. Unfortunately it is not in a cubic structure but he material is hard, stable, pure and inexpensive, so should still be an attractive Faraday oscillator. We recently received a patent on this material. We also synthesized K 2 Tb(Ge 2 O 7 ) containing stable octahedral Tb 4+ ions, which appears to be the first example of a well-characterized Tb 4+ complex. Given that Tb 4+ has been proposed as a benign surrogate for more treacherous tetravalent ions such as Cf 4+ and Bk 4+ , we think that Tb 4+ silicates can be a particularly useful study for actinide immobilization and related work. We also began reaction studies with rare earths and both ReO 2 and RuO 2 . These resulted in large single crystals of species like RE 5 Ru 2 O 12 , RE 4 Re 2 O 11 , REReO 4 and RE 2 ReO 5 . Several of these samples have already been sent to ORNL for magnetic and neutron diffraction studies.

36 MATERIALS SCIENCE↗

Evaluation and beneficiation of high sulfur and high alkali fly ashes for use as supplementary cementitious materials in concrete

Coal-fuel power plants with semi-dry or dry flue gas desulfurization (FGD) systems produce high sulfur and/or high alkali fly ashes due to comingling of fly ash with FGD products. Such fly ashes do not meet the SO 3 content limit (5.0% max.) of ASTM C618 or are unable to mitigate the alkali-silica reaction. The mineralogy of the sulfur present in these ashes can vary significantly (e.g., CaSO 4 , CaSO 3 , Na2SO 4 ) based on the FGD technology used and this affects the performance of these fly ashes in concrete. Thus, the single SO 3 % limit of ASTM C618 is unable to capture the complexity and performance of fly ash, and this results in elimination of potentially viable pozzolans for concrete. Here this study performs a systematic investigation of the effect of SO 3 type and content in fly ash on various performance parameters of cement-fly ash pastes and mortars, including workability (flow and flow retention), pore fluid pH, setting time, strength development, and potential for deleterious expansion. To better quantify and understand these effects, the study considers both real and doped fly ashes (i.e., a blend of specification-compliant fly ash with target sulfur compounds). The poor performance observed in the case of setting time and pore fluid pH was successfully mitigated using chemical admixtures.

36 MATERIALS SCIENCE↗

The interplay between vapour, liquid, and solid phases in laser powder bed fusion

The capability of producing complex, high performance metal parts on demand has established laser powder bed fusion (LPBF) as a promising additive manufacturing technology, yet deeper understanding of the laser-material interaction is crucial to exploit the potential of the process. By simultaneous in-situ synchrotron x-ray and schlieren imaging, we probe directly the interconnected fluid dynamics of the vapour jet formed by the laser and the depression it produces in the melt pool. The combined imaging shows the formation of a stable plume over stable surface depressions, which becomes chaotic following transition to a full keyhole. We quantify process instability across several parameter sets by analysing keyhole and plume morphologies, and identify a previously unreported threshold of the energy input required for stable line scans. The effect of the powder layer and its impact on process stability is explored. These high-speed visualisations of the fluid mechanics governing LPBF enable us to identify unfavourable process dynamics associated with unwanted porosity, aiding the design of process windows at higher power and speed, and providing the potential for in-process monitoring of process stability.

36 MATERIALS SCIENCE↗

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗

The Art of Automation: Translating Electron Microscopy Workflows Into Automated Processes

Acquiring data using a scanning transmission electron microscope (STEM) is a complex, multi-step process. The intricacy of the process depends on the type of sample, composition of the material, desired results of the experiment, resolution requirement and other experimental factors. Each experiment presents unique complications, such as sample drift and contamination, that the microscopist must consider when acquiring data. All these challenges are handled fluidly and expertly by experienced microscopists, but to reach new levels of innovation in material development, including greater reproducibility, throughput, and precision, the automation of these workflows is essential. The initial phase of this work involved translating intuition-based workflows into discrete, programmable steps. Some common key stages in STEM workflows are the initial tuning, scanning the sample for areas of interest, and then acquiring the data. Each stage can be broken further into specific parameter adjustments, such as aberration correction and dwell time optimization, depending on the experiment. When deconstructing various experiments each step was assessed for automation feasibility based on the amount of real time operator decisions. There are steps that lend themselves to automation more readily than others, such as course focusing and sample screening, but there is potential for full automation of all stages with time. As an initial step, an automated montage routine was developed, allowing for the efficient acquisition of large portions of the sample without requiring continuous intervention from the operator. The automation of this small process of the procedure demonstrates the value of this capability. A major challenge in automation arises from discrepancies between commanded, reported and actual stage movements. Using systematic tests, stage movement was quantified. This error can be corrected algorithmically for more accurate workflows in the future. Expanding automation capabilities would result in larger, more efficient data acquisition which allows for more robust statistical analysis. Additionally, this work lays the groundwork for a closed loop system where machine learning algorithms would intake automatically acquired data and make real time decisions. By progressively automating this instrument, this work establishes the foundation for fully automated experimentation in transmission electron microscopy.

97 MATHEMATICS AND COMPUTING↗

Harvesting Energy from Wastewater by Converting Sewage

This project aims were to develop and demonstrate a scalable, integrated process to convert sewage sludge into renewable natural gas (RNG), enabling wastewater treatment plants (WWTPs) to become net energy producers. The system proposal integrates autothermal hydrothermal liquefaction (AT-HTL), supercritical salt precipitation (SCSP), and hydrothermal gasification (HTG), collectively forming the Supercritical Sludge-to-Gas (SC-S2G) platform. Initially, batch hydrothermal liquefaction reactions were used to screen sewage sludge using AT-HTL (later termed RI-HTL) conversion to biocrude, aqueous and char phases compared to hydrothermal liquefaction (HTL). Significant improvement in biocrude yield using peroxide addition at O:C ratio of 0.05 and under conditions of 300°C for 10 minutes gave 57% biocrude yield and 85% fluid carbon yield (biocrude plus aqueous), while minimizing the loss of carbon to char solids (~7%). Hence, RI-HTL was shown to be effective for conversion of real sewage sludge. The corrosion of the alloy reactor tubes or vessels is an important factor when developing a process that includes an oxidant and a chemically complex feed like sewage sludge. We investigated the corrosion rates on metal alloys at 350°C for 240 hours. Corrosion rates of 0.21 and 0.26 mpy for 304L and 316L stainless steel were measured respectively. The corrosion information obtained in this investigation was utilized by PNNL for design, materials sourcing and construction of the pilot scale continuous flow system.

09 BIOMASS FUELS↗

Geothermal Heat Recovery Complex: Large-Scale, Deep Direct-Use System in a Low-Temperature Sedimentary Basin (Final Report)

A feasibility study of using deep direct-use (DDU) geothermal energy to heat agricultural research facilities (ARFs) was conducted at the University of Illinois at Urbana-Champaign (U of IL) and its similar application to military facilities in the Illinois Basin (ILB). The geothermal energy system (GES) investigated utilizes low-temperature (30–90°C; 90–190°F) geothermal fluid (i.e., brine) from an extraction well that is part of a deep, two-well (doublet) system that extends to the bottom of the ILB. The geothermal reservoir modeled, the Mt. Simon Sandstone (MSS), is about 1,280 m (4,200 feet) deep and 457 m thick (1,500 feet) beneath the U of IL. The DDU GES surface infrastructure includes heat exchangers connected in-parallel to pipelines carrying the geothermal fluid and fresh cold and hot water. Analysis of the GES indicated that the MSS can provide a baseload of 2 MMBtu/hr to heat the ARFs by extracting 954 m³/d (6,000 barrels/day [bbl/d]) of geothermal fluid that has a temperature of 44–46 °C (111–115 °F). In addition to analyzing the levelized cost of heat (LCOH) and life cycle costs, the environmental effects of the DDU GES were evaluated, including reduced greenhouse gas (GHG) emissions and water consumption. Multiple system designs were evaluated and then ranked based on their maximum heating performance, energy efficiency, and cost recovery. This feasibility study identified the key components of the fully-integrated DDU technology that can be implemented, both technically and economically. The results and information from this study provides end-users and policy makers with guidance for additional research on the specific components of DDU technology such that its widespread use can provide an uninterruptible energy source, increase resilience from extreme weather conditions, reduce U.S. dependency on fossil fuels, and reduce greenhouse gas (GHG) emissions. The site-specific part of this study gives U of IL administrators a realistic and pragmatic assessment of the financial resources necessary to add a DDU GES in the MSS to the campus’ energy portfolio.

15 GEOTHERMAL ENERGY↗

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗

Using defects as a ‘fossil record’ to help interpret complex processes during additive manufacturing: as applied to raster-scanned electron beam powder bed additively manufactured Ti–6Al–4V

Abstract Defects in parts produced by additive manufacturing, instead of simply being perceived as deleterious, can act as important sources of information associated with the complex physical processes that occur during materials deposition and subsequent thermal cycles. Indeed, they act as materials-state ‘fossil’ records of the dynamic AM process. The approach of using defects as epoch-like records of prior history has been developed while studying additively manufactured Ti–6Al–4V and has given new insights into processes that may otherwise remain either obscured or unquantified. Analogous to ‘epochs,’ the evolution of these defects often is characterized by physics that span across a temporal length scale. To demonstrate this approach, a broad range of analyses including optical and electron microscopy, X-ray computed tomography, energy-dispersive spectroscopy, and electron backscatter diffraction have been used to characterize a raster-scanned electron beam Ti–6Al–4V sample. These analysis techniques provide key characteristics of defects such as their morphology, location within the part, complex compositional fields interacting with the defects, and structures on the free surfaces of defects. Observed defects have been classified as banding, spherical porosity, and lack of fusion. Banding is directly related to preferential evaporation of Al, which has an influence on mechanical properties. Lack-of-fusion defects can be used to understand columnar grain growth, fluid flow of melt pools, humping, and spattering events. Graphical abstract

36 MATERIALS SCIENCE↗

Thermo-Mechanical Distortion of Tungsten-Coated Steel During High Heat Flux Testing Using Plasma Arc Lamps

An experimental setup and a test section were designed and fabricated for high heat flux testing (HHFT) of neutron-irradiated specimens using water-wall plasma arc lamps. Because of the radiological considerations and limitations of reactor irradiation, the size of the test articles was limited to disks less than 10 mm in diameter. The specimen was clamped onto an actively cooled block, and clamping allowed the insertion of several thermocouples on the back surface of the specimen through a copper (Cu) block. Five vacuum plasma sprayed tungsten (W)–coated F82H steel specimens were subjected to HHFT. Surface profilometry measurements, which were conducted after HHFT, revealed central bowing of the top W surface. This type of residual distortion occurred for all of the specimens, and the larger the specimens were, the larger was the distortion.In an attempt to understand specimen distortion and address the science questions related to the testing of subsize specimens during HHFT, a simplified thermo-mechanical model was developed. By using a measured temperature in the Cu as an isothermal boundary condition, the model eliminated the need for coupling cooling fluid flow models with stress models, greatly simplifying the analysis. The main variable in the proposed model is hC, i.e., the thermal contact conductance between the F82H and the Cu washer. Inelastic properties, including hardening properties, were considered for F82H steel and Cu. Numerical simulation results demonstrated a buildup of residual deformation during HHFT and a very complex state of stress and deformation during typical heat flux (HF) cycling. Additionally, hoop stress evolution during a high heat flux cycle reveals that F82H at an interface with W would be mainly in compression during HF application and experienced a transition to a tension state during cooldown. Also, specimen distortion evolves during each HF cycle, as the specimen bows downward during HF application and upward during the cooldown period between HF cycles. The final specimen distortion, i.e., upward bowing of the specimen center, was qualitatively predicted for hC values of 4000 to 5000 W/(m 2 ·K). This hC range of values, for which bulging is obtained, is at the lower spectrum of the range of values for hC, consistent with the low thermal contact conductance expected from the unpolished F82H surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fabricating Single Crystal Quantum Dot Solids

The central goal of this research program was to establish fundamental synthesis and processing principles to enable the fabrication of single crystal quantum dot solids (QDS) with programmable structure (i.e., hexagonal and square) and composition. Our approach towards that goal leveraged access to and experience with unique in-situ, multi-probe characterization techniques to understand and ultimately control the fundamental relationship between processing conditions and nucleation and growth of QDS. The program integrated synthesis and fabrication (Hanrath) with in-situ TEM analysis (Kourkoutis) and computational modeling (Clancy) to gain atomic-level insights into the underlying physical phenomena governing assembly and attachment and to guide the development of optimized processing methods. We have established a foundational understanding of physicochemical interactions of self-assembly at a fluid interface, superlattice structure transformation pathways, the critical role of disorder during the initial dimerization of colloidal quantum dot monomers, and the residual strain distribution within the inter-dot epitaxial bridge which hampers the formation of novel electronics states of the quantum dot solids. Collectively, these insights have contributed towards advancing the mechanistic understanding of the complex choreography of assembly and attachment as well as understanding what currently limits further advances in high-fidelity quantum dot solids and tiles.

36 MATERIALS SCIENCE↗

Effects of oxide surface chemistry on diffusioosmosis

Diffusioosmosis is the movement of fluid induced by gradients in solute concentration. Recent studies suggest that in low permeability rocks, it may be the dominant mode of reactant transport and thus control rates of diagenesis, which cannot be adequately explained by pressure-driven flow alone. In this paper we investigate how the equilibrium between an oxide mineral and the surrounding fluid phase influences the diffusioosmotic velocity. We have developed a theory for the case of a thin double layer, where the Debye length is smaller than the characteristic pore size. Several factors contribute to the total fluid velocity: the chemical structure of the mineral surface, the electrolyte type and concentration gradient, and the solution pH. Individual factors can act in concert or in opposition, leading to widely varying magnitudes and directions of the velocity. The numerical results are within the range of the limited experimental data. Our results highlight how surface charging and surface complexation impact the flow, and how they depend on pH.

Diagenesis↗

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FY22 Progress on Multicontinuum Methods in Containment

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as containment transport and/or containment of gases in the subsurface. To estimate rock fluid flow properties and subsequently predict physically realistic processes such as patterns and timing of water, gas, and energy (e.g., heat) movement in the subsurface, laboratory spontaneous water imbibition with simultaneous temperature measurement and numerical modeling methods are presented in the FY22 progress report. A multiple-overlapping-continua conceptual model is used to explain and predict observed complex multi-phenomenological laboratory test behavior during spontaneous imbibition experiments. This report primarily addresses two complexities that arise during the experiments: 1) capturing the late-time behavior of spontaneous imbibition tests with dual porosity; and 2) understanding the thermal perturbation observed at or ahead of the imbibing wetting front, which are associated with adsorption of water in initially dry samples. We use numerical approaches to explore some of these issues, but also lay out a plan for further laboratory experimentation and modeling to best understand and leverage these unique observations.

58 GEOSCIENCES↗