Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code optimization”

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 307 records · Page 17

A numerical study of the effect of various reactions, pressure and gas mixture ratio on the density distribution of etchant species (H, Br, Br{sup +}, and HBr{sup +}) in HBr/He plasma

In this study, a fluid model has been used to study the effect of gas mixing ratio and pressure on the density distribution of important etchant species, i.e., hydrogen (H), bromine (Br), Br{sup +}, and HBr{sup +} in HBr/He plasma. Our simulation results show that the densities of active etchant species H, Br, and HBr{sup +} increase with the increase in pressure as well as the HBr fraction in HBr/He mixture. On the contrary, the density of Br{sup +} decreases with the increase in He percentage in HBr/He mixture and with the increase in the pressure. Time averaged reaction rates (of the reactions involved in the production and consumption of these species) have been calculated to study the effect of these reactions on the density distribution of these species. The spatial distribution of these species is explained with the help of the time averaged reaction rates. Important reactions have been identified that contribute considerably to the production and consumption of these active species. The code has been optimized by identifying 26 reactions (out of 40 reactions which contribute in the production and consumption of these species) that have insignificant effect on the densities of H, Br, Br{sup +}, and HBr{sup +}. This shows that out of 40 reactions, only 14 reactions can be used to calculate the density and distribution of the important species in HBr/He plasma discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GronOR: Massively Parallel and GPU-Accelerated Non-Orthogonal Configuration Interaction for Large Molecular Systems

GronOR is a program package for non-orthogonal configuration interaction calculations for an electronic wave function built in terms of anti-symmetrized products of multi-configuration molecular fragment wave functions. The two-electron integrals that have to be processed may be expressed in terms of atomic orbitals or in terms of an orbital basis determined from the molecular orbitals of the fragments. The code has been specifically designed for execution on distributed memory massively parallel and Graphics Processing Unit (GPU)-accelerated computer architectures, using an MPI+OpenACC/OpenMP programming approach. The task-based execution model used in the implementation allows for linear scaling with the number of nodes on the largest pre-exascale architectures available, provides hardware fault resiliency, and enables effective execution on systems with distinct central processing unit-only and GPU-accelerated partitions. The code interfaces with existing multi-configuration electronic structure codes that provide optimized molecular fragment orbitals, configuration interaction coefficients, and the required integrals. Algorithm and implementation details, parallel and accelerated performance benchmarks, and an analysis of the sensitivity of the accuracy of results and computational performance to thresholds used in the calculations are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An ALMA survey of sub-millimeter galaxies in the extended Chandra deep field south: Physical properties derived from ultraviolet-to-radio modeling

The ALESS survey has followed up on a sample of 122 sub-millimeter sources in the Extended Chandra Deep Field South at 870 μm with the Atacama Large Millimeter Array (ALMA), allowing us to pinpoint the positions of sub-millimeter galaxies (SMGs) to ∼0.3 arcsec and to find their precise counterparts at different wavelengths. This enabled the first compilation of the multi-wavelength spectral energy distributions (SEDs) of a statistically reliable survey of SMGs. In this paper, we present a new calibration of the magphys SED modeling code that is optimized to fit these ultraviolet-to-radio SEDs of z>1 star-forming galaxies using an energy balance technique to connect the emission from stellar populations, dust attenuation, and dust emission in a physically consistent way. We derive statistically and physically robust estimates of the photometric redshifts and physical parameters (such as stellar masses, dust attenuation, star formation rates (SFRs), and dust masses) for the ALESS SMGs. We find that the ALESS SMGs have median stellar mass M{sub ∗}=8.9±0.1×10{sup 10} M{sub ⊙}, median SFR = 280±70 M{sub ⊙} yr{sup −1}, median overall V-band dust attenuation A{sub V}=1.9±0.2 mag, median dust mass M{sub dust}=(5.6±1.0)×10{sup 8} M{sub ⊙}, and median average dust temperature T{sub dust}≃40 K. We find that the average intrinsic SED of the ALESS SMGs resembles that of local ultra-luminous infrared galaxies in the infrared range, but the stellar emission of our average SMG is brighter and bluer, indicating lower dust attenuation, possibly because they are more extended. We explore how the average SEDs vary with different parameters (redshift, sub-millimeter flux, dust attenuation, and total infrared luminosity), and we provide a new set of SMG templates that can be used to interpret other SMG observations. To put the ALESS SMGs into context, we compare their stellar masses and SFRs with those of less actively star-forming galaxies at the same redshifts. We find that at z≃2, about half of the SMGs lie above the star-forming main sequence (with SFRs three times larger than normal galaxies of the same stellar mass), while half are consistent with being at the high-mass end of the main sequence. At higher redshifts (z≃3.5), the SMGs tend to have higher SFRs and stellar masses, but the fraction of SMGs that lie significantly above the main sequence decreases to less than a third.

