Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SCALE”

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 919 records · Page 51

Competition between roughness and strength for scale-dependent surfaces

Rocks famously have scale-dependent strength, yet the actual dependence is notoriously hard to measure or incorporate into any theoretical framework. Natural rough surfaces present an opportunity to solve the problem. Surfaces sliding in shear evolve as protrusions collide. These asperities can deform or break, thus creating a new surface shape. In particular, natural surfaces have roughness at all scales as well as scale-dependent strength. Based on a scaling analysis, we have previously suggested that the scale-dependent aspect ratio of steady-state surfaces should be proportional to the scale-dependent shear strain at yield. If true, scale-dependent strength could easily be inferred from natural surfaces. Thus, moving beyond the scaling argument to a rigorous treatment of scale-dependent strength for multiscale rough surfaces in shear is important. However, analytic frameworks for analyzing multiscale problems are challenging, as conventional continuum mechanics typically involves a single value for a material property across scales. Here, in this work, we build on the formalism of Persson (2001) that presents a method to compute contact area for rough surfaces with a prescribed topographic spectrum using a stochastic differential equation. The Persson formalism allows for plastic yield under normal loading of otherwise elastic materials and leaves open the possibility of scale-dependent yield stress. In this study, we pursue this route to develop a theory and numerical results for the yielding of a rough, elastoplastic surface with scale-dependent yield stress. Here, we examine surfaces for which the power spectrum of the topography 𝐶 and yield stress 𝑌 follow power laws as a function of scale 𝜆, such that 𝐶∼𝜆 −𝑚 and 𝑌∼𝜆 −𝑛 , respectively. In this formal treatment of the problem, we focus on surfaces in contact and the resulting yield and do not impose shear. Numerical solutions show that the deviation from the elastic scaling solution is bounded as expected by the prior 1D heuristic scaling argument that anticipates the Hurst exponent as 1−𝑛. We also show that the plasticity is expected to erode the contacts if 𝑚 is lower than 𝑛−3, which corresponds to a Hurst exponent lower than 1−𝑛/2. This result is rigorously sound for 2D, i.e., realistic surfaces, and quantitatively different than the prior scaling argument. The theory now permits a correspondingly quantitative approach to interpreting natural surfaces.

elasticity↗

Diffusional growth of cloud droplets in homogeneous isotropic turbulence: DNS, scaled-up DNS, and stochastic model

This paper presents a novel methodology to use direct numerical simulation (DNS) to study the impact of isotropic homogeneous turbulence on the condensational growth of cloud droplets. As shown by previous DNS studies, the impact of turbulence increases with the computational domain size, that is, with the Reynolds number, because larger eddies generate higher and longer-lasting supersaturation fluctuations that affect growth of individual cloud droplets. The traditional DNS can only simulate a limited range of scales because of the excessive computational cost that comes from resolving all scales involved, that is, from large scales at which the turbulent kinetic energy (TKE) is introduced down to the Kolmogorov microscale, and from following every single droplet. The novel approach is referred to as the “scaled-up DNS”. The scaling up is done in two parts, first by increasing both the computational domain and the Kolmogorov microscale and second by using super-droplets instead of real droplets. To ensure proper dissipation of TKE and scalar variance at small scales, molecular transport coefficients are appropriately scaled up with the grid length. For the scaled-up domains, say, meters and tens of meters, one needs to follow billions of real droplets. This is not computationally feasible, and so-called super-droplets are applied in scaled-up DNS simulations. Each super-droplet represents an ensemble of identical real droplets, and the number of real droplets represented by a super-droplet is referred to as the multiplicity attribute. After simple tests showing the validity of the methodology, scaled-up DNS simulations are conducted for five domains, the largest of 64 3 m 3 volume using a DNS of 256 3 grid points and various multiplicities. All simulations are carried out with vanishing mean vertical velocity and with no mean supersaturation, similarly to past DNS studies. As expected, the supersaturation fluctuations as well as the spread in droplet size distribution increase with the domain size, with the droplet radius variance increasing in time t as t 1∕2 as identified in previous DNS studies. Scaled-up simulations with different multiplicities document numerical convergence of the scaled-up solutions. Finally, we compare the scaled-up DNS results with a simple stochastic model that calculates supersaturation fluctuations based on the vertical velocity fluctuations updated using the Langevin equation. Overall, the results document similar scaling to previous small-domain DNS simulations and support the notion that the stochastic subgrid-scale model is a valuable tool for the multi-scale simulation of droplet spectral evolution applying a large-eddy simulation model.

