Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generator parameter estimation”

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

Estimation of the time for steam generator trip due to cyber intrusions

The time required to trip a pressurized water reactor (PWR) by inserting malicious signals into its steam generator (SG) control system has been studied using the Generic PWR (GPWR) Simulator. A semi-analytical model is developed to approximately reproduce the simulator response and understand the dynamics of the control unit. A series of two proportional-integral controllers determines control action according to preset constants, the readings from the feedwater level sensor, and those from feedwater and steam flowrate transmitters. It is observed that the most important factor that determines whether a trip will occur is how much additional water is added to or withheld from the SG over time compared to normal operating conditions. In order to determine the effects of control action on the SG, changes in mass inventory are considered. This approach models the SG water level as a function of mass inventory and has a backward temporal memory. A Python interface is developed for the GPWR framework to automatically simulate different spoofing scenarios and post-process the related data. We observe that the trip times predominantly depend on flow mismatch and/or level errors. Controller parameters, including the integral time and gain constants, either speed up or slow down the rate of progression to a trip setpoint but do not cause a trip by themselves. The reactor can trip on a high-level signal when the reading crosses above 78%, increased from its reference level of 57%, or a low-level reading when it is below 25%. The present results show roughly how long the operators would have to respond to an attack, given a specific set of spoofing signals within the issue space analyzed. Furthermore, we have generated a simple surface by fitting a combination of exponential functions to the data obtained from the GPWR Simulator. In general, trips on a low level have been observed to occur faster than those on a high level.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Life cycle energy and environmental impacts of hydroprocessed renewable jet fuel production from pennycress

There is a growing interest in making the overall aviation industry more environmentally-friendly by developing advanced biofuels. In this context, the main driver is establishing feedstocks with lower environmental loads and higher energy use efficiency. The objective of this study was to evaluate the energy requirement and environmental impacts of hydroprocessed renewable jet fuel (HRJ) production from pennycress. The system boundary included pennycress production, logistics, oil extraction and conversion to HRJ, wastewater treatment, and electricity and steam generation. Energy allocation method was used for the analysis. The functional unit was considered to be one GJ of HRJ, and the biorefinery size was considered to be 18.9 million l/yr (5 million gal/yr) HRJ. Energy use efficiency of pennycress-based HRJ production was estimated to be 3.9–5.5 (90% central range), which was higher than that of HRJ from other oilseeds, such as canola. Global warming potential of pennycress-based HRJ (35–49 kgCO 2 eq/GJ of HRJ (90% central range)) was also lower than those of HRJ from similar oilseeds and the commercial jet fuel baseline. Pennycress production was the main contributor to the total energy use and environmental impacts of HRJ production, mainly due to the use of nitrogen fertilizer for production, and fuel combustion for farm operations. Pennycress yield was identified as the most sensitive parameter affecting the energy use efficiency and environmental impacts of HRJ production. Furthermore, the outcomes of this study are useful for identifying the performance targets towards improving energy use efficiency and environmental performance of pennycress-based HRJ production at the commercial scale.

09 BIOMASS FUELS↗

Estimates of the wavenumber wavelet power spectrum of magnetic fluctuations during magnetic reconnection

Fluctuation analyses of experimental observations generally lack high temporal resolution and are in frequency-space f, contrary to theoretical efforts in wavenumber-space k. This is due to the inherent limits of the Fourier transform, though it is prominent due to the ease of diagnostic implementation. Advances in wavelet-based analysis have provided relief due to its temporal resolution, but in its common use, is still hard to compare to theoretical models. By using the two-point correlation technique in conjunction with large data sets, a wavelet power spectrum in wavenumber-space can be created. Dubbed the wavenumber wavelet power spectrum, this spectrum relates wavenumber to power in time. Further, this analysis technique more closely connects characterizations of experimentally observed fluctuations with other system parameters and theoretical predictions. In this article, we develop the wavenumber wavelet power spectrum using magnetic fluctuations caused by tearing instability driven magnetic reconnection in reproducible, high temperature laboratory plasmas. These dynamic magnetic fluctuations generated in reversed field pinch plasmas are broadband, ranging from the low frequency, 10's of kHz, up to the ion gyroradii frequencies, 100's of kHz. The dominant fluctuations have poloidal and toroidal mode numbers (m,n)=(1,6−10) and can grow to 2%–3% of the mean magnetic field. During these reconnection events, ions, and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases, all on a semi-periodic basis. The newly developed spectrum provides better temporal resolution of spectrum characteristics to correlate with these particle energization phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗

