Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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 703 records · Page 39

Coexistence of single-particle and collective excitation in Ni 61

Here, the high-spin states in 61 Ni have been studied using the fusion evaporation reaction, 50 Ti( 14 C, 3n) 61 Ni at an incident beam energy of 40 MeV. A Compton suppressed multi-HPGe detector setup, consisting of six clover detectors and three single-crystal HPGe detectors, was used to detect the de-exciting $\gamma$ rays from the excited states. The level scheme has been extended up to an excitation energy of 12.8 MeV and a tentative J π =35/2 + . The low-lying negative parity levels are found to be generated by single-particle excitation within the fp shell and also excitations to the g 9/2 orbitals as explained well with shell model calculations using the GXPF1Br+V MU (modified) interaction. Two rotational structures of regular E2 sequences with small to moderate axial deformation have been established at higher excitation energy. Most interestingly, two sequences of M1 transitions are reported for the first time and described as magnetic rotational bands. The shears mechanism for both the bands can be described satisfactorily by the geometrical model. The shell model calculation involving the cross shell excitation beyond the fp shell well reproduce the M1 and E2 sequences. The shell model predicted B(M1) values for the magnetic rotational band B1 show the decreasing trend with spin as expected with closing of the shears.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cooling Design and Thermal Analysis for Thermal Shields of a Cryocooler-Cooled Superconducting ECR Ion Source MARS-D Magnet

Here, a demonstrative NbTi based Mixed Axial and Radial field System (MARS-D) is being developed for a next-Generation Electron Cyclotron Resonance Ion Source (ECRIS) at Lawrence Berkeley National Laboratory (LBL), which employs a novel closed-loop coil design scheme that more efficiently utilizes conductor fields and extend the application of NbTi for high frequency (up to 45 GHz) ECR operation. The NbTi MARS-D magnet consists of a single hexagonally shaped closed-loop coil and a set of auxiliary solenoids. A cryostat for cooling the MARS-D magnet is under design at LBL. The MARS-D magnet working around 4.2 K will be bath-cooled in liquid helium using multiple two-stage cryocoolers. An intermediate temperature thermal radiation shield is adopted to reduce the heat leakage imposed on 4.2 K coil cold mass from room temperature. The thermal shield is conduction-cooled by the first-stage cold heads of four two-stage cryocoolers and the cold head of a single-stage cryocooler shared with nine binary leads. The temperature in the area of the shield that warm ends of HTS leads are mounted on is expected no higher than 60 K, which is limited by maximum allowable working temperature of HTS leads. The paper presents thermal analysis on the thermal radiation shield including heat loads and effects of eddy current induced during quench on its material selection.

ECRIS↗

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

The predominance of Kohn–Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high-performance computing (HPC). With recent trends in HPC leading toward increasing reliance on heterogeneous accelerator-based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high levels of performance that have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn–Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.

97 MATHEMATICS AND COMPUTING↗

Combining LANDSAT MSS, aerial photographs and ground measurements to estimate rangeland productivity

The production of a vegetation map of over 2.2 million acres with detail down to the plant community level, and the production of estimates of rangeland productivity (pounds of usable forage per acre for cattle) for a 500,000 acre subset of area with a design goal for accuracy and precision of + or - 20% at the 80% confidence level, are considered. The data consist of five groups: maps of area, LANDSAT data, digital terrain data, large scale aerial photography, and ground plots. An outline of the data acquisition and data reduction schemes are presented.

Gialdini, M. J.↗

Active control of sound transmission/radiation from elastic plates by vibration inputs. II - Experiments

Actively controlled harmonic force inputs were applied experimentally to reduce the sound transmitted through an elastic circular plate. The control implementation used a time domain least mean square adaptive algorithm with two error sensors. The control forces were applied directly to the plate by point force vibration inputs, while the error information and performance were measured in the radiated acoustic field by microphones. Test cases were also performed in which the error sensors were accelerometers mounted on the plate. When accelerometers were used as error sensors, the controller performance was degraded; leading to the conclusion that minimizing plate motion does not necessarily lead to an associated decrease in radiated sound levels. In contrast, the results show excellent attenuation of the transmitted sound field when microphone error sensors were used. This result was consistent over a range of frequencies. In addition, the experimental results are compared to previously derived analytical results and the effect of using a point or global minimization scheme is discussed.

Metcalf, V. L.↗

An approximate Riemann solver for magnetohydrodynamics (that works in more than one dimension)