54 ENVIRONMENTAL SCIENCES↗

Structural similitude and design of scaled down laminated models

The excellent mechanical properties of laminated composite structures make them prime candidates for wide variety of applications in aerospace, mechanical and other branches of engineering. The enormous design flexibility of advanced composites is obtained at the cost of large number of design parameters. Due to complexity of the systems and lack of complete design based informations, designers tend to be conservative in their design. Furthermore, any new design is extensively evaluated experimentally until it achieves the necessary reliability, performance and safety. However, the experimental evaluation of composite structures are costly and time consuming. Consequently, it is extremely useful if a full-scale structure can be replaced by a similar scaled-down model which is much easier to work with. Furthermore, a dramatic reduction in cost and time can be achieved, if available experimental data of a specific structure can be used to predict the behavior of a group of similar systems. This study investigates problems associated with the design of scaled models. Such study is important since it provides the necessary scaling laws, and the factors which affect the accuracy of the scale models. Similitude theory is employed to develop the necessary similarity conditions (scaling laws). Scaling laws provide relationship between a full-scale structure and its scale model, and can be used to extrapolate the experimental data of a small, inexpensive, and testable model into design information for a large prototype. Due to large number of design parameters, the identification of the principal scaling laws by conventional method (dimensional analysis) is tedious. Similitude theory based on governing equations of the structural system is more direct and simpler in execution. The difficulty of making completely similar scale models often leads to accept certain type of distortion from exact duplication of the prototype (partial similarity). Both complete and partial similarity are discussed. The procedure consists of systematically observing the effect of each parameter and corresponding scaling laws. Then acceptable intervals and limitations for these parameters and scaling laws are discussed. In each case, a set of valid scaling factors and corresponding response scaling laws that accurately predict the response of prototypes from experimental models is introduced. The examples used include rectangular laminated plates under destabilizing loads, applied individually, vibrational characteristics of same plates, as well as cylindrical bending of beam-plates.

Simitses, G. J.↗

Scaling Input Distributions for Probabilistic Models - 19472

It is important to address scaling issues when building models to support risk-based decisions for radioactive waste disposal sites. Currently, probabilistic risk assessment and probabilistic performance assessment models are run for multiple realizations where each realization uses one value from each input distribution over the entire modeled time. The input distributions represent the current state of knowledge in the mean of the input variable over the spatial region of interest and over the time period of the model run. This paper focuses on the time factor in input distribution scaling because it presents more serious conceptual and distributional challenges than the spatial factor. One of the challenges in temporal scaling is the effect of non-linearity on propagation of distribution uncertainty through to uncertainty in results. Radioactive waste disposal models tend to have both linear and non-linear components in their models. Dose assessment calculations are typically represented by linear models, while contaminant transport models have linear and non-linear components. In principle, temporal scaling of input parameter distributions for a linear model is straightforward, but scaling of input distributions for non-linear models is not. This paper uses examples to investigate the effects of input distribution temporal scaling on simple linear and non-linear models. Different scaling scenarios were investigated by comparing three analysis tools: analytical mathematical representations, statistical software programming, and GoldSim Monte Carlo simulation software. R statistical software is an open source statistical programming language often used for statistical analyses and graphical representation of data. GoldSim is a robust probabilistic dynamic modeling platform used in performance assessment modeling of radioactive waste sites. A simple linear model, a quadratic model and a multiplicative model were the cases explored for each of these tools. For the three analysis tools considered and with each of the scaling cases, 1,000 realizations were implemented from distributions of the random variables of interest. This is a sufficient number of realizations for simple models to capture the model space. An annual temporal resolution of a 100-year simulation period was chosen for all models and analysis tools. The effects of temporal scaling were demonstrated by investigating two different options for time scaling. Option 1 was set up so that the input distributions represented the expected value on an annual time scale, and a new value from that distribution was chosen every year (at every time step which equaled a year) for 100 years. In Option 2, the input distributions were scaled to represent the expected mean value of the input variable over the 100 years of simulation, and one value from that distribution was chosen for the entire elapsed time. Results show agreement across the three analysis tools. As expected, the linear models have identical results for the two different scaling Options, and the non-linear models diverge in the results for different time scale implementation. The implications of this work point to the need for careful consideration of input distribution temporal scaling. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Scaling impacts on environmental controls and spatial heterogeneity of soil organic carbon stocks

