SEARCH · Search NASA
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.
The investigation reported is concerned with the construction of the resolvent for any given kernel function. In problems with ill-behaved inhomogeneous terms as, for instance, in the aerodynamic problem of flow over a flapped airfoil, direct numerical methods become very difficult. A description is presented of a solution method by resolvent which can be employed in such problems.
QMCPACK is open source scientific software designed to perform Quantum Monte Carlo (QMC) simulation, a first-principles method for describing many-body systems. The evaluation of each Monte Carlo move requires finding the determinant of a dense matrix in the wave functions. This calculation forms a key computational kernel in QMCPACK. After each accepted event, the wave function matrix undergoes a rank-one update to represent a single particle move within the system. The matrix inverse is updated via the Sherman-Morrison formula; occasionally, the explicit inverse must be recomputed to maintain numerical stability.; We propose an alternate approach to this kernel that utilizes QR factorization to maintain stability without refactorization. In addition, algorithms based on a novel delayed update scheme are explored in this effort. This strategy enables probability evaluation for multiple successive Monte Carlo moves, with application of accepted moves to the wave function matrix delayed until an event is denied, or a predetermined limit on acceptances p is reached. Accepted events grouped in this manner are then applied to the matrix en bloc with enhanced arithmetic intensity and computational efficiency. Delayed update Sherman-Morrison probability evaluation algorithms with CPU and GPU-accelerated implementations are tested, profiled, and analyzed. A QR based delayed update algorithm remains in development. Results are evaluated against existing methods for numerical stability and efficiency; emphasis is placed on large systems, for which acceleration is critical.
We compare various radiative feedbacks in response to short-term climate variations estimated from reanalysis, satellite, and global climate model datasets using the combined Kernel-Gregory approach over the Arctic (60-90°N). Except for the shortwave cloud feedback, all feedbacks are positive over the Arctic. Lapse rate feedback is comparable with surface albedo feedback. Relative to the tropics (30°S-30°N), the temperature feedback (lapse rate feedback plus Planck feedback deviation from its global mean) is the largest contributor to Arctic amplification among all feedbacks, followed by surface albedo feedback. The choice of meteorological fields has a small impact on the cloud feedback estimate, but different top-of-the-atmosphere clear-sky fluxes datasets can result in a wide spread of estimated cloud feedbacks. The net cloud feedback over the Arctic is positive due to a larger positive in longwave than the negative shortwave feedback, but has the largest uncertainty. Relative to the tropics, the longwave cloud feedback largely contributes to the Arctic amplification.
Real-time computing has traditionally been considered largely in the context of single-processor and embedded systems, and indeed, the terms real-time computing, embedded systems, and control systems are often mentioned in closely related contexts. However, real-time computing in the context of multinode systems, specifically high-performance, cluster-computing systems, remains relatively unexplored. Imposing real-time constraints on a parallel (cluster) computing environment introduces a variety of challenges with respect to the formal verification of the system's timing properties. In this paper, we give a motivating example to demonstrate the need for such a system--- an application to estimate the electromechanical states of the power grid--- and we introduce a formal method for performing verification of certain temporal properties within a system of parallel processes. We describe our work towards a full real-time implementation of the target application--- namely, our progress towards extracting a key mathematical kernel from the application, the formal process by which we analyze the intricate timing behavior of the processes on the cluster, as well as timing measurements taken on our test cluster to demonstrate use of these concepts.
The U.S. Department of Energy Office of Nuclear Energy (DOE NE) and the Idaho National Laboratory (INL) Advanced Reactor Technologies (ART) Advanced Gas Reactor (AGR) Fuel Development and Qualification program (referred to as AGR Fuel program hereafter) are pursuing qualification of tristructural isotropic (TRISO) coated particle fuel for use in high temperature gas cooled reactors (HTGRs). The AGR Fuel program was established to provide a fuel qualification data set in support of the licensing and operation of an HTGR. BWX Technologies Nuclear Operations Group (BWXT-NOG) was subcontracted to fabricate the fuel for the AGR program. Several investments and innovations were realized in preparation to fabricate fuel for the AGR-5/6/7 irradiation experiments that brought fuel fabrication fully out of the laboratory and into engineering-scale operations. These included: • Increased the kernel fabrication line capacity and uniformity • Upgraded ancillary support equipment and processes for the tristructural isotropic (TRISO) coating furnace • Demonstrated efficient production of the matrix precursor powder by dry jet milling of co mingled components • Demonstrated an engineering-scale method for quick and efficient overcoating TRISO particles with the matrix precursor • Demonstrated an automated, multi cavity compacting system with a volumetric feed system • Demonstrated a combined-cycle thermal treatment furnace These changes increased production rates of some of these processes by an order of magnitude or more while eliminating the use of flammable solvents, multiple grinding and sorting operations, and the weighing out of individual die charges. Three fuel kernel lots were fabricated for production of the fuel for AGR-5/6/7. The initial lot (J52R-16-39316) was certified to fuel specifications but was not used because the kernels had a high fraction of internal fissures that caused an unacceptable fraction of the kernels to fragment when charged to the coating furnace where the TRISO coating would be deposited. Fragmented kernels increased the dispersed uranium in the particles and produced a worrisome fraction of dimpled particles with an elevated probability of in-pile failure. After some efforts to identify the cause of the fissure formation, two additional lots were produced with much lower fissure fractions, J52R-16-69317 and 69318. The latter kernel lot was a backup to the first and was not needed. Multiple kernel batches were composited to form each of the lots so as to simulate a commercial-scale operation where kernel batches would also be composited. Multiple TRISO coating runs were performed and the product characterized so that several could be composited into a TRISO lot. TRISO lot J52R-16-98005 conformed to all fuel specifications except the mean outer pyrocarbon (OPyC) layer thickness was thinner than specified. Furthermore, it was determined that the TRISO lot had a dispersed uranium fraction (DUF) that might result in the compacts not meeting the DUF specification. A review of the role of the OPyC layer and consequences of the DUF by the Technical Coordination Team and INL concluded that the fuel was acceptable for use in the AGR-5/6/7 irradiation experiment. The TRISO particles were overcoated with the matrix precursor that had been produced in a jet-milling operation. The overcoating was performed in equipment originally designed to coat pharmaceuticals. The overcoating process performed well; producing highly spherical and uniform overcoats requiring little upgrading and no recycle or rework. TRISO particles were overcoated with the matrix precursor to achieve nominal volumetric packing fractions (PFs) of TRISO particles of 25% and 40% for the irradiation experiments. The 40% PF compacts occupy the first and fifth test capsule in the test train while the inner three capsules are loaded with 25% PF compacts. The resinated graphite matrix precursor powder was a derivative of the German A3-27 matrix formulation, which differs from previous AGR irradiation campaigns that used an A3-3 formulation. Jet milling of the matrix powder precursor produced a finer mean graphite particle size than the milling operations used for the A3-3 matrix powder precursor. Changes made in the matrix formula and equipment yielded compacts with significantly higher matrix density than was attained in previous AGR irradiation campaigns. The changes in the matrix formulation and the means of milling the powders also complicated resolution of the three fuel compact defect fractions, DUF, exposed kernel fraction (EKF), and the silicon carbide defect fraction (SDF). Characterization data from BWXT-NOG had some anomalous results, so samples of the fuel compacts and overcoated TRISO particles were also analyzed by Oak Ridge National Laboratory (ORNL) to ensure that the defect fractions were accurately characterized.
Simplified treatment of resonance elastic scattering model considering thermal motion of heavy nuclides and the energy dependence of the resonance cross section was implemented into NJOY [1]. In order to solve deterministic slowing down equation considering the effect of up-scattering without iterative calculations, scattering kernel for heavy nuclides is pre-calculated by the formula derived by Ouisloumen and Sanchez [2], and neutron spectrum in up-scattering term is expressed by NR approximation. To check the verification of the simplified treatment, the treatment is applied to U-238 for the energy range from 4 eV to 200 eV. Calculated multi-group capture cross section of U-238 is greater than that of conventional method and the increase of the capture cross sections is remarkable as the temperature becomes high. Therefore Doppler coefficient calculated in UO{sub 2} fuel pin is calculated more negative value than that on conventional method. The impact on Doppler coefficient is equivalent to the results of exact treatment of resonance elastic scattering reported in previous studies [2-7]. The agreement supports the validation of the simplified treatment and therefore this treatment is applied for other heavy nuclide to evaluate the Doppler coefficient in MOX fuel. The result shows that the impact of considering thermal agitation in resonance scattering in Doppler coefficient comes mainly from U-238 and that of other heavy nuclides such as Pu-239, 240 etc. is not comparable in MOX fuel. (authors)
With this study, our goal is compression of massive-scale grid-structured data, such as the multi-terabyte output of a high-fidelity computational simulation. For such data sets, we have developed a new software package called TuckerMPI, a parallel C++/MPI software package for compressing distributed data. The approach is based on treating the data as a tensor, i.e., a multidimensional array, and computing its truncated Tucker decomposition, a higher-order analogue to the truncated singular value decomposition of a matrix. The result is a low-rank approximation of the original tensor-structured data. Compression efficiency is achieved by detecting latent global structure within the data, which we contrast to most compression methods that are focused on local structure. In this work, we describe TuckerMPI, our implementation of the truncated Tucker decomposition, including details of the data distribution and in-memory layouts, the parallel and serial implementations of the key kernels, and analysis of the storage, communication, and computational costs. We test the software on 4.5 and 6.7 terabyte data sets distributed across 100 s of nodes (1,000 s of MPI processes), achieving compression ratios between 100 and 200,000×, which equates to 99--99.999% compression (depending on the desired accuracy) in substantially less time than it would take to even read the same dataset from a parallel file system. Moreover, we show that our method also allows for reconstruction of partial or down-sampled data on a single node, without a parallel computer so long as the reconstructed portion is small enough to fit on a single machine, e.g., in the instance of reconstructing/visualizing a single down-sampled time step or computing summary statistics. The code is available at https://gitlab.com/tensors/TuckerMPI.
Purpose: Commissioning data for intensity modulated proton therapy (IMPT) must be post-processed by fits to ad-hoc functions to derive the dose calculation kernel parameters in a treatment planning system (TPS). Whether from experimental measurement or Monte Carlo simulation, the limited and noisy nature of such data makes this task very challenging. We present a method to improve the modeling of the lateral dose distribution of clinical energy proton beams in water to commission an in-house IMPT dose calculation engine. Methods: A linear sum of three Gaussian distribution functions was fitted to the lateral dose data in logarithmic scale. Starting values of fitting solutions were determined from the Generalized Highland Approximation. We exhaustively optimized the combinations of data weights with upper bounds of the fitting solutions to minimize confidence intervals of the fitting solutions while maintaining the coefficient of determination (R{sup 2}). Results: Across all energies, average confidence bounds improved 72.88% [Max: 88.28%, Min: 55.05%] for small angle coulomb scattering, 114.25% [409.13%, 66.72%,] for nuclear scattering, and 68.66% [141.09%, 33.27%] for large angle coulomb scattering, while the coefficients of determination of the fits (R{sup 2}) remained comparable. On average R {sup 2} only changed 0.18% and were very close to 1 (approx. 0.999). Wilcoxon signed rank tests comparing unweighted/unbounded fits with weighted/bounded fits averaged 0.0146 (Max: 0.177, Min: 7.05×10−{sup 7}) for small angle Coulomb, 0.0903 (0.945, 7.05×10−{sup 7}) for nuclear, and 0.254 (0.871, 1.86×10−{sup 6}) for large angle Coulomb scattering. This allows rejection of the null hypothesis for small angle Coulomb scattering at the 0.015 level and nuclear interaction at the 0.1 level. Conclusion: Optimal weights assigned to IMPT lateral dose data minimized fitting to stochastic noise in the tail region. Optimizing the upper bounds of fitting parameters improved the robustness of fitting solutions. Both led to a more physically meaningful dose calculation kernel for the TPS. NIH/NCI K25CA168984, Eagles Cancer Research Career Development, The Lawrence W. and Marilyn W. Matteson Fund for Cancer Research, Mayo ASU Seed Grant, and The Kemper Marley Foundation.
We compute the energy loss of heavy fermions moving in a plasma, taking into account the modification of the photon collective modes induced by collisions using a Bhatnagar-Gross–Krook collisional kernel. We include contributions from both hard and soft scatterings of the heavy fermion using a collisionally modified hard-thermal-loop resummed propagator. Using this method, one does not need to introduce a separation scale between hard- and soft-momentum exchanges. To place our calculation in context, we review other theoretical approaches to computing the collisional energy loss of fermions and discuss the systematics and results obtained in each approach compared to using a resummed propagator for both hard and soft momentum exchanges. Our final results indicate that self-consistently including the effect of collisions in the self-energies of the resummed propagator results in an increased energy loss compared to using collisionless hard-thermal-loop propagators. The effect becomes larger as the magnitude of the coupling constant and the velocity of the fermion increase. Published by the American Physical Society 2024
A brief review of Stokes' problem for the ellipsoid as a reference surface is given. Another solution of the problem using an ellipsoidal kernel, which represents an iterative form of Stokes' integral, is suggested with a relative error of the order of the flattening. On studying of Rapp's method in detail the procedures of improving its convergence are discussed.
We introduce a new tool to planetary geology for quantifying the spatial arrangement of vent fields and volcanic provinces using non parametric kernel density estimation. Unlike parametricmethods where spatial density, and thus the spatial arrangement of volcanic vents, is simplified to fit a standard statistical distribution, non parametric methods offer more objective and data driven techniques to characterize volcanic vent fields. This method is applied to Syria Planum volcanic vent catalog data as well as catalog data for a vent field south of Pavonis Mons. The spatial densities are compared to terrestrial volcanic fields.
WEC-Sim is a midfidelity numerical tool for modeling wave energy conversion devices. The code uses the MATLAB SimMechanics package to solve multibody dynamics and models wave interactions using hydrodynamic coefficients derived from frequency-domain boundary-element methods. This paper presents the new modeling features introduced in the latest release of WEC-Sim. The first feature discussed conversion of the fluid memory kernel to a state-space form. This enhancement offers a substantial computational benefit after the hydrodynamic body-to-body coefficients are introduced and the number of interactions increases exponentially with each additional body. Additional features include the ability to calculate the wave-excitation forces based on the instantaneous incident wave angle, allowing the device to weathervane, as well as import a user-defined wave elevation time series. A review of the hydrodynamic theory for each feature is provided and the successful implementation is verified using test cases.
Analysis and interpretation of time series data is easiest when the data values occur at uniform intervals in time, but actual data may have differing data sampling frequencies, such as monthly and daily readings. Applying data analysis techniques, such as smoothing, to such a data set may not give a representative result between time segments. The ability to automatically distinguish time segments of differing data frequency would provide a means for applying data analysis independently to each segment, though a suitable blending at segment boundaries would be required. A method for detecting frequency changes was developed and applied to Gaussian and median smoothing of hydraulic head data from groundwater wells at the U.S. Department of Energy Hanford Site in southeastern Washington state. The process identifies time segments of high-frequency (daily) or low-frequency (greater than daily) data using adjusted-bandwidth Gaussian kernel density estimation and a threshold value, which are further refined to address small blocks of low-frequency data within larger blocks of high-frequency data. User-selectable levels of smoothing are then applied independently to the time segments prior to combining the segment results for a single smoothed data set. This time segment identification approach provides effective low- and high-frequency data separation, which provides a method to apply data analysis independently to each time segment.
Purpose: To derive and introduce a new correction factor kV, the “volume effect correction factor”, that accounts for not only the dose averaging over the detector's sensitive volume but also the secondary electron generation and transport inclusive of the disturbance of the field of secondary electrons within the detector. Materials and Methods: Mathematical convolutions and Fourier's convolution theorem have been used. Monte Carlo simulations of photon pencil beams were performed using EGSnrc. Detector constructions were adapted from manufacturers' information. Results: For the calculation of kV, the three basic convolution kernels have to be taken into account: the dose deposition kernel KD(x) (fluence to dose), the photon fluence response kernel KM(x) (photon fluence to detector signal) and the “dose response kernel” K(x) (dose to detector signal). K(x) is calculated from FT[K(x)] = [1/sqrt(2”)]FT[KM(x)]/FT[KD(x)], where the magnitude of kV(x) can be thereby calculated for arbitrary photon beam profiles and the areanormalized K(x). Conclusions: n order to take into account for the dimensions of dosimetric detectors in narrow photon beams, the “volume effect correction factor” kV has been introduced into the fundamental equation of probe-type dosimetry, and the convolution method has proven to be a method for the derivation of its numerical values. For narrow photon beams, whose width is comparable to the secondary electron ranges, kV can reach very high values, but it can be shown that the signals of small diamond detectors are well representing the absorbed dose to water averaged over the detector volume.
Highlights: > We study crossover ensembles of Jacobi family of random matrices. > We consider correlations for orthogonal-unitary and symplectic-unitary crossovers. > We use the method of skew-orthogonal polynomials and quaternion determinants. > We prove universality of spectral correlations in crossover ensembles. > We discuss applications to quantum conductance and communication theory problems. - Abstract: In a recent paper (S. Kumar, A. Pandey, Phys. Rev. E, 79, 2009, p. 026211) we considered Jacobi family (including Laguerre and Gaussian cases) of random matrix ensembles and reported exact solutions of crossover problems involving time-reversal symmetry breaking. In the present paper we give details of the work. We start with Dyson's Brownian motion description of random matrix ensembles and obtain universal hierarchic relations among the unfolded correlation functions. For arbitrary dimensions we derive the joint probability density (jpd) of eigenvalues for all transitions leading to unitary ensembles as equilibrium ensembles. We focus on the orthogonal-unitary and symplectic-unitary crossovers and give generic expressions for jpd of eigenvalues, two-point kernels and n-level correlation functions. This involves generalization of the theory of skew-orthogonal polynomials to crossover ensembles. We also consider crossovers in the circular ensembles to show the generality of our method. In the large dimensionality limit, correlations in spectra with arbitrary initial density are shown to be universal when expressed in terms of a rescaled symmetry breaking parameter. Applications of our crossover results to communication theory and quantum conductance problems are also briefly discussed.
In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + integral(0 to 1) g(x,t) F(t, y(t)) dt, 0 less than or equal to x less than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integrals equations arise, e.g., when one applies Green's function techniques to nonlinear two-point boundary value problems of the form U''(x) = f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and g(l) = y(sub 1), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trapezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.
In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.
The cosmic structure formed from baryon acoustic oscillations (BAO) in the early universe is imprinted in the galaxy distribution observable in large-scale surveys and is used as a standard ruler in contemporary cosmology. Typically, BAOs are detected as a preferential length scale in two-point statistics, which gives little information about the location of the BAO structures in real space. The aim of the algorithm described in this paper is to find probable centers of BAOs in the cosmic matter distribution. The algorithm convolves the three-dimensional distribution of matter density with a spherical shell kernel of variable radius placed at different locations. The locations that correspond to the highest values of the convolution correspond to the probable centers of BAOs. This method is realized in an open-source, computationally efficient algorithm. Here, we describe the algorithm and present the results of applying it to the SDSS DR9 CMASS survey and associated mock catalogs. A detailed performance study demonstrates the ability of the algorithm to locate BAO centers and in doing so presents a novel detection of the BAO scale in galaxy surveys.