Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sum reduction”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Unitarization from geometry

We study the perturbative unitarity of scattering amplitudes in general dimensional reductions of Yang–Mills theories and general relativity on closed internal manifolds. For the tree amplitudes of the dimensionally reduced theory to have the expected high-energy behavior of the higher-dimensional theory, the masses and cubic couplings of the Kaluza–Klein states must satisfy certain sum rules that ensure there are nontrivial cancellations between Feynman diagrams. These sum rules give constraints on the spectra and triple overlap integrals of eigenfunctions of Laplacian operators on the internal manifold and can be proven directly using Hodge and eigenfunction decompositions. One consequence of these constraints is that there is an upper bound on the ratio of consecutive eigenvalues of the scalar Laplacian on closed Ricci-flat manifolds with special holonomy. This gives a sharp bound on the allowed gaps between Kaluza–Klein excitations of the graviton that also applies to Calabi–Yau compactifications of string theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Modification of ion-temperature-gradient turbulence by impurities in stellarator plasmas

Recent nonlinear gyrokinetic simulations have shown that impurities can strongly modify the turbulent heat flux in stellarator plasmas. Here, the ion-temperature-gradient (ITG) dispersion relation in a plasma containing impurities is analytically solved in certain limits and an expression for the modification of the ITG growth rate by impurities is derived. The analytical expression is the sum of three terms corresponding to three different physical causes (impurity density gradient, impurity temperature gradient and dilution) of the change in the growth rate. The scalings predicted analytically for the modification of the growth rate are shown to be reproduced by linear gyrokinetic simulations. The conditions for reduction or increase of the ITG growth by impurities are also correctly predicted by the analytical solution to the dispersion relation. Finally, a remarkable correlation is found between the analytical expression for the modification of the growth rate and the modification of the turbulent heat flux obtained from nonlinear gyrokinetic simulations.

