Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sum of squares”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Constraints on Λ CDM extensions from the SPT-3G 2018 $EE$ and $TE$ power spectra

Here, we present constraints on extensions to the Λ CDM cosmological model from measurements of the E-mode polarization autopower spectrum and the temperature-E-mode cross-power spectrum of the cosmic microwave background (CMB) made using 2018 SPT-3G data. The extensions considered vary the primordial helium abundance, the effective number of relativistic degrees of freedom, the sum of neutrino masses, the relativistic energy density and mass of a sterile neutrino, and the mean spatial curvature. We do not find clear evidence for any of these extensions, from either the SPT-3G 2018 dataset alone or in combination with baryon acoustic oscillation and Planck data. None of these model extensions significantly relax the tension between Hubble-constant, H 0 , constraints from the CMB and from distance-ladder measurements using Cepheids and supernovae. The addition of the SPT-3G 2018 data to Planck reduces the square-root of the determinants of the parameter covariance matrices by factors of 1.3–2.0 across these models, signaling a substantial reduction in the allowed parameter volume. We also explore CMB-based constraints on H 0 from combined SPT, Planck, and ACT DR4 datasets. While individual experiments see some indications of different H 0 values between the TT, TE, and EE spectra, the combined H 0 constraints are consistent between the three spectra. For the full combined datasets, we report H 0 = 67.49 ± 0.53 km s -1 Mpc -1 , which is the tightest constraint on H0 from CMB power spectra to date and in 4.1σ tension with the most precise distance-ladder-based measurement of H 0 . The SPT-3G survey is planned to continue through at least 2023, with existing maps of combined 2019 and 2020 data already having ~ 3.5 x lower noise than the maps used in this analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

Influence of Vibrotactile Feedback on Controlling Tilt Motion After Spaceflight

We hypothesize that adaptive changes in how inertial cues from the vestibular system are integrated with other sensory information leads to perceptual disturbances and impaired manual control following transitions between gravity environments. The primary goals of this ongoing post-flight investigation are to quantify decrements in manual control of tilt motion following short-duration spaceflight and to evaluate vibrotactile feedback of tilt as a sensorimotor countermeasure. METHODS. Data is currently being collected on 9 astronaut subjects during 3 preflight sessions and during the first 8 days after Shuttle landings. Variable radius centrifugation (216 deg/s, <20 cm radius) in a darkened room is utilized to elicit otolith reflexes in the lateral plane without concordant canal or visual cues. A Tilt-Translation Sled (TTS) is capable of synchronizing pitch tilt with fore-aft translation to align the resultant gravitoinertial vector with the longitudinal body axis, thereby eliciting canal reflexes without concordant otolith or visual cues. A simple 4 tactor system was implemented to provide feedback when tilt position exceeded predetermined levels in either device. Closed-loop nulling tasks are performed during random tilt steps or sum-of-sines (TTS only) with and without vibrotactile feedback of chair position. RESULTS. On landing day the manual control performance without vibrotactile feedback was reduced by >30% based on the gain or the amount of tilt disturbance successfully nulled. Manual control performance tended to return to baseline levels within 1-2 days following landing. Root-mean-square position error and tilt velocity were significantly reduced with vibrotactile feedback. CONCLUSIONS. These preliminary results are consistent with our hypothesis that adaptive changes in vestibular processing corresponds to reduced manual control performance following G-transitions. A simple vibrotactile prosthesis improves the ability to null out tilt motion within a limited range of motion disturbances.

Wood, S. J.↗
DOE OSTI2026

Impact of new physics on the JUNO-long-baseline synergy in the neutrino mass ordering determination

The determination of the neutrino mass ordering is one of the flagship goals in particle physics. A well-known and powerful synergy emerges when combining high-precision measurements of the effective atmospheric mass-squared splitting from electron antineutrino disappearance in reactor experiments with that from muon (anti)neutrino disappearance in accelerator-based long-baseline experiments. To fully exploit this synergy, percent-level precision in the atmospheric mass splitting is required—a target that JUNO is expected to achieve within a few months of data taking. This motivated the formulation of a mass ordering sum rule for neutrino disappearance channels, which shows that by combining data from T2K and NOvA with JUNO after one year of operation, the neutrino mass ordering can be determined at the 3⁢𝜎 confidence level. Since JUNO has recently started taking data, it is timely to ask whether this sum rule remains robust in the presence of new physics. We identify the necessary conditions for new physics to affect the sum rule and demonstrate that, in some cases, such effects could lead to an incorrect inference of the mass ordering. As concrete examples, we consider scalar nonstandard interactions (SNSI) and neutrinos coupled to an ultralight scalar field. We find that, for SNSI, current constraints render any modification of the sum rule negligible, whereas in the latter case, the inference of the ordering requires caution. Nevertheless, these effects can be disentangled, illustrating how the sum rule can also be used to search for new physics.

