SEARCH · Engineering Papers

Results for “Normal Distribution”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74

All sources NASA NTRS DOE OSTI

DOE OSTI2016

WE-H-BRA-01: BEST IN PHYSICS (THERAPY): Nano-Dosimetric Kinetic Model for Variable Relative Biological Effectiveness of Proton and Ion Beams

Purpose: Recent clonogenic cell survival and γH2AX studies suggest proton relative biological effectiveness (RBE) may be a non-linear function of linear energy transfer (LET) in the distal edge of the Bragg peak and beyond. We sought to develop a multiscale model to account for non-linear response phenomena to aid in the optimization of intensity-modulated proton therapy. Methods: The model is based on first-principle simulations of proton track structures, including secondary ions, and an analytical derivation of the dependence on particle LET of the linear-quadratic (LQ) model parameters α and β. The derived formulas are an extension of the microdosimetric kinetic (MK) model that captures dissipative track structures and non-Poissonian distribution of DNA damage at the distal edge of the Bragg peak and beyond. Monte Carlo simulations were performed to confirm the non-linear dose-response characteristics arising from the non-Poisson distribution of initial DNA damage. Results: In contrast to low LET segments of the proton depth dose, from the beam entrance to the Bragg peak, strong deviations from non-dissipative track structures and Poisson distribution in the ionization events in the Bragg peak distal edge govern the non-linear cell response and result in the transformation α=(1+c-1 L) α-x+2(c-0 L+c-2 L^2 )(1+c-1 L) β-x and β=(1+c-1 L)^2 β-x. Here L is the charged particle LET, and c-0,c-1, and c-2 are functions of microscopic parameters and can be served as fitting parameters to the cell-survival data. In the low LET limit c-1, and c-2 are negligible hence the linear model proposed and used by Wilkins-Oelfke for the proton treatment planning system can be retrieved. The present model fits well the recent clonogenic survival data measured recently in our group in MDACC. Conclusion: The present hybrid method provides higher accuracy in calculating the RBE-weighted dose in the target and normal tissues.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2014

Dynamics of Confined Flexible and Unentangled Polymer Melts in Highly Adsorbing Cylindrical Pores

Inspired by the recent neutron spin echo experiments (NSE) designed to address the dynamic phenomena in polymer melts that are induced by interactions with a confining surface, we performedcoarse-grained molecular dynamics simulations to replicate the experimental results in order to provide new molecular insight for the observations. The results show excellent agreement in the values obtained for the normalized coherent single chain dynamic structure factor, S(Q;Dt)/S(Q;0) , between experiments and simulations. The simulations indicate that using different chain molecular weights that are used to achieve scattering contrast results in the uneven distribution of scatterers in the radial direction because lower molecular weight chains prefer to adsorb in the confining surface. In the bulk configuration, both simulations and experiments confirm that the polymer chains follow Rouse dynamics. However, under confinement, the Rouse modes are suppressed. The mean-squared radius of gyration, R2g, and the average relative shape anisotropy k2 of the conformation of the polymer chains indicate a pancake-like conformation near the surface and a bulk-like conformation near the center of the confining cylinder. Direct visualization of the polymers in the simulation confirm the pancake-like conformation of the adsorbed chains and the presence of trains, loops and tails in the region between the adsorbed chains and the chains not in contact with the surface. Despite the presence of these different conformations, the average form factor of the confined chains still follows the Debye function which describes linear ideal chains, which is in agreement with small angle neutron scattering (SANS) experiments. The experimentally inaccessible mean squared displacement of the confined monomers, calculated as a function of radial distance from the pore surface, was obtained in the simulations. The simulation shows a gradual increase of the MSD from the adsorbed, but mobile layer, to that similar to the bulk far away from the surface.

Carrillo, Jan-Michael Y.↗

DOE OSTI2023

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 1: Experimental protocols and surface changes

