Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analysis of Variance”

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 1,639 records · Page 91

Tidal and atmospheric forcing of the upper ocean in the Gulf of California. I - Sea surface temperature variability

SST variability in the northern Gulf of California is examined on the basis of findings of two years of satellite infrared imagery (1984-1986). Empirical orthogonal functions of the temporal and spatial SST variance for 20 monthly mean images show that the dominant SST patterns are generated by spatially varying tidal mixing in the presence of seasonal heating and cooling. Atmospheric forcing of the northern gulf appears to occur over large spatial scales. Area-averaged SSTs for the Guaymas Basin, island region, and northern basin exhibit significant fluctuations which are highly correlated. These fluctuations in SST correspond to similar fluctuations in the air temperature which are related to synoptic weather events over the gulf. A regression analysis of the SST relative to the fortnightly tidal range shows that tidal mixing occurs over the sills in the island region as well as on the shallow northern shelf. Mixing over the sills occurs as a result of large breaking internal waves of internal hydraulic jumps which mix over water in the upper 300-500 m.

Paden, Cynthia A.↗

Multimodel Analysis of Future Land Use and Climate Change Impacts on Ecosystem Functioning

Land use and climate changes both affect terrestrial ecosystems. Here, we used three combinations of Shared Socioeconomic Pathways and Representative Concentration Pathways (SSP1xRCP26, SSP3xRCP60, and SSP5xRCP85) as input to three dynamic global vegetation models to assess the impacts and associated uncertainty on several ecosystem functions: terrestrial carbon storage and fluxes, evapotranspiration, surface albedo, and runoff. We also performed sensitivity simulations in which we kept either land use or climate (including atmospheric CO2) constant from year 2015 on to calculate the isolated land use versus climate effects. By the 2080–2099 period, carbon storage increases by up to 87 ± 47 Gt (SSP1xRCP26) compared to present day, with large spatial variance across scenarios and models. Most of the carbon uptake is attributed to drivers beyond future land use and climate change, particularly the lagged effects of historic environmental changes. Future climate change typically increases carbon stocks in vegetation but not soils, while future land use change causes carbon losses, even for net agricultural abandonment (SSP1xRCP26). Evapotranspiration changes are highly variable across scenarios, and models do not agree on the magnitude or even sign of change of the individual effects. A calculated decrease in January and July surface albedo (up to −0.021 ± 0.007 and −0.004 ± 0.004 for SSP5xRCP85) and increase in runoff (+67 ± 6 mm/year) is largely driven by climate change. Overall, our results show that future land use and climate change will both have substantial impacts on ecosystem functioning. However, future changes can often not be fully explained by these two drivers and legacy effects have to be considered.

Krause, A.↗

Multimodel Analysis of Future Land Use and Climate Change Impacts on Ecosystem Functioning

Land use and climate changes both affect terrestrial ecosystems. Here, we used three combinations of Shared Socioeconomic Pathways and Representative Concentration Pathways (SSP1xRCP26, SSP3xRCP60, and SSP5xRCP85) as input to three dynamic global vegetation models to assess the impacts and associated uncertainty on several ecosystem functions: terrestrial carbon storage and fluxes, evapotranspiration, surface albedo, and runoff. We also performed sensitivity simulations in which we kept either land use or climate (including atmospheric CO2) constant from year 2015 on to calculate the isolated land use versus climate effects. By the 2080–2099 period, carbon storage increases by up to 87 ± 47 Gt (SSP1xRCP26) compared to present day, with large spatial variance across scenarios and models. Most of the carbon uptake is attributed to drivers beyond future land use and climate change, particularly the lagged effects of historic environmental changes. Future climate change typically increases carbon stocks in vegetation but not soils, while future land use change causes carbon losses, even for net agricultural abandonment (SSP1xRCP26). Evapotranspiration changes are highly variable across scenarios, and models do not agree on the magnitude or even sign of change of the individual effects. A calculated decrease in January and July surface albedo (up to −0.021 ± 0.007 and −0.004 ± 0.004 for SSP5xRCP85) and increase in runoff (+67 ± 6 mm/year) is largely driven by climate change. Overall, our results show that future land use and climate change will both have substantial impacts on ecosystem functioning. However, future changes can often not be fully explained by these two drivers and legacy effects have to be considered.

