Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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,279 records · Page 71

Final report on "Modeling Diurnal Variations of California Land Biosphere CO2 Fluxes"

In Mediterranean climates, the season of water availability (winter) is out of phase with the season of light availability and atmospheric demand for moisture (summer). Multi-year half-hourly observations of sap flow velocities in 26 evergreen trees in a small watershed in Northern California show that different species of evergreen trees have different seasonalities of transpiration: Douglas-firs respond immediately to the first winter rain, while Pacific madrones have peak transpiration in the dry summer. Using these observations, we have derived species-specific parameterization of normalized sap flow velocities in terms of insolation, vapor pressure deficit and near-surface soil moisture. A simple 1-D boundary layer model showed that afternoon temperatures may be higher by 1 degree Celsius in an area with Douglas-firs than with Pacific madrones. The results point to the need to develop a new representation of subsurface moisture, in particular pools beneath the organic soil mantle and the vadose zone. Our ongoing and future work includes coupling our new parameterization of transpiration with new representation of sub-surface moisture in saprolite and weathered bedrock. The results will be implemented in a regional climate model to explore vegetation-climate feedbacks, especially in the dry season.

58 GEOSCIENCES↗

Framework for Probabilistic Projections of Energy-Relevant Streamflow Indicators under Climate Change Scenarios for the U.S.

This project focuses on uncertainty in streamflow forecasting under climate change conditions. The objective is to develop easy to use methodologies that can be applied across a range of river basins to estimate changes in water availability for realistic projections of climate change. There are three major components to the project: Empirical downscaling of regional climate change projections from a range of Global Climate Models; Developing a methodology to use present day information on the climate controls on the parameterizations in streamflow models to adjust the parameterizations under future climate conditions (a trading-space-for-time approach); and Demonstrating a bottom-up approach to establishing streamflow vulnerabilities to climate change. The results reinforce the need for downscaling of climate data for regional applications, and further demonstrates the challenges of using raw GCM data to make local projections. In addition, it reinforces the need to make projections across a range of global climate models. The project demonstrates the potential for improving streamflow forecasts by using model parameters that are adjusted for future climate conditions, but suggests that even with improved streamflow models and reduced climate uncertainty through the use of downscaled data, there is still large uncertainty is the streamflow projections. The most useful output from the project is the bottom-up vulnerability driven approach to examining possible climate and land use change impacts on streamflow. Here, we demonstrate an inexpensive and easy to apply methodology that uses Classification and Regression Trees (CART) to define the climate and environmental parameters space that can produce vulnerabilities in the system, and then feeds in the downscaled projections to determine the probability top transitioning to a vulnerable sate. Vulnerabilities, in this case, are defined by the end user.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the uncertainties of aerosol indirect effects and impacts on decadal-scale climate variability in NCAR CAM5 and CESM1

The main goal of this project is to systematically quantify the major uncertainties of aerosol indirect effects due to the treatment of moist turbulent processes that drive aerosol activation, cloud macrophysics and microphysics in response to anthropogenic aerosol perturbations using the CAM5/CESM1. To achieve this goal, the P.I. hired a postdoctoral research scientist (Dr. Anna Fitch) who started her work from the Nov.1st.2012. In order to achieve the project goal, the first task that the Postdoc. and the P.I. did was to quantify the role of subgrid vertical velocity variance on the activation and nucleation of cloud liquid droplets and ice crystals and its impact on the aerosol indirect effect in CAM5. First, we analyzed various LES cases (from dry stable to cloud-topped PBL) to check whether this isotropic turbulence assumption used in CAM5 is really valid. It turned out that this isotropic turbulence assumption is not universally valid. Consequently, from the analysis of LES, we derived an empirical formulation relaxing the isotropic turbulence assumption used for the CAM5 aerosol activation and ice nucleation, and implemented the empirical formulation into CAM5/CESM1, and tested in the single-column and global simulation modes, and examined how it changed aerosol indirect effects in the CAM5/CESM1. These results were reported in the poster section in the 18th Annual CESM workshop held in Breckenridge, CO during Jun.17-20.2013. While we derived an empirical formulation from the analysis of couple of LES from the first task, the general applicability of that empirical formulation was questionable, because it was obtained from the limited number of LES simulations. The second task we did was to derive a more fundamental analytical formulation relating vertical velocity variance to TKE using other information starting from basic physical principles. This was a somewhat challenging subject, but if this could be done in a successful way, it could be directly implemented into the CAM5 as a practical parameterization, and substantially contributes to achieving the project goal. Through an intensive research for about one year, we found appropriate mathematical formulation and tried to implement it into the CAM5 PBL and activation routine as a practical parameterized numerical code. During these processes, however, the Postdoc applied for another position in Sweden, Europe, and accepted a job offer there, and left NCAR in August 2014. In Sweden, Dr. Anna Fitch is still working on this subject in a part time, planning to finalize the research and to write the paper in a near future.