The spatial heterogeneity of land surfaces affects energy, moisture, and greenhouse gas exchanges with the atmosphere. However, representing the heterogeneity of terrestrial hydrological and biogeochemical processes in Earth system models (ESMs) remains a critical scientific challenge. We report the impact of spatial scaling on environmental controls, spatial structure, and statistical properties of soil organic carbon (SOC) stocks across the US state of Alaska. We used soil profile observations and environmental factors such as topography, climate, land cover types, and surficial geology to predict the SOC stocks at a 50 m spatial scale. These spatially heterogeneous estimates provide a data set with reasonable fidelity to the observations at a sufficiently high resolution to examine the environmental controls on the spatial structure of SOC stocks. We upscaled both the predicted SOC stocks and environmental variables from finer to coarser spatial scales ( s = 100, 200, and 500 m and 1, 2, 5, and 10 km) and generated various statistical properties of SOC stock estimates. We found different environmental factors to be statistically significant predictors at different spatial scales. Only elevation, temperature, potential evapotranspiration, and scrub land cover types were significant predictors at all scales. The strengths of control (the median value of geographically weighted regression coefficients) of these four environmental variables on SOC stocks decreased with increasing scale and were accurately represented using mathematical functions ( R 2 = 0.83–0.97). The spatial structure of SOC stocks across Alaska changed with spatial scale. Although the variance (sill) and unstructured variability (nugget) of the calculated variograms of SOC stocks decreased exponentially with scale, the correlation length (range) remained relatively constant across scale. The variance of predicted SOC stocks decreased with spatial scale over the range of 50 m to ~ 500 m, and remained constant beyond this scale. The fitted exponential function accounted for 98 % of variability in the variance of SOC stocks. We found moderately accurate linear relationships between mean and higher-order moments of predicted SOC stocks ( R 2 ∼ 0.55–0.63). Current ESMs operate at coarse spatial scales (50–100 km), and are therefore unable to represent environmental controllers and spatial heterogeneity of high-latitude SOC stocks consistent with observations. We conclude that improved understanding of the scaling behavior of environmental controls and statistical properties of SOC stocks could improve ESM land model benchmarking and perhaps allow representation of spatial heterogeneity of biogeochemistry at scales finer than those currently resolved by ESMs.

54 ENVIRONMENTAL SCIENCES↗

Fragmentation Prediction for a Brittle Ceramic: An Integrated Two-Scale Model Approach

A new implementation of a multi-scale fragmentation model is introduced that combines existing independent macro and micro-scale models into a single material model within the Library of Advanced Materials for Engineering (LAME) that can be used directly by the finite element structural mechanics code SIERRA/SM, developed at Sandia National Laboratories. The implementation not only joins existing independent fragmentation models into a single framework, but also enables further development by virtue of being integrated into SIERRA/SM. Integration with a finite element structural mechanics code for a true two-way coupling of the fragmentation models that were previously independent. Hereafter, the prior implementation will be referred to as the sequentially applied fragmentation model. In contrast, the newly developed model will be called the integrated two-scale model. The sequentially applied fragmentation model consists of two independent models applied sequentially to account for the macro-scale and microscale fragmentation. The macro-scale fragmentation behavior is calculated using an elastic gradient damage explicit material model developed to model fracture in brittle materials. The micro-scale fragmentation model is an implementation of the 1D cohesive-zone model developed by Zhou, et al. The micro-scale model is implemented by post-processing the macro-scale data that provides initial conditions for the 1D analysis. The micro-scale model is necessary to predict the resulting fragment size distribution at length scales inaccessible to the finite element based macro-scale model. The full range of fragmentation behavior is then predicted by combining the results from the two independent models. The sequentially applied fragmentation model was used to predict the fragments that result from impact loading of a brittle material, which in that case was a ceramic cylinder of commercial UO{sub 2} nuclear fuel. Understanding the fragmentation behavior of such materials is important for evaluating the risk of airborne particulate release under impact loading scenarios. This work will verify the implementation of the integrated two-scale model using the same experiment modeled by the existing sequentially applied model. The purpose of this manuscript is to introduce the new integrated two-scale fragmentation model and validate the implementation against the existing sequentially applied model for a specific case of impact loading of a brittle material. This validation case will not fully utilize the benefits of having the model integrated into the finite element code as a material model. Rather, this study will demonstrate consistency between the two fragmentation models to ensure that the success of the sequentially applied model will be duplicated by the integrated two-scale approach. Once validated, further development of the model can be made by utilizing the integration into the finite element code to create a fully two-way coupled fragmentation model that will be extended to include higher fidelity models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The renormalization scale-setting problem in QCD

