Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spreading models”

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 595 records · Page 33

SYNMAG PHOTOMETRY: A FAST TOOL FOR CATALOG-LEVEL MATCHED COLORS OF EXTENDED SOURCES

Obtaining reliable, matched photometry for galaxies imaged by different observatories represents a key challenge in the era of wide-field surveys spanning more than several hundred square degrees. Methods such as flux fitting, profile fitting, and PSF homogenization followed by matched-aperture photometry are all computationally expensive. We present an alternative solution called 'synthetic aperture photometry' that exploits galaxy profile fits in one band to efficiently model the observed, point-spread-function-convolved light profile in other bands and predict the flux in arbitrarily sized apertures. Because aperture magnitudes are the most widely tabulated flux measurements in survey catalogs, producing synthetic aperture magnitudes (SYNMAGs) enables very fast matched photometry at the catalog level, without reprocessing imaging data. We make our code public and apply it to obtain matched photometry between Sloan Digital Sky Survey ugriz and UKIDSS YJHK imaging, recovering red-sequence colors and photometric redshifts with a scatter and accuracy as good as if not better than FWHM-homogenized photometry from the GAMA Survey. Finally, we list some specific measurements that upcoming surveys could make available to facilitate and ease the use of SYNMAGs.

79 ASTRONOMY AND ASTROPHYSICS↗

How robust are the size measurements of high-redshift compact galaxies?

Massive quiescent galaxies at z ≈ 2 are apparently much more compact than galaxies of comparable mass today. How robust are these size measurements? We perform comprehensive simulations to determine possible biases and uncertainties in fitting single-component light distributions to real galaxies. In particular, we examine the robustness of the measurements of the luminosity, size, and other structural parameters. We devise simulations with increasing realism to systematically disentangle effects due to the technique (specifically using GALFIT) and the intrinsic structures of the galaxies. By accurately capturing the detailed substructures of nearby elliptical galaxies and then rescaling their sizes and signal-to-noise to mimic galaxies at different redshifts, we confirm that the massive quiescent galaxies at z ≈ 2 are significantly more compact intrinsically than their local counterparts. Their observed compactness is not a result of missing faint outer light due to systematic errors in modeling. In fact, we find that fitting multi-component galaxies with a single Sérsic profile, the procedure most commonly adopted in the literature, biases the inferred sizes higher by up to 10%-20%, which accentuates the amount of size evolution required. If the sky estimation has been done robustly and the model for the point-spread function is fairly accurate, GALFIT can retrieve the properties of single-component galaxies over a wide range of signal-to-noise ratios without introducing any systematic errors.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmark calculations of pure neutron matter with realistic nucleon-nucleon interactions

Here, we report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner–Bethe–Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne v 18 interaction and two of its simplified versions, and four of the new Norfolk Δ -full chiral effective field theory potentials. Primarily because of the advancements in the auxiliary-field diffusion Monte Carlo, we observe good agreement among the three many-body techniques up to nuclear saturation density—the maximum difference in the energy per particle is within 1.5 MeV for all the potentials we consider. At higher densities, the divergences become more important, and are mainly driven by the Fermi hypernetted chain/single-operator calculations. We also study the connection between nucleon-nucleon scattering data and the energy per particle of pure neutron matter. Our results suggest that fitting to higher-energy nucleon-nucleon scattering helps reduce the spread of energies among the models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Guided Mode Evolution and Ionization Injection in Meter-Scale Multi-GeV Laser Wakefield Accelerators

Here, we show that multi-GeV laser wakefield electron accelerators in meter-scale, low density hydrodynamic plasma waveguides operate in a new nonlinear propagation regime dominated by sustained beating of lowest order modes of the ponderomotively modified channel; this occurs whether or not the injected pulse is linearly matched to the guide. For a continuously doped gas jet, this emergent mode beating effect leads to axially modulated enhancement of ionization injection and a multi-GeV energy spectrum of multiple quasimonoenergetic peaks; the same process in a locally doped jet produces single multi-GeV peaks with <10% energy spread. A three-stage model of drive laser pulse evolution and ionization injection characterizes the beating effect and explains our experimental results.