58 GEOSCIENCES↗

Continuous Evaluation of Fast Processes in Climate Models Using Arm Measurements

Under the support of this grant, we investigated the fast process of interaction of clouds, shallow convection, and boundary layer turbulence and their parameterizations. Main accomplishments involve two things. One is the understanding of the physical mechanisms of low cloud feedbacks; the second is the development and evaluation of convection and cloud parameterizations in climate models.

54 ENVIRONMENTAL SCIENCES↗

Final Report

We participated in a FASTER SCM intercomparison for which we ran our SCM for 3 years at the SGP to analyze statistics of the precipitation field (Song et al., 2013). An important feature of these simulations was the use of relaxation forcing to observed T, q, which decouples the model convection from the forcing and allows precipitation errors to emerge. Because the GISS cumulus parameterization includes a trigger that prevents convection until sufficient lifting is present, and because convection at the SGP is usually triggered by mesoscale motions that are not represented in the forcing when relaxation is applied, the duration of SCM precipitation is shorter than observed (Del Genio and Wolf, 2012) and thus its mean precipitation less than observed. However, its diurnal cycle phase is correct, and it is the only operational cumulus parameterization in the intercomparison that does not produce excessive warm season precipitation under weak large-scale forcing conditions.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of an Incremental Ventilation Energy Model for Estimating Impacts of Air Sealing and Mechanical Ventilation

Changing the rate of airflow through a home affects the annual thermal conditioning energy. Large-scale changes to airflow rates of the housing stock can significantly alter the energy consumption of the residential energy sector. However, the complexity of existing residential energy models hampers the ability to estimate the impact of policy changes on a state or nationwide level. The Incremental Ventilation Energy (IVE) model developed in this study was designed to combine the output of simple airflow models and a limited set of home characteristics to estimate the associated change in energy demand of homes. The IVE model was designed specifically to enable modelers to use existing databases of home characteristics to determine the impact of policy on ventilation at a population scale. In this report, we describe the IVE model and demonstrate that its estimates of energy change are comparable to the estimates of a wellvalidated, complex residential energy model when applied to homes with limited parameterization. Homes with extensive parameterization would be more accurately characterized by complex residential energy models. The demonstration included a range of home types, climates, and ventilation systems that cover a large fraction of the residential housing sector.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Constraints on cloud feedback from analysis of ARM observations and models

Over the performance period, the PI and his team performed innovative analyses of surface remote sensing data from the Atmospheric Radiation Measurement (ARM) Program Southern Great Plains, Tropical Western Pacific, and North Slope of Alaska sites to elucidate some of the physical processes underlying cloud feedbacks on anthropogenic climate change. ARM data strongly influenced parameterization evaluation and development in the GISS (Goddard Institute for Space Studies) General Circulation Model (GCM) over this time. In this report we discuss how ARM funding was utilized for data analyses, model evaluation, and cumulus and stratiform cloud parameterization improvements, highlighting especially the areas of low-cloud feedbacks, cloud phase, and convective entrainment and downdrafts.

54 ENVIRONMENTAL SCIENCES↗

Determining Reduced Order Models for Optimal Stochastic Reduced Order Models

