Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bulk techniques”

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 793 records · Page 44

Hydrogen Refueling Reference Station Lot Size Analysis for Urban Sites

Hydrogen Fueling Infrastructure Research and Station Technology (H2FIRST) is a project initiated by the DOE in 2015 and executed by Sandia National Laboratories and the National Renewable Energy Laboratory to address R&D barriers to the deployment of hydrogen fueling infrastructure. One key barrier to the deployment of fueling stations is the land area they require (i.e. "footprint"). Space is particularly a constraint in dense urban areas where hydrogen demand is high but space for fueling stations is limited. This work presents current fire code requirements that inform station footprint, then identifies and quantifies opportunities to reduce footprint without altering the safety profile of fueling stations. Opportunities analyzed include potential new methods of hydrogen delivery, as well as alternative placements of station technologies (i.e. rooftop/underground fuel storage). As interest in heavy-duty fueling stations and other markets for hydrogen grows, this study can inform techniques to reduce the footprint of heavy-duty stations as well. This work characterizes generic designs for stations with a capacity of 600 kg/day hydrogen dispensed and 4 dispenser hoses. Three base case designs (delivered gas, delivered liquid, and on-site electrolysis production) have been modified in 5 different ways to study the impacts of recently released fire code changes, colocation with gasoline refueling, alternate delivery assumptions, underground storage of hydrogen, and rooftop storage of hydrogen, resulting in a total of 32 different station designs. The footprints of the base case stations range from 13,000 to 21,000 ft 2 . A significant focus of this study is the NFPA 2 requirements, especially the prescribed setback distances for bulk gaseous or liquid hydrogen storage. While the prescribed distances are large in some cases, these setback distances are found to have a nuanced impact on station lot size; considerations of the delivery truck path, traffic flow, parking, and convenience store location are also important. Station designs that utilize underground and rooftop storage can reduce footprint but may not be practical or economical. For example, burying hydrogen storage tanks underground can reduce footprint, but the cost savings they enable depend on the cost of burial and the cost land. Siting and economic analysis of station lot sizes illustrate the benefit of smaller station footprints in the flexibility and cost savings they can provide. This study can be used as a reference that provides examples of the key design differences that fueling stations can incorporate, the approximate sizes of generic station lots, and considerations that might be unique to particular designs.

08 HYDROGEN↗

Removal of Metal Oxides from Reactor Metals/Alloys by Laser Ablation

Storage of low and intermediate level radioactive waste in a deep geological repository (L and ILW DGR): - Large volumes of metallic waste = increased production of hydrogen gas. - Large fraction of the metallic wastes are out-of-core reactor components. - Carbon steel feeder pipes. - Radioactive contamination is concentrated within the surface oxides. - Removal of surface contaminants would lead to volume reduction of metallic waste that is stored in L and ILW DGR. - Current decontamination techniques generate large volumes of secondary waste. Objectives: Oxide removal from carbon steel using a ytterbium fiber laser (1064 nm) with a particular focus on identifying the key laser parameters that affect the removal efficiency. Advantages of laser ablation: Generate low volumes of waste, Flexibility of deployment, Remote operation. Grow an oxide film on carbon steel with the incorporation of cobalt. Optimize laser cleaning technology for the removal of surface oxides. Incorporate an oxide collection system for the ablated particles. Conclusions: Objective 1 - Simulated inactive oxides on carbon steel, using Co{sup 2+} as a surrogate for Co radionuclide. Objective 2 - Optimized laser parameters (200 MW/cm{sup 2}, 3 scans, 15 μm overlap); Complete removal of cobalt from surface and bulk metal with optimized cleaning settings. Objective 3 - Glass fiber media H14 HEPA filters had highest collection yield; High oxide collection yields (∼100%). Future Work: Further testing of the liquid trap will verify the efficacy of the HEPA filter and determine the amount of ablated oxides that pass through. Testing of the laser ablation system needs to be done on other reactor materials along with the incorporation of different transition metal cations. Improvements to the oxide collection system are required to improve collection yields.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Picosecond carrier dynamics in InAs and GaAs revealed by ultrafast electron microscopy