A key problem in making precise perturbative QCD predictions is to set the proper renormalization scale of the running coupling. The conventional scale-setting procedure assigns an arbitrary range and an arbitrary systematic error to fixed-order pQCD predictions. In fact, this ad hoc procedure gives results which depend on the choice of the renormalization scheme, and it is in conflict with the standard scale-setting procedure used in QED. Predictions for physical results should be independent of the choice of the scheme or other theoretical conventions. We review current ideas and points of view on how to deal with the renormalization scale ambiguity and show how to obtain renormalization scheme- and scale-independent estimates. We begin by introducing the renormalization group (RG) equation and an extended version, which expresses the invariance of physical observables under both the renormalization scheme and scale-parameter transformations. The RG equation provides a convenient way for estimating the scheme- and scale-dependence of a physical process. We then discuss self-consistency requirements of the RG equations, such as reflexivity, symmetry, and transitivity, which must be satisfied by a scale-setting method. Four typical scale setting methods suggested in the literature, i.e., the Fastest Apparent Convergence (FAC) criterion, the Principle of Minimum Sensitivity (PMS), the Brodsky–Lepage–Mackenzie method (BLM), and the Principle of Maximum Conformality (PMC), are introduced. Basic properties and their applications are discussed. We pay particular attention to the PMC, which satisfies all of the requirements of RG invariance. Using the PMC, all non-conformal terms associated with the β-function in the perturbative series are summed into the running coupling, and one obtains a unique, scale-fixed, scheme-independent prediction at any finite order. The PMC provides the principle underlying the BLM method, since it gives the general rule for extending BLM up to any perturbative order; in fact, they are equivalent to each other through the PMC–BLM correspondence principle. Thus, all the features previously observed in the BLM literature are also adaptable to the PMC. The PMC scales and the resulting finite-order PMC predictions are to high accuracy independent of the choice of the initial renormalization scale, and thus consistent with RG invariance. The PMC is also consistent with the renormalization scale-setting procedure for QED in the zero-color limit. The use of the PMC thus eliminates a serious systematic scale error in perturbative QCD predictions, greatly improving the precision of empirical tests of the Standard Model and their sensitivity to new physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Airborne Measurements of Scale‐Dependent Latent Heat Flux Impacted by Water Vapor and Vertical Velocity Over Heterogeneous Land Surfaces During the CHEESEHEAD19 Campaign

Abstract The water vapor transport associated with latent heat flux (LE) in the planetary boundary layer (PBL) is critical for the atmospheric hydrological cycle, radiation balance, and cloud formation. The spatiotemporal variability of LE and water vapor mixing ratio ( r v ) are poorly understood due to the scale‐dependent and nonlinear atmospheric transport responses to land surface heterogeneity. Here, airborne in situ measurements with the wavelet technique are utilized to investigate scale‐dependent relationships among LE, vertical velocity ( w ) variance (), and r v variance () over a heterogeneous surface during the Chequamegon Heterogeneous Ecosystem Energy‐balance Study Enabled by a High‐density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign. Our findings reveal distinct scale distributions of LE, , and at 100 m height, with a majority scale range of 120 m–4 km in LE, 32 m–2 km in , and 200 m–8 km in . The scales are classified into three scale ranges, the turbulent scale (8–200 m), large‐eddy scale (200 m–2 km), and mesoscale (2–8 km) to evaluate scale‐resolved LE contributed by and . The large‐eddy scale in PBL contributes over 70% of the monthly mean total LE with equal parts (50%) of contributions from and . The monthly temporal variations mainly come from the first two major contributing classified scales in LE, , and . These results confirm the dominant role of the large‐eddy scale in the PBL in the vertical moisture transport from the surface to the PBL, while the mesoscale is shown to contribute an additional ∼20%. This analysis complements published scale‐dependent LE variations, which lack detailed scale‐dependent vertical velocity and moisture information.

54 ENVIRONMENTAL SCIENCES↗

On the linkage between the k –5/3 spectral and k –7/3 cospectral scaling in high-Reynolds number turbulent boundary layers