A wall interference assessment/correction system

A Wall Signature method, the Hackett method, has been selected to be adapted for the 12-ft Wind Tunnel wall interference assessment/correction (WIAC) system in the present phase. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in turn for estimating wall interferences at the model location. The Wall Signature method will be formulated for application to the unique geometry of the 12-ft Tunnel. The development and implementation of a working prototype will be completed, delivered and documented with a software manual. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free-air and confined wind-tunnel flow fields for each of the test articles over a range of test configurations. Specifically, the pressure signature at the test section wall will be computed for the tunnel case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method-Wall Signature method. The performance of the WIAC method then may be evaluated by comparing the corrected parameters with those for the free-air simulation. Each set of wind tunnel/test article numerical simulations provides data to validate the WIAC method. A numerical wind tunnel test simulation is initiated to validate the WIAC methods developed in the project. In the present reported period, the blockage correction has been developed and implemented for a rectangular tunnel as well as the 12-ft Pressure Tunnel. An improved wall interference assessment and correction method for three-dimensional wind tunnel testing is presented in the appendix.

Lo, Ching F.↗

A probabilistic analysis of cumulative carbon emissions and long-term planetary warming

Efforts to mitigate and adapt to long-term climate change could benefit greatly from probabilistic estimates of cumulative carbon emissions due to fossil fuel burning and resulting CO 2 -induced planetary warming. Here we demonstrate the use of a reduced-form model to project these variables. We performed simulations using a large-ensemble framework with parametric uncertainty sampled to produce distributions of future cumulative emissions and consequent planetary warming. A hind-cast ensemble of simulations captured 1980–2012 historical CO 2 emissions trends and an ensemble of future projection simulations generated a distribution of emission scenarios that qualitatively resembled the suite of Representative and Extended Concentration Pathways. The resulting cumulative carbon emission and temperature change distributions are characterized by 5–95th percentile ranges of 0.96–4.9 teratonnes C (Tt C) and 1.4 °C–8.5 °C, respectively, with 50th percentiles at 3.1 Tt C and 4.7 °C. Within the wide range of policy-related parameter combinations that produced these distributions, we found that low-emission simulations were characterized by both high carbon prices and low costs of non-fossil fuel energy sources, suggesting the importance of these two policy levers in particular for avoiding dangerous levels of climate warming. With this analysis we demonstrate a probabilistic approach to the challenge of identifying strategies for limiting cumulative carbon emissions and assessing likelihoods of surpassing dangerous temperature thresholds.

54 ENVIRONMENTAL SCIENCES↗

Effects of Temperature on Development and Voltinism of Chaetodactylus krombeini (Acari: Chaetodactylidae): Implications for Climate Change Impacts

Temperature plays an important role in the growth and development of arthropods, and thus the current trend of climate change will alter their biology and species distribution. We used Chaetodactylus krombeini (Acari: Chaetodactylidae), a cleptoparasitic mite associated with Osmia bees (Hymenoptera: Megachilidae), as a model organism to investigate how temperature affects the development and voltinism of C. krombeini in the eastern United States. The effects of temperature on the stage-specific development of C. krombeini were determined at seven constant temperatures (16.1, 20.2, 24.1, 27.5, 30.0, 32.4 and 37.8°C). Parameters for stage-specific development, such as threshold temperatures and thermal constant, were determined by using empirical models. Results of this study showed that C. krombeini eggs developed successfully to adult at all temperatures tested except 37.8°C. The nonlinear and linear empirical models were applied to describe quantitatively the relationship between temperature and development of each C. krombeini stage. The nonlinear Lactin model estimated optimal temperatures as 31.4, 32.9, 32.6 and 32.5°C for egg, larva, nymph, and egg to adult, respectively. In the linear model, the lower threshold temperatures were estimated to be 9.9, 14.7, 13.0 and 12.4°C for egg, larva, nymph, and egg to adult, respectively. The thermal constant for each stage completion were 61.5, 28.1, 64.8 and 171.1 degree days for egg, larva, nymph, and egg to adult, respectively. Under the future climate scenarios, the number of generations (i.e., voltinism) would increase more likely by 1.5 to 2.0 times by the year of 2100 according to simulation. Lastly, the findings herein firstly provided comprehensive data on thermal development of C. krombeini and implications for the management of C. krombeini populations under global warming were discussed.