Alves, Gustavo F. S. [Fermi National Accelerator L↗
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 NTRS1999

Gaussian Decomposition of Laser Altimeter Waveforms

We develop a method to decompose a laser altimeter return waveform into its Gaussian components assuming that the position of each Gaussian within the waveform can be used to calculate the mean elevation of a specific reflecting surface within the laser footprint. We estimate the number of Gaussian components from the number of inflection points of a smoothed copy of the laser waveform, and obtain initial estimates of the Gaussian half-widths and positions from the positions of its consecutive inflection points. Initial amplitude estimates are obtained using a non-negative least-squares method. To reduce the likelihood of fitting the background noise within the waveform and to minimize the number of Gaussians needed in the approximation, we rank the "importance" of each Gaussian in the decomposition using its initial half-width and amplitude estimates. The initial parameter estimates of all Gaussians ranked "important" are optimized using the Levenburg-Marquardt method. If the sum of the Gaussians does not approximate the return waveform to a prescribed accuracy, then additional Gaussians are included in the optimization procedure. The Gaussian decomposition method is demonstrated on data collected by the airborne Laser Vegetation Imaging Sensor (LVIS) in October 1997 over the Sequoia National Forest, California.

Hofton, Michelle A.↗
DOE OSTI2026

Impact of New Physics on the JUNO-Long-Baseline Synergy in Neutrino Mass Ordering Determination

The determination of the neutrino mass ordering is one of the flagship goals in particle physics. A well-known and powerful synergy emerges when combining high-precision measurements of the effective atmospheric mass-squared splitting from electron antineutrino disappearance in reactor experiments with that from muon (anti)neutrino disappearance in accelerator-based long-baseline experiments. To fully exploit this synergy, percent-level precision in the atmospheric mass splitting is required-a target that JUNO is expected to achieve within a few months of data taking. This motivated the formulation of a mass ordering sum rule for neutrino disappearance channels, which shows that by combining data from T2K and NOvA with JUNO after one year of operation, the neutrino mass ordering can be determined at the $3σ$ confidence level. Since JUNO has recently started taking data, it is timely to ask whether this sum rule remains robust in the presence of new physics. We identify the necessary conditions for new physics to affect the sum rule and demonstrate that, in some cases, such effects could lead to an incorrect inference of the mass ordering. As concrete examples, we consider Scalar Non-Standard Interactions (SNSI) and neutrinos coupled to an ultralight scalar field. We find that, for SNSI, current constraints render any modification of the sum rule negligible, whereas in the latter case, the inference of the ordering requires caution. Nevertheless, these effects can be disentangled, illustrating how the sum rule can also be used to search for new physics.

Alves, Gustavo F.S. [Fermilab; Northwestern U.; Sa↗
DOE OSTI2025

Exploring quantum statistics for massive Dirac and Majorana neutrinos using spinor-helicity techniques

Recently, there has been interest in the applicability of quantum statistics to distinguish Dirac from Majorana neutrinos in multineutrino final states. In particular, debate has arisen over the validity of the Dirac-Majorana confusion theorem in these processes, i.e., that any distinction between the Dirac and Majorana processes goes to zero as the neutrino mass goes to zero. Here we approach this problem equipped with spinor-helicity methods generalized for massive Dirac and Majorana fermions. We explicitly calculate all helicity amplitudes, and their squares, for the decay of a light scalar particle to two neutrinos and two oppositely charged leptons. This allows us to pinpoint the crucial steps which could lead to claims of a violation of the confusion theorem. We show that, if the correct antisymmetrization of Dirac to Majorana amplitudes is used, identification of which is clear in this framework, and all relevant contributions are appropriately summed, a scalar decay into two charged leptons and two neutrinos satisfies the Dirac-Majorana confusion theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

On the mass-luminosity relation.

The results of a least-squares study of the mass-luminosity relation for eclipsing and visual binary stars consisting of main sequence components are presented. Two methods are discussed. First, the values of the coefficients A and B in the relation log M = A + BM sub BOL are determined. Then a technique which permits the determination of alpha and beta in the relation M = alpha L beta, when only the sum of the masses, and not the individual masses of each component, is known. The results and a comparison of the two methods are discussed. It is found that the following mass-luminosity relation represents the observational data satisfactorily: log M = 0.504 - 0.103 M sub BOL, -8 less than or equal to M sub BOL less than or equal to +10.5. A discussion of the data and of the possibility that separate mass-luminosity relations may exist for visual and eclipsing binaries is given. The possibility that more than one mass-luminosity relation is required in the range -8 less than or equal to M sub BOL less than or equal to +13 is also discussed.

Mccluskey, G. E., Jr.↗
DOE OSTI2023

Number of sites-based solver for determining coverages from steady-state mean-field micro-kinetic models

Kinetic models parameterized by ab-initio calculations have led to significant improvements in understanding chemical reactions in heterogeneous catalysis. These studies have been facilitated by implementations which determine steady-state coverages and rates of mean-field micro-kinetic models. As implemented in the open-source kinetic modeling program, CatMAP, the conventional solution strategy is to use a root-finding algorithm to determine the coverage of all intermediates through the steady-state expressions, constraining all coverages to be non-negative and to properly sum to unity. Though intuitive, this root-finding strategy causes issues with convergence to solution due to these imposed constraints. In this work, we avoid explicitly imposing these constraints, solving the mean-field steady-state micro-kinetic model in the space of number of sites instead of solving it in the space of coverages. Here, we transform the constrained root-finding problem to an unconstrained least-squares minimization problem, leading to significantly improved convergence in solving micro-kinetic models and thus enabling the efficient study of more complex catalytic reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Modeling radiation forces acting on TOPEX/Poseidon for precision orbit determination

Geodetic satellites such as GEOSAT, SPOT, ERS-1, and TOPEX/Poseidon require accurate orbital computations to support the scientific data they collect. Until recently, gravity field mismodeling was the major source of error in precise orbit definition. However, albedo and infrared re-radiation, and spacecraft thermal imbalances produce in combination no more than a 6-cm radial root-mean-square (RMS) error over a 10-day period. This requires the development of nonconservative force models that take the satellite's complex geometry, attitude, and surface properties into account. For TOPEX/Poseidon, a 'box-wing' satellite form was investigated that models the satellite as a combination of flat plates arranged in a box shape with a connected solar array. The nonconservative forces acting on each of the eight surfaces are computed independently, yielding vector accelerations which are summed to compute the total aggregate effect on the satellite center-of-mass. In order to test the validity of this concept, 'micro-models' based on finite element analysis of TOPEX/Poseidon were used to generate acceleration histories in a wide variety of orbit orientations. These profiles are then compared to the box-wing model. The results of these simulations and their implication on the ability to precisely model the TOPEX/Poseidon orbit are discussed.

Marshall, J. A.↗
NASA NTRS2010

Rapid Corner Detection Using FPGAs

In order to perform precision landings for space missions, a control system must be accurate to within ten meters. Feature detection applied against images taken during descent and correlated against the provided base image is computationally expensive and requires tens of seconds of processing time to do just one image while the goal is to process multiple images per second. To solve this problem, this algorithm takes that processing load from the central processing unit (CPU) and gives it to a reconfigurable field programmable gate array (FPGA), which is able to compute data in parallel at very high clock speeds. The workload of the processor then becomes simpler; to read an image from a camera, it is transferred into the FPGA, and the results are read back from the FPGA. The Harris Corner Detector uses the determinant and trace to find a corner score, with each step of the computation occurring on independent clock cycles. Essentially, the image is converted into an x and y derivative map. Once three lines of pixel information have been queued up, valid pixel derivatives are clocked into the product and averaging phase of the pipeline. Each x and y derivative is squared against itself, as well as the product of the ix and iy derivative, and each value is stored in a WxN size buffer, where W represents the size of the integration window and N is the width of the image. In this particular case, a window size of 5 was chosen, and the image is 640 480. Over a WxN size window, an equidistance Gaussian is applied (to bring out the stronger corners), and then each value in the entire window is summed and stored. The required components of the equation are in place, and it is just a matter of taking the determinant and trace. It should be noted that the trace is being weighted by a constant k, a value that is found empirically to be within 0.04 to 0.15 (and in this implementation is 0.05). The constant k determines the number of corners available to be compared against a threshold sigma to mark a valid corner. After a fixed delay from when the first pixel is clocked in (to fill the pipeline), a score is achieved after each successive clock. This score corresponds with an (x,y) location within the image. If the score is higher than the predetermined threshold sigma, then a flag is set high and the location is recorded.

Morfopoulos, Arin C.↗
NASA NTRS1986

Efficient decentralized consensus protocols

Decentralized consensus protocols are characterized by successive rounds of message interchanges. Protocols which achieve a consensus in one round of message interchange require O(N-squared) messages, where N is the number of participants. In this paper, a communication scheme, based on finite projective planes, which requires only O(N sq rt N) messages for each round is presented. Using this communication scheme, decentralized consensus protocols which achieve a consensus within two rounds of message interchange are developed. The protocols are symmetric, and the communication scheme does not impose any hierarchical structure. The scheme is illustrated using blocking and nonblocking commit protocols, decentralized extrema finding, and computation of the sum function.

Lakshman, T. V.↗
DOE OSTI2021

Fission Product Yield Data Adjustment in a Prototype Version of TSURFER

The TSURFER (Tool for Sensitivity/Uncertainty analysis of Response Functionals using Experimental Results) module of Oak Ridge National Laboratory’s (ORNL’s) SCALE code system has been updated to perform nuclear data adjustments for fixed-source irradiation/depletion problems. TSURFER uses a generalized linear least squares (GLLS) approach to consolidate a prior set of measured responses and corresponding calculated values to create the most self-consistent set of nuclear data. Traditionally, TSURFER adjustments have been performed for multigroup nuclear data such as reaction cross sections. In this work, TSURFER is expanded to perform adjustments to independent fission product yields and branching ratios that need equality constraints. To preserve equality constraints after the data adjustment procedure, an updated GLLS formulation includes a new Lagrange multiplier that forces data adjustment to sum to 0 for a given fission yield/branching ratio parent. A test problem illustrates that the newly updated TSURFER module satisfies the required constraint that adjustments for fission yield data sum to 0.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

The radiated noise from isotropic turbulence revisited

The noise radiated from isotropic turbulence at low Mach numbers and high Reynolds numbers, as derived by Proudman (1952), was the first application of Lighthill's Theory of Aerodynamic Noise to a complete flow field. The theory presented by Proudman involves the assumption of the neglect of retarded time differences and so replaces the second-order retarded-time and space covariance of Lighthill's stress tensor, Tij, and in particular its second time derivative, by the equivalent simultaneous covariance. This assumption is a valid approximation in the derivation of the second partial derivative of Tij/derivative of t exp 2 covariance at low Mach numbers, but is not justified when that covariance is reduced to the sum of products of the time derivatives of equivalent second-order velocity covariances as required when Gaussian statistics are assumed. The present paper removes these assumptions and finds that although the changes in the analysis are substantial, the change in the numerical result for the total acoustic power is small. The present paper also considers an alternative analysis which does not neglect retarded times. It makes use of the Lighthill relationship, whereby the fourth-order Tij retarded-time covariance is evaluated from the square of similar second order covariance, which is assumed known. In this derivation, no statistical assumptions are involved. This result, using distributions for the second-order space-time velocity squared covariance based on the Direct Numerical Simulation (DNS) results of both Sarkar and Hussaini(1993) and Dubois(1993), is compared with the re-evaluation of Proudman's original model. These results are then compared with the sound power derived from a phenomenological model based on simple approximations to the retarded-time/space covariance of Txx. Finally, the recent numerical solutions of Sarkar and Hussaini(1993) for the acoustic power are compared with the results obtained from the analytic solutions.

Lilley, Geoffrey M.↗
NASA NTRS2004

Fresnel Lenses for Wide-Aperture Optical Receivers

Wide-aperture receivers for freespace optical communication systems would utilize Fresnel lenses instead of conventional telescope lenses, according to a proposal. Fresnel lenses weigh and cost much less than conventional lenses having equal aperture widths. Plastic Fresnel lenses are commercially available in diameters up to 5 m large enough to satisfy requirements for aperture widths of the order of meters for collecting sufficient light in typical long-distance free-space optical communication systems. Fresnel lenses are not yet suitable for high-quality diffraction-limited imaging, especially in polychromatic light. However, optical communication systems utilize monochromatic light, and there is no requirement for high-quality imaging; instead, the basic requirement for an optical receiver is to collect the incoming monochromatic light over a wide aperture and concentrate the light onto a photodetector. Because of lens aberrations and diffraction, the light passing through any lens is focused to a blur circle rather than to a point. Calculations for some representative cases of wide-aperture non-diffraction-limited Fresnel lenses have shown that it should be possible to attain blur-circle diameters of less than 2 mm. Preferably, the blur-circle diameter should match the width of the photodetector. For most high-bandwidth communication applications, the required photodetector diameters would be about 1 mm. In a less-preferable case in which the blur circle was wider than a single photodetector, it would be possible to occupy the blur circle with an array of photodetectors. As an alternative to using a single large Fresnel lens, one could use an array of somewhat smaller lenses to synthesize the equivalent aperture area. Such a configuration might be preferable in a case in which a single Fresnel lens of the requisite large size would be impractical to manufacture, and the blur circle could not be made small enough. For example one could construct a square array of four 5-m-diameter Fresnel lenses to obtain the same light-collecting area as that of a single 10-m-diameter lens. In that case (see figure), the light collected by each Fresnel lens could be collimated, the collimated beams from the four Fresnel lenses could be reflected onto a common offaxis paraboloidal reflector, and the paraboloidal reflector would focus the four beams onto a single photodetector. Alternatively, detected signal from each detector behind each lens would be digitized before summing the signals.

Hemmati, Hamid↗
NASA NTRS1979

Free convection about a rectangular prismatic crystal growing from a solution

Sodium chlorate crystals growing from unstirred aqueous solutions were observed by Schlieren techniques. Growth rate data were best correlated by a Sherwood number with finite interface kinetics and a characteristic length equal to the crystal height plus the inverse reciprocal sum of the horizontal dimensions. The growth rate was slightly higher when the plume was laminar than when it was irregular. The horizontal growth rate was greater than the vertical growth rate, with the ratio tending to decrease with increasing plume instability. The threshold Grashof number for plume instability was greatly decreased with decreasing height for crystals with the same horizontal dimensions. The velocity of the eddies approached a constant value with increasing distance above the crystal surface, the terminal value being roughly proportional to the ratio of the square root of the Rayleigh number to the crystal characteristic length.

Chen, P.-S.↗
DOE OSTI2025

Optimized Auxiliary Functions for Robust Mitigation of Finite-Size Errors in Periodic Hybrid Density Functional Theory

When calculating properties of periodic systems at the thermodynamic limit (TDL), the dominant source of finite size error (FSE) arises from the long-range Coulomb interaction, and can manifest as a slowly converging quadrature error when approximating an integral in the reciprocal space by a finite sum. The singularity subtraction (SS) method offers a systematic approach for reducing this quadrature error and thus the FSE. Here, in this work, we first investigate the performance of the SS method in the simplest setting, aiming at reducing the FSE in exact exchange calculations by subtracting the Coulomb contribution with a single, adjustable Gaussian auxiliary function. We demonstrate that a simple fitting method can robustly estimate the optimal Gaussian width and leads to rapid convergence toward the TDL. Furthermore, we suggest new forms of the auxiliary function, whose optimal parameters could also be determined through least-squares fitting. For a range of semiconductors and insulators, the proposed auxiliary functions achieve robust, millihartree-level accuracy in hybrid density functional theory calculations, including cases with sparse k-meshes and large basis sets.

Quiton, Stephen Jon [University of California, Ber↗
DOE OSTI2022

Analysis of coherent Thomson scattering from a low temperature plasma

In this work, the spectrum of coherent Thomson scattering (CTS) induced by a periodic ponderomotive perturbation in a low-density low temperature plasma is considered. The analysis is performed for the case when the period of the resulting optical lattice is less than the Debye screening length in the plasma by solving an electron Boltzmann equation, where the total force is the sum of the periodic force due to the optical lattice and the electrostatic force due to self-consistent electric field in the plasma. An analogy between the CTS spectra calculated here and coherent Rayleigh scattering spectra in a neutral gas is established. For relatively low intensity for the optical lattice, the calculated CTS spectra are nearly Gaussian with widths slightly wider than the incoherent Thomson widths. We demonstrate that at higher intensities the line shape narrows and saturates to a width approximately half of that found at low lattice intensities. The proportionality of the spectral width to the square root of the electron temperature allows one to extract the electron temperature from the saturated spectra. Possible application of CTS for remote measuring the electron temperature in plasma is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