There is now substantial literature on climate model studies of equatorial or tropical stratospheric SO 2 injections that aim to counteract the surface warming produced by rising concentrations of greenhouse gases. Here we present the results from the first systematic intercomparison of climate responses in three Earth system models wherein the injection of SO 2 occurs at different latitudes in the lower stratosphere: CESM2-WACCM6, UKESM1.0 and GISS-E2.1-G. The first two use a modal aerosol microphysics scheme, while two versions of GISS-E2.1-G use a bulk aerosol (One-Moment Aerosol, OMA) and a two-moment (Multiconfiguration Aerosol TRacker of mIXing state, MATRIX) microphysics approach, respectively. Our aim in this work is to determine commonalities and differences between the climate model responses in terms of the distribution of the optically reflective sulfate aerosols produced from the oxidation of SO 2 and in terms of the surface response to the resulting reduction in solar radiation. A focus on understanding the contribution of characteristics of models transport alongside their microphysical and chemical schemes, and on evaluating the resulting stratospheric responses in different models, is given in the companion paper (Bednarz et al., 2023). The goal of this exercise is not to evaluate these single-point injection simulations as stand-alone proposed strategies to counteract global warming; instead we determine sources and areas of agreement and uncertainty in the simulated responses and, ultimately, the possibility of designing a comprehensive intervention strategy capable of managing multiple simultaneous climate goals through the combination of different injection locations. We find large disagreements between GISS-E2.1-G and the CESM2-WACCM6 and UKESM1.0 models regarding the magnitude of cooling per unit of aerosol optical depth (AOD) produced, which varies from 4.7 K per unit of AOD in CESM2-WACCM6 to 16.7 K in the GISS-E2.1-G version with two-moment aerosol microphysics. By normalizing the results with the global mean response in each of the models and thus assuming that the amount of SO 2 injected is a free parameter that can be managed independently, we highlight some commonalities in the overall distributions of the aerosols, in the inter-hemispheric surface temperature response and in shifts to the Intertropical Convergence Zone, as well as some areas of disagreement, such as the extent of the aerosol confinement in the equatorial region and the efficiency of the transport to polar latitudes. In conclusion, we demonstrate that it is possible to use these simulations to produce more comprehensive injection strategies in multiple climate models. However, large differences in the injection magnitudes can be expected, potentially increasing inter-model spreads in some stratospheric quantities (such as aerosol distribution) while reducing the spread in the surface response in terms of temperature and precipitation; furthermore, the selection of the injection locations may be dependent on the models' specific stratospheric transport.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

An Analysis of Post Attack Impacts and Effects of Learning Parameters on Vulnerability Assessment of Power Grid

Due to the increasing number of heterogeneous devices connected to electric power grid, the attack surface increases the threat actors. Game theory and machine learning are being used to study the power system failures caused by external manipulation. Most of existing works in the literature focus on one-shot process of attacks and fail to show the dynamic evolution of the defense strategy. In this paper, we focus on an adversarial multistage sequential game between the adversaries of the smart electric power transmission and distribution system. We study the impact of exploration rate and convergence of the attack strategies (sequences of action that creates large scale blackout based on the system capacity) based on the reinforcement learning approach. We also illustrate how the learned attack actions disrupt the normal operation of the grid by creating transmission line outages, bus voltage violations, and generation loss. This simulation studies are conducted on IEEE 9 and 39 bus systems. The results show the improvement of the defense strategy through the learning process. The results also prove the feasibility of the learned attack actions by replicating the disturbances created in simulated power system.

attack impacts↗

DOE OSTI2023

Field validation of isotropic analytical models for simulating fabric shades