The use of parameterized reduced order models(PROMs) within the stochastic reduced order model (SROM) framework is a logical progression for both methods. In this report, five different parameterized reduced order models are selected and critiqued against the other models along with truth model for the example of the Brake-Reuss beam. The models are: a Taylor series using finite difference, a proper orthogonal decomposition of the the output, a Craig-Bampton representation of the model, a method that uses Hyper-Dual numbers to determine the sensitivities, and a Meta-Model method that uses the Hyper-Dual results and constructs a polynomial curve to better represent the output data. The methods are compared against a parameter sweep and a distribution propagation where the first four statistical moments are used as a comparison. Each method produces very accurate results with the Craig-Bampton reduction having the least accurate results. The models are also compared based on time requirements for the evaluation of each model where the Meta- Model requires the least amount of time for computation by a significant amount. Each of the five models provided accurate results in a reasonable time frame. The determination of which model to use is dependent on the availability of the high-fidelity model and how many evaluations can be performed. Analysis of the output distribution is examined by using a large Monte-Carlo simulation along with a reduced simulation using Latin Hypercube and the stochastic reduced order model sampling technique. Both techniques produced accurate results. The stochastic reduced order modeling technique produced less error when compared to an exhaustive sampling for the majority of methods.

97 MATHEMATICS AND COMPUTING↗

Uniaxial Creep as a Control on Mercury Intrusion Capillary Pressure in Consolidating Rock Salt

The nature of geologic disposal of nuclear waste in salt formations requires validated and verified two-phase flow models of transport of brine and gas through intact, damaged, and consolidating crushed salt. Such models exist in other realms of subsurface engineering for other lithologic classes (oil and gas, carbon sequestration etc. for clastics and carbonates) but have never been experimentally validated and parameterized for salt repository scenarios or performance assessment. Models for waste release scenarios in salt back-fill require phenomenological expressions for capillary pressure and relative permeability that are expected to change with degree of consolidation, and require experimental measurement to parameterize and validate. This report describes a preliminary assessment of the influence of consolidation (i.e. volume strain or porosity) on capillary entry pressure in two phase systems using mercury injection capillary pressure (MICP). This is to both determine the potential usefulness of the mercury intrusion porosimetry method, but also to enable a better experimental design for these tests. Salt consolidation experiments are performed using novel titanium oedometers, or uniaxial compression cells often used in soil mechanics, using sieved run-of-mine salt from the Waste Isolation Pilot Plant (WIPP) as starting material. Twelve tests are performed with various starting amounts of brine pore saturation, with axial stresses up to 6.2 MPa (~900 psi) and temperatures to 90°C. This corresponds to UFD Work Package 15SN08180211 milestone “FY:15 Transport Properties of Run-of-Mine Salt Backfill – Unconsolidated to Consolidated”. Samples exposed to uniaxial compression undergo time-dependent consolidation, or creep, to various degrees. Creep volume strain-time relations obey simple log-time behavior through the range of porosities (~50 to 2% as measured); creep strain rate increases with temperature and applied stress as expected. Mercury porosimetry is used to determine characteristic capillary pressure curves from a series of consolidation tests and show characteristic saturation-capillary pressure curves that follow the common van Genuchten (1978, 1980) formulation at low stresses. Higher capillary pressure data are suspect due to the large potential for sample damage, including fluid inclusion decrepitation and pore collapse. Data are supportive of use of the Leverett “J” function (Leverett, 1941) to use for scaling characteristic curves at different degrees of consolidation, but better permeability determinations are needed to support this hypothesis. Recommendations for further and refined testing are made with the goal of developing a self- consistent set of constitutive laws for granular salt consolidation and multiphase (brine-air) flow.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

HARP targets pion production cross section and yield measurements. Implications for MiniBooNE neutrino flux

The prediction of the muon neutrino flux from a 71.0 cm long beryllium target for the MiniBooNE experiment is based on a measured pion production cross section which was taken from a short beryllium target (2.0 cm thick - 5% nuclear interaction length) in the Hadron Production (HARP) experiment at CERN. To verify the extrapolation to our longer target, HARP also measured the pion production from 20.0 cm and 40.0 cm beryllium targets. The measured production yields, d 2 N π± (p; θ )=dpd Ω, on targets of 50% and 100% nuclear interaction lengths in the kinematic rage of momentum from 0.75 GeV/c to 6.5 GeV/c and the range of angle from 30 mrad to 210 mrad are presented along with an update of the short target cross sections. The best fitted extended Sanford-Wang (SW) model parameterization for updated short beryllium target π + production cross section is presented. Yield measurements for all three targets are also compared with that from the Monte Carlo predictions in the MiniBooNE experiment for different SW parameterization. The comparisons of v μ flux predictions for updated SW model is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Determining Best Estimates and Uncertainties in Cloud Microphysical Parameters from ARM Field Data: Implications for Models, Retrieval Schemes and Aerosol-Cloud-Radiation Interactions