43 PARTICLE ACCELERATORS↗

A General N -Moment Normalization Method for Deriving Raindrop Size Distribution Scaling Relationships

A general drop size distribution (DSD) normalization method is formulated in terms of generalized power series relating any DSD moment to any number and combination of reference moments. This provides a consistent framework for comparing the variability of normalized DSD moments using different sets of reference moments, with no explicit assumptions about the DSD functional form (e.g., gamma). It also provides a method to derive any unknown moment plus an estimate of its uncertainty from one or more known moments, which is relevant to remote sensing retrievals and bulk microphysics schemes in weather and climate models. Herein, the approach is applied to a large dataset of disdrometer-observed and bin microphysics-modeled DSDs. As expected, the spread of normalized moments decreases as the number of reference moments is increased, quantified by the logarithmic standard deviation of the normalized moments, σ. Averaging σ for all combinations of reference moments and normalized moments of integer order 0–10, 42.9%, 81.3%, 93.7%, and 96.9% of spread are accounted for applying one-, two-, three-, and four-moment normalizations, respectively. Thus, DSDs can be well characterized overall using three reference moments, whereas adding a fourth reference moment contributes little independent information. The spread of disdrometer-observed DSD moments from uncertainty associated with drop count statistics generally lies between values of σ using two- and three-moment normalizations. However, this uncertainty has little impact on the derived DSD scaling relationships or σ when considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sandia's Research in Support of COVID-19 Pandemic Response: Computing and Information Sciences

This report summarizes the goals and findings of eight research projects conducted under the Computing and Information Sciences (CIS) Research Foundation and related to the COVID- 19 pandemic. The projects were all formulated in response to Sandia's call for proposals for rapid-response research with the potential to have a positive impact on the global health emergency. Six of the projects in the CIS portfolio focused on modeling various facets of disease spread, resource requirements, testing programs, and economic impact. The two remaining projects examined the use of web-crawlers and text analytics to allow rapid identification of articles relevant to specific technical questions, and categorization of the reliability of content. The portfolio has collectively produced methods and findings that are being applied by a range of state, regional, and national entities to support enhanced understanding and prediction of the pandemic's spread and its impacts.

59 BASIC BIOLOGICAL SCIENCES↗

Stack plume characterization and model assessment with lidar data

A mobile lidar system was used to obtain plume rise and dispersion measurements for validating plume models. Predictions of rise and spread from the two-thirds law were found to agree well with observations. These predictions were used in combination with Lamb's results for neutrally buoyant particle diffusion in the mixed layer to estimate dispersion when convective turbulence became important.

Weil, J. C.↗

Sm-Nd constraints on early lunar differentiation and the evolution of KREEP

The internal Sm-Nd systematics of lunar KREEP basalt 15386 and its implications for the evolutionary history of KREEP on the moon are reported. The age of the 15386 rocklet indicated by the measured Sm-Nd mineral isochron is consistent with the reported Rb-Sr and (Ar-39)-(Ar-40) ages of other KREEP-rich samples, which is interpreted as the time 15386 crystallized from a liquid on or near the lunar surface. This suggests that the eruption of these incompatible element-rich liquids was related to deep impact events during the final bombardment phase of the surface of the moon. It is concluded that the extremely narrow spread of Sm-Nd model ages for these samples around 4.36 AE, and the compatibility of this age with that indicated by the U-Pb and Rb-Sr systems, indicate that the source of later KREEP volcanism was produced in the closing stages of an early global scale lunar differentiation episode.

Carlson, R. W.↗

The Nazca-South American convergence rate and the recurrence of the great 1960 Chilean earthquake