Fabric roller shades are common shading materials used in commercial and residential buildings. Accurately characterizing and modeling shades helps practitioners select the appropriate product and its control strategy based on climate and occupants' priorities, such as visual comfort and view to outdoors. Previous studies established a generalized method for modeling complex fenestration systems using data-driven tabulated bidirectional scattering distribution functions. However, deploying such a method at scale to all fabric shading products on the market is too costly and time-consuming. Analytical models that are based on a limited set of measurements (e.g., normal-normal and normal-hemispherical visible transmittance and reflectance, and directional cut-off angles) can be used to model the wide variety of shading products on the market. This study evaluates the performance of two isotropic analytical models, Roos-Wienold and Modified-Kotey, for modeling fabric roller shades, with a focus on the model's ability to predict occupant visual comfort. The performance evaluation was conducted through laboratory and field measurements and simulations. The results showed that both models are sufficient for predicting vertical illuminance at seated eye-level. Roos-Wienold model was able to predict binary visual comfort classification (glare/no-glare) under a wide range of luminance conditions, while Modified-Kotey model did not perform as well under high-contrast low-adaptation conditions. Both models are insufficient in predicting visual comfort at a four-point scale (e.g., imperceptible, perceptible, disturbing, intolerable). The two isotropic models become less accurate when the fabric exhibits high anisotropy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2018

Turbulent transport stabilization by ICRH minority fast ions in low rotating JET ILW L-mode plasmas

The first experimental demonstration that fast ion induced stabilization of thermal turbulent transport takes place also at low values of plasma toroidal rotation has been obtained in JET ILW (ITER-like wall) L-mode plasmas with high ( 3 He)-D ICRH (ion cyclotron resonance heating) power. A reduction of the gyro-Bohm normalized ion heat flux and higher values of the normalized ion temperature gradient have been observed at high ICRH power and low NBI (neutral beam injection) power and plasma rotation. Gyrokinetic simulations indicate that ITG (ion temperature gradient) turbulence stabilization induced by the presence of high-energetic 3 He ions is the key mechanism in order to explain the experimental observations. Two main mechanisms have been identified to be responsible for the turbulence stabilization: a linear electrostatic wave-fast particle resonance mechanism and a nonlinear electromagnetic mechanism. The dependence of the stabilization on the 3 He distribution function has also been studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Plasma physical parameters along CME-driven shocks. II. Observation–simulation comparison

In this work, we compare the spatial distribution of the plasma parameters along the 1999 June 11 coronal mass ejection (CME)-driven shock front with the results obtained from a CME-like event simulated with the FLIPMHD3D code, based on the FLIP-MHD particle-in-cell method. The observational data are retrieved from the combination of white-light coronagraphic data (for the upstream values) and the application of the Rankine–Hugoniot equations (for the downstream values). The comparison shows a higher compression ratio X and Alfvénic Mach number M{sub A} at the shock nose, and a stronger magnetic field deflection d toward the flanks, in agreement with observations. Then, we compare the spatial distribution of M{sub A} with the profiles obtained from the solutions of the shock adiabatic equation relating M{sub A}, X, and θ{sub Bn} (the angle between the upstream magnetic field and the shock front normal) for the special cases of parallel and perpendicular shock, and with a semi-empirical expression for a generically oblique shock. The semi-empirical curve approximates the actual values of M{sub A} very well, if the effects of a non-negligible shock thickness δ{sub sh} and plasma-to magnetic pressure ratio β{sub u} are taken into account throughout the computation. Moreover, the simulated shock turns out to be supercritical at the nose and sub-critical at the flanks. Finally, we develop a new one-dimensional Lagrangian ideal MHD method based on the GrAALE code, to simulate the ion-electron temperature decoupling due to the shock transit. Two models are used, a simple solar wind model and a variable-γ model. Both produce results in agreement with observations, the second one being capable of introducing the physics responsible for the additional electron heating due to secondary effects (collisions, Alfvén waves, etc.).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Kinetic transverse dispersion relation for relativistic magnetized electron-positron plasmas with Maxwell-Jüttner velocity distribution functions

