Search NASA⌕ Search

SEARCH · Search NASA

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,153 records · Page 64

Variance-Reduced Accelerated First-Order Methods: Central Limit Theorems and Confidence Statements

In this paper, we consider a strongly convex stochastic optimization problem and propose three classes of variable sample-size stochastic first-order methods: (i) the standard stochastic gradient descent method, (ii) its accelerated variant, and (iii) the stochastic heavy-ball method. In each scheme, the exact gradients are approximated by averaging across an increasing batch size of sampled gradients. We prove that when the sample size increases at a geometric rate, the generated estimates converge in mean to the optimal solution at an analogous geometric rate for schemes (i)–(iii). Based on this result, we provide central limit statements, whereby it is shown that the rescaled estimation errors converge in distribution to a normal distribution with the associated covariance matrix dependent on the Hessian matrix, the covariance of the gradient noise, and the step length. If the sample size increases at a polynomial rate, we show that the estimation errors decay at a corresponding polynomial rate and establish the associated central limit theorems (CLTs). Under certain conditions, we discuss how both the algorithms and the associated limit theorems may be extended to constrained and nonsmooth regimes. As a result, we provide an avenue to construct confidence regions for the optimal solution based on the established CLTs and test the theoretical findings on a stochastic parameter estimation problem.

Lei, Jinlong↗

Evaluating the Applicability of Data on Total Solar-Radiation Intensity Derived from Various Sources of Actinometric Information

In designing generating solar-engineering devices, the first priority is the predesign study of the potential of local resources of solar energy at the places where they may be located. Currently, actinometric information is obtained in several ways. They include direct measurement of solar-radiation characteristics at meteorological stations or from satellites and analytical methods for calculating intensity of solar radiation. Despite the significant number of sources of actinometric information, it seems impossible to evaluate the potential of incoming solar energy in the region, to substantiate the composition and parameters of solar-power-plant equipment, and to predict its operation at the initial stage. Using the example of the town of Kamyshin in Volgograd oblast, the sufficiency of required actinometric information is analyzed for further solar-engineering calculations and the data accuracy and limits of their applicability are estimated.

14 SOLAR ENERGY↗

The SMAP Level-4 ECO Project: Linking the Terrestrial Water and Carbon Cycles

The SMAP (Soil Moisture Active Passive) Level-4 projects aims to develop a fully coupled hydrology-vegetation data assimilation algorithm to generate improved estimates of modeled hydrological fields and carbon fluxes. This includes using the new NASA Catchment-CN (Catchment-Carbon-Nitrogen) model, which combines the Catchment land surface hydrology model with dynamic vegetation components from the Community Land Model version 4 (CLM4). As such, Catchment-CN allows a more realistic, fully coupled feedback between the land hydrology and the biosphere. The L4 ECO project further aims to inform the model through the assimilation of Soil Moisture Active Passive (SMAP) brightness temperature observations as well as observations of Moderate Resolution Imaging Spectroradiometer (MODIS) fraction of absorbed photosynthetically active radiation (FPAR). Preliminary results show that the assimilation of SMAP observations leads to consistent improvements in the model soil moisture skill. An evaluation of the Catchment-CN modeled vegetation characteristics showed that a calibration of the model's vegetation parameters is required before an assimilation of MODIS FPAR observations is feasible.

MODIS↗

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.↗

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.↗

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↗