The seismic slip rate along the Chile Trench estimated from the slip in the great 1960 earthquake and the recurrence history of major earthquakes has been interpreted as consistent with the subduction rate of the Nazca plate beneath South America. The convergence rate, estimated from global relative plate motion models, depends significantly on closure of the Nazca - Antarctica - South America circuit. NUVEL-1, a new plate motion model which incorporates recently determined spreading rates on the Chile Rise, shows that the average convergence rate over the last three million years is slower than previously estimated. If this time-averaged convergence rate provides an appropriate upper bound for the seismic slip rate, either the characteristic Chilean subduction earthquake is smaller than the 1960 event, the average recurrence interval is greater than observed in the last 400 years, or both. These observations bear out the nonuniformity of plate motions on various time scales, the variability in characteristic subduction zone earthquake size, and the limitations of recurrence time estimates.

Stein, S.↗

Low-frequency acoustic propagation in high winds

The effect of refraction due to wind gradients on low-frequency acoustic signals over large propagation distances is considered. The results of an investigation of low-frequency propagation in high winds are presented and compared to normal mode theory for the downwind propagation case. In the downwind direction, the low-frequency signal was channeled close to the ground resulting in a divergence of 3 dB per doubling of propagation distance. Predictions of both ray tracing and normal mode theoretical models supported the downwind cylindrical spreading. In effect, the low-frequency noise signals propagate further in the presence of wind in the downwind direction. No low-frequency acoustic shadow zone was observed in the upwind direction.

Willshire, William L., Jr.↗

Asymptotic integration algorithms for nonhomogeneous, nonlinear, first order, ordinary differential equations

New methods for integrating systems of stiff, nonlinear, first order, ordinary differential equations are developed by casting the differential equations into integral form. Nonlinear recursive relations are obtained that allow the solution to a system of equations at time t plus delta t to be obtained in terms of the solution at time t in explicit and implicit forms. Examples of accuracy obtained with the new technique are given by considering systems of nonlinear, first order equations which arise in the study of unified models of viscoplastic behaviors, the spread of the AIDS virus, and predator-prey populations. In general, the new implicit algorithm is unconditionally stable, and has a Jacobian of smaller dimension than that which is acquired by current implicit methods, such as the Euler backward difference algorithm; yet, it gives superior accuracy. The asymptotic explicit and implicit algorithms are suitable for solutions that are of the growing and decaying exponential kinds, respectively, whilst the implicit Euler-Maclaurin algorithm is superior when the solution oscillates, i.e., when there are regions in which both growing and decaying exponential solutions exist.

Walker, K. P.↗

A comparative study of several compressibility corrections to turbulence models applied to high-speed shear layers

Several recently published compressibility corrections to the standard k-epsilon turbulence model are used with the Navier-Stokes equations to compute the mixing region of a large variety of high speed flows. These corrections, specifically developed to address the weakness of higher order turbulence models to accurately predict the spread rate of compressible free shear flows, are applied to two stream flows of the same gas mixing under a large variety of free stream conditions. Results are presented for two types of flows: unconfined streams with either (1) matched total temperatures and static pressures, or (2) matched static temperatures and pressures, and a confined stream.

Viegas, John R.↗

Strontium isotopic variations of Neoproterozoic seawater - Implications for crustal evolution

High-precision Sr isotopic data were obtained on carbonate samples from the Neoproterozoic Shaler Group, Victoria Island (Canada). Results indicate that, between ca. 790 and 850 Ma, the Sr-87/Sr-86 ratio of seawater varied betweeen 0.70676 and 0.70561, with the minimum value at about 830 Ma. A curve of the Sr-87/Sr-86 seawater ratio vs. age showed that the new data substantially improve the existing isotopic record of Sr in seawater for the period 790-850 Ma. The Sr isotopic system data were coupled with data for the Nd isotopic system to model changes in the seafloor spreading rates (hydrothermal flux) and the continental erosion for the period 500-900 Ma. Results indicate that hydrothermal flux reached a maximum value at ca. 830 Ma, while a maximum in erosion rate occurred at ca. 570 Ma. These peaks are considered to be related to the developments in the Pan-African and related orogenic events.

Asmerom, Yemane↗

Automatic array alignment in data-parallel programs

