Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ADVECTION”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Thermal stability of a thin disk with magnetically driven winds

The absence of thermal instability in the high/soft state of black hole X-ray binaries, in disagreement with the standard thin disk theory, has been a long-standing riddle for theoretical astronomers. We have tried to resolve this question by studying the thermal stability of a thin disk with magnetically driven winds in the M-dot −Σ plane. It is found that disk winds can greatly decrease the disk temperature and thus help the disk become more stable at a given accretion rate. The critical accretion rate, M-dot {sub crit}, corresponding to the thermal instability threshold, is significantly increased in the presence of disk winds. For α = 0.01 and B {sub φ} = 10B {sub p}, the disk is quite stable even for a very weak initial poloidal magnetic field [β{sub p,0}∼2000,β{sub p}=(P{sub gas}+P{sub rad})/(B{sub p}{sup 2}/8π)]. However, when B {sub φ} = B {sub p} or B {sub φ} = 0.1B {sub p}, a somewhat stronger (but still weak) field (β{sub p,} {sub 0} ∼ 200 or β{sub p,} {sub 0} ∼ 20) is required to make the disk stable. Nevertheless, despite the great increase of M-dot {sub crit}, the luminosity threshold, corresponding to instability, remains almost constant or decreases slowly with increasing M-dot {sub crit} due to decreased gas temperature. The advection and diffusion timescales of the large-scale magnetic field threading the disk are also investigated in this work. We find that the advection timescale can be smaller than the diffusion timescale in a disk with winds, because the disk winds take away most of the gravitational energy released in the disk, resulting in the decrease of the magnetic diffusivity η and the increase of the diffusion timescale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE DISCOVERY OF A NEW INSTABILITY IN A HYPERACCRETION FLOW AND ITS IMPLICATION FOR GAMMA-RAY BURSTS

A hyperaccretion flow around a stellar mass black hole is thought to be the most plausible engine that powers gamma-ray bursts (GRBs). The flow efficiently cools via neutrino emission at ∼> 0.003-0.01 M {sub ☉} s{sup –1} (corresponding to a luminosity of ∼10{sup 50} erg s{sup –1}), while neither neutrino nor photon emission is efficient below this rate, so the flow should be advection-dominated. We carefully solve how a transition occurs from the advection-dominated to the neutrino-dominated branches, and find that the slope of the thermal equilibrium curve is negative in the surface density-accretion rate (Σ- M-dot ) plane, a condition for viscous instability, at radii smaller than ∼12 R {sub g} (with R {sub g} being the gravitational radius). We also confirm that the flow is thermally stable. The consequence of this instability is the formation of a clumpy structure in the flow. This is because the larger (respectively smaller) surface density is, the smaller (respectively larger) the mass accretion rate from the region in question becomes, leading to growth of the density contrast. The timescale for clump formation is estimated to be shorter than 0.1 s. The observational implication is discussed in the context of GRBs. We suggest that this might explain the origin of the large variability observed in the prompt emission of GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Broadening of Modeled Cloud Droplet Spectra Using Bin Microphysics in an Eulerian Spatial Domain

This study investigates droplet size distribution (DSD) characteristics from condensational growth and transport in Eulerian dynamical models with bin microphysics. A hierarchy of modeling frameworks is utilized, including parcel, one-dimensional (1D), and three-dimensional large-eddy simulation (LES). The bin DSDs from the 1D model, which includes only vertical advection and condensational growth, are nearly as broad as those from the LES and in line with observed DSD widths for stratocumulus clouds. These DSDs are much broader than those from Lagrangian microphysical calculations within a parcel framework that serve as a numerical benchmark for the 1D tests. In contrast, the bin-modeled DSDs are similar to the Lagrangian microphysical benchmark for a rising parcel in which Eulerian transport is not considered. These results indicate that numerical diffusion associated with vertical advection is a key contributor to broadening DSDs in the 1D model and LES. This DSD broadening from vertical numerical diffusion is unphysical, in contrast to the physical mixing processes that previous studies have indicated broaden DSDs in real clouds. It is proposed that artificial DSD broadening from vertical numerical diffusion compensates for underrepresented horizontal variability and mixing of different droplet populations in typical LES configurations with bin microphysics, or the neglect of other mechanisms that broaden DSDs such as growth of giant cloud condensation nuclei. These results call into question the ability of Eulerian dynamical models with bin microphysics to investigate the physical mechanisms for DSD broadening, even though they may reasonably simulate overall DSD characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Observing Profiles of Derived Kinematic Field Quantities Using a Network of Profiling Sites