54 ENVIRONMENTAL SCIENCES↗

Why DA and DB white dwarfs do not show coronal activity and p-mode oscillations

The problems of nonradiative heating of outer atmospheric layers and p-mode oscillations in white dwarfs caused by acoustic waves generated in convective zones are discussed. These effects have been studied by calculating the cutoff periods for adiabatic and isothermal waves propagating in atmospheres of DA and DB stars with Teff greater than or equal 20,000 K and log g = 6-9. The obtained cutoff periods are approximately bounded by 0.01 and 40 sec for high- and low-gravity white dwarfs, respectively. Expected amplitudes of p-mode oscillations corresponding to trapped acoustic waves with small angular wave numbers are estimated, indicating that the amplitudes could be observed as Doppler shifts of spectral lines which might be detectable if adequate spectral resolution were available. The luminosity variations corresponding to these amplitudes are unlikely to be observable when all damping processes are accounted for. Results also indicate that the present theory of convection predicts some irregularities in the behavior of physical parameters.

Musielak, Z. E.↗

Electron Neutrino Charged-Current Quasielastic Scattering in the MINERvA Experiment

The electron-neutrino charged-current quasielastic (CCQE) cross section on nuclei is an important input parameter for electron neutrino appearance oscillation experiments. Current experiments typically begin with the muon neutrino cross section and apply theoretical corrections to obtain a prediction for the electron neutrino cross section. However, at present no experimental verification of the estimates for this channel at an energy scale appropriate to such experiments exists. We present the cross sections for a CCQE-like process determined using the MINERvA detector, which are the first measurements of any exclusive reaction in few-GeV electron neutrino interactions. The result is given as differential cross-sections vs. the electron energy, electron angle, and square of the four-momentum transferred to the nucleus, $Q^{2}$. We also compute the ratio to a muon neutrino cross-section in $Q^{2}$ from MINERvA. We find satisfactory agreement between these measurements and the predictions of the GENIE generator. We furthermore report on a photon-like background unpredicted by the generator which we interpret as neutral-coherent diffractive scattering from hydrogen.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗

Quantification and assessment of the atmospheric boundary layer height measured during the AWAKEN experiment by a scanning LiDAR

The atmospheric boundary layer (ABL) height plays a key role in many atmospheric processes as one of the dominant flow length scales. However, a systematic quantification of the ABL height over the entire range of scales (i.e., with periods ranging from one minute to one year) is still lacking in literature. In this work, the ABL height is quantified based on high-resolution measurements collected by a scanning pulsed Doppler LiDAR during the recent American WAKE experimeNt (AWAKEN) campaign. The high availability of ABL height estimates (≈2200 collected over one year and each of them based on 10-min averaged statistics) allows to robustly assess five different ABL height models, i.e., one for convective thermal conditions and four for stable conditions. Thermal condition is quantified by a stability parameter spanning three orders of magnitude and probed by near-ground 3D sonic anemometry. The free-atmosphere stability, quantified by the Brunt–Väisälä frequency, is both calculated from simultaneous radiosonde measurements and obtained from the best fit of two of the chosen ABL height models. Good agreement is found between the data and three of the chosen models, quantified by mean absolute errors on the ABL height between 281 and 585 m. Furthermore, the seasonal variability of the convective ABL height model parameters (−15% to +23% with respect to the year baseline) agrees with the variability of buoyancy-generated turbulence caused by the variation in solar radiation throughout the year.

17 WIND ENERGY↗

A New Analytical Expression for Estimating the Adiabatic Temperature Rise in Lithium-Ion Batteries During High-Power Pulses