FORTRAN 90 and other data-parallel languages express parallelism in the form of operations on data aggregates such as arrays. Misalignment of the operands of an array operation can reduce program performance on a distributed-memory parallel machine by requiring nonlocal data accesses. Determining array alignments that reduce communication is therefore a key issue in compiling such languages. We present a framework for the automatic determination of array alignments in array-based, data-parallel languages. Our language model handles array sectioning, reductions, spreads, transpositions, and masked operations. We decompose alignment functions into three constituents: axis, stride, and offset. For each of these subproblems, we show how to solve the alignment problem for a basic block of code, possibly containing common subexpressions. Alignments are generated for all array objects in the code, both named program variables and intermediate results. We assign computation to processors by virtue of explicit alignment of all temporaries; the resulting work assignment is in general better than that provided by the 'owner-computes' rule. Finally, we present some ideas for dealing with control flow, replication, and dynamic alignments that depend on loop induction variables.

Chatterjee, Siddhartha↗

Downward Link of Solar Activity Variations Through Wave Driven Equatorial Oscillations (QBO and SAO)

Signatures of the 11-year solar activity/irradiance cycle are observed in the Quasi Biennial Oscillation (QBO) of the lower stratosphere. At these altitudes, the QBO is understood to be the result of "downward control" exerted by the wave mean flow interactions that drive the phenomenon. It is reasonable then to speculate that the QBO is a natural conduit to lower altitudes of solar activity variations in radiance (SAV). To test this hypothesis, we conducted experiments with a 2D version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). To account for the SAV, we change the solar heating rate on a logarithmic scale from 0.1% at the surface to 1% at 50 kin to 10% at 100 km. With the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the solar activity irradiance effect (SAE) on the zonal circulation at low latitudes. The numerical results obtained show that, under certain conditions, the SAE is significant in the zonal circulation and does extend to lower altitudes where the SAV is small. The differences in the wind velocities can be as large as 5 m/s at 20 kin. We carried out two numerical experiments with integrations over more than 20 years: 1) With the QBO period "tuned" to be 30 months, of academic interest but instructive, the seasonal cycle in the solar forcing [through the Semi-annual Oscillation (SAO)] acts as a strong pacemaker to produce a firm lock on the period and phase of the QBO. The SAE then shows up primarily as a distinct but relatively weak amplitude modulation. 2) With the QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak compared with (1). The SAV in the seasonal cycle then causes variations in the QBO period and phase, and this amplifies the SAE to produce relatively large variations in the wind field. We conclude that, under realistic conditions as in (2), the solar seasonal forcing, with activity variations in radiance, causes the phase and period of the QBO to change to produce a relatively large solar activity effect in the zonal circulation at low latitudes and low altitudes.

Mengel, J. G.↗

Thermal Non-equilibrium Consistent with Widespread Cooling

Time correlation analysis has been used to show widespread cooling in the solar corona; this cooling has been interpreted as a result of impulsive (nanoflare) heating. In this work, we investigate wide-spread cooling using a 3D model for a solar active region which has been heated with highly stratified heating. This type of heating drives thermal non-equilibrium solutions, meaning that though the heating is effectively steady, the density and temperature in the solution are not. We simulate the expected observations in narrowband EUV images and apply the time correlation analysis. We find that the results of this analysis are qualitatively similar to the observed data. We discuss additional diagnostics that may be applied to differentiate between these two heating scenarios.

Winebarger, A.↗

Thermal Non-equilibrium Consistent with Widespread Cooling

Time correlation analysis has been used to show widespread cooling in the solar corona; this cooling has been interpreted as a result of impulsive (nanoflare) heating. In this work, we investigate wide-spread cooling using a 3D model for a solar active region which has been heated with highly stratified heating. This type of heating drives thermal non-equilibrium solutions, meaning that though the heating is effectively steady, the density and temperature in the solution are not. We simulate the expected observations in narrowband EUV images and apply the time correlation analysis. We find that the results of this analysis are qualitatively similar to the observed data. We discuss additional diagnostics that may be applied to differentiate between these two heating scenarios.

Winebarger, A.↗