Abstract Observations of thermodynamic and kinematic parameters associated with derivatives of the thermodynamics and wind fields, namely, advection, vorticity, divergence, and deformation, can be obtained by applying Green’s theorem to a network of observing sites. The five nodes that comprise the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) profiling network, spaced 50–80 km apart, are used to obtain measurements of these parameters over a finite region. To demonstrate the applicability of this technique at this location, it is first applied to gridded model output from the High-Resolution Rapid Refresh (HRRR) numerical weather prediction model, using profiles from the locations of ARM network sites, so that values calculated from this method can be directly compared to finite difference calculations. Good agreement is found between both approaches as well as between the model and values calculated from the observations. Uncertainties for the observations are obtained via a Monte Carlo process in which the profiles are randomly perturbed in accordance with their known error characteristics. The existing size of the ARM network is well suited to capturing these parameters, with strong correlations to model values and smaller uncertainties than a more closely spaced network, yet it is small enough that it avoids the tendency for advection to go to zero over a large area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Sorption Modeling and Verification for Off-Gas Treatment

The project has made progress toward developing a comprehensive modeling capability for the capture of target species in off gas evolved during the reprocessing of nuclear fuel. The effort has integrated experimentation, model development, and computer code development for adsorption and absorption processes. For adsorption, a modeling library has been initiated to include (a) equilibrium models for uptake of off-gas components by adsorbents, (b) mass transfer models to describe mass transfer to a particle, diffusion through the pores of the particle and adsorption on the active sites of the particle, and (c) interconnection of these models to fixed bed adsorption modeling which includes advection through the bed. For single-component equilibria, a Generalized Statistical Thermodynamic Adsorption (GSTA) code was developed to represent experimental data from a broad range of isotherm types; this is equivalent to a Langmuir isotherm in the two-parameter case, and was demonstrated for Kr on INL-engineered sorbent HZ PAN, water sorption on molecular sieve A sorbent material (MS3A), and Kr and Xe capture on metal-organic framework (MOF) materials. The GSTA isotherm was extended to multicomponent systems through application of a modified spreading pressure surface activity model and generalized predictive adsorbed solution theory; the result is the capability to estimate multicomponent adsorption equilibria from single-component isotherms. This advance, which enhances the capability to simulate systems related to off-gas treatment, has been demonstrated for a range of real-gas systems in the literature and is ready for testing with data currently being collected for multicomponent systems of interest, including iodine and water on MS3A. A diffusion kinetic model for sorbent pellets involving pore and surface diffusion as well as external mass transfer has been established, and a methodology was developed for determining unknown diffusivity parameters from transient uptake data. Two parallel approaches have been explored for integrating the kernels described above into a mass-transport model for adsorption in fixed beds. In one, the GSTA isotherm kernel has been incorporated into the MOOSE framework; in the other approach, a focused finite-difference framework and PDE kernels have been developed. Issues, including oscillatory behavior in MOOSE solutions to advection-diffusion problems, and opportunities have been identified for each approach, and a path forward has been identified toward developing a stronger modeling platform. Experimental systems were established for collection of microscopic kinetics and equilibria data for single and multicomponent uptake of gaseous species on solid sorbents. The systems, which can operate at ambient temperature to 250°C and dew points from -69 to 17°C, are useful for collecting data needed for modeling performance of sorbents of interest. Experiments were conducted to determine applicable models and parameters for isotherms and mass transfer for water and/or iodine adsorption on MS3A. Validation experiments were also conducted for water adsorption on fixed beds of MS3A. For absorption, work involved modeling with supportive experimentation. A dynamic model was developed to simulate CO 2 absorption with chemical reaction using high alkaline content water solutions. A computer code was developed to implement the model based upon transient mass and energy balances. Experiments were conducted in a laboratory-scale column to determine model parameters. The influence of geometric parameters and operating variables on CO 2 absorption was studied over a wide range of conditions. This project has resulted in 7 publications, with 3 manuscripts in preparation. Also, 15 presentations were given at national meetings of ANS and AIChE and at Material Recovery and Waste Forms Campaign Working Group meetings.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

TGLF Recalibration for ITER Standard Case Parameters FY2015: Theory and Simulation Performance Target Final Report

