Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scaling”

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 217 records · Page 12

Residual stress accumulation in large-scale Ti-6Al-4V wire-arc additive manufacturing

Large-scale additive manufacturing (AM) is of interest in the manufacturing industry to produce components of dimensions on the order of hundreds of millimeters to meters in scale. Wire-arc AM (WAAM) of Ti-6Al-4V (Ti64) is an attractive technique for large-scale AM in the aerospace industry, give the high strength-to-weight ratio of the material and high deposition rate of the process. However, due to the large scale, significant distortion and residual stresses are developed in the material during deposition and cooling, potentially leading to part failure. WAAM of a prototypical large-scale Ti64 machine tool component is studied in this work via finite element analysis (FEA). Element activation/deactivation technique is employed to simulate deposition and the resulting distortion and residual stress (RS) predictions are analyzed to evaluate the possibility of crack formation. Significant distortion (∼10 mm) and RS (>1300 MPa) is predicted, suggesting the formation of cracks and possible crack propagation into the build region. Incorporation of fillets significantly reduces RS concentration regions, thus reducing the possibility of part failure. Material deposition sequence also affects the RS pattern in the build.

Mathews, Ritin↗

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scaling of Transition State Vibrational Frequencies and Application of d -Band Theory to the Brønsted–Evans–Polanyi Relationship on Surfaces

Semiempirical energy relations provide a means of estimating thermodynamic properties. Specifically, linear scaling relationships (LSRs) and Brønsted–Evans–Polanyi (BEP) relationships correlate adsorption energies between adsorbates across surfaces and reaction energies with activation barriers, respectively. Although vibrational scaling relations (VSRs) exist between adsorbates at identical sites, scaling between vibrational frequencies of adsorbed local minima and transition states is lacking. Here, we present density functional theory calculations for AH X (A = C, N, O) diffusions on transition metal surfaces and reveal linear scaling between frequencies of local minima and the transition state between those minima. Using d-band theory and linear muffin tin orbital theory (LMTO), we derive the slopes of these transition state vibrational scaling relations (TSVSRs) and, in so doing, provide a rigorous theory extending the original BEP relations developed for solution chemistry to surface chemistry. Furthermore, with a single reference DFT calculation, we predict the slopes and quantify uncertainty in the predictions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Soil Respiration Variability and Correlation Across a Wide Range of Temporal Scales

Abstract The high temporal variability of the soil‐to‐atmosphere CO 2 flux (soil respiration, R S ) has been studied at hourly to multiannual time scales but remains less well understood than R S spatial variability. How R S fluxes vary and are autocorrelated at various time lags has practical implications for sampling and more fundamentally for our understanding of its abiotic and biotic underlying mechanisms. We examined the variability, correlation, and sampling requirements of R S over a wide range of temporal scales in a temperate deciduous forest in eastern Maryland, USA, using both automated (temporally continuous, N = 30,036 over 10 months) and survey (spatially diverse, temporally sparse, N = 1,912 over 17 months) data. Data from a global R S database were also used to examine interannual variability in comparable forests. The coefficient of variability of successive measurements generally varied from the minute (median coefficient of variation 16%) to hourly and daily (11–12%) time scales. Successive R S values measured at a given collar exhibited a strong hour‐to‐hour correlation ( r = 0.931) and a moderate correlation at a 2‐hr lag (0.289); day‐to‐day (i.e., 24 hr lag) hourly observations were uncorrelated. Daily R S means were well correlated at a 1‐day lag ( r = 0.856) but not at any further time lag. In a linear mixed‐effects model predicting R S , soil temperature and moisture exerted consistently strong effects regardless of time scale, and model coefficient of variability was generally high (>80%). These results provide new opportunities to explore the drivers and variability of R S fluxes, quantify sampling requirements, and improve error propagation.

Bond‐Lamberty, Ben↗

Stabilization of electron-scale turbulence by electron density gradient in national spherical torus experiment