Connections between the “–5/3” spectral and “–7/3” cospectral scaling exponents characterizing the inertial subranges of the wall-normal energy spectrum and the turbulent momentum flux cospectrum are explored in the equilibrium layer of high-Reynolds number turbulent boundary layers. Previous laboratory experiments and field measurements featured here in the atmospheric boundary layer show that the “–7/3” scaling in the momentum flux cospectrum F uw (k) commences at lower wavenumbers (around kz = 3) than the “–5/3” scaling in the wall-normal energy spectrum E ww (k) (around kz = 6), where k is the streamwise wavenumber and z is the distance from the surface. A satisfactory explanation as to why F uw (k) attains its “–7/3” inertial subrange scaling earlier than E ww (k) in wavenumber space remains elusive. A cospectral budget (CSB) model subject to several simplifications and closure schemes offers one viewpoint. In its simplest form, the CSB model assumes a balance at all k between the production term and a Rotta-like pressure decorrelation term with a prescribed wavenumber-dependent relaxation time scale. It predicts the “–7/3” scaling for F uw (k) from the “–5/3” scaling in E ww (k), thereby recovering earlier results derived from dimensional considerations. A finite flux transfer term was previously proposed to explain anomalous deviations from the “–7/3” cospectral scaling in the inertial subrange using a simplified spectral diffusion closure. However, this explanation is not compatible with an earlier commencement of the “–7/3” scaling in F uw (k). An alternative explanation that does not require a finite flux transfer is explored here. By linking the relaxation time scale in the slow-component of the Rotta model to the turbulent kinetic energy (TKE) spectrum, the earlier onset of the “–7/3” scaling in Fuw(k) is recovered without attainment of a “–5/3” scaling in E ww (k). In conclusion, the early onset of the “–7/3” scaling at smaller k is related to a slower than k–5/3 decay in the TKE spectrum at the crossover from production to inertial scales.

58 GEOSCIENCES↗

Servo-aero-gravo-elastic (SAGE) scaling and its application to a 13-MW downwind turbine

Reduced scale wind turbines can be extremely cost-effective to test new rotor concepts since prototype costs are heavily dependent on the rotor diameter. Ideally, the scaled model would have the same non-dimensional deflections, dynamics, and control behavior as the full-scale model. This would provide a high-fidelity demonstration of the full-scale performance, which is ideal if the full-scale turbine has significant aeroelastic interactions. To this end, servo-aero-gravo-elastic (SAGE) scaling is developed and applied to a 13-MW turbine that is scaled to a 20% scale model. The scaling preserves the tip-speed ratio, the rotor speed normalized by the flapping frequency, and the tip deflections normalized by the blade length. In addition, the controller employs the same control structure (gain-scheduled pitch control and variable speed torque control) and is scaled dynamically (e.g., matching non-dimensional time constant of the pitch angle, etc.). Furthermore, the thrust, gravity, and centrifugal moments are scaled such that the load angles are preserved as a function of a non-dimensional wind speed. However, the environmental scaling must consider differences in Reynolds number (since this parameter cannot be held constant) and subsequent changes in the axial induction factor. While the presented results showcase these differences during operational conditions, the non-dimensional tip deflections remain comparable through all wind speed ranges, indicating the viability of the SAGE scaling method in matching full-scale aeroelastic responses.

Kaminski, Meghan (ORCID:0000000322334784)↗

A NEW METHOD TO DIRECTLY MEASURE THE JEANS SCALE OF THE INTERGALACTIC MEDIUM USING CLOSE QUASAR PAIRS