land use change↗

Revisiting the Global Patterns of Seasonal Cycle in Sea Surface Salinity

Argo profiling floats and L-band passive microwave remote sensing have significantly improved the global sampling of sea surface salinity (SSS) in the past 15 years, allowing the study of the range of SSS seasonal variability using concurrent satellite and in situ platforms. Here, harmonic analysis was applied to four 0.25° satellite products and two 1° in situ products between 2016 and 2018 to determine seasonal harmonic patterns. The 0.25° World Ocean Atlas (WOA) version 2018 was referenced to help assess the harmonic patterns from a long-term perspective based on the 3-year period. The results show that annual harmonic is the most characteristic signal of the seasonal cycle, and semiannual harmonic is important in regions influenced by monsoon and major rivers. The percentage of the observed variance that can be explained by harmonic modes varies with products, with values ranging between 50% and 72% for annual harmonic and between 15% and 19% for semiannual harmonic. The large spread in the explained variance by the annual harmonic reflects the large disparity in nonseasonal variance (or noise) in the different products. Satellite products are capable of capturing sharp SSS features on meso- and frontal scales and the patterns agree well with the WOA 2018. These products are, however, subject to the impacts of radiometric noises and are algorithm dependent. The coarser-resolution in situ products may underrepresent the full range of high-frequency small scale SSS variability when data record is short, which may have enlarged the explained SSS variance by the annual harmonic.

L-band↗

The Impacts of Freshwater Input and Surface Wind Velocity on the Strength and Extent of a Large High Latitude River Plume

Arctic Ocean physical and biogeochemical properties are strongly influenced by freshwater input from land and through the Bering Strait, where the mean currents transport water northward from the Bering Sea. The Yukon River is one of the largest rivers in North America and the Arctic, contributing large quantities of freshwater and terrigenous material to the coastal ocean in the northern Bering Sea. However, a detailed analysis of the coastal hydrodynamics at the outflow of the river has not been conducted in this remote but regionally important river. A three-dimensional hydrodynamic model was built to represent the lower Yukon River and coastal ocean for the ice-free months in 7 years. On average, a large anticyclonic eddy persisted at the main outflow of the Yukon that recirculates water back toward the coast where the currents converge to form a mean northward transport along the delta. Interannual spatial variance in salinity was relatively small, while there was substantial variance in u and v current velocity. u velocity spatial variance was correlated to the volume of freshwater discharge across years, while v velocity spatial variance was correlated to the N–S wind velocity. During strong wind events, plume structure was substantially altered: southerly winds deepened the plume and enhanced northward transport, while northerly winds shoaled and strengthened the pycnocline, and reversed the flow toward the south. The variability in plume dispersion on short time scales due to wind forcing has implications for where terrigenous material is processed in and settles out of the water column.

river plume↗

Material Response Analysis of a Titan Entry Heatshield

Accurate calculation of thermal protection material response is critical to the vehicle design for missions to the Saturn moon Titan. In this study, Icarus, a three-dimensional, unstructured, finite-volume material response solver under active development at NASA Ames Research Center, is used to compute the in-depth material response of the Huygens spacecraft along its November 11 entry trajectory. The heatshield analyzed in this study consists of a five-layer stack-up of Phenolic Impregnated Carbon Ablator (PICA), aluminum honeycomb, adhesive, and face sheetmaterials. During planetary entry, the PICA outer layer is expected to undergo pyrolysis. A surface energy balance boundary condition that captures both time- and spatial-variance of surface properties during entry is used in the simulation.

material response↗

MC21 analysis of the MIT PWR benchmark: Hot zero power results

MC21 Monte Carlo results have been compared with hot zero power measurements from an operating pressurized water reactor (PWR), as specified in a new full core PWR performance benchmark from the MIT Computational Reactor Physics Group. Included in the comparisons are axially integrated full core detector measurements, axial detector profiles, control rod bank worths, and temperature coefficients. Power depressions from grid spacers are seen clearly in the MC21 results. Application of Coarse Mesh Finite Difference (CMFD) acceleration within MC21 has been accomplished, resulting in a significant reduction of inactive batches necessary to converge the fission source. CMFD acceleration has also been shown to work seamlessly with the Uniform Fission Site (UFS) variance reduction method. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon

Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.