We use a kinetic treatment to study the linear transverse dispersion relation for a magnetized isotropic relativistic electron-positron plasma with finite relativistic temperature. The explicit linear dispersion relation for electromagnetic waves propagating along a constant background magnetic field is presented, including an analytical continuation to the whole complex frequency plane for the case of Maxwell-Jüttner velocity distribution functions. This dispersion relation is studied numerically for various temperatures. For left-handed solutions, the system presents two branches, the electromagnetic ordinary mode and the Alfvén mode. In the low frequency regime, the Alfvén branch has two dispersive zones, the normal zone (where ∂ω/∂k > 0) and an anomalous zone (where ∂ω/∂k < 0). We find that in the anomalous zone of the Alfvén branch, the electromagnetic waves are damped, and there is a maximum wave number for which the Alfvén branch is suppressed. We also study the dependence of the Alfvén velocity and effective plasma frequency with the temperature. We complemented the analytical and numerical approaches with relativistic full particle simulations, which consistently agree with the analytical results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Effects of chemical segregation on ductility-anisotropy in high strength Fe-Mn-Al-C lightweight austenitic steels

Ductile Fe-Mn-Al-C lightweight steels offer great potential as high specific strength materials for weight critical applications. In the present work, very high yield (~1 GPa) and ultimate tensile strength levels (~1.5 GPa) were obtained along the rolling and transverse directions of Fe-30Mn-8.5Al-0.9Si-0.9C-0.5Mo (wt.%), with high tensile elongation to failure of ~35%. This excellent mechanical behavior was attributed to the nanoscale, ordered κ-carbide precipitates formed during aging after rolling. However, tensile ductility was not isotropic through all directions of the rolled plate. In particular, the ductility was consistently lower than 5% true strain along the normal direction. Upon detailed microstructural investigations, micro-scale lamellar chemical micro-segregation bands were detected, which caused the observed anisotropic embrittlement along certain plate directions. These bands resulted in low failure strains along the normal direction since it is oriented perpendicular to these chemical segregation bands, resulting in a failure type reminiscent of delamination-induced cleavage fracture. Chemical segregation causes an inhomogeneous κ-carbide precipitate distribution upon aging, producing more pronounced mechanical anisotropy after peak-aging as compared to the solution heat treated samples. Surprisingly, deformation twinning was observed after tensile deformation of solution-heat treated specimens via electron-backscatter diffraction experiments, which is anomalous for a high stacking faulty energy (SFE) material. Further, it was proposed that chemical micro-segregation promotes locally different SFE regions, leading to deformation twinning in this nominally high SFE material. Overall, this work demonstrates the effects of chemical segregation on mechanical properties and microstructural evolution in a high performance, lightweight Fe-Mn-Al-C steel.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Tolerance of the Brachial Plexus to High-Dose Reirradiation

Purpose: To study the tolerance of the brachial plexus to high doses of radiation exceeding historically accepted limits by analyzing human subjects treated with reirradiation for recurrent tumors of the head and neck. Methods and Materials: Data from 43 patients who were confirmed to have received overlapping dose to the brachial plexus after review of radiation treatment plans from the initial and reirradiation courses were used to model the tolerance of this normal tissue structure. A standardized instrument for symptoms of neuropathy believed to be related to brachial plexus injury was utilized to screen for toxicity. Cumulative dose was calculated by fusing the initial dose distributions onto the reirradiation plan, thereby creating a composite plan via deformable image registration. The median elapsed time from the initial course of radiation therapy to reirradiation was 24 months (range, 3-144 months). Results: The dominant complaints among patients with symptoms were ipsilateral pain (54%), numbness/tingling (31%), and motor weakness and/or difficulty with manual dexterity (15%). The cumulative maximum dose (Dmax) received by the brachial plexus ranged from 60.5 Gy to 150.1 Gy (median, 95.0 Gy). The cumulative mean (Dmean) dose ranged from 20.2 Gy to 111.5 Gy (median, 63.8 Gy). The 1-year freedom from brachial plexus–related neuropathy was 67% and 86% for subjects with a cumulative Dmax greater than and less than 95.0 Gy, respectively (P=.05). The 1-year complication-free rate was 66% and 87%, for those reirradiated within and after 2 years from the initial course, respectively (P=.06). Conclusion: The development of brachial plexus–related symptoms was less than expected owing to repair kinetics and to the relatively short survival of the subject population. Time-dose factors were demonstrated to be predictive of complications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2016