We proposed to analyze in-situ cloud data collected during ARM/ASR field campaigns to create databases of cloud microphysical properties and their uncertainties as needed for the development of improved cloud parameterizations for models and remote sensing retrievals, and for evaluation of model simulations and retrievals. In particular, we proposed to analyze data collected over the Southern Great Plains (SGP) during the Mid-latitude Continental Convective Clouds Experiment (MC3E), the Storm Peak Laboratory Cloud Property Validation Experiment (STORMVEX), the Small Particles in Cirrus (SPARTICUS) Experiment and the Routine AAF Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations (RACORO) field campaign, over the North Slope of Alaska during the Indirect and Semi-Direct Aerosol Campaign (ISDAC) and the Mixed-Phase Arctic Cloud Experiment (M-PACE), and over the Tropical Western Pacific (TWP) during The Tropical Warm Pool International Cloud Experiment (TWP-ICE), to meet the following 3 objectives; derive statistical databases of single ice particle properties (aspect ratio AR, dominant habit, mass, projected area) and distributions of ice crystals (size distributions SDs, mass-dimension m-D, area-dimension A-D relations, mass-weighted fall speeds, single-scattering properties, total concentrations N, ice mass contents IWC), complete with uncertainty estimates; assess processes by which aerosols modulate cloud properties in arctic stratus and mid-latitude cumuli, and quantify aerosol’s influence in context of varying meteorological and surface conditions; and determine how ice cloud microphysical, single-scattering and fall-out properties and contributions of small ice crystals to such properties vary according to location, environment, surface, meteorological and aerosol conditions, and develop parameterizations of such effects.In this report we describe the accomplishments that we made on all 3 research objectives.

54 ENVIRONMENTAL SCIENCES↗

Topics in Theoretical Physics