Understanding the limits of spatiotemporal carrier dynamics, especially in III-V semiconductors, is key to designing ultrafast and ultrasmall optoelectronic components. However, identifying such limits and the properties controlling them has been elusive. Here, using scanning ultrafast electron microscopy, in bulk n-GaAs and p-InAs, we simultaneously measure picosecond carrier dynamics along with three related quantities: subsurface band bending, above-surface vacuum potentials, and surface trap densities. We make two unexpected observations. First, we uncover a negative-time contrast in secondary electrons resulting from an interplay among these quantities. Second, despite dopant concentrations and surface state densities differing by many orders of magnitude between the two materials, their carrier dynamics, measured by photoexcited band bending and filling of surface states, occur at a seemingly common timescale of about 100 ps. This observation may indicate fundamental kinetic limits tied to a multitude of material and surface properties of optoelectronic III-V semiconductors and highlights the need for techniques that simultaneously measure electro-optical kinetic properties.

42 ENGINEERING↗

Structure of Molten Alkali Chlorides at Charged Interfaces and the Prediction and Interpretation of Their X-ray Reflectivity

The fundamental properties of molten salts have been the subject of research that spans a century. Yet, in the past few years, there has been an unprecedented surge in interest for these systems in the bulk and under confinement by walls and interfaces including under applied potentials. This is driven by the prospect of exciting and very practical energy technologies, including those in the solar and nuclear fields. This article sets to answer two simple but fundamental questions. How does the liquid structure of alkali chlorides change at a real interface when it is charged? Also, how would such changes on the liquid side of the interface be detected in X-ray reflectivity experiments? We use an interface mimicking conductive diamond, which because of its lattice spacing, is an excellent choice for reflectivity experiments. The reason for our interest in X-ray reflectivity is that, as opposed to electrochemical measurements alone, this is likely the only technique in which atomic level information at the liquid side of the interface can be gained under the extreme temperature environments of molten salts. As it will become apparent, the interpretation of reflectivity results in terms of atomic positions is complex when multiple species with different X-ray contrasts on the liquid side are considered. A theoretical scheme termed “the peaks and antipeaks analysis of reflectivity” originally introduced in our prior work (J. Phys. Chem. C2019, 123 (8), 4914–4925) is expanded to interpret the structural changes at the interface as a function of applied electrical bias.

36 MATERIALS SCIENCE↗

Investigating the mechanical stability of flexible metal–organic frameworks

As we continue to develop metal–organic frameworks (MOFs) for potential industrial applications, it becomes increasingly imperative to understand their mechanical stability. Notably, amongst flexible MOFs, structure-property relationships regarding their compressibility under pressure remain unclear. In this work, we conducted in situ variable pressure powder X-ray diffraction (PXRD) measurements up to moderate pressures (<1 GPa) using a synchrotron source on two families of flexible MOFs: (i) NU-1400 and NU-1401, and (ii) MIL-88B, MIL-88B-(CH 3 ) 2 , and MIL-88B-(CH 3 ) 4 . In this project scope, we found a positive correlation between bulk moduli and degree of flexibility, where increased rigidity (e.g., smaller swelling or breathing amplitude) arising from steric hindrance was deleterious, and observed reversibility in the unit cell compression of these MOFs. This study serves as a primer for the community to begin to untangle the factors that engender flexible frameworks with mechanical resilience.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal Probes of Phonon-Coupled Kitaev Spin Liquids: From Accurate Extraction of Quantized Edge Transport to Anyon Interferometry

Recent theoretical studies inspired by experiments on the Kitaev magnet α-RuCl 3 highlight the nontrivial impact of phonons on the thermal Hall conductivity of chiral topological phases. Here, we introduce mixed mesoscopic-macroscopic devices that allow refined thermal-transport probes of non-Abelian spin liquids with Ising topological order. These devices feature a quantum-coherent region with quantized or negligible phonon conductance, flanked by macroscopic lobes that facilitate efficient thermalization between chiral Majorana edge modes and bulk phonons. We show that our devices enable (i) accurate determination of the quantized thermal Hall conductivity, (ii) identification of non-Abelian Ising anyons via the temperature dependence of the thermal conductance, and, most interestingly, (iii) single-anyon detection through heat-based anyon interferometry. Analogous results apply broadly to phonon-coupled chiral topological orders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experimental Study of the Fundamental Properties of Warm Dense Mixtures