Determination of grain boundary mobility during recrystallization by statistical evaluation of electron backscatter diffraction measurements

One of the key aspects influencing microstructural design pathways in metallic systems is grain boundary motion. The present work introduces a method by means of which direct measurement of grain boundary mobility vs. misorientation dependence is made possible. The technique utilizes datasets acquired by means of serial electron backscatter diffraction (EBSD) measurements. The experimental EBSD measurements are collectively analyzed, whereby datasets were used to obtain grain boundary mobility and grain aspect ratio with respect to grain boundary misorientation. The proposed method is further validated using cellular automata (CA) simulations. Single crystal aluminium was cold rolled and scratched in order to nucleate random orientations. Subsequent annealing at 300 °C resulted in grains growing, in the direction normal to the scratch, into a single deformed orientation. Growth selection was observed, wherein the boundaries with misorientations close to Σ7 CSL orientation relationship (38° 〈111〉) migrated considerably faster. The obtained boundary mobility distribution exhibited a non-monotonic behavior with a maximum corresponding to misorientation of 38° ± 2° about 〈111〉 axes ± 4°, which was 10–100 times higher than the mobility values of random high angle boundaries. Correlation with the grain aspect ratio values indicated a strong growth anisotropy displayed by the fast growing grains. The observations have been discussed in terms of the influence of grain boundary character on grain boundary motion during recrystallization. - Highlights: • Statistical microstructure method to measure grain boundary mobility during recrystallization • Method implementation independent of material or crystal structure • Mobility of the Σ7 boundaries in 5N Al was calculated as 4.7 × 10{sup –8} m{sup 4}/J ⋅ s. • Pronounced growth selection in the recrystallizing nuclei in Al • Boundary mobility values during recrystallization 2–3 orders of magnitude larger vis-à-vis curvature driven motion.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Sensitivity and uncertainty of the IFR-1 BISON benchmark

The fuel performance code BISON is being used to evaluate metallic fuel for a new fast-spectrum test reactor called the Versatile Test Reactor, which is being considered for adoption by the US Department of Energy. To quantify the accuracy of BISON predictions, researchers at Oak Ridge National Laboratory have been developing a series of benchmarks based on legacy metallic fuel experiments. As part of this effort, the sensitivity of BISON predictions to variations in model inputs and the uncertainties associated with BISON predictions must be established. This paper summarizes efforts to perform a comprehensive sensitivity analysis (SA) and uncertainty quantification (UQ) on a benchmark based on the IFR-1 experiment.For the SA, at least one input was chosen from every BISON model and physics module used in the benchmark. The inputs were varied individually in a series of BISON simulations. Here, the resulting variations in benchmark predictions were normalized to calculate sensitivities. These sensitivities were then used to inform input selections for the UQ.The UQ was performed using the Monte Carlo UQ method. A literature review was conducted to estimate uncertainty distributions for the selected inputs, and values were sampled randomly from each distribution in a series of BISON simulations. Variations in the benchmark predictions were used to estimate uncertainty distributions and confidence intervals. It was found that nearly 100% of the benchmark predictions matched the corresponding legacy values within the confidence intervals. However, this is at least partially because the confidence intervals associated with benchmark predictions were wide. The uncertainty contributions of assumptions in the benchmark, experimental uncertainties, and BISON models were quantified. Some analysis was performed to identify inputs that contributed to the uncertainties. Finally, recommendations are made for future benchmark and future BISON development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

DFT-Based Calculation of Molecular Hyperpolarizability and SFG Intensity of Symmetric and Asymmetric Stretch Modes of Alkyl Groups