An approximate Riemann solver is developed for the governing equations of ideal magnetohydrodynamics (MHD). The Riemann solver has an eight-wave structure, where seven of the waves are those used in previous work on upwind schemes for MHD, and the eighth wave is related to the divergence of the magnetic field. The structure of the eighth wave is not immediately obvious from the governing equations as they are usually written, but arises from a modification of the equations that is presented in this paper. The addition of the eighth wave allows multidimensional MHD problems to be solved without the use of staggered grids or a projection scheme, one or the other of which was necessary in previous work on upwind schemes for MHD. A test problem made up of a shock tube with rotated initial conditions is solved to show that the two-dimensional code yields answers consistent with the one-dimensional methods developed previously.

Powell, Kenneth G.↗

Exploring Sustained Food Production in the Veggie Vegetable Production System

Since 2014, the Veggie vegetable production system on the International Space Station has intermittently supplied astronaut crews with fresh produce. To assess the potential for continuous crop production in Veggie and develop a baseline for future space crop production systems, a 120-day study was conducted to determine methodology for inputs, optimal yield, food safety standards, and crew involvement. ‘Amara’ mustard and ‘Red Russian’ kale were grown as initial crops, followed by ‘Extra Dwarf’ pakchoi and shungiku as final crops. Previous grow-outs in Veggie have included harvests at 28-35 days after initiation. In this study, a 56-day grow-out with multiple harvests from the same plants was compared to the conventional, single harvest Veggie schedule. Unlike previous Veggie studies which grew all plants simultaneously, this test staggered initiation and harvest, aiming for consistent and increased production. Plant pillows were initiated in pairs weekly and positioned to reduce shading. Completed pillows were immediately replaced with fresh ones. The multi-harvest scheme used fewer pillows, totaling 46% less pillow mass than the single harvest method. Yield varied by crop and harvest scheme, when compared as fresh edible biomass production across 56-day time increments. ‘Red Russian’ kale yielded similarly across harvest schemes. In the multiple-harvest schedule, ‘Amara’ mustard and shungiku yielded 23% to 25% higher, respectively, while ‘Extra Dwarf’ pakchoi had 43% lower yield. Microbial analysis of the plants indicated no culturable human pathogens. Microbial load of a given plant appeared to depend more on system age than plant age; across harvest methods, aerobic plate counts from final crops were higher than those of initial crops. This project also considered the complexity of crew involvement in a continuous production scenario. New crew procedures that periodically remove plant material from the Veggie root mat would be needed under continuous production to prevent potential pathogens and unpleasant odors. This study supports future space crop production scenarios and was funded by NASA’s Human Research Program.

Jess M Bunchek↗

Sustained Veggie: A Continuous Food Production Comparison

The International Space Station’s Veggie system intermittently supplies the crew with fresh produce. To assess the potential for continuous crop production in Veggie and develop a baseline for future space crop production systems, a 120-day study was conducted to determine methodology for inputs, optimal yield, and crew involvement. ‘Amara’ mustard and ‘Red Russian’ kale were grown as initial crops, followed by ‘Extra Dwarf’ pak choi and shungiku as final crops. Previous grow-outs in Veggie have included harvests at 28-35days after initiation. In this study, a 56-day grow-out with multiple harvests from the same plants was compared to the conventional, single harvest Veggie schedule. Unlike previous Veggie studies which grew all plants simultaneously, this test staggered initiation and harvest, aiming for consistent and increased production. Plant pillows were initiated in pairs weekly and positioned to reduce shading. Completed pillows were immediately replaced with fresh ones. The multi-harvest scheme used fewer pillows, totaling 46% less pillow mass than the single harvest method. Scaled to 56-day increments, yield varied by crop and harvest scheme. ‘Red Russian’ kale yielded similarly across harvest schemes. In the multiple-harvest schedule, ‘Amara’ mustard and shungiku yielded 23% to 25% higher, respectively, while ‘Extra Dwarf’ pak choi had43% lower yield. Microbial analysis of the plants indicated no culturable human pathogens. Microbial load of a given plant appears to depend more on system age than plant age. Across harvest methods, aerobic plate counts from final crops were higher than those of initial crops. This project also considered the complexity of crew involvement in a continuous production scenario. New crew procedures that periodically remove plant material from the Veggie root mat are needed under continuous production to prevent potential pathogens and unpleasant odors. This study supports future space crop production scenarios and was funded by NASA’s Human Research Program.

Jess M Bunchek↗

Updates to Implicit Edge-Based Gradient Methods