The aim of the proposed research was to provide a multi-scale study of the properties of warm dense hydrocarbons, by studying the thermodynamic properties through the equation of state and the microscopic properties by x-ray scattering. Understanding the fundamental properties of warm dense mixtures is an intellectual challenge due to the complexity of the system. Unlike liquids or gases, where constituent particle interactions occur through collisions or bonding between valence electrons, in these strongly coupled systems the atoms will be partially ionized and compressed so tightly together that interactions between the inner core electrons can play a role in the systems chemistry. Advances in computational capabilities and development of new theoretical models have been used to predict the properties of mixtures but there is currently no experimental data of the fundamental interaction between the particles in mixtures to test these predictions against. The proposed research was to experimentally investigate the interaction of the different species in hydrocarbons by measuring the compressibility of substances with different carbon and hydrogen ratios and the complexity of the microscopic interactions through elastic and inelastic x-ray scattering. To study the bulk properties of warm dense hydrocarbons we established the Warm Dense Matter Research Laboratory (WDMRL) in the Institute for Shock Physics at Washington State University. The goal was for experiments in the WDMRL to determine shock loading conditions of interest in the hydrocarbon mixtures. We planned on using Hugoniot EOS measurements and a range of carbon and hydrogen concentrations to determine conditions when the EOS of the mixture varied significantly from that of the classical mixing model. Even though shock transit measurements through aluminum foils suggested pressures upto 400GPa, the experiments in the WDMRL were unsuccessful in getting usable shockwave compression data above 100GPa in polystyrene which was below the pressure of interest for hydrocarbon mixtures (>200GPa). To complete the project, we tested a technique using x-ray phase contrast imaging to map the location of tracer layers in a test sample of polycarbonate to record the material motion in dynamically compressed samples. A technique that will be useful for future warm dense matter experiments. These experiments used <200nm gold layers in polycarbonate samples to measure the material velocity and shock speed using x-ray phase contrast imaging. These results were compared to continuum surface measurements performed at ISP and show that the tracer layer technique can measure hydrodynamic properties accurately in dynamically compressed materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Iodination of PEGylated Polymers Counteracts the Inhibition of Fibrinogen Adsorption by PEG

Poly(ethylene glycol), PEG, known to inhibit protein adsorption, is widely used on the surfaces of biomedical devices when biofilm formation is undesirable. Poly(desaminotyrosyl-tyrosine ethyl ester carbonate), PDTEC, PC for short, has been a promising coating polymer for insertion devices, and it has been anticipated that PEG plays a similar role if it is copolymerized with PC. Earlier studies show that no fibrinogen (Fg) is adsorbed onto PC polymers with PEG beyond the threshold weight percentage. This is attributed to the phase separation of PEG. Further, iodination of the PC units in the PC polymer, (I 2 PC), has been found to counteract this Fg-repulsive effect by PEG. Here, in this study, we employ surface-sensitive X-ray techniques to demonstrate the surface affinity of Fg toward the air–water interface, particularly in the presence of self-assembled PC-based film, in which its constituent polymer units are assumed to be much more mobile as a free-standing film. Fg is found to form a Gibbs monolayer with its long axis parallel to the aqueous surface, thus maximizing its interactions with hydrophobic interfaces. It influences the amount of insoluble, surface-bound I 2 PC likely due to the desorption of the formed Fg–I 2 PC complex and/or the penetration of Fg onto the I 2 PC film. The results show that the phase behavior at the liquid–polymer interface shall be taken into account for the surface behavior of bulk polymers surrounded by tissue. The ability of PEG units rearranging into a protein-blocking layer, rather than its mere presence in the polymer, is the key to antifouling characteristics desired for polymeric coating on insertion devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Picosecond Carrier Dynamics in InAs and GaAs Revealed by Ultrafast Electron Microscopy

Understanding the limits of spatio-temporal carrier dynamics, especially in III-V semiconductors, is key to designing ultrafast and ultra-small optoelectronic components. However, identifying such limits and the properties controlling them has been elusive. Here, using scanning ultrafast electron microscopy (SUEM), in bulk n-GaAs and p-InAs, we simultaneously measure picosecond carrier dynamics along with three related quantities: sub-surface band bending, above-surface vacuum potentials, and surface trap densities. We make two surprising observations. First, we uncover a negative-time contrast in secondary electrons resulting from an interplay among these quantities. Second, despite dopant concentrations and surface state densities differing by many orders of magnitude between the two materials, their carrier dynamics, measured by photo-excited band bending and filling of surface states, occur at a seemingly common timescale of about 100 ps. This observation may indicate fundamental kinetic limits tied to a multitude of material and surface properties of optoelectronic III-V semiconductors, and highlights the need for techniques that simultaneously measure electro-optical kinetic properties.

36 MATERIALS SCIENCE↗

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Start-up of Silo 130 Waste Retrieval at La Hague: first Step Towards Reducing Legacy Inventory - 20020