A coupled, thermal-electrochemical model is used to explain why the Joule heating assumption (i.e., I 2 R) does not provide a good representation of the temperature rise during high power pulses in lithium-ion batteries, even in cases where the reversible heat generation can be neglected. The poor agreement occurs because the internal resistance changes during the pulse due to the opposing effects of mass transport in the electrolyte, which raises the resistance, and heat generation (temperature rise), which lowers the resistance. These insights are used to propose a new analytical expression for predicting the temperature rise during adiabatic pulses with limited experimental and physical data. The expression accounts for thermal effects using a Taylor series expansion of an Arrhenius-type equation. It accounts for transport effects using an approximate solution to the one-dimensional diffusion equation. The new expression is shown to accurately estimate the simulated, adiabatic temperature rise across a range of loadings (1 to 4 mAh cm -2 ) and C-rates (1C to 10C) for cells containing LiNi 0.5 Mn 0.3 Co 0.2 O 2 positive electrodes and graphite negative electrodes. In conclusion, it is also shown to accurately estimate the experimental, adiabatic temperature rise measured for cells with LiFePO 4 positive electrodes and mesocarbon negative electrodes with minimal changes to the fitted parameters.

25 ENERGY STORAGE↗

Quantitative Ceilometer-Radar Studies of Clouds Field Campaign Report

The Long Island Solar Farm (LISF) is a 32-megawatt solar photovoltaic power plant built through a collaboration including BP Solar, the Long Island Power Authority (LIPA), and the U.S. Department of Energy (DOE). The LISF, located on the Brookhaven National Laboratory site, began delivering power to the LIPA grid in November 2011, and is currently one of the largest solar photovoltaic power plants in the eastern United States. It is generating enough renewable energy to power approximately 4,500 homes and is helping New York State meet its clean energy and carbon reduction goals. Brookhaven National Laboratory (BN)L has the LISF instrumented and is capturing solar insolation and power data that can be used for research purposes. Additional information on LISF is available here: https://www.bnl.gov/lisf/. A network of nine high-definition (HD cameras) and 32 pyranometers are deployed within the LISF for the purpose of monitoring the location and characteristics of clouds and available global horizontal irradiance at the surface across the region. Information collected by this network every 30 s is currently used in BNL’s solar NowCasting algorithm, which forecasts near-term solar energy availability accounting for the behavior of clouds. The most important parameters for forecasting the near-term solar energy availability are accurate estimations of their horizontal extend, horizontal motion, and cloud base height. Currently, images from nearly located cameras are used to estimate the cloud base height. In order to evaluate the potential of this method for estimating the cloud base height, we decided to deploy a Vaisala ceilometer at the LISF.

14 SOLAR ENERGY↗

Impact of aerosols on ice crystal size

The interactions between aerosols and ice clouds represent one of the largest uncertainties in global radiative forcing from pre-industrial time to the present. In particular, the impact of aerosols on ice crystal effective radius (R ei ), which is a key parameter determining ice clouds' net radiative effect, is highly uncertain due to limited and conflicting observational evidence. Here we investigate the effects of aerosols on R ei under different meteorological conditions using 9-yearsatellite observations. We find that the responses of R ei to aerosol loadings are modulated by water vapor amount in conjunction with several other meteorological parameters. While there is a significant negative correlation between Rei and aerosol loading in moist conditions, consistent with the “Twomey effect” for liquid clouds, a strong positive correlation between the two occurs in dry conditions. Simulations based on a cloud parcel model suggest that water vapor modulates the relative importance of different ice nucleation modes, leading to the opposite aerosol impacts between moist and dry conditions. When ice clouds are decomposed into those generated from deep convection and formed in situ, the water vapor modulation remains in effect for both ice cloud types, although the sensitivities of R ei to aerosols differ noticeably between them due to distinct formation mechanisms. The water vapor modulation can largely explain the difference in the responses of R ei to aerosol loadings in various seasons. A proper representation of the water vapor modulation is essential for an accurate estimate of aerosol–cloud radiative forcing produced by ice clouds.

54 ENVIRONMENTAL SCIENCES↗

Foreground Bias from Parametric Models of Far-IR Dust Emission