This work was motivated by the observation, as early as 2008, that GYRO simulations of some ITER operating scenarios exhibited nonlinear zonal-flow generation large enough to effectively quench turbulence inside r /a ~ 0.5. This observation of flow-dominated, low-transport states persisted even as more accurate and comprehensive predictions of ITER profiles were made using the state-of-the-art TGLF transport model. This core stabilization is in stark contrast to GYRO-TGLF comparisons for modern-day tokamaks, for which GYRO and TGLF are typically in very close agreement. So, we began to suspect that TGLF needed to be generalized to include the effect of zonal-flow stabilization in order to be more accurate for the conditions of reactor simulations. While the precise cause of the GYRO-TGLF discrepancy for ITER parameters was not known, it was speculated that closeness to threshold in the absence of driven rotation, as well as electromagnetic stabilization, created conditions more sensitive the self-generated zonal-flow stabilization than in modern tokamaks. Need for nonlinear zonal-flow stabilization: To explore the inclusion of a zonal-flow stabilization mechanism in TGLF, we started with a nominal ITER profile predicted by TGLF, and then performed linear and nonlinear GYRO simulations to characterize the behavior at and slightly above the nominal temperature gradients for finite levels of energy transport. Then, we ran TGLF on these cases to see where the discrepancies were largest. The predicted ITER profiles were indeed near to the TGLF threshold over most of the plasma core in the hybrid discharge studied (weak magnetic shear, q > 1). Scanning temperature gradients above the TGLF power balance values also showed that TGLF overpredicted the electron energy transport in the low-collisionality ITER plasma. At first (in Q3), a model of only the zonal-flow stabilization (Dimits shift) was attempted. Although we were able to construct an ad hoc model of the zonal flows that fit the GYRO simulations, the parameters of the model had to be tuned to each case. A physics basis for the zonal flow model was lacking. Electron energy transport at short wavelength: A secondary issue – the high-k electron energy flux – was initially assumed to be independent of the zonal flow effect. However, detailed studies of the fluctuation spectra from recent multiscale (electron and ion scale) GYRO simulations provided a critical new insight into the role of zonal flows. The multiscale simulations suggested that advection by the zonal flows strongly suppressed electron-scale turbulence. Radial shear of the zonal E×B fluctuation could not compete with the large electron-scale linear growth rate, but the k x -mixing rate of the E×B advection could. This insight led to a preliminary new model for the way zonal flows saturate both electron- and ion-scale turbulence. It was also discovered that the strength of the zonal E×B velocity could be computed from the linear growth rate spectrum. The new saturation model (SAT1), which replaces the original model (SAT0), was fit to the multiscale GYRO simulations as well as the ion-scale GYRO simulations used to calibrate the original SAT0 model. Thus, SAT1 captures the physics of both multiscale electron transport and zonal-flow stabilization. In future work, the SAT1 model will require significant further testing and (expensive) calibration with nonlinear multiscale gyrokinetic simulations over a wider variety of plasma conditions – certainly more than the small set of scans about a single C-Mod L-mode discharge. We believe the SAT1 model holds great promise as a physics-based model of the multiscale turbulent transport in fusion devices. Correction to ITER performance predictions: Finally, the impact of the SAT1model on the ITER hybrid case is mixed. Without the electron-scale contribution to the fluxes, the Dimits shift makes a significant improvement in the predicted fusion power as originally posited. Alas, including the high-k electron transport reduces the improvement, yielding a modest net increase in predicted fusion power compared to the TGLF prediction with the original SAT0 model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Transient Kinetic Experiments within the High Conversion Domain: The Case of Ammonia Decomposition