Vibrational sum frequency generation (SFG) spectroscopy has been extensively used for obtaining structural information of molecular functional groups at two-dimensional (2D) interfaces buried in the gas or liquid medium. Although the SFG experiment can be done elegantly, interpreting the measured intensity in terms of molecular orientation with respect to the lab coordinate is quite complicated. One of the main reasons is the difficulty of determining the hyperpolarizability tensors of even simple molecules that govern their SFG responses. The single-bond polarizability derivative (SBPD) model has been proposed to estimate the relative magnitude of SFG-active hyperpolarizability by assuming that the perturbation associated to each vibration is negligible. In this study, density functional theory (DFT) was used to calculate the polarizability and dipole derivative tensors of the CH 3 stretch mode of CH 3 I, CH 3 CH 2 I, CH 3 OH, CH 3 CH 2 OH. Then, the hyperpolarizability tensors of symmetric and asymmetric vibration modes were calculated considering the Boltzmann distribution of representative conformers, which allowed theoretical calculation of their SFG intensities at all polarization combinations as a function of the tilt angle of the CH 3 group with respect to the surface normal direction. Then the ratios of the calculated SFG intensities for the CH 3 peaks used in experimental studies for the CH 3 tilt angle determination were compared. This comparison clearly showed the effect of vibrational coupling among neighboring functional groups. As a result, this study presents new parameters that can be used in determining the average tilt angle of the CH 3 group at the 2D interface with SFG measurements as well as limitations of the method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Supernovae and their host galaxies – VII. The diversity of Type Ia supernova progenitors

ABSTRACT We present an analysis of the light curve (LC) decline rates (Δm15) of 407 normal and peculiar supernovae (SNe) Ia and global parameters of their host galaxies. As previously known, there is a significant correlation between the Δm15 of normal SNe Ia and global ages (morphologies, colours, and masses) of their hosts. On average, those normal SNe Ia that are in galaxies from the Red Sequence (early-type, massive, old hosts) have faster declining LCs in comparison with those from the Blue Cloud (late-type, less massive, younger hosts) of the colour–mass diagram. The observed correlations between the Δm15 of normal SNe Ia and hosts’ parameters appear to be due to the superposition of at least two distinct populations of faster and slower declining normal SNe Ia from older and younger stellar components. We show, for the first time, that the Δm15 of 91bg- and 91T-like SNe is independent of host morphology and colour. The distribution of hosts on the colour–mass diagram confirms the known tendency for 91bg-like SNe to occur in globally red/old galaxies, while 91T-like events prefer blue/younger hosts. On average, the youngest global ages of 02cx-like SNe hosts and their positions in the colour–mass diagram hint that these events likely originate from young population, but they differ from 91T-like events in the LC decline rate. Finally, we discuss the possible explosion channels and present our favoured SN Ia models that have the potential to explain the observed SN–host relations.

Hakobyan, A. A.↗

DOE OSTI2020

Analytic Thermal Model of an Optical Fiber Based Gamma Thermometer and its Application in a University Research Reactor

This paper describes and validates, by comparison with numerical modeling results, an analytical model of thermal transport in an optical fiber based gamma thermometer (OFBGT) that is appropriate for use in university research reactors. The maximum temperature difference between the thermal mass and the outer sheath ( Δ T ) for the OFBGT design that we have considered is approximately 50 ∘ C , for the OFBGT in the Central Irradiation Facility of the Ohio State University Research Reactor (OSURR) with the reactor operating at full power (450 kW). The maximum value of Δ T that is predicted by the analytic model for the OFBGT design is smaller by approximately 1.1 °C than the value of Δ T which is predicted by the numerical model, for the same OFBGT design, but including all the details of the design. We have used the analytical model of thermal transport in an OFBGT to determine a normalized Modulation Transfer Function M T F ′ ( k ) ) for the OFBGT. We conclude that MTF ′ ( k OSU ) > 0.99 , where k OSU is the spatial frequency for the axial dependence of the reactor power distribution in the

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Towards Integrating Data Quality Assessments and Radiometer Uncertainty for Determining the Expanded Uncertainty of Three-Component Solar Radiation Measurements