Although the baryons in the intergalactic medium (IGM) trace dark matter fluctuations on megaparsec scales, on smaller scales ∼100 kpc, fluctuations are suppressed because the finite temperature gas is pressure supported against gravity, analogous to the classical Jeans argument. This Jeans filtering scale, which quantifies the small-scale structure of the IGM, has fundamental cosmological implications. First, it provides a thermal record of heat injected by ultraviolet photons during cosmic reionization events, and thus constrains the thermal and reionization history of the universe. Second, the Jeans scale determines the clumpiness of the IGM, a critical ingredient in models of cosmic reionization. Third, it sets the minimum mass scale for gravitational collapse from the IGM, and hence plays a pivotal role in galaxy formation. Unfortunately, it is extremely challenging to measure the Jeans scale via the standard technique of analyzing purely longitudinal Lyα forest spectra, because the thermal Doppler broadening of absorption lines along the line-of-sight, is highly degenerate with Jeans smoothing. In this work, we show that the Jeans filtering scale can be directly measured by characterizing the coherence of correlated Lyα forest absorption in close quasar pairs, with separations small enough ∼100 kpc to resolve it. We present a novel technique for this purpose, based on the probability density function (PDF) of phase angle differences of homologous longitudinal Fourier modes in close quasar pair spectra. A Bayesian formalism is introduced based on the phase angle PDF, and Markov Chain Monte Carlo techniques are used to characterize the precision of a hypothetical Jeans scale measurement, and explore degeneracies with other thermal parameters governing the IGM. A semi-analytical model of the Lyα forest is used to generate a large grid (500) of thermal models from a dark matter only simulation. Our full parameter study indicates that a realistic sample of only 20 close quasar pair spectra can pinpoint the Jeans scale to ≅ 5% precision, independent of the amplitude T{sub 0} and slope γ of the temperature-density relation of the IGM T=T{sub 0}(ρ/ ρ-bar ){sup γ-1}. This exquisite sensitivity arises because even long-wavelength one-dimensional Fourier modes ∼10 Mpc, i.e., two orders of magnitude larger than the Jeans scale, are nevertheless dominated by projected small-scale three-dimensional (3D) power. Hence phase angle differences between all modes of quasar pair spectra actually probe the shape of the 3D power spectrum on scales comparable to the pair separation. We show that this new method for measuring the Jeans scale is unbiased and is insensitive to a battery of systematics that typically plague Lyα forest measurements, such as continuum fitting errors, imprecise knowledge of the noise level and/or spectral resolution, and metal-line absorption.

79 ASTRONOMY AND ASTROPHYSICS↗

Bridging the length scales on mechanical property evaluation (Final Report)

The development of small-scale mechanical testing in combination with microstructural investigation is of great interest to the nuclear materials community for both materials development and monitoring applications. Dramatically reducing the sample sizes to reduce radioactivity and obtaining mechanical properties of irradiated samples is truly intriguing. Moreover, such studies promise a range of benefits including cost reduction, fundamental insight in structure-property relationships, increased statistics on less sample material, and reinvestigation of prior irradiated and tested reactor samples while simultaneously enabling the generation of mechanical test data on ion beam irradiated materials with limited penetration depths. Small scale materials testing on sub-sized samples has been studied for several decades, though it has only been after the development of micro-testing based on Focused Ion Beam (FIB) sample manufacturing in that orders of magnitude smaller samples could really be investigated in a quantitative manner. In recent years, small scale mechanical testing techniques at a number of length scales has been developed for both unirradiated and irradiated (ion and neutron) materials. Technological advances made in this field have enabled ex-situ and in-situ transmission electron microscopy (TEM) and scanning electron microscopy (SEM) examination, thus leading to more accurate measurements as well as additional mechanistic information. A recent review of the benefits of these techniques show that these techniques are at a stage to tackle multi-scale ranges of materials investigations and can be utilized to obtain fundamental science-based understanding of nuclear materials. Considering the tremendous advances made, one can see how small-scale mechanical testing techniques combined with modeling can enable true small scale to bulk scale mechanical property correlations. However, for the engineering community to adapt this approach fully, one needs to demonstrate that a) that these techniques can produce results with high fidelity and reproducibility, b) generate engineering stress-strain data that one can utilize to understand bulk behavior, and c) generate new insight into relevant phenomena fostering the true understanding of radiation damage and microstructure in materials for nuclear applications. It is the objective of this proposal to bridge the length scale between macro- and micro- scale mechanical testing of unirradiated and irradiated materials. This involves the development and demonstration of procedures for multi-scale mechanical testing that enable high fidelity reproducibility of data and the generation bulk property data from small-scale mechanical tests. Through this, the proposal aims to enhance the confidence in the obtained data at the smaller length scales and enhance the insight provided from these techniques for bulk scale applications on both unirradiated and irradiated materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scaling impacts on environmental controls and spatial heterogeneity of soil organic carbon stocks