On June 25, 2019, the first grapple of waste and debris was lifted from Silo 130, one of the most challenging legacy waste storage silos at La Hague. This was the achievement of more than ten years of studies and efforts focused on solving some very unique technical and safety challenges, and on designing an innovative, efficient solution. The waste is from de-cladding materials from past Gas-Cooled Reactor fuel reprocessing operations. It is an unsorted mixture of magnesium cladding pieces, graphite debris, and various other metallic waste types. These materials can be chemically reactive and therefore need specific dispositions for handling and storing. The building is more than 40 years old and does not meet modern standards. French safety authorities have required the silo to be emptied. After about three years of operation, the bulk of this waste will be removed and placed in safe storage, pending the selection of a final conditioning solution. Removal of bottom debris will then be possible. When finished, this operation will bring a conclusion to a significant legacy-related hazard on the site. The silo will then be available for final cleaning and de-activation. This paper will describe the technical, safety, and economical challenges, and the various decisions and optimizations that led to the definition of the process and technology to be implemented. Several innovative features were developed and implemented, in all sections of the project: project organization, progressive approach for the waste disposition strategy, light building design, customized and robust retrieval technology, processing technologies involving the most modern robotic and sorting techniques, storage method... Several of these first-of-a-kind innovative features have been developed and qualified specifically for this application, using the detailed qualification procedure specified for all La Hague projects. The facility was built in less than five years, overcoming various difficulties along the line. Commissioning is now finished with the combined efforts of the project, engineering, start-up, and operating teams assisted by dedicated suppliers. The facility is now ramping-up. The paper will also report on the commissioning and start-up of the facility and describe the first operational results. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

ATR-SEIRAS Method to Measure Interfacial pH during Electrocatalytic Nitrate Reduction on Cu

This study reports the accuracy and applications of an attenuated total reflectance–surface-enhanced infrared absorption spectroscopy (ATR–SEIRAS) technique to indirectly measure the interfacial pH of the electrolyte within 10 nm of the electrocatalyst surface. This technique can be used in situ to study aqueous electrochemical reactions with a calibration range from pH 1–13, time resolution down to 4 s, and an average 95% confidence interval of 14% that varies depending on the pH region (acidic, neutral, or basic). The method is applied here to electrochemical nitrate reduction at a copper cathode to demonstrate its capabilities, but is broadly applicable to any aqueous electrochemical reaction (such as hydrogen evolution, carbon dioxide reduction, or oxygen evolution) and the electrocatalyst may be any SEIRAS-active thin film (e.g., silver, gold, or copper). The time-resolved results show a dramatic increase in the interfacial pH from pH 2–7 in the first minute of operation during both constant current and pulsed current experiments where the bulk pH is unchanged. Attempts to control the pH polarization at the surface by altering the electrochemical operating conditions—lowering the current or increasing the pulse frequency—showed no significant change, demonstrating the challenge of controlling the interfacial pH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Milestone 1.2.15: Feasibility of In Situ Accelerated Aluminum Coupon Radiolysis by Synchrotron X-Rays

Extended dry storage of aluminum-clad spent nuclear fuel (ASNF) requires an assessment of the radiolytic generation of molecular hydrogen gas (H2) from the ASNF’s corrosion layers. H2 accumulation can potentially lead to storage canister embrittlement/rupture and the accumulation of flammable gas mixtures in accident scenarios. To date, computational models have been developed for the prediction of radiolytic H2 generation from samples exposed to lower absorbed dose regimes (up to ~3 MGy). However, greater accuracy is needed for higher absorbed doses (> 25 MGy) where the concentration of H2 ultimately reaches a steady-state. Data in this area is limited due to the amount of time taken (months) to accumulate such high gamma doses. Furthermore, a recent study observed radiation-induced damage to the surface of pre-corroded aluminum alloy coupons at high absorbed gamma doses. Despite this observation and its implications, there is currently no computational connection between the radiolytic formation of H2 and changes in the composition and morphology of the cladding’s corrosion layers with absorbed dose. To develop a better understanding of the processes and effects described above, the present study aimed to evaluate the feasibility of leveraging synchrotron x-ray capabilities for coupled accelerated irradiation and in situ surface characterization of ASNF alloys. To that end, National Synchrotron Light Source II (NSLS-II) beamtime was secured and a series of aluminum alloy 1100 (AA1100) wires, prepared under a variety of conditions, were interrogated using ex situ x-ray diffraction (XRD) and scanning electron microscopy (SEM) techniques at Idaho National Laboratory (INL) and ex situ and in situ XRD capabilities at the NSLS-II x-ray powder diffraction (XPD) beamline. Results indicated that performing synchrotron XRD using the available XPD beamline configuration can provide sufficient radiation dose but cannot discern changes in the composition and morphology of sample corrosion layers. Additionally, ex situ XRD and SEM data from Idaho National Laboratory (INL) indicated that the pre-corrosion of AA1100 wires did not generate an appreciably thick corrosion layer, as compared with previous aluminum coupon studies. These thinner surface corrosion layers were not sufficient for detection by both ex situ and in situ synchrotron XRD, specifically under consideration of the NSLS-II XPD beamline’s configuration primarily capturing information from the sample’s bulk material, i.e., aluminum metal. In summary, the use of the NSLS-II XPD beamline for promoting accelerated radiolysis and in situ surface characterization is not appropriate for this program’s goals.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Nano-imaging of strain-tuned stripe textures in a Mott crystal