79 ASTRONOMY AND ASTROPHYSICS↗

Spatial noise correlations in a Si/SiGe two-qubit device from Bell state coherences

Here, we study spatial noise correlations in a Si/SiGe two-qubit device with integrated micromagnets. Our method relies on the concept of decoherence-free subspaces, whereby we measure the coherence time for two different Bell states, designed to be sensitive only to either correlated or anticorrelated noise, respectively. From these measurements we find weak correlations in low-frequency noise acting on the two qubits, while no correlations could be detected in high-frequency noise. We expect nuclear spin noise to have an uncorrelated nature. A theoretical model and numerical simulations give further insight into the additive effect of multiple independent (anti)correlated noise sources with an asymmetric effect on the two qubits as can result from charge noise. Such a scenario in combination with nuclear spins is plausible given the data and the known decoherence mechanisms. This work is highly relevant for the design of optimized quantum error correction codes for spin qubits in quantum dot arrays, as well as for optimizing the design of future quantum dot arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Encoding an Oscillator into Many Oscillators

An outstanding challenge for quantum information processing using bosonic systems is Gaussian errors such as excitation loss and added thermal noise errors. Thus, bosonic quantum error correction is essential. Most bosonic quantum error correction schemes encode a finite-dimensional logical qubit or qudit into noisy bosonic oscillator modes. In this case, however, the infinite-dimensional bosonic nature of the physical system is lost at the error-corrected logical level. On the other hand, there are several proposals for encoding an oscillator mode into many noisy oscillator modes. However, these oscillator-into-oscillators encoding schemes are in the class of Gaussian quantum error correction. Therefore, these codes cannot correct practically relevant Gaussian errors due to the established no-go theorems that state that Gaussian errors cannot be corrected by using only Gaussian resources. Here, we circumvent these no-go results and show that it is possible to correct Gaussian errors by using Gottesman-Kitaev-Preskill (GKP) states as non-Gaussian resources. In particular, we propose a non-Gaussian oscillator-into-oscillators code, namely the GKP two-mode squeezing code, and demonstrate that it can quadratically suppress additive Gaussian noise errors in both the position and momentum quadratures except for a small sublogarithmic correction. Furthermore, we demonstrate that our GKP two-mode squeezing code is near optimal in the weak noise limit by proving via quantum information theoretic tools that quadratic noise suppression is optimal when we use two physical oscillator modes. Lastly, in this work we show that our non-Gaussian oscillator encoding scheme can also be used to correct excitation loss and thermal noise errors, which are dominant error sources in many realistic bosonic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

WinPSSP : a revamp of the computer program PSSP and its performance solving the crystal structures of small organic compounds and solids of biological and pharmaceutical interest