The spatial heterogeneity of land surfaces affects energy, moisture, and greenhouse gas exchanges with the atmosphere. However, representing heterogeneity of terrestrial hydrological and biogeochemical processes in earth system models (ESMs) remains a critical scientific challenge. We report the impact of spatial scaling on environmental controls, spatial structure, and statistical properties of soil organic carbon (SOC) stocks across the US state of Alaska. We used soil profile observations and environmental factors such as topography, climate, land cover types, and surficial geology to predict the SOC stocks at a 50 m spatial scale. These spatially heterogeneous estimates provide a dataset with reasonable fidelity to the observations at a sufficiently high resolution to examine the environmental controls on the spatial structure of SOC stocks. We upscaled both the predicted SOC stocks and environmental variables from finer to coarser spatial scales ( s = 100, 200, 500 m, 1, 2, 5, 10 km) and generated various statistical properties of SOC stock estimates. We found different environmental factors to be statistically significant predictors at different spatial scales. Only elevation, temperature, potential evapotranspiration, and scrub land cover types were significant predictors at all scales. The strengths of control (the median value of geographically weighted regression coefficients) of these four environmental variables on SOC stocks decreased with increasing scale and were accurately represented using mathematical functions ( R 2 = 0.83–0.97). The spatial structure of SOC stocks across Alaska changed with spatial scale. Although the variance (sill) and unstructured variability (nugget) of the calculated variograms of SOC stocks decreased exponentially with scale, the correlation length (range) remained relatively constant across scale. The variance of predicted SOC stocks decreased with spatial scale over the range of 50 to ~ 500 m, and remained constant beyond this scale. The fitted exponential function accounted for 98% of variability in the variance of SOC stocks. We found moderately-accurate linear relationships between mean and higher-order moments of predicted SOC stocks ( R 2 ~ 0.55–0.63). Current ESMs operate at coarse spatial scales (50–100 km), and are therefore unable to represent environmental controllers and spatial heterogeneity of high-latitude SOC stocks consistent with observations. We conclude that improved understanding of the scaling behavior of environmental controls and statistical properties of SOC stocks can improve ESM land model benchmarking and perhaps allow representation of spatial heterogeneity of biogeochemistry at scales finer than those currently resolved by ESMs.

54 ENVIRONMENTAL SCIENCES↗

Scaling effects in the static and dynamic response of graphite-epoxy beam-columns

Scale model technology represents one method of investigating the behavior of advanced, weight-efficient composite structures under a variety of loading conditions. It is necessary, however, to understand the limitations involved in testing scale model structures before the technique can be fully utilized. These limitations, or scaling effects, are characterized. in the large deflection response and failure of composite beams. Scale model beams were loaded with an eccentric axial compressive load designed to produce large bending deflections and global failure. A dimensional analysis was performed on the composite beam-column loading configuration to determine a model law governing the system response. An experimental program was developed to validate the model law under both static and dynamic loading conditions. Laminate stacking sequences including unidirectional, angle ply, cross ply, and quasi-isotropic were tested to examine a diversity of composite response and failure modes. The model beams were loaded under scaled test conditions until catastrophic failure. A large deflection beam solution was developed to compare with the static experimental results and to analyze beam failure. Also, the finite element code DYCAST (DYnamic Crash Analysis of STructure) was used to model both the static and impulsive beam response. Static test results indicate that the unidirectional and cross ply beam responses scale as predicted by the model law, even under severe deformations. In general, failure modes were consistent between scale models within a laminate family; however, a significant scale effect was observed in strength. The scale effect in strength which was evident in the static tests was also observed in the dynamic tests. Scaling of load and strain time histories between the scale model beams and the prototypes was excellent for the unidirectional beams, but inconsistent results were obtained for the angle ply, cross ply, and quasi-isotropic beams. Results show that valuable information can be obtained from testing on scale model composite structures, especially in the linear elastic response region. However, due to scaling effects in the strength behavior of composite laminates, caution must be used in extrapolating data taken from a scale model test when that test involves failure of the structure.

Jackson, Karen E.↗

Advances in time-scale algorithms

The term clock is usually used to refer to a device that counts a nearly periodic signal. A group of clocks, called an ensemble, is often used for time keeping in mission critical applications that cannot tolerate loss of time due to the failure of a single clock. The time generated by the ensemble of clocks is called a time scale. The question arises how to combine the times of the individual clocks to form the time scale. One might naively be tempted to suggest the expedient of averaging the times of the individual clocks, but a simple thought experiment demonstrates the inadequacy of this approach. Suppose a time scale is composed of two noiseless clocks having equal and opposite frequencies. The mean time scale has zero frequency. However if either clock fails, the time-scale frequency immediately changes to the frequency of the remaining clock. This performance is generally unacceptable and simple mean time scales are not used. First, previous time-scale developments are reviewed and then some new methods that result in enhanced performance are presented. The historical perspective is based upon several time scales: the AT1 and TA time scales of the National Institute of Standards and Technology (NIST), the A.1(MEAN) time scale of the US Naval observatory (USNO), the TAI time scale of the Bureau International des Poids et Measures (BIPM), and the KAS-1 time scale of the Naval Research laboratory (NRL). The new method was incorporated in the KAS-2 time scale recently developed by Timing Solutions Corporation. The goal is to present time-scale concepts in a nonmathematical form with as few equations as possible. Many other papers and texts discuss the details of the optimal estimation techniques that may be used to implement these concepts.