Calvo, Iván [Research Centre for Energy, Environme↗
DOE OSTI2017

Extreme-Scale Algorithms & Software Resilience (EASIR) Architecture-Aware Algorithms for Scalable Performance and Resilience on Heterogeneous Architectures

This project addresses both communication-avoiding algorithms, and reproducible floating-point computation. Communication, i.e. moving data, either between levels of memory or processors over a network, is much more expensive per operation than arithmetic (measured in time or energy), so we seek algorithms that greatly reduce communication. We developed many new algorithms for both dense and sparse, and both direct and iterative linear algebra, attaining new communication lower bounds, and getting large speedups in many cases. We also extended this work in several ways: (1) We minimize writes separately from reads, since writes may be much more expensive than reads on emerging memory technologies, like Flash, sometimes doing asymptotically fewer writes than reads. (2) We extend the lower bounds and optimal algorithms to arbitrary algorithms that may be expressed as perfectly nested loops accessing arrays, where the array subscripts may be arbitrary affine functions of the loop indices (eg A(i), B(i,j+k, k+3*m-7, …) etc.). (3) We extend our communication-avoiding approach to some machine learning algorithms, such as support vector machines. This work has won a number of awards. We also address reproducible floating-point computation. We define reproducibility to mean getting bitwise identical results from multiple runs of the same program, perhaps with different hardware resources or other changes that should ideally not change the answer. Many users depend on reproducibility for debugging or correctness. However, dynamic scheduling of parallel computing resources, combined with nonassociativity of floating point addition, makes attaining reproducibility a challenge even for simple operations like summing a vector of numbers, or more complicated operations like the Basic Linear Algebra Subprograms (BLAS). We describe an algorithm that computes a reproducible sum of floating point numbers, independent of the order of summation. The algorithm depends only on a subset of the IEEE Floating Point Standard 754-2008, uses just 6 words to represent a “reproducible accumulator,” and requires just one read-only pass over the data, or one reduction in parallel. New instructions based on this work are being considered for inclusion in the future IEEE 754-2018 floating-point standard, and new reproducible BLAS are being considered for the next version of the BLAS standard.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Time Resolved Spectroscopy of Eclipsing Polars

No changes have been made since the last annual progress report was submitted in conjunction with a unilateral NCX. Dr. Barrett was affected by an STScI Reduction in Force (RIF). He is now employed by the Johns Hopkins University and plans to continue his research there. No expenses have been charged to this grant, however the FUSE data for the eclipsing polar V1432 Aql has been received and processed using CALFWSE v3.0.6. The resulting summed spectrum has been used for a preliminary analysis of the interstellar absorption towards V1432 Aql. We find a hydrogen column density of less than 1.5e21 cm^-2. We have used this result in the paper "X-Ray Emission and Optical Polarization of V1432 Aquilae: An Asynchronous Polar" to fix the hydrogen column density in the soft (<0.5 keV) X-ray band when analyzing the XMM-Newton spectra of this polar. This has enabled us to find an accurate temperature for the blackbody component of 88+/-2 eV, which is significantly higher than that of other polars (20 - 40 eV). We hope to complete our analysis of the phase-resolved emission line spectra of V1432 Aql and to prepare the results for publication in a refereed journal. We hope to begin work on this star within the next few months.

Barrett, Paul↗
DOE OSTI2017

Simultaneous 3D MR elastography of the in vivo mouse brain

The feasibility of sample interval modulation (SLIM) magnetic resonance elastography (MRE) for the in vivo mouse brain is assessed, and an alternative SLIM-MRE encoding method is introduced. In SLIMMRE, the phase accumulation for each motion direction is encoded simultaneously by varying either the start time of the motion encoding gradient (MEG), SLIM-phase constant (SLIM-PC), or the initial phase of the MEG, SLIM-phase varying (SLIM-PV). SLIM-PC provides gradient moment nulling, but the mutual gradient shift necessitates increased echo time (TE). SLIM-PV requires no increased TE, but exhibits nonuniform flow compensation. Comparison was to conventional MRE using six C57BL/6 mice. For SLIMPC, the Spearman’s rank correlation to conventional MRE for the shear storage and loss modulus images were 80% and 76%, respectively, and likewise for SLIM-PV, 73% and 69%, respectively. The results of the Wilcoxon rank sum test showed that there were no statistically significant differences between the spatially averaged shear moduli derived from conventional-MRE, SLIM-PC, and SLIM-PV acquisitions. Both SLIM approaches were comparable to conventional MRE scans with Spearman’s rank correlation of 69%-80% and with 3 times reduction in scan time. As a result, the SLIM-PC method had the best correlation, and SLIM-PV may be a useful tool in experimental conditions, where both measurement time and T2 relaxation is critical.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1987

Algorithm Sorts Groups Of Data

For efficient sorting, algorithm finds set containing minimum or maximum most significant data. Sets of data sorted as desired. Sorting process simplified by reduction of each multielement set of data to single representative number. First, each set of data expressed as polynomial with suitably chosen base, using elements of set as coefficients. Most significant element placed in term containing largest exponent. Base selected by examining range in value of data elements. Resulting series summed to yield single representative number. Numbers easily sorted, and each such number converted back to original set of data by successive division. Program written in BASIC.

Evans, J. D.↗
NASA NTRS1993

A Spectral Algorithm for Envelope Reduction of Sparse Matrices

The problem of reordering a sparse symmetric matrix to reduce its envelope size is considered. A new spectral algorithm for computing an envelope-reducing reordering is obtained by associating a Laplacian matrix with the given matrix and then sorting the components of a specified eigenvector of the Laplacian. This Laplacian eigenvector solves a continuous relaxation of a discrete problem related to envelope minimization called the minimum 2-sum problem. The permutation vector computed by the spectral algorithm is a closest permutation vector to the specified Laplacian eigenvector. Numerical results show that the new reordering algorithm usually computes smaller envelope sizes than those obtained from the current standard algorithms such as Gibbs-Poole-Stockmeyer (GPS) or SPARSPAK reverse Cuthill-McKee (RCM), in some cases reducing the envelope by more than a factor of two.

Barnard, Stephen T.↗
DOE OSTI2014

MO-G-17A-07: Improved Image Quality in Brain F-18 FDG PET Using Penalized-Likelihood Image Reconstruction Via a Generalized Preconditioned Alternating Projection Algorithm: The First Patient Results

Purpose: To investigate the performance of a new penalized-likelihood PET image reconstruction algorithm using the 1{sub 1}-norm total-variation (TV) sum of the 1st through 4th-order gradients as the penalty. Simulated and brain patient data sets were analyzed. Methods: This work represents an extension of the preconditioned alternating projection algorithm (PAPA) for emission-computed tomography. In this new generalized algorithm (GPAPA), the penalty term is expanded to allow multiple components, in this case the sum of the 1st to 4th order gradients, to reduce artificial piece-wise constant regions (“staircase” artifacts typical for TV) seen in PAPA images penalized with only the 1st order gradient. Simulated data were used to test for “staircase” artifacts and to optimize the penalty hyper-parameter in the root-mean-squared error (RMSE) sense. Patient FDG brain scans were acquired on a GE D690 PET/CT (370 MBq at 1-hour post-injection for 10 minutes) in time-of-flight mode and in all cases were reconstructed using resolution recovery projectors. GPAPA images were compared PAPA and RMSE-optimally filtered OSEM (fully converged) in simulations and to clinical OSEM reconstructions (3 iterations, 32 subsets) with 2.6 mm XYGaussian and standard 3-point axial smoothing post-filters. Results: The results from the simulated data show a significant reduction in the 'staircase' artifact for GPAPA compared to PAPA and lower RMSE (up to 35%) compared to optimally filtered OSEM. A simple power-law relationship between the RMSE-optimal hyper-parameters and the noise equivalent counts (NEC) per voxel is revealed. Qualitatively, the patient images appear much sharper and with less noise than standard clinical images. The convergence rate is similar to OSEM. Conclusions: GPAPA reconstructions using the 1{sub 1}-norm total-variation sum of the 1st through 4th-order gradients as the penalty show great promise for the improvement of image quality over that currently achieved with clinical OSEM reconstructions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Towards standardized nuclear reactors: Seismic isolation and the cost impact of the earthquake load case

We report nuclear energy has a key role to play in global decarbonization. Impediments to the widespread deployment of reactors are their projected high capital cost and levelized cost of energy, and time required to analyze, design, license, construct, and commission them. The earthquake load case is a key cost driver for a new build nuclear plant, because near-surface soils and seismic hazard are different at each site, requiring site-specific analysis, design, engineering, qualification, licensing, and regulatory review, essentially making every design First-of-a-Kind (FoaK). To enable deployment at the scale needed for deep decarbonization, the cost and time impact of the seismic load case must be significantly mitigated, and plants must be standardized. Seismic base isolation has been proven to considerably reduce the earthquake response of structures and equipment but has yet to be applied to a nuclear power plant in the United States, in part because the financial impacts, positive or negative, are not known. Because there are no recent non-proprietary data to characterize the influence of the seismic load case on capital cost, it is difficult to confidently quantify the financial benefits of seismic isolation. Scheme-level designs of two fundamentally different advanced reactor buildings were developed to assemble cost data on the influence of the seismic load case. Both buildings were equipped with three bespoke pieces of safety-related equipment and analyzed for incremented levels of earthquake shaking to quantify the seismic penalty on equipment, in terms of vessel weights and horizontal accelerations. Using analysis results, a questionnaire was developed and transmitted to nuclear utilities, reactor developers, engineers, and equipment suppliers to collect cost data on engineering and fabrication costs for these unique pieces of safety-class equipment. Synthesis of the cost data showed that the seismic load case significantly affects the capital cost (sum of engineering and fabrication cost) of safety-class equipment, with engineering costs being comparable to fabrication costs. Standardization of safety-class equipment is made possible by seismic isolation, that is, equipment designed for minimal seismic robustness can resist earthquake shaking at a site of much higher seismic hazard. The average reduction in the capital cost of the safety-related equipment, enabled by seismic isolation, is a factor of two for FoaK equipment and a factor of five for standardized equipment.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Physically Realistic Estimates of Electric Water Heater Demand Response Resource

Multiple grid integration studies have examined the value of demand response, including from water heaters, to bulk power systems. Common shortcomings of these studies include physically unrealistic assumptions. For example, studies often assume that all electric water heater load can be shed for up to an hour to provide capacity or contingency service, that all electric water heater load is shiftable, and that round-trip efficiencies for shifting service are 100%. Another issue is that simply summing up estimates of individual water heaters' flexibility bounds, an attractive idea for constructing the MW-scale resources necessary for direct inclusion in grid investment and operational models, can result in significant overestimates of actual resource, because, e.g., a water heater with the ability to reduce load may not be able sustain the reduction for as long as the aggregate model suggests is possible. This presentation demonstrates the impact of estimating the flexibility of electric water heaters with physically realistic models suitable for grid service analysis. Electric water heater flexibility representations are constructed based on ResStock simulations of New England. Specifically, we examine contingency resource (related to capacity, shed, and contingency reserve services) and energy shifting resource from electric resistance and heat pump water heaters. Different estimation methods are compared; and the surrogate models' event responses are validated against events directly simulated in EnergyPlus. We also demonstrate the impact of using electric water heaters for contingency reserves and energy shifting in detailed models of possible future ISO New England power systems.

building energy modeling↗
DOE OSTI2020

Impact of Drying on Meso- and Nanoscale Structures of Citrus Fiber: A Study by SFG, ATR-IR, XRD, and DLS

Citrus fibers are a sidestream of the pectin extraction process from citrus peel. This sidestream can be converted into a functional ingredient through a shear-induced homogenization process. One technical challenge with this material is that dehydration and subsequent rehydration result in reduction of viscosity compared to the original product. In this study, various drying methods were compared with never-dried fibers to investigate the structural changes underlying the viscosity loss. Infrared and X-ray diffraction analyses confirmed no changes in chemical composition and crystalline structure of citrus fibers. Here, the dynamic light scattering and sum frequency generation analyses of citrus fiber suspension showed that the rehydration process could not fully disperse aggregated fibers, which appears to be the main cause for the viscosity loss.

09 BIOMASS FUELS↗
DOE OSTI2023

Optimal Sizing of Movable Energy Resources for Enhanced Resilience in Distribution Systems: A Techno-Economic Analysis

This article introduces a techno-economic analysis aimed at identifying the optimal total size of movable energy resources (MERs) to enhance the resilience of electric power supply. The core focus of this approach is to determine the total size of MERs required within the distribution network to expedite restoration after extreme events. Leveraging distribution line fragility curves, the proposed methodology generates numerous line outage scenarios, with scenario reduction techniques employed to minimize computational burden. For each reduced multiple line outage scenario, a systematic reconfiguration of the distribution network, represented as a graph, is executed using tie-switches within the system. To evaluate each locational combination of MERs for a specific number of these resources, the expected load curtailment (ELC) is calculated by summing the load curtailment within microgrids formed due to multiple line outages. This process is repeated for all possible locational combinations of MERs to determine minimal ELC for each MER total size. For every MER total size, the minimal ELCs are determined. Finally, a techno-economic analysis is performed using power outage cost and investment cost of MERs to pinpoint an optimal total size of MERs for the distribution system. To demonstrate the effectiveness of the proposed approach, case studies are conducted on the 33-node and the modified IEEE 123-node distribution test systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1972

Phase-lock demodulation of a PM signal contaminated with incidental AM

Signals from phase-modulated satellite transmitters usually exhibit some degree of incidental amplitude modulation. The effects of incidental AM are analyzed when this type of signal is demodulated by a phase-lock receiver which does not employ a limiter preceding the loop phase detector. The presence of incidental AM causes a reduction in the receiver output signal-to-noise ratio. The tolerable level of AM decreases in proportion to the phase modulation index Beta. For a square-wave modulating signal, a 1 db reduction results at the receiver PM channel output when Beta = 1 radian and the percentage of AM = 23, Beta = 1.2 radians and the percentage of AM = 16, or Beta = 1.5 radians and the percentage of AM = 4. Although only the PM channel of the receiver is used ordinarily, utilizing both the AM and PM channel by summing offers an improvement in S/N relative to the S/N ratio of the PM channel if the percentage of incidental AM is greater than fifteen.

Robinson, G. B.↗
NASA NTRS1982

Chromospheric, transition layer and coronal emission of metal deficient stars

It is shown that while MgII k line emission decreases for metal deficient stars, the Ly alpha emission increases. The sum of chromospheric hydrogen and metallic emission appears to be independent of metal abundances. The total chromospheric energy loss is estimated to be 0.0004 F sub bol. The chromospheric energy input does not seem to decrease for increasing age. The transition layer emission is reduced for metal deficient stars, but it is not known whether the reduction is larger than can be explained by curve of growth effects only. Coronal X-ray emission was measured for 4 metal deficient stars. Within a 12 limit it could still be consistent with the emission of solar abundance stars.

Boehm-Vitense, E.↗
NASA NTRS1994

Atmospheric pressure loading parameters from very long baseline interferometry observations

Atmospheric mass loading produces a primarily vertical displacement of the Earth's crust. This displacement is correlated with surface pressure and is large enough to be detected by very long baseline interferometry (VLBI) measurements. Using the measured surface pressure at VLBI stations, we have estimated the atmospheric loading term for each station location directly from VLBI data acquired from 1979 to 1992. Our estimates of the vertical sensitivity to change in pressure range from 0 to -0.6 mm/mbar depending on the station. These estimates agree with inverted barometer model calculations (Manabe et al., 1991; vanDam and Herring, 1994) of the vertical displacement sensitivity computed by convolving actual pressure distributions with loading Green's functions. The pressure sensitivity tends to be smaller for stations near the coast, which is consistent with the inverted barometer hypothesis. Applying this estimated pressure loading correction in standard VLBI geodetic analysis improves the repeatability of estimated lengths of 25 out of 37 baselines that were measured at least 50 times. In a root-sum-square (rss) sense, the improvement generally increases with baseline length at a rate of about 0.3 to 0.6 ppb depending on whether the baseline stations are close to the coast. For the 5998-km baseline from Westford, Massachusetts, to Wettzell, Germany, the rss improvement is about 3.6 mm out of 11.0 mm. The average rss reduction of the vertical scatter for inland stations ranges from 2.7 to 5.4 mm.

Macmillan, D. S.↗
DOE OSTI2015

Fuel assembly design for APR1400 with low CBC

APR 1400 is a PWR (Pressurized Water Reactor) with rated power of 3983 MWth and 241 assemblies. Recently, demand for extremely longer cycle up to 24 months is increasing with challenge of higher critical boron concentration (CBC). In this paper, assembly design method of selecting Gd-rods is introduced to reduce CBC. The purpose of the method is to lower the critical boron concentration of the preliminary core loading pattern (PLP), and consequently to achieve more negative or less positive moderator temperature coefficient (MTC). In this method, both the ratio of the number of low-Gd rod to the number of high-Gd rod (r) and assembly average Gd wt% (w) are the decision variables. The target function is the amount of soluble boron concentration reduction, which can be converted to Δk{sub TARGET}. A set of new designed fuel assembly satisfies an objective function, min [f=∑{sub i}(Δk{sub FA}−Δk{sub i})], and enables a final loading pattern to reach a target CBC. The constraints required to determine a set of Δk are physically realizable pair, (r,w), and the sum of Δk of new designed assemblies as close to Δk{sub TARGET} as possible. New Gd-bearing assemblies selected based on valid pairs of (r,w) are replaced with existing assemblies in a PLP. This design methodology is applied to Shin-Kori Unit 3 Cycle 1 used as a reference model. CASMO-3/MASTER code is used for depletion calculation. CASMO-3/MASTER calculations with new designed assemblies produce lower CBC than the expected CBC, proving that the proposed method works successful.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Quantify Plasma Response to Non-Axisymmetric (3D) Magnetic Fields in Tokamaks, Final Report for FES (Fusion Energy Sciences) FY2014 Joint Research Target

The goal of the 2014 Joint Research Target (JRT) has been to conduct experiments and analysis to investigate and quantify the response of tokamak plasmas to non-axisymmetric (3D) magnetic fields. Although tokamaks are conceptually axisymmetric devices, small asymmetries often result from inaccuracies in the manufacture and assembly of the magnet coils, or from nearby magnetized objects. In addition, non-axisymmetric fields may be deliberately applied for various purposes. Even at small amplitudes of order 10 -4 of the main axisymmetric field, such “3D” fields can have profound impacts on the plasma performance. The effects are often detrimental (reduction of stabilizing plasma rotation, degradation of energy confinement, localized heat flux to the divertor, or excitation of instabilities) but may in some case be beneficial (maintenance of rotation, or suppression of instabilities). In general, the magnetic response of the plasma alters the 3D field, so that the magnetic field configuration within the plasma is not simply the sum of the external 3D field and the original axisymmetric field. Typically the plasma response consists of a mixture of local screening of the external field by currents induced at resonant surfaces in the plasma, and amplification of the external field by stable kink modes. Thus, validated magnetohydrodynamic (MHD) models of the plasma response to 3D fields are crucial to the interpretation of existing experiments and the prediction of plasma performance in future devices. The non-axisymmetric coil sets available at each facility allow well-controlled studies of the response to external 3D fields. The work performed in support of the 2014 Joint Research Target has included joint modeling and analysis of existing experimental data, and collaboration on new experiments designed to address the goals of the JRT. A major focus of the work was validation of numerical models through quantitative comparison to experimental data, in order to gain confidence in their ability to predict the performance of ITER and other burning plasmas. An important element of this effort was also to verify experimentally the limits of validity of the models. The rest of this overview will summarize the key results of the work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