Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sharp features”

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.

504 records · Page 28

Synthesis and morphology of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals with tungsten bronze structure in RE{sub 2}O{sub 3}-BaO-Nb{sub 2}O{sub 5}-B{sub 2}O{sub 3} glasses (RE: Sm, Eu, Gd, Dy, Er)

Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals with a tetragonal tungsten bronze (TTB) structure are synthesized using a conventional glass crystallization technique in 2.3RE{sub 2}O{sub 3}-27.4BaO-34.3Nb{sub 2}O{sub 5}-36B{sub 2}O{sub 3} (mol%) (RE=Sm, Eu, Gd, Dy, and Er) glasses. One sharp crystallization peak is observed at {approx}670 Degree-Sign C in both powdered and bulk glasses, and the formation of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals with unit cell parameters of a{approx}1.24 nm and c{approx}0.39 nm was confirmed. It is found from high resolution transmission electron microscope observations that the morphology of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals is ellipsoidal. Their average particle size is in the range of 15-60 nm and decreases with decreasing ionic radius of RE{sup 3+} being present in the precursor glasses. The optical transparent crystallized glass (bulk) shows the total photoluminescence (PL) quantum yield of 53% in the visible region of Eu{sup 3+} ions, suggesting a high potential of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals as PL materials. - Graphical abstract: This figure shows a TEM photograph for the heat-treated (667 Degree-Sign C, 3 h) sample of 2.3Dy{sub 2}O{sub 3}-27.4BaO-34.3Nb{sub 2}O{sub 5}-36B{sub 2}O{sub 3}. An ellipsoidal-shaped Ba{sub 1-x}Dy{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystal with diameters of 17 and 28 nm is observed. The ellipsoidal morphology is a common feature in Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals synthesized by the crystallization of 2.3RE{sub 2}O{sub 3}-27.4BaO-34.3Nb{sub 2}O{sub 5}-36B{sub 2}O{sub 3} glasses. Highlights: Black-Right-Pointing-Pointer Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals with a tetragonal tungsten bronze structure are synthesized. Black-Right-Pointing-Pointer A glass crystallization technique was applied. Black-Right-Pointing-Pointer The morphology of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals is ellipsoidal. Black-Right-Pointing-Pointer The average particle size with 15-60 nm decreases with decreasing ionic radius of RE{sup 3+}. Black-Right-Pointing-Pointer Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals have a high potential as photoluminescence materials.

36 MATERIALS SCIENCE↗

Coherent Vibrational Wavepacket Dynamics in Platinum(II) Dimers and Their Implications

Vibrational coherence in the metal-metal-to-ligand-charge transfer (MMLCT) excited state of cyclometalated platinum dimers with a pseudo C 2 symmetry was investigated where two nearly degenerate transitions from HOMO (metal-metal σ* orbital) to higher energy ligand π* orbitals could be simultaneously induced. We observed oscillatory features in femtosecond degenerate transient absorption (TA) signals from complexes [(ppy)Pt(μ- t Bu 2 pz)] 2 (1) and anti-[(ppy)Pt(μ-pyt)] 2 (2), which are attributed to coherent nuclear motions that modulate the HOMO (antibonding σ*) energy level, and hence the energy for the MMLCT transition. The characteristics of such coherent nuclear motions, such as the oscillatory frequency and the dephasing time, differ between 1 and 2 and are explained by mainly two structural factors that could influence the vibrational coherence, the Pt-Pt distance (2.97 Å for 1 vs. 2.85 Å for 2) and molecular shape (1 in an “A” frame vs. 2 in an “H” frame). Because the electronic coupling between the two Pt atoms determines the energy splitting of the bonding σ and antibonding σ* orbital, the Pt-Pt stretching mode coupled to the MMLCT transition changes the inter Pt distance from that of the ground state. Interestingly, while 1 shows a single PtPt stretching frequency of 120 cm-1 in the MMLCT state, 2 exhibits multiple downshifted frequencies (80 and 105 cm-1) in the MMLCT state along with a shorter vibrational dephasing time than 1. Based on the ground state optimized structures and Raman calculations, the changes evident in the vibrational wavepaket dynamics in 2 are closely correlated with the “H” framed geometry in 2 compared to the “A” frame in 1, leading to the sharp increase in - interaction between ppy ligands. Although the TA experiments do not directly reveal the ultrafast intersystem crossing (ISC) because of a strong coherent spike at early timescales, the dependence of the vibrational wavepacket dynamics on molecular geometry can be understood based on previously proposed potential energy surfaces (PES) as function of PtPt distance, suggesting that the interaction between the cyclometalating ligands can be a key factor in determining the PtPt shortening and the related energy relaxation dynamics in the Pt(II) dimers. Further experiments using femtosecond broadband TA and X-ray scattering spectroscopy are planned to investigate directly the ISC and PtPt contraction, to support the relationship between ground state molecular geometry and photoinduced structural changes in the Pt(II) dimers.

Kim, Pyosang↗

SU-C-BRA-04: Automated Segmentation of Head-And-Neck CT Images for Radiotherapy Treatment Planning Via Multi-Atlas Machine Learning (MAML)

Purpose: Accurate image segmentation is a crucial step during image guided radiation therapy. This work proposes multi-atlas machine learning (MAML) algorithm for automated segmentation of head-and-neck CT images. Methods: As the first step, the algorithm utilizes normalized mutual information as similarity metric, affine registration combined with multiresolution B-Spline registration, and then fuses together using the label fusion strategy via Plastimatch. As the second step, the following feature selection strategy is proposed to extract five feature components from reference or atlas images: intensity (I), distance map (D), box (B), center of gravity (C) and stable point (S). The box feature B is novel. It describes a relative position from each point to minimum inscribed rectangle of ROI. The center-of-gravity feature C is the 3D Euclidean distance from a sample point to the ROI center of gravity, and then S is the distance of the sample point to the landmarks. Then, we adopt random forest (RF) in Scikit-learn, a Python module integrating a wide range of state-of-the-art machine learning algorithms as classifier. Different feature and atlas strategies are used for different ROIs for improved performance, such as multi-atlas strategy with reference box for brainstem, and single-atlas strategy with reference landmark for optic chiasm. Results: The algorithm was validated on a set of 33 CT images with manual contours using a leave-one-out cross-validation strategy. Dice similarity coefficients between manual contours and automated contours were calculated: the proposed MAML method had an improvement from 0.79 to 0.83 for brainstem and 0.11 to 0.52 for optic chiasm with respect to multi-atlas segmentation method (MA). Conclusion: A MAML method has been proposed for automated segmentation of head-and-neck CT images with improved performance. It provides the comparable result in brainstem and the improved result in optic chiasm compared with MA. Xuhua Ren and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000), and the Shanghai Pujiang Talent Program (#14PJ1404500).

60 APPLIED LIFE SCIENCES↗

Onset of normal and inverse homoclinic bifurcation in a double plasma system near a plasma fireball

Plasma fireballs are generated due to a localized discharge and appear as a luminous glow with a sharp boundary, which suggests the presence of a localized electric field such as electrical sheath or double layer structure. The present work reports the observation of normal and inverse homoclinic bifurcation phenomena in plasma oscillations that are excited in the presence of fireball in a double plasma device. The controlling parameters for these observations are the ratio of target to source chamber (n{sub T}/n{sub S}) densities and applied electrode voltage. Homoclinic bifurcation is noticed in the plasma potential fluctuations as the system evolves from narrow to long time period oscillations and vice versa with the change of control parameter. The dynamical transition in plasma fireball is demonstrated by spectral analysis, recurrence quantification analysis (RQA), and statistical measures, viz., skewness and kurtosis. The increasing trend of normalized variance reflects that enhancing n{sub T}/n{sub S} induces irregularity in plasma dynamics. The exponential growth of the time period is strongly indicative of homoclinic bifurcation in the system. The gradual decrease of skewness and increase of kurtosis with the increase of n{sub T}/n{sub S} also reflect growing complexity in the system. The visual change of recurrence plot and gradual enhancement of RQA variables DET, L{sub max}, and ENT reflects the bifurcation behavior in the dynamics. The combination of RQA and spectral analysis is a clear evidence that homoclinic bifurcation occurs due to the presence of plasma fireball with different density ratios. However, inverse bifurcation takes place due to the change of fireball voltage. Some of the features observed in the experiment are consistent with a model that describes the dynamics of ionization instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Signatures of tunable superconductivity in a trilayer graphene moiré superlattice

Understanding the mechanism of high-transition-temperature (high-T c ) superconductivity is a central problem in condensed matter physics. It is often speculated that high-Tc superconductivity arises in a doped Mott insulator1 as described by the Hubbard model 2,3,4 . An exact solution of the Hubbard model, however, is extremely challenging owing to the strong electron–electron correlation in Mott insulators. Therefore, it is highly desirable to study a tunable Hubbard system, in which systematic investigations of the unconventional superconductivity and its evolution with the Hubbard parameters can deepen our understanding of the Hubbard model. Here we report signatures of tunable superconductivity in an ABC-trilayer graphene (TLG) and hexagonal boron nitride (hBN) moiré superlattice. Unlike in ‘magic angle’ twisted bilayer graphene, theoretical calculations show that under a vertical displacement field, the ABC-TLG/hBN heterostructure features an isolated flat valence miniband associated with a Hubbard model on a triangular superlattice 5,6 where the bandwidth can be tuned continuously with the vertical displacement field. Upon applying such a displacement field we find experimentally that the ABC-TLG/hBN superlattice displays Mott insulating states below 20 kelvin at one-quarter and one-half fillings of the states, corresponding to one and two holes per unit cell, respectively. Upon further cooling, signatures of superconductivity (‘domes’) emerge below 1 kelvin for the electron- and hole-doped sides of the one-quarter-filling Mott state. The electronic behaviour in the ABC-TLG/hBN superlattice is expected to depend sensitively on the interplay between the electron–electron interaction and the miniband bandwidth. By varying the vertical displacement field, we demonstrate transitions from the candidate superconductor to Mott insulator and metallic phases. Our study shows that ABC-TLG/hBN heterostructures offer attractive model systems in which to explore rich correlated behaviour emerging in the tunable triangular Hubbard model.

36 MATERIALS SCIENCE↗

Sharp detection of low-dimensional structure in probability measures via dimensional logarithmic Sobolev inequalities

Identifying low-dimensional structure in high-dimensional probability measures is an essential pre-processing step for efficient sampling. To identify this structure, we approximate the target measure as a perturbation of an arbitrary reference measure along a few directions in $\mathbb{R}^{d}$. These directions are determined by minimizing an upper bound on the Kullback–Leibler (KL) divergence between the target and its approximation. Our contribution improves upon previous works by leveraging dimensional logarithmic Sobolev inequalities to refine the bound on the KL divergence. These inequalities lead to a uniformly tighter bound on the KL divergence, thereby enhancing the identification of the most significant perturbation directions. In particular, when the target and reference are both Gaussian, minimizing the resulting bound is equivalent to minimizing the KL divergence. We further demonstrate the applicability of this analysis to the squared Hellinger distance, where analogous reasoning shows that the dimensional Poincaré inequality offers improved bounds.

Bayesian inference↗

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗

Tunable correlated Chern insulator and ferromagnetism in a moiré superlattice

Studies of two-dimensional electron systems in a strong magnetic field revealed the quantum Hall effect, a topological state of matter featuring a finite Chern number C and chiral edge states. Haldane later theorized that Chern insulators with integer quantum Hall effects could appear in lattice models with complex hopping parameters even at zero magnetic field. The ABC-trilayer graphene/hexagonal boron nitride (ABC-TLG/hBN) moiré superlattice provides an attractive platform with which to explore Chern insulators because it features nearly flat moiré minibands with a valley-dependent, electrically tunable Chern number. In this work, we report the experimental observation of a correlated Chern insulator in an ABC-TLG/hBN moiré superlattice. We show that reversing the direction of the applied vertical electric field switches the moiré minibands of ABC-TLG/hBN between zero and finite Chern numbers, as revealed by large changes in magneto-transport behaviour. For topological hole minibands tuned to have a finite Chern number, we focus on quarter filling, corresponding to one hole per moiré unit cell. The Hall resistance is well quantized at h/2e 2 (where h is Planck’s constant and e is the charge on the electron), which implies C = 2, for a magnetic field exceeding 0.4 tesla. The correlated Chern insulator is ferromagnetic, exhibiting substantial magnetic hysteresis and a large anomalous Hall signal at zero magnetic field. Finally, our discovery of a C = 2 Chern insulator at zero magnetic field should open up opportunities for discovering correlated topological states, possibly with topological excitations, in nearly flat and topologically nontrivial moiré minibands.

36 MATERIALS SCIENCE↗

Models for predicting objective function weights in prostate cancer IMRT

Purpose: To develop and evaluate the clinical applicability of advanced machine learning models that simultaneously predict multiple optimization objective function weights from patient geometry for intensity-modulated radiation therapy of prostate cancer. Methods: A previously developed inverse optimization method was applied retrospectively to determine optimal objective function weights for 315 treated patients. The authors used an overlap volume ratio (OV) of bladder and rectum for different PTV expansions and overlap volume histogram slopes (OVSR and OVSB for the rectum and bladder, respectively) as explanatory variables that quantify patient geometry. Using the optimal weights as ground truth, the authors trained and applied three prediction models: logistic regression (LR), multinomial logistic regression (MLR), and weighted K-nearest neighbor (KNN). The population average of the optimal objective function weights was also calculated. Results: The OV at 0.4 cm and OVSR at 0.1 cm features were found to be the most predictive of the weights. The authors observed comparable performance (i.e., no statistically significant difference) between LR, MLR, and KNN methodologies, with LR appearing to perform the best. All three machine learning models outperformed the population average by a statistically significant amount over a range of clinical metrics including bladder/rectum V53Gy, bladder/rectum V70Gy, and dose to the bladder, rectum, CTV, and PTV. When comparing the weights directly, the LR model predicted bladder and rectum weights that had, on average, a 73% and 74% relative improvement over the population average weights, respectively. The treatment plans resulting from the LR weights had, on average, a rectum V70Gy that was 35% closer to the clinical plan and a bladder V70Gy that was 29% closer, compared to the population average weights. Similar results were observed for all other clinical metrics. Conclusions: The authors demonstrated that the KNN and MLR weight prediction methodologies perform comparably to the LR model and can produce clinical quality treatment plans by simultaneously predicting multiple weights that capture trade-offs associated with sparing multiple OARs.

60 APPLIED LIFE SCIENCES↗

Reactor Design and Decommissioning - An Overview of International Activities in Post Fukushima Era1 - 12396

Accidents at the Fukushima Dai-ichi reactors as a result of the devastating earthquake and tsunami of March 11, 2011 have not only dampened the nuclear renaissance but have also initiated a re-examination of the design and safety features for the existing and planned nuclear reactors. Even though failures of some of the key site features at Fukushima can be attributed to events that in the past would have been considered as beyond the design basis, the industry as well as the regulatory authorities are analyzing what features, especially passive features, should be designed into the new reactor designs to minimize the potential for catastrophic failures. It is also recognized that since the design of the Fukushima BWR reactors which were commissioned in 1971, many advanced safety features are now a part of the newer reactor designs. As the recovery efforts at the Fukushima site are still underway, decisions with respect to the dismantlement and decommissioning of the damaged reactors and structures have not yet been finalized. As it was with Three Mile Island, it could take several decades for dismantlement, decommissioning and clean up, and the project poses especially tough challenges. Near-term assessments have been issued by several organizations, including the IAEA, the USNRC and others. Results of such investigations will lead to additional improvements in system and site design measures including strengthening of the anti-tsunami defenses, more defense-in-depth features in reactor design, and better response planning and preparation involving reactor sites. The question also arises what would the effect be on the decommissioning scene worldwide, and what would the effect be on the new reactors when they are eventually retired and dismantled. This paper provides an overview of the US and international activities related to recovery and decommissioning including the decommissioning features in the reactor design process and examines these from a new perspective in the post Fukushima -accident era. Accidents at the Fukushima Daiichi reactors in the aftermath of the devastating earthquake and tsunami of March 11, 2011 have slowed down the nuclear renaissance world-wide and may have accelerated decommissioning either because some countries have decided to halt or reduce nuclear, or because the new safety requirements may reduce life-time extensions. Even in countries such as the UK and France that favor nuclear energy production existing nuclear sites are more likely to be chosen as sites for future NPPs. Even as the site recovery efforts continue at Fukushima and any decommissioning decisions are farther into the future, the accidents have focused attention on the reactor designs in general and specifically on the Fukushima type BWRs. The regulatory authorities in many countries have initiated a re-examination of the design of the systems, structures and components and considerations of the capability of the station to cope with beyond-design basis events. Enhancements to SSCs and site features for the existing reactors and the reactors that will be built will also impact the decommissioning phase activities. The newer reactor designs of today not only have enhanced safety features but also take into consideration the features that will facilitate future decommissioning. Lessons learned from past management and operation of reactors as well as the lessons from decommissioning are incorporated into the new designs. However, in the post-Fukushima era, the emphasis on beyond-design-basis capability may lead to significant changes in SSCs, which eventually will also have impact on the decommissioning phase. Additionally, where some countries decide to phase out the nuclear power, many reactors may enter the decommissioning phase in the coming decade. While the formal updating and expanding of existing guidance documents for accident cleanup and decommissioning would benefit by waiting until the Fukushima project has progressed sufficiently for that experience to be reliably interpreted, the development of structured on-line sharing of information and especially the creation of an on-line compendium of methods, tools, and techniques by which damaged fuel and other unique situations have been addressed can be addressed sooner and maintained as new problems and solutions arise and are resolved. The IAEA's new 'WEB 2.0 tool' CONNECT is expected to play a significant role in this and related information-sharing activities. The trend in some countries such as the United States has been to re-license the existing reactors for additional twenty years, beyond the original design life. Given the advances in technology over the past four decades, and considering that the newer designs incorporate significant improvements in safety systems, it may not be economical or technically feasible to retrofit enhancements into some of the older reactors. In such cases, the reactors may be retired from service and decommissioned. Overall, the energy demand in the world continues to rise, with sharp increases in the Asian countries, and nuclear power's role in the world's energy supply is expected to continue. Events at Fukushima have led to a re-examination on many fronts, including reactor design and regulatory requirements. Further changes may occur in these areas in the post-Fukushima era. These changes in turn will also impact the world-wide decommissioning scene and the decommissioning phase of the future reactors. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SU-E-T-616: Comparison of Plan Dose Accuracy for Anterior Vs. Lateral Fields in Proton Therapy of Prostate Cancer

Purpose: With the anticipated introduction of in vivo range verification methods, the use of anterior fields for proton therapy of prostate cancer may become an attractive treatment option, and improve upon the dose distributions achievable with conventional lateral-opposed fields. This study aimed to evaluate and compare the planned dose accuracy for lateral versus anterior oblique field arrangements. Methods: Four patients with low/intermediate risk prostate cancer, participating in a clinical trial at our institution, were selected for this study. All patients were treated using lateral-opposed fields (LAT). The clinical target volume (CTV) received a total dose of 79.2 Gy in 44 fractions. Anterior oblique research plans (ANT) were created using the clinical planning system, and featured beams with ±35-degree gantry angle, 1.2 cm aperture margins, 3-mm range compensator smearing and no range uncertainty margins. Monte Carlo (MC) simulations were performed for both beam arrangements using TOPAS. Dose volume histograms were analyzed and compared for planned and MC dose distributions. Differences between MC and planned DVH parameters were computed as a percentage of the total prescribed dose. Results: For all patients, CTV dose was systematically lower (∼2–2.5%) for MC than the plan. This discrepancy was slightly larger (∼0.5%) for LAT compared to ANT plans for all cases. Although the dose differences for bladder and anterior rectal wall remained within 0.7% for all LAT cases, they were slightly larger for ANT plans, especially for case 3 due to larger patient size and MC-plan range difference. The EUD difference for femoral heads was within 0.6% for both LAT and ANT cases. Conclusion: The dose calculated by the treatment planning system using pencil beam algorithm agrees with MC to within 2.5% and is comparable for lateral and anterior scenarios. The dose agreement in the anterior rectal wall is range- and hence, patient-dependent for ANT treatments.

60 APPLIED LIFE SCIENCES↗

Polymorphous band structure model of gapping in the antiferromagnetic and paramagnetic phases of the Mott insulators MnO, FeO, CoO, and NiO

The existence of band gaps in both the antiferromagnetic (AFM) and paramagnetic (PM) phases of the classic NaCl-structure Mott insulators MnO, FeO, CoO, and NiO is traditionally viewed and taught as a manifestation of strong correlation whereby insulation results from electrons moving across the lattice forming states with doubly occupied d orbitals on certain atomic sites and empty d orbitals on other sites. Within such theories, the gap of the AFM and PM phases of these oxides emerges even in the absence of spatial symmetry breaking. The need for such a correlated picture is partially based on the known failure of the commonly used band models for the PM phase that assume for such a spin disordered state the macroscopically averaged NaCl structure, where all transition metal (TM) sites are symmetry-equivalent (a monomorphous description), producing a gapless PM state with zero magnetic moments, in sharp conflict with experiment. Here, we seek to understand the minimum theoretical description needed to capture the leading descriptors of ground state Mott insulation in the classic, 3d monoxide Mott systems—gapping and moment formation in the AFM and PM phase. As noted by previous authors, the spin-ordered AFM phase in these materials already shows in band theory a significant band gap when one doubles the NaCl unit cell by permitting different potentials for transition-metal atoms with different spins. For the spin-disordered PM phase, we allow analogously larger NaCl-type supercells where each TM site can have different spin direction and local bonding environments (i.e., disordered), yet the total spin is zero. Such a polymorphous description has the flexibility to acquire symmetry-breaking energy-lowering patterns that can lift the degeneracy of the d orbitals and develop large on-site magnetic moments without violating the global, averaged NaCl symmetry. Electrons are exchanged between spin-up and spin-down bands to create closed-shell insulating configurations that lend themselves to a single determinental description. It turns out that such a polymorphous description of the structure within the single-determinant, mean-field, Bloch periodic band structure approach (based on DFT+U) allows large on-site magnetic moments to develop spontaneously, leading to significant (1–3 eV) band gaps and large local moments in the AFM and PM phases of the classic NaCl-structure Mott insulators MnO, FeO, CoO, and NiO in agreement with experiment. We adapt to the spin disordered polymorphous configurations the “special quasirandom structure” (SQS) construct known from the theory of disordered substitutional alloys whereby supercell approximants which represent the best random configuration average (not individual snapshots) for finite (64, 216 atoms, or larger) supercells of a given lattice symmetry are constructed. We conclude that the basic features of these paradigmatic Mott insulators can be approximated by the physics included in energy-lowering symmetry broken DFT.

36 MATERIALS SCIENCE↗

Mixing with applications to inertial-confinement-fusion implosions

Approximate one-dimensional (1D) as well as 2D and 3D simulations are playing an important supporting role in the design and analysis of future experiments at National Ignition Facility. This work is mainly concerned with 1D simulations, used extensively in design and optimization. We couple a 1D buoyancy-drag mix model for the mixing zone edges with a 1D inertial confinement fusion simulation code. This analysis predicts that National Ignition Campaign (NIC) designs are located close to a performance cliff, so modeling errors, design features (fill tube and tent) and additional, unmodeled instabilities could lead to significant levels of mix. The performance cliff we identify is associated with multimode plastic ablator (CH) mix into the hot-spot deuterium and tritium (DT). The buoyancy-drag mix model is mode number independent and selects implicitly a range of maximum growth modes. Our main conclusion is that single effect instabilities are predicted not to lead to hot-spot mix, while combined mode mixing effects are predicted to affect hot-spot thermodynamics and possibly hot-spot mix. Combined with the stagnation Rayleigh-Taylor instability, we find the potential for mix effects in combination with the ice-to-gas DT boundary, numerical effects of Eulerian species CH concentration diffusion, and ablation-driven instabilities. With the help of a convenient package of plasma transport parameters developed here, we give an approximate determination of these quantities in the regime relevant to the NIC experiments, while ruling out a variety of mix possibilities. Plasma transport parameters affect the 1D buoyancy-drag mix model primarily through its phenomenological drag coefficient as well as the 1D hydro model to which the buoyancy-drag equation is coupled.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparing Velocity Profiles Along the Rod Length of a Helical Coil Steam Generator Model

A model of a specific type of tube and shell heat exchanger, the helical coil steam generator (HCSG), was designed and constructed to study multiple fluidic properties that attribute to the efficiency of the heat transfer it provides. Tube and shell heat exchangers are currently used throughout industry in chemical, power and petroleum engineering. Most commonly, a primary flow is circulated through a shell, while the secondary flow is run through several small tubes that are bundled together. Methods to increase the heat transfer between the two flows include adding baffles to force shell side flow over more tubes, fins added in-between layers of tubes and passing tubes multiple times through the shell to increase heat exchanger effectiveness. The HCSG incorporates many of the proposed methods but also relies on a complex geometry to increase the heat transfer between tube and shell side flow. In an area where cost effectiveness is favorable, such as inside a nuclear reactor, helically banked steam generators makes extensive use of the room with no sharp tube bends. This increases the pressure loss on the secondary flow and can solve potential temperature manufacturing issues at the welds of junctions. Additionally, they have exhibited pure cross flow behaviors with high heat transfer coefficients and axial temperature profiles that correspond to pure counter-flow with high mean temperature differences. Early work showed that helical banks can have a significant increase of 16 to 43% higher heat transfer coefficient than straight pipe heat exchangers. While many previous studies have looked at the heat transfer, pressure drop and flow fields of a tube and shell heat exchanger, they all focus on tube bundles that are either in-line or staggered. The early flow visualization studies on these types of bundle arrangements reported that attached eddies were formed alternatively behind each of the rods and were occasionally dragged into the main stream of shell side flow. Later studies suggested that observable periodic vortex shedding was dependent on geometry. Experimental data taken using techniques such as laser Doppler anemometry (LDV) by Simonin and Barcouda (1986) as mentioned in Hassan (2002), hot-wire anemometry with Tsuchida (1982) and particle image velocimetry (PIV) conducted by researchers as late as 2007 have been used to characterize turbulent flow and distinguish flow field structures within in-line and staggered tube bundles. By studying these flow field structures, their contribution to fluidic properties such as heat transfer efficiency and pressure drop can be characterized effectively. Literature review lead to a design that incorporated five layers of coils wrapped around one another. The outermost and middle layers of tubes coiled at a specific pitch in a clockwise direction while the other two layers coiled counterclockwise. The proposed design incorporates a constantly changing lateral pitch ratio, a, between adjacent rods. The constructed test section isolates a single interface between two adjacent coiling rod bundles of twelve rods. Three distinct cross sectional planes along the rod bundles were looked at to study the effect of the local geometry to flow fields, Section A, B and C. Particle Image Velocimetry (PIV) frames of 896 x 848 pixels were captured at each section across a specific height. Profiles of mean velocity fields for each section are compared and discussed as possible features of the changing lateral pitch ratio.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY2012 Report on Fast Reactor Toolset Work (PROTEUS-Fast)

As part of the NEAMS activity in DOE, a fast reactor simulation program was launched in April 2007, termed SHARP. The initial goal of SHARP was to develop a suite of modern simulation tools specifically for the analysis and design of sodium cooled fast reactors, and this goal has since been extended to consider all reactor types of interest. As part of that work, a high‐fidelity deterministic neutron transport code named UNIC was started, providing a common framework for fast reactor neutronics tools. Focusing on UNIC, spatial de‐homogenization increased the typical neutronics problem with a few hundred million degrees of freedom (DOF) into ones with hundreds of trillions of DOF. Such problems are beyond the capabilities of existing supercomputers. With the advent of future computer hardware, these high‐fidelity problems can be solved but only after a multiyear development program to adapt codes such as UNIC to that technology. UNIC can presently be used on smaller problems (single assembly, multiple assemblies) with today’s hardware, but it is also necessary to provide tools to solve the larger problems of interest, albeit without the aforementioned level of spatial de‐homogenization. To address the desire to demonstrate the coupling capability of NEAMS, part of the FY12 funding was focused on preparing the legacy codes, termed PROTEUS‐Fast, which are fully capable of solving full‐scale fast reactor problems. A bulk of the work this year was focused on preparing and distributing the PROTEUS‐Fast codes for use by all NEAMS users. The demonstration of a full coupling capability will require more follow on work, but significant progress was made this year preparing the interface coupling routines and preparing PROTEUS‐Fast for use. The primary targets of interest in PROTEUS‐Fast are DIF3D, DIF3D‐K (kinetics), and VARI3D (reactivity coefficients), with REBUS (fuel cycle) being an auxiliary component updated for convenience. All codes in PROTEUS‐Fast were verified and added to the BuildBot regression test. Unfortunately, all of the older codes suffered significant problems from the deprecated coding practices used within them and resulting conflicts with modern compilers. After some work, all of these issues were resolved and we feel confident in releasing them for use by others NEAMS members. With regard to specific upgrades, DIF3D‐VARIANT was upgraded with a few months of effort to allow the large scale problems of interest to be executed. After significant study, we choose not to modify VARI3D to include a transport option, but to develop a new code, PERSENT (partially funded by an external project) which uses the DIF3D‐VARIANT transport solution option. The outstanding features missing from MC 2 ‐3 were added, and MC 2 ‐3 has been released to RSICC. All of the remaining PROTEUS‐Fast codes are either finished with the copyright and licensing procedure or are actively going through it. While that process will delay the targeted August 2012 release to RSICC for a few more months, the codes are ready for immediate use and can be connected to the framework either by us or by others as the interface readers and documentation will be released with the PROTEUS‐Fast codes.

Wolters, E. R.↗

Implications of high-pressure oxygen hydrates on radiolytic oxygen in Jovian icy moons

Various icy moons, such as Europa and Ganymede, have thin oxygen atmospheres and exhibit spectral features attributed to oxygen held in their surface ices. The oxygen forms from the radiolysis of water. The interiors of these bodies are subject to high pressures and it is not known how deep into icy moons oxygen-bearing ices can penetrate, or the structures formed by the oxygen–water system at high pressure. Here, we show that oxygen hydrates are stable to 2.6 GPa, allowing them to penetrate deep into icy moons, both above and below proposed sub-surface liquid-water oceans. Similarities between oxygen and hydrogen hydrates indicate potentially enhanced recombination rates, transforming them back into water and offering a resolution to the discrepancy between predicted and measured radiolysis rates. In addition to the low-pressure CS-II clathrate, our results find three high-pressure phases in the oxygen–water system: an ST clathrate, a C 0 hydrate, and a filled ice isomorphous with methane hydrate III. This shows a vast storage potential for molecular oxygen in icy moons and indicates that Europa could still be absorbing oxygen into its crustal ice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The HITRAN2016 molecular spectroscopic database

This paper describes the contents of the 2016 edition of the HITRAN molecular spectroscopic compilation. The new edition replaces the previous HITRAN edition of 2012 and its updates during the intervening years. The HITRAN molecular absorption compilation is comprised of five major components: the traditional line-by-line spectroscopic parameters required for high-resolution radiative-transfer codes, infrared absorption cross-sections for molecules not yet amenable to representation in a line-by-line form, collision-induced absorption data, aerosol indices of refraction, and general tables such as partition sums that apply globally to the data. The new HITRAN is greatly extended in terms of accuracy, spectral coverage, additional absorption phenomena, added line-shape formalisms, and validity. Moreover, molecules, isotopologues, and perturbing gases have been added that address the issues of atmospheres beyond the Earth. Of considerable note, experimental IR cross-sections for almost 200 additional significant molecules have been added to the database.

Gordon, I. E.↗

Tropospheric Ozone Assessment Report: Database and Metrics Data of Global Surface Ozone Observations

In support of the first Tropospheric Ozone Assessment Report (TOAR) a relational database of global surface ozone observations has been developed and populated with hourly measurement data and enhanced metadata. A comprehensive suite of ozone data products including standard statistics, health and vegetation impact metrics, and trend information, are made available through a common data portal and a web interface. These data form the basis of the TOAR analyses focusing on human health, vegetation, and climate relevant ozone issues, which are part of this special feature. Cooperation among many data centers and individual researchers worldwide made it possible to build the world's largest collection of in-situ hourly surface ozone data covering the period from 1970 to 2015. By combining the data from almost 10,000 measurement sites around the world with global metadata information, new analyses of surface ozone have become possible, such as the first globally consistent characterisations of measurement sites as either urban or rural/remote. Exploitation of these global metadata allows for new insights into the global distribution, and seasonal and long-term changes of tropospheric ozone and they enable TOAR to perform the first, globally consistent analysis of present-day ozone concentrations and recent ozone changes with relevance to health, agriculture, and climate. Considerable effort was made to harmonize and synthesize data formats and metadata information from various networks and individual data submissions. Extensive quality control was applied to identify questionable and erroneous data, including changes in apparent instrument offsets or calibrations. Such data were excluded from TOAR data products. Limitations of a posteriori data quality assurance are discussed. As a result of the work presented here, global coverage of surface ozone data for scientific analysis has been significantly extended. Yet, large gaps remain in the surface observation network both in terms of regions without monitoring, and in terms of regions that have monitoring programs but no public access to the data archive. Therefore future improvements to the database will require not only improved data harmonization, but also expanded data sharing and increased monitoring in data-sparse regions.

54 ENVIRONMENTAL SCIENCES↗