Theory and experiments have shown that electron temperature gradient (ETG) turbulence on the electron gyro-scale, k{sub ⊥}ρ{sub e} ≲ 1, can be responsible for anomalous electron thermal transport in NSTX. Electron scale (high-k) turbulence is diagnosed in NSTX with a high-k microwave scattering system [D. R. Smith et al., Rev. Sci. Instrum. 79, 123501 (2008)]. Here we report on stabilization effects of the electron density gradient on electron-scale density fluctuations in a set of neutral beam injection heated H-mode plasmas. We found that the absence of high-k density fluctuations from measurements is correlated with large equilibrium density gradient, which is shown to be consistent with linear stabilization of ETG modes due to the density gradient using the analytical ETG linear threshold in F. Jenko et al. [Phys. Plasmas 8, 4096 (2001)] and linear gyrokinetic simulations with GS2 [M. Kotschenreuther et al., Comput. Phys. Commun. 88, 128 (1995)]. We also found that the observed power of electron-scale turbulence (when it exists) is anti-correlated with the equilibrium density gradient, suggesting density gradient as a nonlinear stabilizing mechanism. Higher density gradients give rise to lower values of the plasma frame frequency, calculated based on the Doppler shift of the measured density fluctuations. Linear gyrokinetic simulations show that higher values of the electron density gradient reduce the value of the real frequency, in agreement with experimental observation. Nonlinear electron-scale gyrokinetic simulations show that high electron density gradient reduces electron heat flux and stiffness, and increases the ETG nonlinear threshold, consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory

Domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)) is a highly efficient local correlation method. It is known to be accurate and robust and can be used in a black box fashion in order to obtain coupled cluster quality total energies for large molecules with several hundred atoms. While previous implementations showed near linear scaling up to a few hundred atoms, several nonlinear scaling steps limited the applicability of the method for very large systems. In this work, these limitations are overcome and a linear scaling DLPNO-CCSD(T) method for closed shell systems is reported. The new implementation is based on the concept of sparse maps that was introduced in Part I of this series [P. Pinski, C. Riplinger, E. F. Valeev, and F. Neese, J. Chem. Phys. 143, 034108 (2015)]. Using the sparse map infrastructure, all essential computational steps (integral transformation and storage, initial guess, pair natural orbital construction, amplitude iterations, triples correction) are achieved in a linear scaling fashion. In addition, a number of additional algorithmic improvements are reported that lead to significant speedups of the method. The new, linear-scaling DLPNO-CCSD(T) implementation typically is 7 times faster than the previous implementation and consumes 4 times less disk space for large three-dimensional systems. For linear systems, the performance gains and memory savings are substantially larger. Calculations with more than 20 000 basis functions and 1000 atoms are reported in this work. In all cases, the time required for the coupled cluster step is comparable to or lower than for the preceding Hartree-Fock calculation, even if this is carried out with the efficient resolution-of-the-identity and chain-of-spheres approximations. The new implementation even reduces the error in absolute correlation energies by about a factor of two, compared to the already accurate previous implementation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying particle-scale 3D granular dynamics during rapid compaction from time-resolved in situ 2D x-ray images

Understanding the particle-scale dynamics of granular materials during rapid compaction and flow is of fundamental importance for manufacturing, planetary science, geology, and defense applications. Time-resolved 2D radiography and static 3D x-ray tomography are powerful in situ tools for studying particle-scale dynamics but provide detail only in 2D or with significant time-scale limitations, respectively. Here, we introduce a new method that uses 2D in situ x-ray imaging for determining time-resolved 3D particle-scale dynamics in rapidly compressed granular materials. The method employs initial particle packing structures obtained from x-ray tomography, a 2D x-ray image generation algorithm, and an optimization algorithm. We first describe and validate the method using finite element simulations. Furthermore, we then apply the technique to x-ray phase-contrast images obtained during rapid compaction of granular materials with varying particle morphology and sample thickness. The depth-resolved particle-scale dynamics reveal complex velocity and porosity fields evolving heterogeneously along and perpendicular to the compaction direction. We characterize these features, their fluctuations near the compaction front, and the compaction front thickness. Our technique can be applied to understanding granular dynamics during rapid compaction events, and rearrangements during slower, but non-quasi-static, flows.