This award supported a broadly based research effort in theoretical particle physics, including research aimed at uncovering the laws of nature at short (subatomic) and long (cosmological) distances. These theoretical developments apply to experiments in laboratories such as CERN, the facility that operates the Large Hadron Collider outside Geneva, as well as to cosmological investigations done using telescopes and satellites. The results reported here apply to physics beyond the so-called Standard Model of particle physics; physics of high energy collisions such as those observed at the Large Hadron Collider; theoretical and mathematical tools and frameworks for describing the laws of nature at short distances; cosmology and astrophysics; and analytic and computational methods to solve theories of short distance physics. Some specific research accomplishments include + Theories of the electroweak interactions, the forces that give rise to many forms of radioactive decay; + Physics of the recently discovered Higgs boson. + Models and phenomenology of dark matter, the mysterious component of the universe, that has so far been detected only by its gravitational effects. + High energy particles in astrophysics and cosmology. + Algorithmic research and Computational methods for physics of and beyond the Standard Model. + Theory and applications of relativity and its possible limitations. + Topological effects in field theory and cosmology. + Conformally invariant systems and AdS/CFT. This award also supported significant training of students and postdoctoral fellows to lead the research effort in particle theory for the coming decades. These students and fellows worked closely with other members of the group as well as theoretical and experimental colleagues throughout the physics community. Many of the research projects funded by this grant arose in response to recently obtained experimental results in the areas of particle physics and cosmology. We describe a few of these below. Relativity is founded on a symmetry property of nature called "Lorentz Invariance". Like all symmetry properties, it is essential to determine precisely how symmetric nature actually is; that is, do the laws of nature fully respect the symmetry or is there room for tiny symmetry violating effects? An important consequence of Lorentz invariance is the existence of a universal limiting velocity for all physical particles. Light travels at this limiting velocity so it is frequently referred to as simply "the speed of light", but relativity requires that ALL particles travel more slowly than this speed. Once the Higgs particle was discovered in 2012 a natural question was whether or not this particle's speed was consistent with relativity. Although the speed of the Higgs particle is not measurable directly, Cohen has shown that, if the maximal speed of the Higgs particle was not precisely the same as the speed of light, then the Higgs would have some unusual properties. In some cases the Higgs would be unstable to some unusual decay modes; in other cases the interactions of the Higgs with other particles would change the properties of these other particles in ways that could be observed in so-called cosmic rays, very energetic particles (such as photons, protons and other atomic nuclei) coming from space. Once these particles hit the upper atmosphere they produce a "shower" of particles that can be seen by ground-based instruments. If the Higgs has a maximal speed that differs even a tiny bit from the speed of light these showers would look quite different from what is observed. In this way Cohen was able to establish that the Higgs travels with a maximal speed that cannot differ from the speed of light by more than one part in a thousand-trillion. This is by far the most precisely determined property of the Higgs particle. Cohen and Schmaltz reviewed evidence from the Large Hadron Collider (LHC), a particle physics experiment operating at the CERN laboratory near Geneva, for a new particle sometimes called a W'. This evidence included certain unexpected by-products in collisions of protons at very high energy. While the evidence was not significant enough to claim a discovery, it was sufficiently intriguing that many particle theorists worked to construct explanations for this signal. Cohen and Schmaltz were able to determine that such explanations are highly constrained by previous experiments involving collisions of very energetic particles. Nevertheless they were able to construct a theory that adequately explains the LHC data and remain consistent with prior experiments. Their explanation predicts the existence of yet another new particle, called a Z', with a mass slightly greater than that of the W'. This additional particle, if it exists, should be seen as more data is collected from the LHC. Amusingly, there is one collision by-product that has already been seen by the CMS experiment at the LHC that supports the existence of this new particle; however, it is not unlikely that this single event is a so-called "background" event, that is a somewhat atypical by-product of a conventional Standard Model process. This theory for the anomalous LHC data will either be confirmed or excluded with further data-taking at the LHC. The ratio of the number of electrons produced in bottom quark decays over the number of muons produced has been measured at the LHC. This ratio is interesting because it can be predicted very precisely from a basic property of the Standard Model: lepton universality. If lepton universality is correct, the ratio of electrons to muons is predicted to be equal to 1. The first measurements of this ratio find a value different from 1 with a statistical significance of about 3 standard deviations. Schmaltz and collaborator proposed a new extension of the Standard Model which can explain the new data. In addition, Schmaltz and collaborators proposed several new measurements of ratios of decay rates which can confirm or rule out the surprising results from the earlier LHC data. The most recent and precise measurements of the cosmic microwave background from the Planck satellite, from a combination of measurements of the dark matter distribution in the universe, and from a measurement of the expansion rate of the universe today show some disagreement when interpreted in terms of the so-called LambdaCDM model. Schmaltz and collaborators proposed an alternative model to LambdaCDM in which the usual cold dark matter is replaced by a new ``dark sector". This sector consists of a cold dark matter particle which interacts with a newly postulated dark radiation component of the universe. The dark radiation can help explain the discrepancy in measurements of the expansion rate, and the dark matter interactions subtly modify the clumping of dark matter at large scales, thus potentially explaining both kinds of tensions in the data. In two publications Schmaltz described the new model and then performed a precision comparison of the predictions of the model with all currently available cosmological data. The results favor the new model at the level of three standard deviations with current data. Quantum Field Theory (QFT) is the language we use to describe quantum systems which are consistent with Einstein’s theory of Special Relativity. In particular, the requirement of Einstein’s theory that signals not travel faster than the speed of light constrains the types of interactions which particles can engage in. One consequence of relativity is that these interactions cannot preserve particle number. The stronger the interactions, the more severe the particle number violation in a given Relativistic QFT. When particle number violation is strong, it becomes very difficult to adequately parameterize the quantum wave function (which characterizes the state of a quantum system). For example, though we can formulate the QFT which describes the strong force as a set of interactions between quark and gluon particles, we have no clear idea how to express the proton state in terms of these quarks and gluons. This is because the proton, though a bound state of quarks and gluons, is not a state of a fixed number of particles due to strong interactions. Yet, understanding the proton state is very important in order to theoretically predict the reaction rates observed at the LHC in Geneva, which is a proton-proton collider. Katz has formulated a new approach to QFT, which among other things offers a way to adequately approximate the quantum wave function of a bound state at strong coupling. The approximation scheme is related to the fact that any sensible QFT (including that of the strong interactions) is at short distances approximately self-similar upon rescaling of space and time. It turns out that keeping track of the response upon this rescaling is important in efficiently parameterizing the state. Katz and collaborators have used this observation to approximate the state of the proton in toy versions of the strong force. In the late 60s Sheldon Glashow, Abdus Salam and Steven Weinberg (1979 Nobel Prize awardees) proposed a theory unifying weak and electromagnetic interaction which assumed the existence of new particles, the W and Z bosons. The W and Z bosons were eventually detected in high-energy collision in a particle accelerator at CERN, and the recent discovery of the Higgs meson at the Large Hadron Collider (LHC), always at CERN, completed the picture. However, deep theoretical considerations indicate that the theory by Glashow, Weinberg and Salam, often referred to as "the standard model" cannot be the whole story: the existence of new particles and new interactions at yet higher energies is widely anticipated. The experiments at the LHC are looking for these, while theorists, like Brower, Rebbi and collaborators, are investigating models for these new interactions. Working in a large national collaboration with access to the most powerful DOE computers Brower, Rebbi and colleagues have been using calculational techniques, similar to those successfully employed for many years to investigate the interactions among quarks in nucleons, to study theories that can describe the expected "beyond the standard model" (BSM) interactions. Their results, which include also a model for dark matter, have been published in several refereed papers in prestigious journals. Various ideas in topologically interesting field theories predict hypothetical objects such as fractional charges and Majorana excitations. However, such fascinating objects have not been seen in particle physics. Nevertheless, these objects demonstrate possible phenomena that quantum field theory can support. Pi used condensed matter physics as a laboratory to study possible realizations and observable effects of these objects predicted by quantum field theory. In recent times there has developed considerable interest among condensed matter field theorists in precisely the same geometrical and topological structures, which were first discovered in particle physics field theories. From particle physicists' point of view, this is an interesting development, since condensed matter provides an arena in which one can concretely realize particle physics ideas. Moreover, particle physicists can learn new ideas from condensed matter physics. Higgs phenomenon is precisely an important particle physics realization of condensed matter ideas. In contrast to the small distance characterizing condensed matter systems, field theory also describes large distance physics characterizing cosmology. Pi worked on various geometrical effects in the standard theory of cosmology, viz general relativity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Target Soil Impact Verification: Experimental Testing and Kayenta Constitutive Modeling.