We use simple toy models of far-IR dust emission to estimate the accuracy to which the polarization of the cosmic microwave background can be recovered using multi-frequency fits, if the parametric form chosen for the fitted dust model differs from the actual dust emission. Commonly used approximations to the far-IR dust spectrum yield CMB residuals comparable to or larger than the sensitivities expected for the next generation of CMB missions, despite fitting the combined CMB plus foreground emission to precision 0.1 percent or better. The Rayleigh-Jeans approximation to the dust spectrum biases the fitted dust spectral index by (Delta)(Beta)(sub d) = 0.2 and the inflationary B-mode amplitude by (Delta)(r) = 0.03. Fitting the dust to a modified blackbody at a single temperature biases the best-fit CMB by (Delta)(r) greater than 0.003 if the true dust spectrum contains multiple temperature components. A 13-parameter model fitting two temperature components reduces this bias by an order of magnitude if the true dust spectrum is in fact a simple superposition of emission at different temperatures, but fails at the level (Delta)(r) = 0.006 for dust whose spectral index varies with frequency. Restricting the observing frequencies to a narrow region near the foreground minimum reduces these biases for some dust spectra but can increase the bias for others. Data at THz frequencies surrounding the peak of the dust emission can mitigate these biases while providing a direct determination of the dust temperature profile.

cosmic microwave background↗

A Polywell Fusion Reactor Designed for Net Power Generation

A brief history of Polywell progress is recounted. The present PIC simulation explains why the most recent Polywell fusion reactor failed to produce fusion energy. Synchronized variations of multiple parameters would require DC power supplies, not available in historic model testing. Even with DC power, the simulation showed that the trapping of cold electrons would ruin plasma stability during start-up. A theoretical solution to this trapping problem was found in Russian literature describing diocotron-pumping of electrons out of a plasma trap at Kharkov Institute. In Polywell, diocotron-pumping required matching the depth of the potential-well to the electron-beam current falling on a special aperture installed in one of the electromagnets. With diocotron-pumping the reactor was simulated to reach steady-state, maximum-power operation in a few milliseconds of simulated time. These improvements, validated in simulating small-scale DD reactors, were scaled up by a factor of 30 to simulate a large, net-power reactor burning p + {sup 11}B fuel. Power-balance was estimated from a textbook formula for fusion power density by numerically integrating the power density. Unity power-balance required the size of the p + {sup 11}B reactor to be somewhat larger than ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preface to the Special Issue on TOUGH Symposium 2015