36 MATERIALS SCIENCE↗

Exploring the parameter space of MagLIF implosions using similarity scaling. III. Rise-time scaling

Magnetized liner inertial fusion (MagLIF) is a z-pinch magneto-inertial-fusion concept studied at the Z Pulsed Power Facility of Sandia National Laboratories. Two important metrics characterizing current delivery to a z-pinch load are the peak current and the current-rise time, which is roughly the time interval to reach the peak current. It is known that, when driving a z-pinch load with a longer current-rise time, the performance of the z-pinch decreases. However, a theory to understand and quantify this effect is still lacking. Here, in this paper, we utilize a framework based on similarity scaling to analytically investigate the variations in the performance of MagLIF loads when varying the current-rise time, or equivalently, the implosion timescale. To maintain similarity between the implosions, we provide scaling prescriptions of experimental input parameters defining a MagLIF load and derive the expected scaling laws for stagnation conditions and for various performance metrics. We compare predictions of the theory to 2D numerical simulations using the radiation, magneto-hydrodynamic code hydra. For several metrics, we find acceptable agreement between the theory and simulations. Our results show that the voltage φload near the MagLIF load follows a weak scaling law φload ∝$t_{φ}^{-0.12}$ with respect to the characteristic timescale t φ of the voltage source, instead of the ideal φload ∝$t_{φ}^{-1}$ scaling. This occurs because the imploding height of the MagLIF load must increase to preserve end losses. As a consequence of the longer imploding liners, the required total laser preheat energy and delivered electric energy increase. Overall, this study helps understand the trade-offs of the MagLIF design space when considering future pulsed-power generators with shorter and longer current-rise times.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

THE FRONTIER BETWEEN SMALL-SCALE BIPOLES AND EPHEMERAL REGIONS IN THE SOLAR PHOTOSPHERE: EMERGENCE AND DECAY OF AN INTERMEDIATE-SCALE BIPOLE OBSERVED WITH SUNRISE/IMaX

We report on the photospheric evolution of an intermediate-scale ( Almost-Equal-To 4 Mm footpoint separation) magnetic bipole, from emergence to decay, observed in the quiet Sun at high spatial (0.''3) and temporal (33 s) resolution. The observations were acquired by the Imaging Magnetograph Experiment imaging magnetograph during the first science flight of the SUNRISE balloon-borne solar observatory. The bipole flux content is 6 Multiplication-Sign 10{sup 17} Mx, representing a structure bridging the gap between granular scale bipoles and the smaller ephemeral regions. Footpoints separate at a speed of 3.5 km s{sup -1} and reach a maximum distance of 4.5 Mm before the field dissolves. The evolution of the bipole is revealed to be very dynamic: we found a proper motion of the bipole axis and detected a change of the azimuth angle of 90 Degree-Sign in 300 s, which may indicate the presence of some writhe in the emerging structure. The overall morphology and behavior are in agreement with previous analyses of bipolar structures emerging at the granular scale, but we also found several similarities with emerging flux structures at larger scales. The flux growth rate is 2.6 Multiplication-Sign 10{sup 15} Mx s{sup -1}, while the mean decay rate is one order of magnitude smaller. We describe in some detail the decay phase of the bipole footpoints that includes break up into smaller structures, and interaction with preexisting fields leading to cancellation, but it appears to be dominated by an as-yet unidentified diffusive process that removes most of the flux with an exponential flux decay curve. The diffusion constant (8 Multiplication-Sign 10{sup 2} km{sup 2} s{sup -1}) associated with this decay is similar to the values used to describe the large-scale diffusion in flux transport models.

79 ASTRONOMY AND ASTROPHYSICS↗

Does jackknife scale really matter for accurate large-scale structure covariances?