The direct-space methods software Powder Structure Solution Program ( PSSP ) [Pagola & Stephens (2010). J. Appl. Cryst. 43 , 370–376] has been migrated to the Windows OS and the code has been optimized for fast runs. WinPSSP is a user-friendly graphical user interface that allows the input of preliminary crystal structure information, integrated intensities of the reflections and FWHM, the definition of structural parameters and a simulated annealing schedule, and the visualization of the calculated and experimental diffraction data overlaid for each individual solution. The solutions are reported as filename.cif files, which can be used to analyze packing motifs and chemical bonding, and to input the atomic coordinates into the Rietveld analysis software GSAS. WinPSSP performance in straightforward crystal structure determinations has been evaluated using 18 molecular solids with 6–20 degrees of freedom. The free-distribution program as well as multimedia tutorials can be accessed at http://users.uoi.gr/nkourkou/winpssp/.

60 APPLIED LIFE SCIENCES↗

Parallel Simulation of Quantum Networks with Distributed Quantum State Management

Quantum network simulators offer the opportunity to cost-efficiently investigate potential avenues for building networks that scale with the number of users, communication distance, and application demands by simulating alternative hardware designs and control protocols. Several quantum network simulators have been recently developed with these goals in mind. As the size of the simulated networks increases, however, sequential execution becomes time-consuming. Parallel execution presents a suitable method for scalable simulations of large-scale quantum networks, but the unique attributes of quantum information create unexpected challenges. In this work, we identify requirements for parallel simulation of quantum networks and develop the first parallel discrete-event quantum network simulator by modifying the existing serial simulator SeQUeNCe. Our contributions include the design and development of a quantum state manager (QSM) that maintains shared quantum information distributed across multiple processes. We also optimize our parallel code by minimizing the overhead of the QSM and decreasing the amount of synchronization needed among processes. Using these techniques, we observe a speedup of 2 to 25 times when simulating a 1,024-node linear network topology using 2 to 128 processes. We also observe an efficiency greater than 0.5 for up to 32 processes in a linear network topology of the same size and with the same workload. We repeat this evaluation with a randomized workload on a caveman network. We also introduce several methods for partitioning networks by mapping them to different parallel simulation processes. We have released the parallel SeQUeNCe simulator as an open source tool alongside the existing sequential version.

97 MATHEMATICS AND COMPUTING↗

Proteus

Proteus enables Just-In-Time (JIT) compilation and optimization of C/C++ code using LLVM.

Georgakoudis, Giorgis↗

A Note on Compiling Fortran

Fortran modules tend to serialize compilation of large Fortran projects, by introducing dependencies among the source files. If file A depends on file B, (A uses a module defined by B), you must finish compiling B before you can begin compiling A. Some Fortran compilers (Intel ifort, GNU gfortran and IBM xlf, at least) offer an option to ‘‘verify syntax’’, with the side effect of also producing any associated Fortran module files. As it happens, this option usually runs much faster than the object code generation and optimization phases. For some projects on some machines, it can be advantageous to compile in two passes: The first pass generates the module files, quickly; the second pass produces the object files, in parallel. We achieve a 3.8× speedup in the case study below.

97 MATHEMATICS AND COMPUTING↗

High-Performance First-Principles Molecular Dynamics for Predictive Theory and Modeling

This project focused on developing high-performance software tools for First-Principles Molecular Dynamics (FPMD) simulations, and applying them in investigations of materials relevant to energy conversion processes. FPMD is an atomistic simulation method that combines a quantum-mechanical description of electronic structure with the statistical description provided by molecular dynamics (MD) simulations. This reliance on fundamental principles allows FPMD simulations to provide a consistent description of structural, dynamical and electronic properties of a material. This is particularly useful in systems for which reliable empirical models are lacking. FPMD simulations are increasingly used as a predictive tool for applications such as batteries, solar energy conversion, light-emitting devices, electro-chemical energy conversion devices and other materials. During the course of the project, several new features were developed and added to the open-source Qbox FPMD code. The code was further optimized for scalable operation of large-scale, Leadership-Class DOE computers. When combined with Many-Body Perturbation Theory (MBPT) calculations, this infrastructure was used to investigate structural and electronic properties of liquid water, ice, aqueous solutions, nanoparticles and solid-liquid interfaces. Computing both ionic trajectories and electronic structure in a consistent manner enabled the simulation of several spectroscopic properties, such as Raman spectra, infrared spectra, and sum-frequency generation spectra. The accuracy of the approximations used allowed for direct comparisons of results with experimental data such as optical spectra, X-ray and neutron diffraction spectra. The software infrastructure developed in this project, as applied to various investigations of solids, liquids and interfaces, demonstrates that FPMD simulations can provide a detailed, atomic-scale picture of structural, vibrational and electronic properties of complex systems relevant to energy conversion devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Performance Portability of the Aeras Atmosphere Model to Next Generation Architectures using Kokkos