In this paper, we report updates to the implicit edge-based gradient methods originally introduced in [H. Nishikawa, AIAA Paper 2020-3048, 2020]. First, we clarify the relationship between gradient accuracy and truncation error and show that the quadratic method involves a free parameter. Then, we provide a complete description and a simplified matrix form of the implicit gradient systems including a consistent boundary treatment, and derive a set of parameters for achieving fourth-order gradient accuracy on regular tetrahedral grids. A stability analysis is performed for a relaxation scheme used to solve the implicit gradient systems, and the result serves as a guide for choosing parameters. Numerical results are shown for accuracy verification and also for realistic inviscid flow problems in three dimensions, including flows with shock waves.

Weighted Least-Squares↗

Advances in Automation of Overset Structured Volume Mesh Generation and Domain Connectivity

Automation of overset structured surface mesh generation has recently been accomplished by the creation of face, edge, and node meshes based on Boundary Representation solids as the geometry input. The current work continues the automation effort in overset volume mesh generation and domain connectivity based on the auto-generated surface meshes. All near-body curvilinear volume meshes are automatically generated using hyperbolic methods. Automation of this step is enabled by appropriate surface grid point distribution, and selection of boundary-splay and smoothing parameters based on concave and convex surface features. The off-body domain is covered by two automatically generated grid systems. The first contains a single Cartesian mesh with a uniform core enclosing all near-body volume meshes and stretched layers that extend to the far field, while the second consists of a set of small stretched Cartesian grids covering pockets of off-body orphan points. With high quality mesh overlap mostly guaranteed by the surface meshing scheme, orphan points that need to be covered by the second Cartesian mesh system are located away from the fine grid spacing region near the wall. Using line-segment and ray-pierce tests against the surface grids, hole-cutting is accomplished on both near and off-body volume grids resulting in appropriate clearances from the wall. The complete mesh generation automation process is demonstrated on five test cases where flow solutions are also computed and compared with solutions obtained using other methods.

TTT↗

Observation of collective modes of excitations in 59 Co, 59 Ni, and 61 Co and the influence of the $g_{9/2}$ orbital

High-spin states in 59 Co (Z = 27), 59 Ni (Z = 28), and 61 Co have been populated by the fusion evaporation reactions, 48 Ti( 14 C, p2n ) 59 Co, 48 Ti( 14 C, 3 n ) 59 Ni, and 50 Ti( 14 C, p2n ) 61 Co. The 9 MV tandem accelerator at the John D. Fox Laboratory, Florida State University (FSU) was used to accelerate the 14 C beam and the deexciting γ rays were detected by the FSU detector array consisting of six high-purity germanium (HPGe) clover detectors and three single crystals. Directional correlation of the $γ$ rays deexciting oriented states (DCO ratios) and polarization asymmetry measurements helped to establish spin and parities of the excited states whenever possible. The level scheme of 59 Co has been expanded with the inclusion of positive-parity states up to 31/2 + at around 11 MeV. The 59 Ni positive-parity states known from previous study were verified with modifications to some of the spins and parities. On the other hand, the negative-parity states were extended to 31/2 at an excitation energy of 12 MeV. No new transition was observed for 61 Co, but one of the major bands has been reassigned as consisting of positive-parity states by reason of this study. Importantly, excitations observed within the f 7/2 , p 3/2 , f 5/2 , and p 1/2 orbitals, and also across the N = 40 subshell closure into the g 9/2 orbital was established by comparison with large-scale shell-model calculations for the three nuclei studied.

59 ≤ A ≤ 89↗

Efficient continuous Energy-Multigroup hybrid depletion scheme using the Shift Monte Carlo code. Part I: Energy condensation sensitivity analysis

Monte Carlo (MC) codes coupled to depletion solvers are increasingly used to provide high fidelity fuel cycle modeling capabilities. Here, these coupled depletion-MC tools produce accurate results in general but can experience nonphysical spatial oscillations when time steps are large or when a system’s dominance ratio approaches unity. Two substepping techniques have been developed previously to remedy and dampen these spatial oscillations without needing to reduce step sizes. The first approach relied on higher-order techniques to account for spectral changes within steps (extrapolation and interpolation techniques). The second approach used the first order perturbation (FOP) theory to account for the change in the one-group spatial flux distribution within steps. This paper develops a hybrid depletion methodology which, in a way, combines how the flux is handled in both substepping techniques. Specifically, the multigroup (MG) MC Shift code is used to update the flux distribution within steps rather than a one-group FOP solver. A fully reflected pincell is investigated, which is not spatially dependent in the MG representation. Thus, the analysis in this paper is an initial demonstration of hybrid depletion. An upcoming companion paper will focus on how the hybrid depletion dampens spatial oscillations. The hybrid depletion approach is verified to be consistent with previous constant extrapolation depletion (CED) methods. This paper finds that the hybrid CED exhibits some error in the eigenvalue and one group constants within macro steps. To address this discrepancy, a simple interpolation scheme (CELI) is investigated. This work found that CELI sufficiently addresses the discrepancy in spectrum for macro steps up to 100 days. Overall, this work demonstrates that the hybrid depletion method can significantly reduce the number of high fidelity MC executions in a MC-coupled depletion with an acceptable eigenvalue error.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A simple parallel prefix algorithm for compact finite-difference schemes