ABSTRACT The jackknife method gives an internal covariance estimate for large-scale structure surveys and allows model-independent errors on cosmological parameters. Using the SDSS-III BOSS CMASS sample, we study how the jackknife size and number of resamplings impact the precision of the covariance estimate on the correlation function multipoles and the error on the inferred baryon acoustic scale. We compare the measurement with the MultiDark Patchy mock galaxy catalogues, and we also validate it against a set of lognormal mocks with the same survey geometry. We build several jackknife configurations that vary in size and number of resamplings. We introduce the Hartlap factor in the covariance estimate that depends on the number of jackknife resamplings. We also find that it is useful to apply the tapering scheme to estimate the precision matrix from a limited number of resamplings. The results from CMASS and mock catalogues show that the error estimate of the baryon acoustic scale does not depend on the jackknife scale. For the shift parameter α, we find an average error of 1.6 per cent, 2.2 per cent and 1.2 per cent, respectively, from CMASS, Patchy, and lognormal jackknife covariances. Despite these uncertainties fluctuate significantly due to some structural limitations of the jackknife method, our α estimates are in reasonable agreement with published pre-reconstruction analyses. Jackknife methods will provide valuable and complementary covariance estimates for future large-scale structure surveys.

Favole, Ginevra (ORCID:000000028218563X)↗

Sunyaev-Zel'dovich Effect and X-ray Scaling Relations from Weak-Lensing Mass Calibration of 32 SPT Selected Galaxy Clusters

Uncertainty in mass-observable scaling relations is currently the limiting factor for galaxy-cluster-based cosmology. Weak gravitational lensing can provide direct mass calibration and reduce the mass uncertainty. We present new ground-based weak lensing observations of 19 South Pole Telescope (SPT) selected clusters at redshifts 0.29 <= z <= 0.61 and combine them with previously reported space-based observations of 13 galaxy clusters at redshifts 0.576 <= z <= 1.132 to constrain the cluster mass scaling relations with the Sunyaev-Zel'dovich effect (SZE), the cluster gas mass M-gas and Y-X, the product of M-gas and X-ray temperature. We extend a previously used framework for the analysis of scaling relations and cosmological constraints obtained from SPT-selected clusters to make use of weak lensing information. We introduce a new approach to estimate the effective average redshift distribution of background galaxies and quantify a number of systematic errors affecting the weak lensing modelling. These errors include a calibration of the bias incurred by fitting a Navarro-Frenk-White profile to the reduced shear using N-body simulations. We blind the analysis to avoid confirmation bias. We are able to limit the systematic uncertainties to 5.6% in cluster mass (68% confidence). Our constraints on the mass-X-ray observable scaling relation parameters are consistent with those obtained by earlier studies and our constraints for the mass-SZE scaling relation are consistent with the simulation-based prior used in the most recent SPT-SZ cosmology analysis. We can now replace the external mass calibration priors used in previous SPT-SZ cosmology studies with a direct, internal calibration obtained for the same clusters.

cosmology: observations↗

A Bloom Filter Based Scalable Data Integrity Check Tool for Large-Scale Dataset

Large scale HPC applications are becoming increasingly data intensive. At Oak Ridge Leadership Computing Facility (OLCF), we are observing the number of files curated under individual project are reaching as high as 200 millions and project data size is exceeding petabytes. These simulation datasets, once validated, often needs to be transferred to archival system for long term storage or shared with the rest of the research community. Ensuring the data integrity of the full dataset at this scale is paramount important but also a daunting task. This is especially true considering that most conventional tools are serial and file-based, unwieldy to use and/or can't scale to meet user's demand.To tackle this particular challenge, this paper presents the design, implementation and evaluation of a scalable parallel checksumming tool, fsum, which we developed at OLCF. It is built upon the principle of parallel tree walk and work-stealing pattern to maximize parallelism and is capable of generating a single, consistent signature for the entire dataset at extreme scale. We also applied a novel bloom-filter based technique in aggregating signatures to overcome the signature ordering requirement. Given the probabilistic nature of bloom filter, we provided a detailed error and trade-off analysis. Using multiple datasets from production environment, we demonstrated that our tool can efficiently handle both very large files as well as many small-file based datasets. Our preliminary test showed that on the same hardware, it outperforms conventional tool by as much as 4×. It also exhibited near-linear scaling properties when provisioned with more compute resources.