The subject of this report is the performance portability of the Aeras global atmosphere dynamical core (implemented within the Albany multi-physics code) to new and emerging architecture machines using the Kokkos library and programming model. We describe the process of refactoring the finite element assembly process for the 3D hydrostatic model in Aeras and highlight common issues associated with development on GPU architectures. After giving detailed build and execute instructions for Aeras with MPI, OpenMP and CUDA on the Shannon cluster at Sandia National Laboratories and the Titan supercomputer at Oak Ridge National Laboratory, we evaluate the per- formance of the code on a canonical test case known as the baroclinic instability problem. We show a speedup of up to 4 times on 8 OpenMP threads, but we were unable to achieve a speedup on the GPU due to memory constraints. We conclude by providing methods for improving the performance of the code for future optimization.

97 MATHEMATICS AND COMPUTING↗

Transient Physical Modelling for the Coupling of a High Temperature Gas Reactor with Thermal Energy Storage

Improving the flexibility of nuclear reactors to meet fluctuating demand on grids with high renewable energy generation penetration is key to improving reactor economics. Economic optimization using the code HERON gives maximum profit of coupling high temperature gas reactors (HTGRs) to thermal energy storage (TES) systems in the configuration shown in Figure 1 when following the demand in Figure 3. We look to model the physics of operation for a HTGR-TES using the HYBRID library in the code Modelica.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational Optimization of A 3D Printed Collimator

This contribution describes the computational methodology behind an optimization procedure for a scattered beam collimator. The workflow includes producing a file that can be manufactured via additive methods. A conical collimator, optimized for neutron diffraction experiments in a high pressure clamp cell, is presented as an example. In such a case the scattering from the sample is much smaller than that of the pressure cell. Monte Carlo Ray tracing in MCViNE was used to model scattering from a Si powder sample and the cell. A collimator was inserted into the simulation and the number and size of channels were optimized to maximize the rejection of the parasitic signal coming from the complex sample environment. Constraints, provided by the additive manufacturing process as well as a specific neutron diffractometer, were also included in the optimization. The source code and the tutorials are available in c3dp (Islam (2019)).

74 ATOMIC AND MOLECULAR PHYSICS↗

SIMULATIONS OF A MUON LINAC FOR A NEUTRINO FACTORY

The Neutrino Factory baseline design involves a complex chain of accelerators including a single-pass linac, two recirculating linacs and an FFAG. The first linac follows the capture and bunching section and accelerates the muons from about 244 to 900 MeV. It must accept a high emittance beam about 30 cm wide with a 10% energy spread. This linac uses counterwound, shielded superconducting solenoids and 201 MHz superconducting cavities. Simulations have been carried out using several codes including Zgoubi, OptiM, GPT, Elegant and G4beamline, both to determine the optics and to estimate the radiation loads on the elements due to beam loss and muon decay.

43 PARTICLE ACCELERATORS↗

BlackjackBench: Portable Hardware Characterization

DARPA's AACE project aimed to develop Architecture Aware Compiler Environments that automatically characterize the hardware and optimize the application codes accordingly. We present the BlackjackBench suite, a collection of portable micro-benchmarks that automate system characterization, plus statistical analysis techniques for interpreting the results. BlackjackBench discovers the effective sizes and speeds of the hardware environment rather than the often unattainable peak values. We aim at hardware features that can be observed by running executables generated by existing compilers from standard C codes. We characterize the memory hierarchy, including cache sharing and NUMA characteristics of the system, properties of the processing cores affecting execution speed, and the length of the OS scheduler time slot. We show how these features of modern multicores can be discovered programmatically. We also show how the features could interfere with each other resulting in incorrect interpretation of the results, and how established classification and statistical analysis techniques reduce experimental noise and aid automatic interpretation of results.