Accurate solar irradiance data are fundamental for determining the design and performance characteristics of photovoltaic systems. The uncertainty of solar irradiance measurements depends on many factors including radiometer design, calibration, installation, maintenance, and operational environment. The key contributors to this uncertainty can be classified as the measurement uncertainty of a particular radiometer and the operational uncertainty determined for the time of measurement. Radiometer measurement uncertainty estimates (U R ) can be based on well-established methods used as part of the radiometer calibration process. Estimates of operational uncertainties (U o ) require consideration of additional site-specific factors that affect data quality. A method is needed for establishing the accuracy of solar irradiance data by integrating an existing data quality process and measurement uncertainty estimates for specific radiometers. An algorithm has been developed to integrate data quality analyses and measurement uncertainty estimates for three-component solar irradiance data: global horizontal (total hemispheric) irradiance, direct normal (beam) irradiance, and diffuse horizontal (sky) irradiance collected at one- to 60-minute intervals. The algorithm has been tested using one-minute irradiance measurements. The goal of the project is to distribute a user-friendly software package based on the new algorithm.

data integrity↗

DOE OSTI2016

The influence of leaf size and shape on leaf thermal dynamics: does theory hold up under natural conditions?

Laboratory studies on artificial leaves suggest that leaf thermal dynamics are strongly influenced by the two-dimensional size and shape of leaves and associated boundary layer thickness. Hot environments are therefore said to favour selection for small, narrow or dissected leaves. Empirical evidence from real leaves under field conditions is scant and traditionally based on point measurements that do not capture spatial variation in heat load. Here in this study, we used thermal imagery under field conditions to measure the leaf thermal time constant (τ) in summer and the leaf-to-air temperature difference (ΔT) and temperature range across laminae (T range ) during winter, autumn and summer for 68 Proteaceae species. We investigated the influence of leaf area and margin complexity relative to effective leaf width (w e ), the latter being a more direct indicator of boundary layer thickness. Normalized difference of margin complexity had no or weak effects on thermal dynamics, but w e strongly predicted τ and ΔT, whereas leaf area influenced T range . Unlike artificial leaves, however, spatial temperature distribution in large leaves appeared to be governed largely by structural variation. Therefore, we agree that small size, specifically we, has adaptive value in hot environments but not with the idea that thermal regulation is the primary evolutionary driver of leaf dissection.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Penumbra: A spatially distributed, mechanistic model for simulating ground-level incident solar energy across heterogeneous landscapes

Landscape solar energy is a significant environmental driver, yet it remains complicated to model well. Several solar radiation models simplify the complexity of light by estimating it at discrete point locations or by averaging values over larger areas. These modeling approaches may be useful in certain cases, but they are unable to provide spatially distributed and temporally dynamic representations of solar energy across entire landscapes. We created a landscape-scale ground-level shade and solar energy model called Penumbra to address this deficiency. Penumbra simulates spatially distributed ground-level shade and incident solar energy at user-defined timescales by modeling local and distant topographic shading and vegetative shading. Spatially resolved inputs of a digital elevation model, a normalized digital surface model, and landscape object transmittance are used to estimate spatial variations in solar energy at user-defined temporal timesteps. The research goals for Penumbra included: 1) simulations of spatiotemporal variations of shade and solar energy caused by both objects and topographic features, 2) minimal user burden and parameterization, 3) flexible user defined temporal parameters, and 4) flexible external model coupling. We test Penumbra’s predictive skill by comparing the model’s predictions with monitored open and forested sites, and achieve calibrated mean errors ranging from -17.3 to 148.1 μmoles/m 2 /s. Penumbra is a dynamic model that can produce spatial and temporal representations of shade percentage and ground-level solar energy. Outputs from Penumbra can be used with other ecological models to better understand the health and resilience of aquatic, near stream terrestrial, and upland ecosystems.

54 ENVIRONMENTAL SCIENCES↗