Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal Distribution”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Analysis of Completion Design Impact on Cluster Efficiency and Pressure-Based Well Communication in HFTS-2 Delaware Basin

The Hydraulic Fracturing Test Site 2 (HFTS-2) is a joint industry project in the Delaware basin to advance hydraulic fracturing understanding and improve productivity in shale reservoirs. The project integrates multi-disciplinary approaches to evaluate different completion designs, well spacing, inter-well communication, and stimulated rock volume, among other factors. This paper focused on two major areas related to hydraulic fracture performance. First, an analysis of different completion designs on cluster efficiency based on near well Distributed Acoustic Sensing (DAS). Second, an evaluation of well and completion designs on inter-well fracture driven interactions (FDIs) based on downhole pressure monitoring. Fluid/sand distribution and cluster efficiency analyses were based on near wellbore DAS data collected from two adjacent horizontal wells completed in two different landing zones in the Wolfcamp formation. These wells had different completion designs aiming to evaluate the effect of normal vs. extended stage lengths, perforation hole tapering and limited entry. Standard deviation from ideal fluid/sand distribution and waterfall plots were used to evaluate cluster efficiency for each design and stage. Inter-well FDIs analysis was conducted among the horizontal wells and a vertical monitor well. One horizontal well served as the monitor well while the other horizontal well was being treated. The vertical well was instrumented with downhole pressure and temperature gauges to aid monitoring fracture height growth. The pressure response during and after fracturing was characterized based on maximum pressure increase value and slope. Pressure response vs. FDIs trigger factors such distance, cluster efficiency and stage fluid volume were also analyzed. Based on the different completion and perforation designs tested, DAS analysis suggests that limited entry design worked best. Extreme limited entry showed the potential of high perforation erosion and reduced cluster efficiency. Here, the limited entry and tapered perforation design demonstrated potential to improve the cluster efficiency for extended stage lengths. Pressure monitoring across formation units proved to be critical to understand fracture interactions and fracture vertical growth. Pressure communication across different formation units during hydraulic fracturing operation indicate fractures grew upwards during Wolfcamp wells fracturing. However, this pressure communication dissipated over time. High intensity FDIs were recorded when the frac stages Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S018R001/2477607/urtec-2021-5289-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5289 were closer to the pressure gauge location in the monitor wells. Some of these stages that produced high intensity FDIs also had high fluid volume per cluster and low cluster efficiency. The multi-disciplinary and high-quality data collected from HFTS-2 helped to further understand why completion approaches such as limited entry and tapered perforation design are successful in improving cluster efficiency. The DAS data combined with downhole high-resolution pressure measurements also helped to quantify the effect of lower cluster efficiency data on the incidence and intensity of FDIs.