In the development of catalytic materials, a set of standard conditions is needed where the kinetic performance of many samples can be compared. This can be challenging when a sample set covers a broad range of activity. Precise kinetic characterization requires uniformity in the gas and catalyst bed composition. This limits the range of convecting devices to low conversion (generally <20%). While steady-state kinetics offer a snapshot of conversion, yield and apparent rates of the slow reaction steps, transient techniques offer much greater detail of rate processes and hence more information as to why certain catalyst compositions offer better performance. In this work, transient experiments in two transport regimes are compared: an advecting differential plug flow reactor (PFR) and a pure-diffusion temporal analysis of products (TAP) reactor. The decomposition of ammonia was used as a model reaction to test three simple materials: polycrystalline iron, cobalt and a bimetallic preparation of the two. These materials presented a wide range of activity and it was not possible to capture transient information in the advecting device for all samples at the same conditions while ensuring uniformity. We push the boundary for the theoretical estimates of uniformity in the TAP device and find reliable kinetic measurement up to 90% conversion. However, what is more advantageous from this technique is the ability to observe the time-dependence of the reaction rate rather than just singular points of conversion and yield. For example, on the iron sample we observed reversible adsorption of ammonia and on cobalt materials we identify two routes for hydrogen production. From the time-dependence of reactants and product, the dynamic accumulation was calculated. This was used to understand the atomic distribution of H and N species regulated by the surface of different materials. When ammonia was pulsed at 550 °C, the surface hydrogen/nitrogen, (H/N), ratios that evolved for Fe, CoFe and Co were 2.4, 0.25 and 0.3 respectively. This indicates that iron will store a mixture of hydrogenated species while materials with cobalt will predominantly store NH and N. While much is already known about iron, cobalt and ammonia decomposition, the goal of this work was to demonstrate new tools for comparing materials over a wider window of conversion and with much greater kinetic detail. As a result, this provides an approach for detailed kinetic discrimination of more complex industrial samples beyond conversion and yield.

36 MATERIALS SCIENCE↗
DOE OSTI2016

RELATIVE CONTRIBUTION OF THE MAGNETIC FIELD BARRIER AND SOLAR WIND SPEED IN ICME-ASSOCIATED FORBUSH DECREASES

We study 50 cosmic-ray Forbush decreases (FDs) from the Oulu neutron monitor data during 1997–2005 that were associated with Earth-directed interplanetary coronal mass ejections (ICMEs). Such events are generally thought to arise due to the shielding of cosmic rays by a propagating diffusive barrier. The main processes at work are the diffusion of cosmic rays across the large-scale magnetic fields carried by the ICME and their advection by the solar wind. In an attempt to better understand the relative importance of these effects, we analyze the relationship between the FD profiles and those of the interplanetary magnetic field (B) and the solar wind speed (V {sub sw}). Over the entire duration of a given FD, we find that the FD profile is generally (anti)correlated with the B and V {sub sw} profiles. This trend holds separately for the FD main and recovery phases too. For the recovery phases, however, the FD profile is highly anti-correlated with the V {sub sw} profile, but not with the B profile. While the total duration of the FD profile is similar to that of the V {sub sw} profile, it is significantly longer than that of the B profile. Using the convection–diffusion model, a significant contribution of advection by solar wind is found during the recovery phases of the FD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A THEORETICAL STUDY OF THE BUILD-UP OF THE SUN’S POLAR MAGNETIC FIELD BY USING A 3D KINEMATIC DYNAMO MODEL

We develop a three-dimensional kinematic self-sustaining model of the solar dynamo in which the poloidal field generation is from tilted bipolar sunspot pairs placed on the solar surface above regions of strong toroidal field by using the SpotMaker algorithm, and then the transport of this poloidal field to the tachocline is primarily caused by turbulent diffusion. We obtain a dipolar solution within a certain range of parameters. We use this model to study the build-up of the polar magnetic field and show that some insights obtained from surface flux transport models have to be revised. We present results obtained by putting a single bipolar sunspot pair in a hemisphere and two symmetrical sunspot pairs in two hemispheres. We find that the polar fields produced by them disappear due to the upward advection of poloidal flux at low latitudes, which emerges as oppositely signed radial flux and which is then advected poleward by the meridional flow. We also study the effect that a large sunspot pair, violating Hale’s polarity law, would have on the polar field. We find that there would be some effect—especially if the anti-Hale pair appears at high latitudes in the mid-phase of the cycle—though the effect is not very dramatic.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

Effects of Mesoscale Eddies in the Active Mixed Layer: Test of the Parametrisation in Eddy Resolving Simulations