A compact scheme is a discretization scheme that is advantageous in obtaining highly accurate solutions. However, the resulting systems from compact schemes are tridiagonal systems that are difficult to solve efficiently on parallel computers. Considering the almost symmetric Toeplitz structure, a parallel algorithm, simple parallel prefix (SPP), is proposed. The SPP algorithm requires less memory than the conventional LU decomposition and is highly efficient on parallel machines. It consists of a prefix communication pattern and AXPY operations. Both the computation and the communication can be truncated without degrading the accuracy when the system is diagonally dominant. A formal accuracy study was conducted to provide a simple truncation formula. Experimental results were measured on a MasPar MP-1 SIMD machine and on a Cray 2 vector machine. Experimental results show that the simple parallel prefix algorithm is a good algorithm for the compact scheme on high-performance computers.

Sun, Xian-He↗

The Tiny Triplet Finder as a Versatile Track Segment Seeding Engine for Trigger Systems

In high energy physics experiment trigger systems, track segment seeding is a resource consuming function and the primary reason is the computing complexity of the segment finding process. As the Moore's Law is reaching its physical limit, reducing computing complexity should be carefully considered, rather than keep piling up silicon resources. The Tiny Triplet Finder is a scheme that reduces the computing complexity of the segment seeding. As a proof of concept, a 3D track segment seeding engine core based on the Tiny Triplet Finder has been implemented and tested in a low-cost FPGA device. The seeding engine is designed to preselect and group hits (stubs) from detector layers to feed subsequent track fitting stage. The seeding engine consists of a Hough transform space for r-z view and a Tiny Triplet Finder for r-phi view to implement 3D constraints. The seeding engine is organized as a pipeline so that each hit is processed in a single clock cycle. Taking advantage of the register-like storage block scheme which enables effectively resetting of a block RAM within a single clock cycle, clearing or refreshing the seeding engine takes only one clock cycles between two events. The Tiny Triplet Finder is also a generic coincidence finding scheme that can be used for many tasks. As a versatility demonstration, track segment finding performances for two distinctive detector geometries are tested in our seeding engine. In a collider barrel-layer geometry, the fake segment rates are studied for 3D (i.e., both r-phi and r-z views) and 2D (i.e., r-phi or r-z view only) configurations for high hit multiplicity events (>4000 hits/layer in the barrel region). Another detector geometry contains strip plane layers with timing information. The numbers of coincidences, both real or fake, with or without timing ("3D" or "2D") information at various hit multiplicities are studied.

43 PARTICLE ACCELERATORS↗

Learning with Delayed Rewards—A Case Study on Inverse Defect Design in 2D Materials

Defect dynamics in materials are of central importance to a broad range of technologies from catalysis to energy storage systems to microelectronics. Material functionality depends strongly on the nature and organization of defects-their arrangements often involve intermediate or transient states that present a high barrier for transformation. The lack of knowledge of these intermediate states and the presence of this energy barrier presents a serious challenge for inverse defect design, especially for gradient-based approaches. Here, we present a reinforcement learning (RL) [Monte Carlo Tree Search (MCTS)] based on delayed rewards that allow for efficient search of the defect configurational space and allows us to identify optimal defect arrangements in low-dimensional materials. Using a representative case of two-dimensional MoS 2 , we demonstrate that the use of delayed rewards allows us to efficiently sample the defect configurational space and overcome the energy barrier for a wide range of defect concentrations (from 1.5 to 8% S vacancies)-the system evolves from an initial randomly distributed S vacancies to one with extended S line defects consistent with previous experimental studies. Detailed analysis in the feature space allows us to identify the optimal pathways for this defect transformation and arrangement. Comparison with other global optimization schemes like genetic algorithms suggests that the MCTS with delayed rewards takes fewer evaluations and arrives at a better quality of the solution. The implications of the various sampled defect configurations on the 2H to 1T phase transitions in MoS 2 are discussed. In this study, we introduce a RL strategy employing delayed rewards that can accelerate the inverse design of defects in materials for achieving targeted functionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparing event generator predictions and ab initio calculations of ν- 12 C neutral-current quasielastic scattering at 1 GeV