The TOUGH Symposium 2015 was held in Berkeley, California, September 28-30, 2015. The TOUGH family of codes, developed at the Energy Geosciences Division of Lawrence Berkeley National Laboratory (LBNL), is a suite of computer programs for the simulation of multiphase and multicomponent fluid and heat flows in porous and fractured media with applications in many geosciences fields, such as geothermal reservoir engineering, nuclear waste disposal, geological carbon sequestration, oil and gas reservoirs, gas hydrate research, vadose zone hydrology and environmental remediation. Since the first release in the 1980s, many modifications and enhancements have been continuously made to TOUGH and its various descendants (iTOUGH2, TOUGH+, TOUGH-MP, TOUGHREACT, TOUGH+HYDRATE, TMVOC...), at LBNL and elsewhere. Today, these codes are used worldwide in academia, government organizations and private companies in problems involving coupled hydrological, thermal, biogeochemical and geomechanical processes. The Symposia, organized every 2-3 years, bring together developers and users for an open exchange on recent code enhancements and applications. In 2015, the Symposium was attended by one hundred participants, representing thirty-four nationalities. This Special Issue in Computers & Geosciences gathers extended versions of selected Symposium proceedings related to (i) recent enhancements to the TOUGH family of codes and (ii) coupled flow and geomechanics processes modeling. In this special issue, a paper by Jung et al. provides a new base version of TOUGH, TOUGH3. This latest version addresses the computational robustness and efficiency needed for the increasing complexity of the simulated processes and the growing size of model domains. Moreover, it consolidates the serial and parallel versions into a single code, hence, inconsistencies between different versions are less likely to occur, and code maintenance and development become more efficient. Finsterle et al. present a review of the capabilities of iTOUGH2 and illustrate the workflow of an integrated simulation-optimization framework. Liu et al. present an approach to reduce the computational cost of performing Bayesian inverse modeling and demonstrate their approach using a synthetic experiment simulated with TOUGH2. Regarding the mesh generation process, the paper by Doughty presents a generalization of the grid-generation algorithm embedded in TOUGH2 (MESHMAKER) to create one-column grids with fractal dimension, useful for representing flow through fractures or highly heterogeneous media with fractal geometry, and also for inverse methods. Sentís and Gable present a module to interface the open source code Lagrit with TOUGH2, that assists in the processes of mesh generation and model setup, and that can be particularly useful for complex geometries and large models. Bonduà et al. present a pre-processing tool (VORO2MESH) that features a full 3D Voronoi approach for structured and unstructured grids, and an enhanced post-processor (TOUGH2Viewer) for the visualization of 3D Voronoi grids. The literature review shows that coupled fluid flow and geomechanics has gained increasing attention in the last years in many geosciences fields. In this context, several codes of the TOUGH family have been successfully coupled with geomechanics software. The paper by Rutqvist presents a review of some of these TOUGH-based geomechanics models, of which the TOUGH-FLAC simulator – presented for the first time in the TOUGH Symposium 2003 – is well-known worldwide and has been used in a large variety of applications, encompassing the application domains of the TOUGH codes. Blanco-Martín et al. present an extension of TOUGH-FLAC to the finite strain framework, motivated by the existence of large strains associated with rock salt deformation in the context of heat-generating nuclear waste disposal. Following a paper published in the Special Issue of the TOUGH Symposium 2012 in this journal, Kim et al. present the current status of the TOUGH-RSBN simulator for coupled hydro-mechanical modeling of hydraulic fracturing using a discrete fracture network approach, and present some validation and demonstration cases. Rinaldi and Nespoli present TOUGH-SEED, an injection-induced seismicity modeling approach using TOUGH2 and a geomechanics stochastic model, and apply it to the Basel EGS project. Finally, Rinaldi et al. use the PEST protocol implemented in iTOUGH2, combined with TOUGH-FLAC as the forward model, to perform inverse modeling for the In Salah CO2 demonstration site in Algeria. Using pore pressures and InSAR displacements as observations, this approach allows for a multi-physics joint inversion and takes advantage of the inverse capabilities of iTOUGH2 (parameter estimation, sensitivity analysis, uncertainty estimation …), applied to coupled fluid flow and geomechanics. On behalf of the members of the TOUGH Symposium 2015 Organizing Committee, I hope that this collection of papers will provide the readers of Computers & Geosciences with a recent view of the capabilities of the TOUGH family of codes, as well as with current modeling research domains in the geosciences. Last but not least, I would like to express my gratitude to all the reviewers who kindly contributed to this Special Issue, as well as the Journal staff who have made this publication possible.

99 GENERAL AND MISCELLANEOUS↗

A Comprehensive View of Global Potential for Hydro-generated Electricity

In this study, we assess global hydropower potential using runoff and stream flow data, along with turbine technology performance, cost assumptions, and environmental considerations. The results provide the first comprehensive quantification of global hydropower potential including gross, technical, economic, and exploitable estimates. Total global potential of gross, technical, economic, and exploitable hydropower are estimated to be approximately 128, 39, 32, and 27 petawatt hours per year, respectively. The economic and exploitable potential of hydropower are calculated at less than 9 cents/kWh. We find that hydropower has the potential to supply a significant portion of the world energy needs, although this potential varies substantially by region. Globally, hydropower can potentially supply about 1.5 times the total electricity demand in 2005. Estimated hydropower resources in a number of countries are sufficient to accommodate their demand for electricity in 2005, e.g., Brazil (5.6 times), Russia (4.6 times), and Canada (3.5 times). A sensitivity analysis indicates that hydropower estimates are not highly sensitive to five key parameters: design flow (varying by -2% to +1% at less than 9 cents/kWh), cost and financing options (by -7% to +6%), turbine efficiency (by -10% to +10%), stream flow (by -10% to +10%), and fixed charge rate (by -6% to 5%). This sensitivity analysis emphasizes the reliable role of hydropower for future energy systems, when compared to other renewable energy resources with larger uncertainty in their future potentials.

Zhou, Yuyu↗