In eddy resolving simulations, we test a mixed layer mesoscale parametrisation, developed recently by Canuto and Dubovikov [Ocean Model., 2011, 39, 200-207]. With no adjustable parameters, the parametrisation yields the horizontal and vertical mesoscale fluxes in terms of coarse-resolution fields and eddy kinetic energy (EKE). We compare terms of the parametrisation diagnosed from coarse-grained fields with the eddy mesoscale fluxes diagnosed directly from the high resolution model. An expression for the EKE in terms of mean fields has also been found to get a closed parametrisation in terms of the mean fields only. In 40 numerical experiments we simulated two types of flows: idealised flows driven by baroclinic instabilities only, and more realistic flows, driven by wind and surface fluxes as well as by inflow-outflow. The diagnosed quasi-instantaneous horizontal and vertical mesoscale buoyancy fluxes (averaged over 1-2 degrees and 10 days) demonstrate a strong scatter typical for turbulent flows, however, the fluxes are positively correlated with the parametrisation with higher (0.5-0.74) correlations at the experiments with larger baroclinic radius Rossby. After being averaged over 3-4 months, diffusivities diagnosed from the eddy resolving simulations are consistent with the parametrisation for a broad range of parameters. Diagnosed vertical mesoscale fluxes restratify mixed layer and are in a good agreement with the parametrisation unless vertical turbulent mixing in the upper layer becomes strong enough in comparison with mesoscale advection. In the latter case, numerical simulations demonstrate that the deviation of the fluxes from the parametrisation is controlled by dimensionless parameter estimating the ratio of vertical turbulent mixing term to mesoscale advection. An analysis using a modified omega-equation reveals that the effects of the vertical mixing of vorticity is responsible for the two-three fold amplification of vertical mesoscale flux. Possible physical mechanisms, responsible for the amplification of vertical mesoscale flux are discussed.

simulation↗
DOE OSTI2020

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Consequences of systematic model drift in DYNAMO MJO hindcasts with SP-CAM and CAM5

Hindcast simulations of MJO events during the dynamics of the MJO (DYNAMO) field campaign are conducted with two models, one with conventional parameterization (CAM5) and a comparable model that utilizes superparameterization (SP–CAM). SP–CAM is shown to produce a qualitatively better reproduction of the fluctuations of precipitation and low–level zonal wind associated with the first two DYNAMO MJO events compared to CAM5. Interestingly, skill metrics using the real–time multivariate MJO index (RMM) suggest the opposite conclusion that CAM5 has more skill than SP–CAM. This inconsistency can be explained by a systematic increase of RMM amplitude with lead time, which results from a drift of the large–scale wind field in SP–CAM that projects strongly onto the RMM index. CAM5 hindcasts exhibit a contraction of the moisture distribution, in which extreme wet and dry conditions become less frequent with lead time. SP–CAM hindcasts better reproduce the observed moisture distribution, but also have stronger drift patterns of moisture budget terms, such as an increase in drying by meridional advection in SP–CAM. This advection tendency in SP–CAM appears to be associated with enhanced off–equatorial synoptic eddy activity with lead time. In conclusion, systematic drift moisture tendencies in SP–CAM are of similar magnitude to intraseasonal moisture tendencies, and therefore are important for understanding MJO prediction skill.

42 ENGINEERING↗
DOE OSTI2014

Functional approach to exploring climatic and landscape controls on runoff generation: 2 Timing of runoff storm response

Abstract Hortonian overland flow, Dunne overland flow, and subsurface stormflow are the three most dominant mechanisms contributing to both the volume and timing of streamflow in headwater catchments. In this paper, guided by the Dunne diagram, we explore the impacts of climate, soil, and topography on estimated probability distributions of the travel times of each of these three runoff components. In each case, these are expressed in terms of the Connected Instantaneous Response Functions (CIRF) and account for the dynamics of their individual partial effective contributing areas that retain the connectivity to the outlet (instead of the whole catchment area). A spatially distributed hydrological model is used to derive the CIRFs numerically under multiple combinations of climate, soil, and topographic properties. The mean travel times and dimensionless forms of the CIRFs (i.e., scaled by their respective mean travel times) are used to examine both advective and dispersive aspects of catchment's runoff routing response. It is found that the CIRFs, upon nondimensionalization, collapsed to common characteristic shapes, which could be explained in terms of the relative contributions of hillslope and channel network flows, and the size of runoff contributing areas. The contributing areas, particularly for the Dunne overland flow, are themselves found to be governed by the competition between drainage of and recharge to the water table, and could be explained by a dimensionless drainage index which quantifies this competition. The study also reveals simple indicators based on landscape properties that can explain the magnitude of travel times in different catchments.

Li, Hong‐Yi↗
DOE OSTI2016

On the variability of the Priestley‐Taylor coefficient over water bodies