The measurement of neutrino oscillations and exotic physics searches are important parts of the physics program in the near future, with new state-of-the-art experiments planned within the next decade. Future and modern experiments in these fields will make use of nuclear targets. Event generators (EGs) are software used in the analysis of neutrino oscillation experiments. EGs are used to predict kinematic observables for a range of neutrino energies. These simulations make use of simple models of nucleon dynamics. As such, they may fail to capture features of more rigorous theoretical calculations. This work compares EG performance to nuclear theory calculations by comparing observables generated in the two frameworks. We provide a common set of definitions between theory and experiment and assess the adequacy of the implementation of two nuclear physics models in EG simulations. Neutral-current quasielastic (NCQE) scattering events for neutrinos and antineutrinos on a 12 C target are simulated with a specific EG, NEUT, used by the T2K experiment for its analysis. The simulated cross sections are compared to analytic calculations from nuclear theory within the factorization scheme. We compare the NEUT implementation of two different models on nuclear spectral functions: the relativistic Fermi gas (RFG) and the correlated basis spectral function (CBF) to analytic calculations of the same models in the factorization scheme. For both nuclear physics models, we compare the appearance of features in the distributions relevant to experimental analyses. The peak of the cross section dσ/(dΩ dω) is consistent in energy transfer, ω, for RFG and CBF simulations. Qualitatively, the shape of the simulated distribution is similar to the one obtained through theory calculations; however, there are some discrepancies between the theory calculations and the NEUT simulation. While the EG simulations and analytic calculations with the same model of nuclear dynamics show similar overall features, there are still differences between the two. These results demonstrate the importance of benchmarking EGs so their nuclear physics implementations can be improved for the analysis of future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Onset of band structure in 70 Ga

Excited states in the odd-odd nucleus 70 Ga were studied using the 62 Ni( 14 C, αpn) fusion-evaporation reaction at the John D. Fox Superconducting Accelerator Laboratory at Florida State University with a 50 MeV beam. Depopulating γ rays were measured in coincidence with a Compton-suppressed Ge array consisting of three Clover detectors and seven single-crystal detectors. Analysis of γ–γ coincidences led to the addition of 16 new transitions to the 70 Ga level scheme, including transitions associated with the onset of a new positive-parity band structure likely based on the πg 9/2 ⊗ νg 9/2 configuration. Spins and parities were assigned using directional correlation of oriented nuclei (DCO) ratios and linear polarization measurements. Shell-model calculations with the JUN45 effective interaction reproduce the excitation energies of the positive-parity states reasonably well but predict negative-parity states that are typically 400–500 keV lower than observed. Furthermore, total Routhian surface calculations for the lowest positive-parity configuration with signature α = 0 indicate significant triaxiality (γ ≈ −20°) and the development of competing non-collective excitations at J = 8, corresponding to the experimentally observed onset of the positive-parity band. Similar features are found in the calculated surfaces for the lowest negative-parity states with α = 0.

59 ≤ A ≤ 89↗

Strictly Enforcing Invertibility and Conservation in CNN-Based Super Resolution for Scientific Datasets

Abstract Recently, deep convolutional neural networks (CNNs) have revolutionized image “super resolution” (SR), dramatically outperforming past methods for enhancing image resolution. They could be a boon for the many scientific fields that involve imaging or any regularly gridded datasets: satellite remote sensing, radar meteorology, medical imaging, numerical modeling, and so on. Unfortunately, while SR-CNNs produce visually compelling results, they do not necessarily conserve physical quantities between their low-resolution inputs and high-resolution outputs when applied to scientific datasets. Here, a method for “downsampling enforcement” in SR-CNNs is proposed. A differentiable operator is derived that, when applied as the final transfer function of a CNN, ensures the high-resolution outputs exactly reproduce the low-resolution inputs under 2D-average downsampling while improving performance of the SR schemes. The method is demonstrated across seven modern CNN-based SR schemes on several benchmark image datasets, and applications to weather radar, satellite imager, and climate model data are shown. The approach improves training time and performance while ensuring physical consistency between the super-resolved and low-resolution data. Significance Statement Recent advancements in using deep learning to increase the resolution of images have substantial potential across the many scientific fields that use images and image-like data. Most image super-resolution research has focused on the visual quality of outputs, however, and is not necessarily well suited for use with scientific data where known physics constraints may need to be enforced. Here, we introduce a method to modify existing deep neural network architectures so that they strictly conserve physical quantities in the input field when “super resolving” scientific data and find that the method can improve performance across a wide range of datasets and neural networks. Integration of known physics and adherence to established physical constraints into deep neural networks will be a critical step before their potential can be fully realized in the physical sciences.

54 ENVIRONMENTAL SCIENCES↗