This report details experimental testing and constitutive modeling of sandy soil deformation under quasi - static conditions. This is driven by the need to understand constitutive response of soil to target/component behavior upon impact . An experimental and constitutive modeling program was followed to determine elastic - plastic properties and a compressional failure envelope of dry soil . One hydrostatic, one unconfined compressive stress (UCS), nine axisymmetric compression (ACS) , and one uniaxial strain (US) test were conducted at room temperature . Elastic moduli, assuming isotropy, are determined from unload/reload loops and final unloading for all tests pre - failure and increase monotonically with mean stress. Very little modulus degradation was discernable from elastic results even when exposed to mean stresses above 200 MPa . The failure envelope and initial yield surface were determined from peak stresses and observed onset of plastic yielding from all test results. Soil elasto - plastic behavior is described using the Brannon et al. (2009) Kayenta constitutive model. As a validation exercise, the ACS - parameterized Kayenta model is used to predict response of the soil material under uniaxial strain loading. The resulting parameterized and validated Kayenta model is of high quality and suitable for modeling sandy soil deformation under a range of conditions, including that for impact prediction.

54 ENVIRONMENTAL SCIENCES↗

Measurements of pore-scale flow through apertures

Pore-scale aperture effects on flow in pore networks was studied in the laboratory to provide a parameterization for use in transport models. Four cases were considered: regular and irregular pillar/pore alignment with and without an aperture. The velocity field of each case was measured and simulated, providing quantitatively comparable results. Two aperture effect parameterizations were considered: permeability and transmission. Permeability values varied by an order of magnitude between the cases with and without apertures. However, transmission did not correlate with permeability. Despite having much greater permeability the regular aperture case permitted less transmission than the regular case. Moreover, both irregular cases had greater transmission than the regular cases, a difference not supported by the permeabilities. Overall, these findings suggest that pore-scale aperture effects on flow though a pore-network may not be adequately captured by properties such as permeability for applications that are interested in determining particle transport volume and timing.

58 GEOSCIENCES↗

Southern Ocean Uptake in the MPAS-Ocean Model