Wang, Chen↗

EFFECT OF MASKED REGIONS ON WEAK-LENSING STATISTICS

Sky masking is unavoidable in wide-field weak-lensing observations. We study how masks affect the measurement of statistics of matter distribution probed by weak gravitational lensing. We first use 1000 cosmological ray-tracing simulations to examine in detail the impact of masked regions on the weak-lensing Minkowski Functionals (MFs). We consider actual sky masks used for a Subaru Suprime-Cam imaging survey. The masks increase the variance of the convergence field and the expected values of the MFs are biased. The bias then compromises the non-Gaussian signals induced by the gravitational growth of structure. We then explore how masks affect cosmological parameter estimation. We calculate the cumulative signal-to-noise ratio (S/N) for masked maps to study the information content of lensing MFs. We show that the degradation of S/N for masked maps is mainly determined by the effective survey area. We also perform simple {chi}{sup 2} analysis to show the impact of lensing MF bias due to masked regions. Finally, we compare ray-tracing simulations with data from a Subaru 2 deg{sup 2} survey in order to address if the observed lensing MFs are consistent with those of the standard cosmology. The resulting {chi}{sup 2}/n{sub dof} = 29.6/30 for three combined MFs, obtained with the mask effects taken into account, suggests that the observational data are indeed consistent with the standard {Lambda}CDM model. We conclude that the lensing MFs are a powerful probe of cosmology only if mask effects are correctly taken into account.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of the Full Shape of the Thermal Sunyaev–Zel’dovich Power Spectrum from the South Pole Telescope and Herschel–SPIRE Observations

We present a measurement of the full shape of the power spectrum of the thermal Sunyaev–Zel’dovich (tSZ) effect down to arcminute scales using cosmic microwave background (CMB) data from the South Pole Telescope (SPT) over a roughly 100 deg 2 field. The analysis incorporates data from the 2019–2020 seasons of the SPT-3G survey in bands centered at 95, 150, and 220 GHz; from the full SPTpol dataset at 150 GHz; and from the Herschel–SPIRE survey in bands centered at 600 and 857 GHz. We combine data from all the above bands using linear combination (LC) techniques to produce a tSZ or Compton-y map. We modify the LC weights to produce multiple versions of the Compton-y map, including minimum-variance (MV) and foreground-minimized (-min) maps. We measure the auto- and cross-power spectra of a subset of these maps in the range ℓ ∈ [500, 5000]. While this power spectrum includes contributions from signals other than tSZ, we present numerous checks to show that the most challenging foreground signal, the cosmic infrared background (CIB), is much lower than the desired tSZ signal in the scales of interest in this work. The final tSZ power spectrum is measured at 9.3σ with both the MV and CIB-min maps. Our results are consistent with those reported in other CMB surveys across the literature. Using the difference in the tSZ power spectrum from the MV and CIB-min maps, we reconstruct the scale-dependent tSZ–CIB cross correlation $ρ^{\textrm{tSZ}}_{ℓ}$ x CIB, finding 3.1σ evidence for a nonzero correlation coefficient that is positive on large scales and approaches zero for ℓ > 2500. This result represents the deepest tSZ maps ever produced and provides new constraints that can help refine astrophysical feedback mechanisms and models of the intracluster medium.