99 GENERAL AND MISCELLANEOUS↗

Dynamic System Scaling Application to Accelerated Nuclear Fuel Testing

The development of nuclear fuel and materials requires a continuous effort to investigate the acute and prolonged effects of irradiation, thermal-material stress, chemical change, or other conceivable damage mechanics acquired during normal operation or accident scenarios throughout its lifetime. As in-core fuel property measurement techniques advance to support in the real-time, non-invasive, and enhanced accuracy realm, it is the future of fuel development to pursue to higher degrees of control, predictability of integral test behavior via separate effects test (SET), and shorter test time intervals. The fuel development life cycle from initial concept to commercial licensing is approximated to be 20 years and current literature suggests by optimizing fuel performance codes with SETs, the process could possibly be compressed to 5 to 10 years. Recently, a research group in the Idaho National Laboratory (INL) is testing reduced scale fuel rods and increased power density to accelerate evolution of fuel phenomena in metallic fuels. In support of nuclear fuel rod development, compressing fuel test process, and accelerating fuel phenomena, it was the purpose of this study to investigate nuclear fuel performance phenomena via literature review and effectively scale the initial conditions, boundary conditions, and geometric properties to describe the time-dependent response including to fuel burnup, thermal and mechanical stress, transmutations and inter-diffusion, and other relevant observed phenomena. The study was based on BISON simulations of historic EBR-II metallic fuel experiments and Dynamical System Scaling (DSS) method are utilized to assess effects of scaling activity including fuel phenomena acceleration and calculations of time-dependent distortions. The research successfully scaled metallic fuel phenomena, accelerated fuel testing, and assessed the distortions for each scaled case derived.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

E-Beam Irradiation Beamline at Jefferson Lab for Wastewater Treatment

Electron beam (e-beam) irradiation is an advanced method for removing the emerging organic contaminants in wastewater, such as 1,4-dioxane, and per- and polyfluoroalkyl sub?stances (PFAS), which have the potential to be regulated in the near future but are difficult to degrade by conventional treatment methods. This dissertation was based on a collabo?rating study between Jefferson Lab (Jlab) and Hampton Roads Sanitation District (HRSD) to evaluate the impacts of e-beam irradiation on the degradation of 1,4-dioxane and PFAS. It presents the detailed design of an e-beam irradiation beamline, the beamline simulation using general particle tracer (GPT), the installation, and the successful commissioning. The beamline is an extended section to the Upgraded Injector Test Facility (UITF) at Jlab, which is a continuous wave (CW), superconducting radio frequency (SRF) linear accelerator with an energy up to 10 MeV. To let the sample volume receive a significant and uniform dose distribution, the sample target was designed by investigating the dose distribution simulated with a Monte-Carlo code, FLUKA. The optimized beam energy was 8 MeV and the beam current for irradiation was ? 100 nA. With the use of a solenoid and a raster, the electron beam was expanded to the optimized size, 1? ? 15 mm, at the sample target. This dis?sertation also presents the irradiation experiments for 64 1,4-dioxane samples and 48 PFAS samples in different water matrices, including real wastewater, which were prepared and an?alyzed by HRSD. The absorbed doses of samples were obtained by designing a new method that combines the FLUKA simulation and dose mapping experiments. The irradiation re?sults showed that around 95% of 1,4-dioxane was removed at a low dose of ? 2 kGy, and PFAS was possibly removed at a lower dose compared to similar research. This work is a first implementation of SRF technology for wastewater irradiation. It pro?vides a valuable new research platform for wastewater remediation. This will be in support of future industrial applications of advanced SRF technologies and better understanding of the effect of e-beam irradiation on wastewater or water treatment.

Li, Xi↗