The research added to the understanding of Southern Ocean Uptake in the following ways. 1) No single process in the model could be “tuned” to improve Southern Ocean uptake of heat, CO 2 and CFCs. 2) More comprehensive studies revealed that multiple aspects of the model are important; in particular the general circulation control of deep (>200m) stratification (Small et al., 2019) and salinity contributions to the shallow (< 200 m) stratification (DuVivier et al, 2018). These two aspects are not strongly connected. The first can be very much improved by increasing horizontal resolution, primarily because of more transport of high salinity water from the Indian Ocean into the Southern Ocean south of Africa, so we now know that “parameterizing” this salinity source will be one necessary step to improving the Southern Ocean uptake at low resolution. 3) However, the second was found to be more than just a resolution issue, which may explain why many studies have found that higher resolution improves only some aspects of model solutions, but degrades others. Even in cases where the surface wind is too strong, the shallow stratification is sufficient to inhibit vertical mixing even though the salinity stabilization is weak compared to observations. The conclusion is that the process responsible for mixing through the stable salinity barrier above 200m is missing from the model. Therefore, a second necessary step is to find the process then incorporate it into the model. The project focused on ocean surface wave driven Langmuir turbulence as the source of additional mixing by studying Large Eddy Simulations (LES) of Southern Ocean boundary layers as they eroded a realistic salinity barrier. The wind, buoyancy and wave driven boundary layer turbulence was found to be well described by Monin-Obukhov similarity by developing similarity functions of the wave driving, and combining these with the established wind and buoyancy functions (Large et al., 2018a). These functions and a novel analysis of the momentum flux and vertical shear fields from the LES revealed two distinct aspects where the boundary layer turbulence was much more organized than expected, and hence much more amenable to implementation in the model (Large et al., 2018b). First, there is a robust non-local across-shear transport of momentum throughout the ocean boundary layer; a hitherto missing process. Furthermore, this transport is related to the across-shear component of the wind stress in the same way that non-local buoyancy transport is related to the surface buoyancy flux, both in the ocean and atmosphere. Even the coefficients are similar in the convective limit. Second, the turbulent fluxes of momentum and buoyancy are distributed through the ocean boundary layer of the LES in a very orderly fashion given by a single shape function that appears to be applicable to both momentum and buoyancy, with or without surface wave driving. Thus, a very technically effective and economic path to incorporating surface wave driving in the vertical mixing scheme of the model is now clear. The high degree of order found for ocean boundary layer turbulence, including Langmuir turbulence, is more amenable to empirical models than to statistical turbulence models, and the computational cost is very much less. Fortuitously, the empirical K-Profile Parameterization is already implemented in the model and the wave-driving similarity function, the non-local momentum transport and the vertical shape function are straightforward enhancements.

58 GEOSCIENCES↗

Advanced Tool for Computer Aided Battery Engineering (SNL)

This project aims to provide experimental support, including support for validation and parameterization, for mechanical failure modeling work as part of DOE's Computer Aided Engineering for Electric-Drive Vehicle Batteries (CAEBAT) program. This work involves mechanical deformation testing on both charged and discharged cells, including abusive mechanical testing leading to battery failure. The mechanical data generated is provided to the modeling teams to provide empirical parameterization as well as validation for newly developed models. Testing fully charged cells and packs (abusive battery testing) also allows for a better understanding of what conditions are most likely to lead to a potentially hazardous thermal runaway event.

25 ENERGY STORAGE↗

Colloid-Facilitated Transport: Studies Related to CFM Project at GTS

Presentation Summary: Insights have been gained as to how to obtain defensible answers to predict radionuclide transport in fractured granites. However, site specific studies still need to be performed to to gain confidence in the prediction; The CFT ladder should be applied to evaluate the potential for enhanced transport with colloids, but most indications are that only very small fractions of strongly-sorbing radionuclides will be capable of CFT over repository time and distance scales; CFT requires very slow desorption from colloids AND very slow filtration of the RN-bearing colloids (relative to time scales of interest); Interrogating such slow processes is a challenge, especially if they are associated with a very small fraction of colloids or very small fraction of sorption sites on colloids (or both); Intuitively, one might expect that stable colloids generated from wasteform degradation that have radionuclides incorporated into their structure (as opposed to asorption association) might pose the biggest risk; Experiments have informed generic modeling approach, including GDSA; Experiments have provided insights into how experimental designs can be tailored and improved to address site-specific and scenario-specific issues; Different host rocks, EBS vs. natural system, and DPC concept can all, in principle, be addressed via different parameterizations of the generic model, with the understanding that parameterizations must be developed through site- and scenario-specific experimental testing

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Major improvements on the longwave radiative interactions between surface and clouds in the polar regions in atmospheric global circulation model