Raghunathan, S. [Univ. of California, Davis, CA (U↗

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗

Techno-Economic Analysis of Recycling Strategies for Catalyst and Acid During Catalytic Graphitization

With the aim of meeting the urgent demand for active anode materials (AAM) in energy storage systems, bio-based graphite (biographite) emerges as an affordable solution to de-risk the turbulent supply chain of critical minerals. Anode grade biographite requires high crystallinity and purity, which can be achieved by catalytic graphitization with iron, followed by acid washing. Therefore, a well-conceived process integration that recycles catalyst can be the starting point to commercialization. This study evaluates closed-loop catalyst recovery, and byproducts valorization scenarios through a technoeconomic framework to help understand the scale-up potential of biographite. For the acid washing, three reactors in series meet the required biographite purity at 99.95%. Iron and acid recovery can reduce material consumption and waste generation by ~95%, albeit at the expense of ~80% increase in capital costs. Recovery scenarios present similar capital and operational expenses, yielding minimum selling prices (MSP) near $6 kg-1 of biographite. Monte Carlo methodology reveals that feedstock price accounts for ~60% of MSP variance, followed by plant capacity ~20%. The likelihood of reaching a competitive profit margin of 30% in the U.S. AAM market sits at 85% average for recovery scenarios, and 103% when iron oxide is sold as byproduct. Additionally, an IRR >= 15% can be achieved for half of Monte Carlo simulations, representing promising early-stage results. Biographite production offers a strategic pathway to stabilize the anode market beyond China by integrating established technologies for a scalable, economically viable, and sustainable process. The role of catalyst recovery and byproducts utilization is critical for advancing the biomaterials industry.

97 MATHEMATICS AND COMPUTING↗

Godiva IV Simulated Radiation Field Characterization and Variance Reduction

Godiva IV is a system comprised of highly enriched uranium alloyed with molybdenum in the form of fuel plate rings. The reactor, along with its predecessors, was designed with the unique ability to satisfy interests in the super-prompt-critical reactor operation space. Originally, the reactor was part of the Los Alamos Critical Experiments Facility (LACEF) at Technical Area-18 (TA-18). The radiation field around Godiva at this facility was well characterized and understood. As a fast neutron system, the neutron spectrum in and around Godiva was close to a Watt Fission spectrum. The Kiva where Godiva IV was located at LACEF was made of thin, sheet metal walls which did not contribute significantly to the neutron spectrum. Following the transition of LACEF to the National Critical Experiments and Research Center (NCERC) in Nevada, Godiva-IV was moved from TA-18 to the Device Assembly Facility (DAF) at the Nevada National Security Site (NNSS). Part of this move brought renewed interest in radiation field characterization. The new facility introduced significant changes to the environment surrounding Godiva, and preliminary foil irradiation results suggested that the room contribution to the neutron spectrum was significant. Unlike at TA-18, a large thermal neutron signature was added to the fast spectrum from Godiva due to significant room return. A primary goal due to the additional complexity that the room return adds to the Godiva IV radiation emission spectrum was the development of an efficient Monte Carlo N-Particle (MCNP) calculation capable of characterizing the neutron spectrum anywhere in the room around Godiva. A campaign of activation foil irradiations and analysis were completed to support the validation of the MCNP model. The modeling of these foils in MCNP can be easily done with a standard volumetric neutron flux tally. However, given the multitude of locations and reaction rates to be modeled, further steps must be taken to increase the efficiency of these calculations in MCNP. During this study, a benchmark model currently under development for Godiva IV was used. A qualitative assessment of the thermal neutron contributors was performed using spatial neutron distribution plots. Additional detail was added to the model based on the qualitative results showing the thermal spectrum’s large sensitivity to hydrogenous material. Neutron energy spectra was evaluated at discrete locations in the room around Godiva to quantify the relative contribution of various components. It was discovered that the concrete walls are the largest contributor to the thermal signature, with minor contributions from plastic components surrounding Godiva. Following these results, two different variance reduction techniques were implemented to improve the problem efficiency in these calculations. In the first approach, an F5 point detector tally was implemented in the standard Godiva IV criticality problem. The second approach involved a weight-window generator implementation with an F5 point detector tally in a fixed source problem. The weight window implementation reduced the runtime from 42739.55 minutes to 1803.34 minutes (computer time), compared to the F5 KCODE implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Low Phosphorus Availability Decreases Susceptibility of Tropical Primary Productivity to Droughts

Large uncertainties in the susceptibility of tropical forest productivity to precipitation changes hamper climate change projection. Interactions between the availabilities of water and phosphorus could theoretically either increase or decrease the susceptibility of tropical gas exchange to variation of precipitation. The inclusion of phosphorus-water interactions in a land surface model reduces the coefficient of variance, a measure of variability, of biweekly gross primary productivity by a factor of 1.5–2.3 at three tropical forest sites in Brazil, bringing it closer to estimates from eddy covariance measurements and remote sensing. Soil drought conditions are attenuated due to 8–30% lower water consumption during wet periods in presence of phosphorus limitation. When soils are dry, plant phosphorus acquisition is impaired by reduced ion mobility, despite an increase in net phosphorus mineralization. We conclude that water-phosphorus interactions cannot be omitted in analysis of the resilience of tropical ecosystems to precipitation changes.

54 ENVIRONMENTAL SCIENCES↗

Joint Planck and WMAP assessment of low CMB multipoles

The remarkable progress in cosmic microwave background (CMB) studies over past decade has led to the era of precision cosmology in striking agreement with the ΛCDM model. However, the lack of power in the CMB temperature anisotropies at large angular scales (low-ℓ), as has been confirmed by the recent Planck data also (up to 0ℓ=4), although statistically not very strong (less than 3σ), is still an open problem. One can avoid to seek an explanation for this problem by attributing the lack of power to cosmic variance or can look for explanations i.e., different inflationary potentials or initial conditions for inflation to begin with, non-trivial topology, ISW effect etc. Features in the primordial power spectrum (PPS) motivated by the early universe physics has been the most common solution to address this problem. In the present work we also follow this approach and consider a set of PPS which have features and constrain the parameters of those using WMAP 9 year and Planck data employing Markov-Chain Monte Carlo (MCMC) analysis. The prominent feature of all the models of PPS that we consider is an infra-red cut off which leads to suppression of power at large angular scales. We consider models of PPS with maximum three extra parameters and use Akaike information criterion (AIC) and Bayesian information criterion (BIC) of model selection to compare the models. For most models, we find good constraints for the cut off scale k{sub c}, however, for other parameters our constraints are not that good. We find that sharp cut off model gives best likelihood value for the WMAP 9 year data, but is as good as power law model according to AIC. For the joint WMAP 9 + Planck data set, Starobinsky model is slightly preferred by AIC which is also able to produce CMB power suppression up to 0ℓ≤3 to some extent. However, using BIC criteria, one finds model(s) with least number of parameters (power law model) are always preferred.

79 ASTRONOMY AND ASTROPHYSICS↗

Information content of higher order galaxy correlation functions

The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance–redshift relationship. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence, all the cosmological information can be extracted from their two-point correlation function. Subsequent non-linear evolution under gravity, as well as halo and then galaxy formation, generates higher order correlation functions. Since the mapping of the initial to the final density field is, on large scales, invertible, it is often claimed that the information content of the initial field’s power spectrum is equal to that of all the higher order functions of the final, non-linear field. This claim implies that reconstruction of the initial density field from the non-linear field renders analysis of higher order correlation functions of the latter superfluous. We show that this claim is false when the N-point functions are used as standard rulers. Constraints available from joint analysis of the two and three-point correlation functions can, in some cases, exceed those offered by the initial power spectrum. We provide a mathematical justification for this claim and demonstrate it using a large suite of N-body simulations. In particular, we show that for the z = 0 real-space matter field in the limit of vanishing shot-noise, taking modes up to k max = 0.2 h Mpc −1 , using the bispectrum alone offers a factor of 2 reduction in the variance on the cosmic distance scale relative to that available from the linear power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗

Modification of a successive corrections objective analysis for improved higher order calculations

The use of objectively analyzed fields of meteorological data for the initialization of numerical prediction models and for complex diagnostic studies places the requirements upon the objective method that derivatives of the gridded fields be accurate and free from interpolation error. A modification was proposed for an objective analysis developed by Barnes that provides improvements in analysis of both the field and its derivatives. Theoretical comparisons, comparisons between analyses of analytical monochromatic waves, and comparisons between analyses of actual weather data are used to show the potential of the new method. The new method restores more of the amplitudes of desired wavelengths while simultaneously filtering more of the amplitudes of undesired wavelengths. These results also hold for the first and second derivatives calculated from the gridded fields. Greatest improvements were for the Laplacian of the height field; the new method reduced the variance of undesirable very short wavelengths by 72 percent. Other improvements were found in the divergence of the gridded wind field and near the boundaries of the field of data.

Achtemeier, Gary L.↗