The 4d transition metal perovskites Ca n+1 Ru n O 3n+1 have attracted interest for their strongly interacting electronic phases showing pronounced sensitivity to controllable stimuli like strain, temperature, and even electrical current. Through multi-messenger low-temperature nano-imaging, we reveal a spontaneous striped texture of coexisting insulating and metallic domains in single crystals of the bilayer ruthenate Ca 3 (Ti x Ru 1-x ) 2 O 7 across its first-order Mott transition at T ≈ 95K. We image on-demand anisotropic nucleation and growth of these domains under in situ applied uniaxial strain rationalized through control of a spontaneous Jahn-Teller distortion. Our scanning nano-susceptibility imaging resolves the detailed susceptibility of coexisting phases to strain and temperature at the transition threshold. Comparing these nano-imaging results to bulk-sensitive elastoresistance measurements, we uncover an emergent “domain susceptibility” sensitive to both the volumetric phase fractions and elasticity of the self-organized domain lattice. Our combined susceptibility probes afford nano-scale insights into strain-mediated control over the insulator-metal transition in 4d transition metal oxides.

36 MATERIALS SCIENCE↗

On the holographic dual of a symmetry operator at finite temperature

Topological symmetry operators of holographic large 𝑁 CFT 𝐷 ’s are dual to dynamical branes in the gravity dual AdS 𝐷+1 . We use this correspondence to establish a dictionary between thermal expectation values of symmetry operators in the Euclidean CFT 𝐷 and the evaluation of gravitational saddles in the presence of a dynamical brane. Expectation values of 0-form symmetry operators in the CFT 𝐷 are then related to branes wrapped on volume minimizing cycles in the bulk, i.e., the Euclidean continuation of a black hole horizon. We illustrate with some representative examples, including gravity in AdS 3 , duality/triality defects in four-dimensional 𝒩 = 4 super Yang-Mills theory, and the dual of R-symmetry operators probing five-dimensional Bogomol’nyi–Prasad–Sommerfield black holes.

Anomalies↗

Calibration-Free Quantitative Analysis of Lithium-Ion Battery (LiB) Electrode Materials Using Laser-Induced Breakdown Spectroscopy (LIBS)

Recent years have seen a surge in the demand for high-performance battery electrode materials for automotive and various electronic device applications. Improving battery performance requires precise knowledge of the structure–composition properties of active electrode materials. To this effect, quantitative and precise estimation of the composition of advanced electrode materials, containing trace amounts of dopants, provide immense value toward developing next-generation high-capacity battery materials. Herein, we demonstrate the application of calibration-free laser-induced breakdown spectroscopy (LIBS) as a powerful analytical tool for rapid and reliable quantitative spectrochemical characterizations of layered Li metal oxide cathodes containing Mo and Cr dopants (<5 atom %). Specifically, we employ the LIBS using an internal calibration methodology to establish the quantitative elemental ratios of major (Ni, Mn, Co) and trace dopant (Cr, Mo) transition metals to bulk Li contents in diverse cathode material samples such as LiNi 0.5 Mn 0.5 O 2 (NM-50/50), LiNi 0.33 Mn 0.33 Co 0.33 O 2 (NMC), LiNi 0.317 Mn 0.317 Co 0.317 Cr 0.05 O 2 (Cr-doped NMC), and LiNi 0.5– x /2 Mn 0.5– x /2 Mo x O 2 ( x = 0.03, 0.04, 0.05) (Mo-doped NM-50/50) that were synthesized via sol–gel routes. The LIBS-estimated elemental compositions are in good agreement with the nominal stoichiometric values. Furthermore, ex situ LIBS characterizations presented here pave the path for future high-efficacy application of calibration-free quantitative LIBS for rapid in situ analyses of elemental composition changes in battery electrode materials under operation that need not resort to off-site analytical techniques requiring cumbersome sample preparations and/or external standards.

25 ENERGY STORAGE↗

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N &#x3C7; according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