Abstract Deviations in the Priestley‐Taylor (PT) coefficient α PT from its accepted 1.26 value are analyzed over large lakes, reservoirs, and wetlands where stomatal or soil controls are minimal or absent. The data sets feature wide variations in water body sizes and climatic conditions. Neither surface temperature nor sensible heat flux variations alone, which proved successful in characterizing α PT variations over some crops, explain measured deviations in α PT over water. It is shown that the relative transport efficiency of turbulent heat and water vapor is key to explaining variations in α PT over water surfaces, thereby offering a new perspective over the concept of minimal advection or entrainment introduced by PT. Methods that allow the determination of α PT based on low‐frequency sampling (i.e., 0.1 Hz) are then developed and tested, which are usable with standard meteorological sensors that filter some but not all turbulent fluctuations. Using approximations to the Gram determinant inequality, the relative transport efficiency is derived as a function of the correlation coefficient between temperature and water vapor concentration fluctuations ( R Tq ). The proposed approach reasonably explains the measured deviations from the conventional α PT = 1.26 value even when R Tq is determined from air temperature and water vapor concentration time series that are Gaussian‐filtered and subsampled to a cutoff frequency of 0.1 Hz. Because over water bodies, R Tq deviations from unity are often associated with advection and/or entrainment, linkages between α PT and R Tq offer both a diagnostic approach to assess their significance and a prognostic approach to correct the 1.26 value when using routine meteorological measurements of temperature and humidity.

Assouline, Shmuel↗
DOE OSTI2017

Entrainment rate diurnal cycle in marine stratiform clouds estimated from geostationary satellite retrievals and a meteorological forecast model: Entrainment Rate in Stratiform Clouds

The mean diurnal cycle of cloud entrainment rate (w e ) over the northeast Pacific region is for the first time computed by combining, in a mixed-layer model framework, the hourly composited GOES-15 satellite-based cloud top height (H T ) tendency, advection, and large-scale vertical velocity (w) during May to September 2013, with horizontal winds and w taken from the European Centre for Medium-Range Weather Forecasts (ECMWF) model. The tendency term dominates the magnitude and phase of the w e diurnal cycle, with a secondary role of w, and a modest advective contribution. The peak and minimum in w e occur between 20:00–22:00 LT and 9:00–11:00 LT, respectively, in close agreement with the diurnal cycle of turbulence driven by cloud top longwave cooling. Uncertainties in H T and ECMWF fields are assessed with in situ observations and three meteorological reanalysis data sets. This study provides the basis for constructing nearly global climatologies of w e by combining a suite of well-calibrated geostationary satellites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Excess Surface Area in Bioelectrochemical Systems Causes ion Transport Limitations

We investigated ion transport limitations on 3D graphite felt electrodes by growing Geobacter sulfurreducens biofilms with advection to eliminate external mass transfer limitations. We characterized ion transport limitations by: (i) showing that serially increasing NaCl concentration up to 200mM increased current linearly up to a total of þ273% vs. 0mM NaCl under advective conditions; (ii) growing the biofilm with a starting concentration of 200mM NaCl, which led to a maximum current increase of 400% vs. current generation without NaCl, and (iii) showing that un-colonized surface area remained even after steadystate current was reached. After accounting for iR effects, we confirmed that the excess surface area existed despite a non-zero overpotential. The fact that the biofilm was constrained from colonizing and producing further current under these conditions confirmed the biofilms under study here were ion transport-limited. Our work demonstrates that the use of high surface area electrodes may not increase current density when the system design allows ion transport limitations to become dominant.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2018

A clarification and extension of our model of regolith formation on hillslopes: A clarification of our model of regolith formation on hillslopes

In this paper we clarify that our weathering model from 2013 did not explicitly describe weathering of soil moving downhill along hillslopes. Additionally, we re-analyze the role of the term that we neglected that describes loss of regolith mass through mineral dissolution. We derive an equation for this term by including lateral flow of water inside the model hill. For the revised hill model, we define a dimensionless parameter that allows estimation of the effect of lateral flow on the steady-state hillslope. This parameter is equal to the ratio of averaged advective flux of dissolved species out of the hill to the rate of total denudation. The parameter also yields a criterion for the existence of a steady state regolith thickness for systems experiencing unidirectional advection at a constant velocity: for a ridge, the rate of downward flow of water (qy) must be less than the rate of upward movement of rock (E) after normalization by a small parameter, α. This parameter is equal to the equilibrium aqueous concentration divided by the concentration of the reacting mineral in the rock. Alternately, a steady state may exist for the case of both vertical and lateral flow in a hill for any value of erosion rate if the Darcy velocities decrease with depth. Subsurface flow systems play an essential role in the existence of both steady-state hillslopes and steady-state regolith thicknesses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