58 GEOSCIENCES↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-zirconium (Nb-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2021

Traveling Wave Relays for Distribution Feeder Protection with High Penetrations of Distributed Energy Resources

Increased penetration of power electronics interfaced Distributed Energy Resources (DER) like PV, electric vehicle and battery storage in the distribution system (both in numbers and sizes) may cause bi-directional power flows under normal operation. In addition, they contribute low fault current under short circuit conditions. These two changes can impact the proper operation of legacy protection systems. This can add additional challenges during weak grid and islanded operation. In this research, traveling wave (high frequency) signature based protection scheme is proposed for future distribution systems. Simulations are performed under different transient scenarios to provide insight into the high frequency signatures to show the traveling wave behavior. Finally, this paper discusses the importance of sensing and digital processing requirements for traveling-wave protection in distribution system.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

Traveling Wave Relays for Distribution Feeder Protection with High Penetrations of Distributed Energy Resources: Preprint

Increased penetration of power electronics interfaced Distributed Energy Resources (DER) like PV, electric vehicle and battery storage in the distribution system (both in numbers and sizes) may cause bi-directional power flows under normal operation. In addition, they contribute low fault current under short circuit conditions. These two changes can impact the proper operation of legacy protection systems. This can add additional challenges during weak grid and islanded operation. In this research, traveling wave (high frequency) signature based protection scheme is proposed for future distribution systems. Simulations are performed under different transient scenarios to provide insight into the high frequency signatures to show the traveling wave behavior. Finally, this paper discusses the importance of sensing and digital processing requirements for traveling-wave protection in distribution system.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

Traveling Wave Relays for Distribution Feeder Protection with High Penetrations of Distributed Energy Resources

Increased penetration of power electronics interfaced Distributed Energy Resources (DER) like PV, electric vehicle and battery storage in the distribution system (both in numbers and sizes) may cause bi-directional power flows under normal operation. In addition, they contribute low fault current under short circuit conditions. These two changes can impact the proper operation of legacy protection systems. This can add additional challenges during weak grid and islanded operation. In this research, traveling wave (high frequency) signature based protection scheme is proposed for future distribution systems. Simulations are performed under different transient scenarios to provide insight into the high frequency signatures to show the traveling wave behavior. Finally, this paper discusses the importance of sensing and digital processing requirements for traveling-wave protection in distribution system.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2016

MO-G-201-04: Knowledge-Based Planning for Single-Isocenter Stereotactic Radiosurgery to Multiple Brain Metastases

Purpose: Single-isocenter, linac-based SRS for multiple brain metastases (multi-mets) can deliver highly conformal radiation doses and reduce overall patient treatment time compared to other therapy techniques. This study aims to quantify the dosimetric benefits of knowledge-based planning (KBP) for multi-met treatments. Methods: Using a previously-published KBP methodology (an artificial neural network (ANN) trained on single-target linac-based SRS plans), 3D dose distribution predictions for multi-met patients were obtained by treating each brain lesion as a solitary target and subsequently combining individual predictions into a single distribution using a dose-weighted geometric averaging to obtain the best results in the inter-target space. 17 previously-treated multi-met plans, with target numbers ranging from N=2–5, were used to validate the ANN predictions and subsequent KBP auto-planning routine. The fully-deliverable KBP plans were developed by converting dose distribution predictions into patient-specific optimization objectives while maintaining identical target normalizations (typically PTV V100%=D98%). Plan quality improvements were quantified by the difference between SRS quality metrics (QMs): δdQM=QM(clinical)-QM(KBP). QMs of interest were: gradient measure (GM), conformity index (CI), brain V10 and V5, brainstem D0.1cc and heterogeneity index (HI). Finally, overall plan quality was judged via blinded plan comparison by SRS-specializing physicians. Results: Two clinical plans were found to be significant outliers wherein plan quality was dramatically worse than KBP. Despite indicating KBP superiority, these were removed from the QM analysis to prevent skewing the results. In the remaining cases, clinical and KBP QMs were nearly identical with modest improvements in the KBP sample: δGM=0.12±0.56mm, δCI=−0.01±0.04, Brain δV10=0.8±2.6cc, brain δV5=6.3 ±10.7cc, brainstem δD0.1cc=0.06±1.19Gy and δHI= −0.04±0.05. Ultimately, 13/17 KBP plans were deemed superior to the manual plans in blinded physician review. Conclusion: The results demonstrate that KBP-driven automated planning in linac-based single-isocenter treatments for multiple brain metastases is indistinguishable from, or even better than, traditional manual planning. J. Hattangadi: Research Grant; Varian Medical Systems; K.L. Moore: Research Grant: Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Effect of choosing the Q{sub α}-values and daughter density distributions on the magic numbers predicted by α decays

We investigated the effects of choosing the α-decay energy (Q{sub α}) and particular density distributions of the daughter nuclei on predicting unknown magic numbers, or confirming some of them. The Q{sub α} values were calculated using three models, namely the finite-range droplet model (FRDM), the Weizsäcker–Skyrme (WS4) mass model, and a recent fitting formula. The α-decay half-lives (T{sub α}) of superheavy nuclei of Z = 118–124 have been calculated within the generalized density-dependent cluster model. The double-folding α+core potential was determined based on the M3Y-Reid nucleon–nucleon interaction. We have used the two-parameter Fermi distribution to represent the nuclear density, with four sets of parameters that were determined based on electron scattering data and different microscopic calculations. Upon investigating the relative stability of the mentioned superheavy isotopic chains, we found that the Q{sub α}-values deduced from FRDM and WS4 models predict different neutron magic numbers. The examined fitting formula failed to predict any magic number larger than N = 162. Under normalization to the proton and neutron (N) numbers, the different choices of the density distributions showed almost no effect on the general behavior of T{sub α} with N, but changed its values up to one order of magnitude. Even though the maxima and minima of T{sub α} clearly follow the characteristic minima and maxima of Q{sub α}, respectively, the results showed that precise density distribution and deformations, in addition to a correct preformation probability, are necessarily required for estimating an accurate value of the half-life.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Robust Matrix Completion State Estimation in Distribution Systems

Due to the insufficient measurements in the distribution system state estimation (DSSE), full observability and redundant measurements are difficult to achieve without using the pseudo measurements. The matrix completion state estimation (MCSE) combines the matrix completion and power system model to estimate voltage by exploring the low-rank characteristics of the matrix. This paper proposes a robust matrix completion state estimation (RMCSE) to estimate the voltage in a distribution system under a low-observability condition. Tradition state estimation weighted least squares (WLS) method requires full observability to calculate the states and needs redundant measurements to proceed a bad data detection. The proposed method improves the robustness of the MCSE to bad data by minimizing the rank of the matrix and measurements residual with different weights. It can estimate the system state in a low-observability system and has robust estimates without the bad data detection process in the face of multiple bad data. The method is numerically evaluated on the IEEE 33-node radial distribution system. The estimation performance and robustness of RMCSE are compared with the WLS with the largest normalized residual bad data identification (WLS-LNR), and the MCSE.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Robust Matrix Completion State Estimation in Distribution Systems: Preprint

Due to the insufficient measurements in the distribution system state estimation (DSSE), full observability and redundant measurements are difficult to achieve without using the pseudo measurements. The matrix completion state estimation (MCSE) combines the matrix completion and power system model to estimate voltage by exploring the low-rank characteristics of the matrix. This paper proposes a robust matrix completion state estimation (RMCSE) to estimate the voltage in a distribution system under a low-observability condition. Tradition state estimation weighted least squares (WLS) method requires full observability to calculate the states and needs redundant measurements to proceed a bad data detection. The proposed method improves the robustness of the MCSE to bad data by minimizing the rank of the matrix and measurements residual with different weights. It can estimate the system state in a low-observability system and has robust estimates without the bad data detection process in the face of multiple bad data. The method is numerically evaluated on the IEEE 33-node radial distribution system. The estimation performance and robustness of RMCSE are compared with the WLS with the largest normalized residual bad data identification (WLS-LNR), and the MCSE.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Numerical Simulations of Particle Acceleration at Interplanetary Quasi-perpendicular Shocks

Using test particle simulations we study particle acceleration at highly perpendicular (θ{sub Bn}⩾75{sup ∘}) shocks under conditions of modeling magnetic turbulence. We adopt a backward-in-time method to solve the Newton–Lorentz equation using the observed shock parameters for quasi-perpendicular interplanetary shocks, and compare the simulation results with ACE/EPAM observations to obtain the injection energy and timescale of particle acceleration. With our modeling and observations, we find that a large upstream speed is responsible for efficient particle acceleration. Our results also show that the quasi-perpendicular shocks are capable of accelerating thermal particles to high energies of the order of MeV for both kappa and Maxwellian upstream distributions, which may originate from the fact that in our model, the local background magnetic field has a component parallel to the shock normal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-E-J-56: Dosimetric Impact of Patient Roll Variability in Peripheral-Lung SBRT

Purpose: Positional variation about the patient longitudinal axis may introduce deleterious effects to clinical treatment plans, the effects of which are exacerbated by hypo-fractionated treatment regimens. Our goal is to evaluate the dosimetric impact of rotational setup error in stereotactic body radiation therapy (SBRT) for peripherally-located lung lesions. Methods: Six patients who had been treated for peripheral lung lesions using varying SBRT techniques were analyzed. Each patient's planning CT was processed using custom in-house software to simulate 1, 3, and 5 degree roll in both directions about their longitudinal axis. The original plan was fused and re-calculated. Agreement between planned and error-simulated dose distributions was quantified using 3%, 3mm criteria for the 3-dimensional gamma metric, and changes in target coverage and adjacent normal anatomy were also quantified. Results: Increased rotational setup error resulted in reduced target coverage, with changes in equivalent uniform dose ranging from −1.24% to −6.11% for simulated roll greater than or equal to 3 degrees, increasing with magnitude of error and skewed for roll directed toward the involved side. Changes to normal tissue metrics are linearly related to simulated roll, with rib D(1cc) changes ranging from −4.23% to 3.27% relative to the original treatment plan. Error-simulated dose distributions demonstrate agreement by the gamma metric in the planned target of 72.56%, 81.02%, and 99.52% for 5, 3, and 1 degree simulated roll scenarios, respectively. Corresponding agreement for the regional dose in excess of 50% of the maximum were 65.71%,70.77%, and 97.88%. Conclusion: This study suggests that the ability to detect and compensate for rotational positioning variations beyond a threshold level may improve care for patients with peripherally-located lesions beyond that which may be accomplished using only translational axes. Positional error of one degree of roll had little influence, though three or more degrees of error introduced significant discrepancy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Rare earth element components in atmospheric particulates in the Bayan Obo mine region

The Bayan Obo mine, located in Inner Mongolia, China, is the largest light rare earth body ever found in the world. The research for rare earth elements (REEs) enrichment in atmospheric particulates caused by mining and ore processing is fairly limited so far. In this paper, atmospheric particulates including total suspended particulate (TSP) matter and particles with an equivalent aerodynamic diameter less than 10 μm (PM{sub 10}) were collected around the Bayan Obo mine region, in August 2012 and March 2013, to analyze the levels and distributions of REEs in particles. The total concentrations of REEs for TSP were 149.8 and 239.6 ng/m{sup 3}, and those for PM{sub 10} were 42.8 and 68.9 ng/m{sup 3}, in August 2012 and March 2013, respectively. Enrichment factor was calculated for all 14 REEs in the TSP and PM{sub 10} and the results indicated that REEs enrichment in atmosphere particulates was caused by anthropogenic sources and influenced by the strong wind in springtime. The spatial distribution of REEs in TSP showed a strong gradient concentration in the prevailing wind direction. REE chondrite normalized patterns of TSP and PM{sub 10} were similar and the normalized curves inclined to the right side, showing the conspicuous fractionation between the light REEs and heavy REE, which supported by the chondrite normalized concentration ratios calculated for selected elements (La{sub N}/Yb{sub N}, La{sub N}/Sm{sub N}, Gd{sub N}/Yb{sub N}). - Highlights: • TSP and PM{sub 10} samples were collected to analyze the levels and distributions of REE. • Enrichment factors indicated that REE enrichment was caused by anthropogenic sources. • The distribution of REEs showed a strong gradient in the prevailing wind direction. • Obvious fractionation between LREEs and HREEs is observed in atmospheric particulates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Measurement of Top Quark Polarization in $t\bar{t}$ Lepton + Jets Final States at D0

We present a measurement of top quark polarization in $t\bar{t}$ pair production in $p\bar{p}$ collisions at $\sqrt{s}=1.96$ TeV using data corresponding to 9.7 fb$^{-1}$ of integrated luminosity recorded with the D0 detector at the Fermilab Tevatron Collider. We consider final states containing a lepton and at least three jets. The polarization is measured through the distribution of lepton angles along three axes: the beam axis, the helicity axis, and the transverse axis normal to the $t\bar{t}$ production plane. This is the first measurement of top quark polarization at Tevatron using lepton+jet final states and the first measurement of the transverse polarization in $t\bar{t}$ production. The observed distributions are consistent with standard model predictions of nearly no polarization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Postoperative radiotherapy following mastectomy for patients with left-sided breast cancer: A comparative dosimetric study

The purposes of this article were to compare the biophysical dosimetry for postmastectomy left-sided breast cancer using 4 different radiotherapy (RT) techniques. In total, 30 patients with left-sided breast cancer were randomly selected for this treatment planning study. They were planned using 4 RT techniques, including the following: (1) 3-dimensional conventional tangential fields (TFs), (2) tangential intensity-modulated therapy (T-IMRT), (3) 4 fields IMRT (4F-IMRT), and (4) single arc volumetric-modulated arc therapy (S-VMAT). The planning target volume (PTV) dose was prescribed 50 Gy, the comparison of target dose distribution, conformity index, homogeneity index, dose to organs at risk (OARs), tumor control probability (TCP), normal tissue complication probability (NTCP), and number of monitor units (MUs) between 4 plans were investigated for their biophysical dosimetric difference. The target conformity and homogeneity of S-VMAT were better than the other 3 kinds of plans, but increased the volume of OARs receiving low dose (V{sub 5}). TCP of PTV and NTCP of the left lung showed no statistically significant difference in 4 plans. 4F-IMRT plan was superior in terms of target coverage and protection of OARs and demonstrated significant advantages in decreasing the NTCP of heart by 0.07, 0.03, and 0.05 compared with TFs, T-IMRT, and S-VMAT plan. Compared with other 3 plans, TFs reduced the average number of MUs. Of the 4 techniques studied, this analysis supports 4F-IMRT as the most appropriate balance of target coverage and normal tissue sparing.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

GALAXY HALO TRUNCATION AND GIANT ARC SURFACE BRIGHTNESS RECONSTRUCTION IN THE CLUSTER MACSJ1206.2-0847

In this work, we analyze the mass distribution of MACSJ1206.2-0847, particularly focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has a significant amount of intracluster light that is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at z = 1.033 into a giant arc and its counterimage. The multiple image positions and the surface brightness (SFB) distribution of the arc, which is bent around several cluster members, are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a Navarro-Frenk-White profile and the galaxy halos with two parameters for the mass normalization and the extent of a reference halo assuming scalings with their observed near-infrared light. We match the multiple image positions at an rms level of 0.''85 and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii-a scale in between those accessible to dynamical and field strong-lensing mass estimates on the one hand and galaxy-galaxy weak-lensing results on the other hand. The velocity dispersion and halo size of a galaxy with m{sub 160W,AB} = 19.2 and M{sub B,Vega} = -20.7 are {sigma} = 150 km s{sup -1} and r Almost-Equal-To 26 {+-} 6 kpc, respectively, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source, which is smaller by a factor of {approx}5.8 in area, demonstrating the increase in morphological information due to lensing. We conclude that this galaxy likely has star-forming spiral arms with a red (older) central component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Fully Ceramic Microencapsulated fuel in prismatic high-temperature gas-cooled reactors: Sensitivity of reactor behavior during design basis accidents to fuel properties and the potential impact of the SiC defect annealing process

We conducted reactor performance calculations to assess the potential design basis accident performance of HTGR fuel designs. Three Fully Ceramic Microencapsulated (FCM) fueled HTGR designs were developed in a previous work (Lu et al., 2018). The maximum fuel temperature in the cores fueled by these three FCM fuels was predicted to be higher than that in the reference 350-MWt mHTGR core in both normal operating conditions and during representative design basis accidents (Lu and Brown, 2019). To better understand the potential safety margins in mHTGR design basis accidents, we performed thermal-hydraulics sensitivity studies to investigate how maximum fuel temperature varies considering various parameters, e.g. thermal properties, within the ranges corresponding to the differences between the FCM-fueled prismatic mHTGR cores and the reference core with conventional fuel compacts. We found that the difference in the steady-state axial power distribution contributed the most to the difference in the maximum fuel temperature, in both normal operation and design basis accidents.Experimental data suggested that the annealing process of irradiation defects in SiC would be rapid at mHTGR relevant fuel temperatures. The bounding potential impact of the SiC annealing on the maximum fuel temperature was analyzed considering both the thermal conductivity recovery and the Wigner energy release due to the annealing of SiC. We found that the defect annealing process in SiC would at most increase the peak maximum fuel temperature of an FCM-fueled core by 40 K in loss of forced cooling accidents and by 10 K in a control rod withdrawal accident. Finally, additional experiments on the SiC defect annealing kinetics and Wigner energy release in more relevant conditions are needed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

An Improved Convective Ice Parameterization for the NASA GISS Global Climate Model and Impacts on Cloud Ice Simulation

Partitioning of convective ice into precipitating and detrained condensate presents a challenge for GCMs since partitioning depends on the strength and microphysics of the convective updraft. It is an important issue because detrainment of ice from updrafts influences the development of stratiform anvils, impacts radiation, and can affect GCM climate sensitivity. Recent studies have shown that the CMIP5 configurations of the Goddard Institute for Space Studies (GISS) GCM simulated upper-tropospheric ice water content (IWC) that exceeded an estimated upper bound by a factor of 2. Partly in response to this bias, a new GCM parameterization of convective cloud ice has been developed that incorporates new ice particle fall speeds and convective outflow particle size distributions (PSDs) from the NASA African Monsoon Multidisciplinary Analyses (NAMMA), NASA Tropical Composition, Cloud and Climate Coupling (TC4), DOE ARM–NASA Midlatitude Continental Convective Clouds Experiment (MC3E), and DOE ARM Small Particles in Cirrus (SPARTICUS) field campaigns. The new parameterization assumes a normalized gamma PSD with two novel developments: no explicit assumption for particle habit in the calculation of mass distributions, and a formulation for translating ice particle fall speeds as a function of maximum diameter into fall speeds as a function of melted-equivalent diameter. Two parameters (particle volume– and projected area–weighted equivalent diameter) are diagnosed as a function of temperature and IWC in the convective plume, and these parameters constrain the shape and scale of the normalized gamma PSD. The diagnosed fall speeds and PSDs are combined with the GCM’s parameterized convective updraft vertical velocity to partition convective updraft condensate into precipitating and detrained components. A 5-yr prescribed sea surface temperature GCM simulation shows a 30%–50% decrease in upper-tropospheric deep convective IWC, bringing the tropical and global mean ice water path into closer agreement with CloudSat best estimates.

54 ENVIRONMENTAL SCIENCES↗