Motivations: One of the most important tasks in the study of physical and socioeconomic aspects of climate change is the simulation of polar climate and how it changes in response to global warming. One import component of the polar energy budget (and thus of polar climate) is the radiative interaction between surface and clouds. As of today, the vast majority of state-of-the-art global climate models (GCMs), including all three US flagship GCMs (NCAR CESM, GFDL CM3, and GISS E2-H/E2-R), still assume constant surface emissivity overall spectral bands in their longwave radiation treatment. Similarly, a majority of state-of-the-art GCMs assume non-scattering clouds in the longwave spectrum, including both the NCAR CESM and GFDL CM3. The issues of ignoring spectral variation in surface emissivity and scattering by clouds manifest themselves in the far-infrared (IR) over the polar continents because (1) there is so little water vapor in such areas that surface far-IR emission can reach clouds; and (2) the imaginary part of the refractive index of ice has a local minimum over 350-550 cm-1, which implies possible strong scattering effects over this spectral region both within and between the snow surface and ice clouds. Our sensitivity studies using CloudSat-retrieved hydrometeor profiles over the high-elevation Antarctic continent show that, in winter, an appropriate inclusion of snow surface spectral emissivity and ice cloud scattering in radiative transfer calculations can noticeably reduce the monthly-mean surface net downward far-IR flux and net atmospheric far-IR emission over the entire region. For the far-IR alone, the magnitude of these effects is ~1Wm-2 or even larger. Projection Objectives: We propose to carry out the following studies: (1) Develop or update schemes that can be used in the CESM to compute the optical properties of the snow surface and ice clouds consistently across the longwave and the shortwave spectrum; (2) Modify the longwave radiation scheme in the CESM such that it can incorporate surface spectral emissivity and scattering of clouds; (3) Validate schemes developed for (1) and (2) against measurements or benchmark models; (4) Evaluate the effects of these new treatments on the simulated polar energy budget and climate by the CESM. Technical Approaches:1.Use a detailed snowpack radiative transfer model to simulate spectrally resolved emissivity of snow surfaces as a function of snow effective radius. The model will be validated against available high spectral-resolution measurements of snow spectral emissivity (ASTER spectra library v2.0). Currently the CESM prognoses vertically-resolved snow effective radius, and snow shortwave albedo is parameterized based on snow effective radius. Our treatment of snow spectral emissivity here will be consistent with the treatment of snow albedo in the shortwave and seamlessly connected with the snow physical module in the CESM. 2. Update the parameterization of ice cloud optical properties in the CESM, for both longwave and shortwave, based on the latest ice-cloud scattering database developed by co-I Prof. Yang. 3. Modify RRTMG_LW, the current longwave radiation scheme in the CESM, to include the surface spectral emissivity and the scattering of clouds in the longwave calculation. Validate modified scheme against a line-by-line radiative transfer model with multiple scattering capability (LBLRTM/CHARTS). 4. Implement products from Steps 1-3 into the CESM and conduct simulations to evaluate their effects on the simulated polar climate and energy budget. These runs will be compared to control runs and evaluated using routine model-observation diagnostics, including data from the ARM NSA (North Slope of Alaska) site. Moreover, a novel aspect in this step is that we will evaluate the simulated longwave band-by-band fluxes and cloud radiative effect against observational counterparts over both polar continents, using a product recently derived by the lead PI Prof. Huang from collocated AIRS and CERES satellite observations from 2002 to 2013. Benefits and outcomes: This project will significantly improve the longwave radiation scheme in the CESM and its successors and provide more realistic treatment of radiative interactions between polar surfaces and clouds. In the future, it will also be straightforward to extend this framework to include spectral emissivity of all surface types. This is responsive to identified DoE Earth System Modeling Science Drivers to improve the performance of the fully coupled climate system modeling and to help address questions on how rapid changes in cryospheric systems interact with the climate system.

54 ENVIRONMENTAL SCIENCES↗