Stein, S. R.↗

Sub-ply level scaling approach investigated for graphite-epoxy composite beam columns

Scale model graphite-epoxy composite specimens were fabricated using the 'sub-ply level' approach and tested as beam-columns under an eccentric axial load to determine the effect of specimen size on flexural response and failure. In the current research project, although the fiber diameters are not scaled, the thickness of the pre-preg material itself has been scaled by adjusting the number of fibers through the thickness of a single ply. Three different grades of graphite-epoxy composite material (AS4/3502) were obtained from Hercules, Inc., in which the number of fibers through the thickness of a single ply was reduced (Grade 190 with 12 to 16 fibers, Grade 95 with 6 to 8 fibers, and Grade 48 with 3 to 4 fibers). Thus, using the sub-ply level approach, a baseline eight ply quasi-isotropic laminate could be fabricated using either the Grade 48 or Grade 95 material and the corresponding full-scale laminate would be constructed from Grade 95 or standard Grade 190 material, respectively. Note that in the sub-ply level approach, the number of ply interfaces is constant for the baseline and full-scale laminates. This is not true for the ply level and sublaminate level scaled specimens. The three grades of graphite-epoxy composite material were used to fabricate scale model beam-column specimens with in-plane dimensions of 0.5*n x 5.75*n, where n=1,2,4 corresponsing to 1/4, 1/2, and full-scale factors. Angle ply, cross ply, and quasi-isotropic laminate stacking sequences were chosen for the investigation and the test matrices for each laminate type are given. Specimens in each laminate family with the same in-plane dimensions but different thicknesses were tested to isolate the influence of the thickness dimension on the flexural response and failure. Also, specific lay-ups were chosen with blocked plies and dispersed plies for each laminate type. Specimens were subjected to an eccentric axial load until failure. The load offset was introduced through a set of hinges which were attached to the platens of a standard load test machine. Three sets of geometrically scaled hinges were used to ensure that scaled loading conditions were applied. This loading condition was chosen because it promotes large flexural deformations and specimens fail at the center of the beam, away from the grip supports. Five channels of data including applied vertical load, end shortening displacement, strain from gages applied back-to-back at the midspan of the beam, and rotation of the hinge from a bubble inclinometer were recorded for each specimen. The beam-column test configuration was used previously to study size effects in ply level scaled composite specimens of the same material system, sizes, and stacking sequences. Thus, a direct comparison between the two scaling approaches is possible. Ply level scaled beam-columns with angle ply, cross ply, and quasi-isotropic lay-ups exhibited no size dependencies in the flexural response, but significant size effects in strength. The reduction in strength with increasing specimen size was not predicted successfully by analysis techniques. It is anticipated that results from this investigation will lead to a better understanding of the strength scale effect in composite structures.

Jackson, Karen E.↗

Effect of Velocity in Icing Scaling Tests

This paper presents additional results of a study first published in 1999 to determine the effect of scale velocity on scaled icing test results. Reference tests were made with a 53.3-cm-chord NACA 0012 airfoil model in the NASA Glenn Icing Research Tunnel at an airspeed of 67 m/s, an MVD of 40 microns, and an LWC of 0.6 g/cu m. Temperature was varied to provide nominal freezing fractions of 0.8, 0.6, and 0.5. Scale tests used both 35.6- and 27.7-cm-chord 0012 models for 2/3- and 1/2-size scaling. Scale test conditions were found using the modified Ruff (AEDC) scaling method with the scale velocity determined in five ways. Four of the scale velocities were found by matching the scale and reference values of water-film thickness, velocity, Weber number, and Reynolds number. The fifth scale velocity was simply the average of those found by matching the Weber and Reynolds numbers. The resulting scale velocities ranged from 85 to 220 percent of the reference velocity. For a freezing fraction of 0.8, the value of the scale velocity had no effect on how well the scale ice shape simulated the reference shape. For nominal freezing fractions of 0.5 and 0.6, the best simulation of the reference shape was achieved when the scale velocity was the average of the constant-Weber-number and the constant-Reynolds-number velocities.

Anderson, David N.↗