97 MATHEMATICS AND COMPUTING↗

Spatial scale drives patterns in soil bacterial diversity: Spatial scale drives soil diversity

Soil microbial communities are essential for ecosystem function, but linking community composition to biogeochemical processes is challenging because of high microbial diversity and large spatial variability of most soil characteristics. We investigated soil bacterial community structure in a switchgrass stand planted on soil with a history of grassland vegetation at high spatial resolution to determine whether biogeographic trends occurred at the centimeter scale. Moreover, we tested whether such heterogeneity, if present, influenced community structure within or among ecosystems. Pronounced heterogeneity was observed at centimeter scales, with abrupt changes in relative abundance of phyla from sample to sample. At the ecosystem scale (> 10 m), however, bacterial community composition and structure were subtly, but significantly, altered by fertilization, with higher alpha diversity in fertilized plots. Moreover, by comparing these data with data from 1772 soils from the Earth Microbiome Project, it was found that 20% diverse globally sourced soil samples, while grassland soils shared approximately 40% of their operational taxonomic units with the current study. By spanning several orders of magnitude, the analysis suggested that extreme patchiness characterized community structure at smaller scales but that coherent patterns emerged at larger length scales.

16S rRNA↗

Measurement of k T splitting scales in W→ℓν events at $\sqrt{s} = 7\ \mathrm{TeV}$ with the ATLAS detector

A measurement of splitting scales, as defined by the k T clustering algorithm, is presented for final states containing a W boson produced in proton–proton collisions at a centre-of-mass energy of 7 TeV. The measurement is based on the full 2010 data sample corresponding to an integrated luminosity of 36 pb -1 which was collected using the ATLAS detector at the CERN Large Hadron Collider. Cluster splitting scales are measured in events containing W bosons decaying to electrons or muons. The measurement comprises the four hardest splitting scales in a k T cluster sequence of the hadronic activity accompanying the W boson, and ratios of these splitting scales. Backgrounds such as multi-jet and top-quark-pair production are subtracted and the results are corrected for detector effects. Predictions from various Monte Carlo event generators at particle level are compared to the data. Overall, reasonable agreement is found with all generators, but larger deviations between the predictions and the data are evident in the soft regions of the splitting scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Scaling Deep Learning on GPU and Knights Landing clusters

The speed of deep neural networks training has become a big bottleneck of deep learning research and development. For example, training GoogleNet by ImageNet dataset on one Nvidia K20 GPU needs 21 days. To speed up the training process, the current deep learning systems heavily rely on the hardware accelerators. However, these accelerators have limited on-chip memory compared with CPUs. To handle large datasets, they need to fetch data from either CPU memory or remote processors. We use both self-hosted Intel Knights Landing (KNL) clusters and multi-GPU clusters as our target platforms. From an algorithm aspect, current distributed machine learning systems are mainly designed for cloud systems. These methods are asynchronous because of the slow network and high fault-tolerance requirement on cloud systems. We focus on Elastic Averaging SGD (EASGD) to design algorithms for HPC clusters. Original EASGD used round-robin method for communication and updating. The communication is ordered by the machine rank ID, which is inefficient on HPC clusters. First, we redesign four efficient algorithms for HPC systems to improve EASGD's poor scaling on clusters. Async EASGD, Async MEASGD, and Hogwild EASGD are faster \textcolor{black}{than} their existing counterparts (Async SGD, Async MSGD, and Hogwild SGD, resp.) in all the comparisons. Finally, we design Sync EASGD, which ties for the best performance among all the methods while being deterministic. In addition to the algorithmic improvements, we use some system-algorithm codesign techniques to scale up the algorithms. By reducing the percentage of communication from 87% to 14%, our Sync EASGD achieves 5.3x speedup over original EASGD on the same platform. We get 91.5% weak scaling efficiency on 4253 KNL cores, which is higher than the state-of-the-art implementation.

60 APPLIED LIFE SCIENCES↗