Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Detection and Characterization of Preferential Flow Paths in the Downstream Area of a Hazardous Landfill

We have used surface-based electrical resistivity tomography to detect and characterize preferential hydraulic pathways in the immediate downstream area of an abandoned, hazardous landfill. The landfill occupies the void left by a former gravel pit and its base is close to the groundwater table and lacking an engineered barrier. As such, this site is remarkably typical of many small- to medium-sized waste deposits throughout the densely populated and heavily industrialized foreland on both sides of the Alpine arc. Outflows of pollutants lastingly contaminated local drinking water supplies and necessitated a partial remediation in the form of a synthetic cover barrier, which is meant to prevent meteoric water from percolating through the waste before reaching the groundwater table. Any future additional isolation of the landfill in the form of lateral barriers thus requires adequate knowledge of potential preferential hydraulic pathways for outflowing contaminants. Our results, inferred from a suite of tomographically inverted surfaced-based electrical resistivity profiles oriented roughly perpendicular to the local hydraulic gradient, indicate that potential contaminant outflows would predominantly occur along an unexploited lateral extension of the original gravel deposit. This finds its expression as a distinct and laterally continuous high-resistivity anomaly in the resistivity tomograms. This interpretation is ground-truthed through a litholog from a nearby well. Since the probed glacio-fluvial deposits are largely devoid of mineralogical clay, the geometry of hydraulic and electrical pathways across the pore space of a given lithological unit can be assumed to be identical, which allows for an order-of-magnitude estimation of the overall permeability structure. These estimates indicate that the permeability of the imaged extension of the gravel body is at least two to three orders-of-magnitude higher than that of its finer-grained embedding matrix. This corroborates the preeminent role of the high-resistivity anomaly as a potential preferential flow path.

Non-constant discounting in finite horizon: The free terminal time case

This paper derives the HJB (Hamilton-Jacobi-Bellman) equation for sophisticated agents in a finite horizon dynamic optimization problem with non-constant discounting in a continuous setting, by using a dynamic programming approach. A simple example is used in order to illustrate the applicability of this HJB equation, by suggesting a method for constructing the subgame perfect equilibrium solution to the problem.Conditions for the observational equivalence with an associated problem with constantdiscounting are analyzed. Special attention is paid to the case of free terminal time. Strotz¿s model (an eating cake problem of a nonrenewable resource with non-constant discounting) is revisited.

Neuromonitoring after major neurosurgical procedures.

Postoperative care of major neurosurgical procedures is aimed at the prevention, detection and treatment of secondary brain injury. This consists of a series of pathological events (i.e. brain edema and intracranial hypertension, cerebral hypoxia/ischemia, brain energy dysfunction, non-convulsive seizures) that occur early after the initial insult and surgical intervention and may add further burden to primary brain injury and thus impact functional recovery. Management of secondary brain injury requires specialized neuroscience intensive care units (ICU) and continuous advanced monitoring of brain physiology. Monitoring of intracranial pressure (ICP) is a mainstay of care and is recommended by international guidelines. However, ICP monitoring alone may be insufficient to detect all episodes of secondary brain insults. Additional invasive (i.e. brain tissue PO2, cerebral microdialysis, regional cerebral blood flow) and non-invasive (i.e. transcranial doppler, near-infrared spectroscopy, EEG) brain monitoring devices might complement ICP monitoring and help clinicians to target therapeutic interventions (e.g. management of cerebral perfusion pressure, blood transfusion, glucose control) to patient-specific pathophysiology. Several independent studies demonstrate such multimodal approach may optimize patient care after major neurosurgical procedures. The aim of this review is to evaluate some of the available monitoring systems and summarize recent important data showing the clinical utility of multimodal neuromonitoring for the management of main acute neurosurgical conditions, including traumatic brain injury, subarachnoid hemorrhage and stroke.

Velocity-jump instabilities in Hele-Shaw flow of associating polymer solutions

We study fracturelike flow instabilities that arise when water is injected into a Hele-Shaw cell filled with aqueous solutions of associating polymers. We explore various polymer architectures, molecular weights, and solution concentrations. Simultaneous measurements of the finger tip velocity and of the pressure at the injection point allow us to describe the dynamics of the finger in terms of the finger mobility, which relates the velocity to the pressure gradient. The flow discontinuities, characterized by jumps in the finger tip velocity, which are observed in experiments with some of the polymer solutions, can be modeled by using a nonmonotonic dependence between a characteristic shear stress and the shear rate at the tip of the finger. A simple model, which is based on a viscosity function containing both a Newtonian and a non-Newtonian component, and which predicts nonmonotonic regions when the non-Newtonian component of the viscosity dominates, is shown to agree with the experimental data.

How to improve flow during cardiopulmonary bypass in an acardia experimental model.

OBJECTIVES In extreme scenarios, such as hyperacute rejection of heart transplant, an urgent heart explantation might be necessary. The aim of this experimental study was to determine the feasibility and to improve the haemodynamics of a venoarterial cardiopulmonary bypass after cardiectomy. METHODS A venoarterial cardiopulmonary bypass was established in seven calves (56.4 ± 7 kg) by the transjugular insertion to the caval axis of a self-expanding cannula, with a carotid artery return. After baseline measurements (A), ventricular fibrillation was induced (B), great arteries were clamped (C), the heart was excised and the right and left atria remnants, containing the pulmonary veins, were sutured together leaving an atrial septal defect over the cannula in the caval axis (D). Measurements were taken with the pulmonary artery clamped and declamped. RESULTS Initial pump flow was 4.16 ± 0.75 l/min dropping to 2.9 ± 0.63 l/min (P(AB )< 0.001) 10 min after induction of ventricular fibrillation. After cardiectomy with the pulmonary artery clamped, the pump flow increased non-significantly to 3.20 ± 0.78 l/min. After declamping, the flow significantly increased close to baseline levels (3.61 ± 0.73 l/min, P(DB )= 0.009, P(DC )= 0.017), supporting the notion that full cardiopulmonary bypass in acardia is feasible only if adequate drainage of pulmonary circulation is assured to avoid pulmonary congestion and loss of volume from the left-to-right shunt of bronchial vessels.

Production Flow Analysis - cases from manufacturing and service industry

Production flow analysis (PFA) is a well-established methodology used for transforming traditional functional layout into product-oriented layout. The method uses part routings to find natural clusters of workstations forming production cells able to complete parts and components swiftly with simplified material flow. Once implemented, the scheduling system is based on period batch control aiming to establish fixed planning, production and delivery cycles for the whole production unit. PFA is traditionally applied to job-shops with functional layouts, and after reorganization within groups lead times reduce, quality improves and motivation among personnel improves. Several papers have documented this, yet no research has studied its application to service operations management. This paper aims to show that PFA can well be applied not only to job-shop and assembly operations, but also to back-office and service processes with real cases. The cases clearly show that PFA reduces non-value adding operations, introduces flow by evening out bottlenecks and diminishes process variability, all of which contribute to efficient operations management.

Volcanic arcs fed by rapid pulsed fluid flow through subducting slabs

At subduction zones, oceanic lithosphere that has interacted with sea water is returned to the mantle, heats up during descent and releases fluids by devolatilization of hydrous minerals. Models for the formation of magmas feeding volcanoes above subduction zones require largescale transport of these fluids into overlying mantle wedges(1-3). Fluid flow also seems to be linked to seismicity in subducting slabs. However, the spatial and temporal scales of this fluid flow remain largely unknown, with suggested timescales ranging from tens to tens of thousands of years(3-5). Here we use the Li-Ca-Sr isotope systems to consider fluid sources and quantitatively constrain the duration of subduction-zone fluid release at similar to 70 km depth within subducting oceanic lithosphere, now exhumed in the Chinese Tianshan Mountains. Using lithium-diffusion modelling, we find that the wall-rock porosity adjacent to the flowpath of the fluids increased ten times above the background level. We show that fluids released by devolatilization travelled through the slab along major conduits in pulses with durations of about similar to 200 years. Thus, although the overall slab dehydration process is continuous over millions of years and over a wide range of pressures and temperatures, we conclude that the fluids produced by dehydration in subducting slabs are mobilized in short-lived, channelized fluid-flow events.

Plant architecture of Paspalum vaginatum schwartz modified by nitrate and ammonium nutrition

Paspalum vaginatum Schwartz plants were grown under greenhouse conditions in a continuous-flow hydroponic culture, containing NO3- or NH4+or NH4NO3 as nitrogen source. After 30 days, the size of aerial biomass and root system decreased significantly when plants were supplied with NH4+as exclusive nitrogen source. Compared to NO3- treatment, reducing and non-reducing sugars were decreasing together with a significant increase in amino acids content. NH4+-nutrition caused tillers to grow toward an orthogravitropic position (average angle of 68° with respect to the horizontal), and with NO3--nutrition, tillers tended to become diagravitropic (average angle of 23°). With NH4NO3 all the parameters measured had values in between those of the other two sources. Thus, the morphologic differences among plants growing in NO3- or NH4+ nutrition confirm the hypothesis that nitrogen source determines the growth habit of tillers in P. vaginatum by modulating the endogenous levels of reducing-non-reducing sugars.

"Suddenly I get into the zone": Examining discontinuities and nonlinear changes in flow experiences at work

Work-related flow is defined as a sudden and enjoyable merging of action and awareness that represents a peak experience in the daily lives of workers. Employees" perceptions of challenge and skill and their subjective experiences in terms of enjoyment, interest and absorption were measured using the experience sampling method, yielding a total of 6981 observations from a sample of 60 employees. Linear and nonlinear approaches were applied in order to model both continuous and sudden changes. According to the R2, AICc and BIC indexes, the nonlinear dynamical systems model (i.e. cusp catastrophe model) fit the data better than the linear and logistic regression models. Likewise, the cusp catastrophe model appears to be especially powerful for modelling those cases of high levels of flow. Overall, flow represents a nonequilibrium condition that combines continuous and abrupt changes across time. Research and intervention efforts concerned with this process should focus on the variable of challenge, which, according to our study, appears to play a key role in the abrupt changes observed in work-related flow.

Debris flows as a factor of hillslope evolution controlled by a continuous or a pulse process?

Flood effectiveness observations imply that two families of processes describe the formation of debris flow volume. One is related to the rainfall?erosion relationship, and can be seen as a gradual process, and one is related to additional geological/geotechnical events, those named hereafter extraordinary events. In order to discuss the hypothesis of coexistence of two modes of volume formation, some methodologies are applied. Firstly, classical approaches consisting in relating volume to catchments characteristics are considered. These approaches raise questions about the quality of the data rather than providing answers concerning the controlling processes. Secondly, we consider statistical approaches (cumulative number of events distribution and cluster analysis) and these suggest the possibility of having two distinct families of processes. However the quantitative evaluation of the threshold differs from the one that could be obtained from the first approach, but they all agree in the sense of the coexistence of two families of events. Thirdly, a conceptual model is built exploring how and why debris flow volume in alpine catchments changes with time. Depending on the initial condition (sediment production), the model shows that large debris flows (i.e. with important volume) are observed in the beginning period, before a steady-state is reached. During this second period debris flow volume such as is observed in the beginning period is not observed again. Integrating the results of the three approaches, two case studies are presented showing: (1) the possibility to observe in a catchment large volumes that will never happen again due to a drastic decrease in the sediment availability, supporting its difference from gradual erosion processes; (2) that following a rejuvenation of the sediment storage (by a rock avalanche) the magnitude?frequency relationship of a torrent can be differentiated into two phases, the beginning one with large and frequent debris flow and a later one with debris flow less intense and frequent, supporting the results of the conceptual model. Although the results obtained cannot identify a clear threshold between the two families of processes, they show that some debris flows can be seen as pulse of sediment differing from that expected from gradual erosion.

Evaluating effects of spectral training data distribution on continuous field mapping performance

Continuous field mapping has to address two conflicting remote sensing requirements when collecting training data. On one hand, continuous field mapping trains fractional land cover and thus favours mixed training pixels. On the other hand, the spectral signature has to be preferably distinct and thus favours pure training pixels. The aim of this study was to evaluate the sensitivity of training data distribution along fractional and spectral gradients on the resulting mapping performance. We derived four continuous fields (tree, shrubherb, bare, water) from aerial photographs as response variables and processed corresponding spectral signatures from multitemporal Landsat 5 TM data as explanatory variables. Subsequent controlled experiments along fractional cover gradients were then based on generalised linear models. Resulting fractional and spectral distribution differed between single continuous fields, but could be satisfactorily trained and mapped. Pixels with fractional or without respective cover were much more critical than pure full cover pixels. Error distribution of continuous field models was non-uniform with respect to horizontal and vertical spatial distribution of target fields. We conclude that a sampling for continuous field training data should be based on extent and densities in the fractional and spectral, rather than the real spatial space. Consequently, adequate training plots are most probably not systematically distributed in the real spatial space, but cover the gradient and covariate structure of the fractional and spectral space well. (C) 2009 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS). Published by Elsevier B.V. All rights reserved.

Portable School Stop Signs and Other Non-Uniform School Stop Control Devices, HR-190, 1978

Portable (roll-out) stop signs are used at school crossings in over 300 cities in Iowa. Their use conforms to the Code of Iowa, although it is not consistent with the provisions of the Manual on Uniform Traffic Control Devices adopted for nationwide application. A survey indicated that most users in Iowa believe that portable stop signs provide effective protection at school crossings, and favor their continued use. Other non-uniform signs that fold or rotate to display a STOP message only during certain hours are used at school crossings in over 60 cities in Iowa. Their use does not conform to either the Code of Iowa or the Manual on Uniform Traffic Control Devices. Users of these devices also tend to favor their continued use. A survey of other states indicated that use of temporary devices similar to those used in Iowa is not generally sanctioned. Some unsanctioned use apparently occurs in several states, however. A different type of portable stop sign for school crossings is authorized and widely used in one state. Portable stop signs similar to those used in Iowa are authorized in another state, although their use is quite limited. A few reports in the literature reviewed for this research discussed the use of portable stop signs. The authors of these reports uniformly recommended against the use of portable or temporary traffic control devices. Various reasons for this recommendation were given, although data to support the recommendation were not offered. As part of this research, field surveys were conducted at 54 locations in 33 communities where temporary stop control devices were in use at school crossings. Research personnel observed the obedience to stop control and measured the vehicular delay incurred. Stopped delay averaged 1.89 seconds/entering vehicle. Only 36.6 percent of the vehicles were observed to come to a complete stop at the study locations controlled by temporary stop control devices. However, this level of obedience does not differ from that observed at intersections controlled by permanent stop signs. Accident experience was compiled for 76 intersections in 33 communities in Iowa where temporary stop signs were used and, for comparative purposes, at 76 comparable intersections having other forms of control or operating without stop control. There were no significant differences in accident experience An economic analysis of vehicle operating costs, delay costs, and other costs indicated that temporary stop control generated costs only about 12 percent as great as permanent stop control for a street having a school crossing. Midblock pedestrian-actuated signals were shown to be cost effective in comparison with temporary stop signs under the conditions of use assumed. Such signals could be used effectively at a number of locations where temporary stop signs are being used. The results of this research do not provide a basis for recommending that use of portable stop signs be prohibited. However, erratic patterns of use of these devices and inadequate designs suggest that improved standards for their use are needed. Accordingly, nine recommendations are presented to enhance the efficiency of vehicular flow at school crossings, without causing a decline in the level of pedestrian protection being afforded.

Non-viral gene therapy for GDNF production in RCS rat: the crucial role of the plasmid dose.

Glial cell line-derived neurotrophic factor (GDNF) is one of the candidate molecules among neurotrophic factors proposed for a potential treatment of retinitis pigmentosa (RP). It must be administered repeatedly or through sustained releasing systems to exert prolonged neuroprotective effects. In the dystrophic Royal College of Surgeon's (RCS) rat model of RP, we found that endogenous GDNF levels dropped during retinal degeneration time course, opening a therapeutic window for GDNF supplementation. We showed that after a single electrotransfer of 30 μg of GDNF-encoding plasmid in the rat ciliary muscle, GDNF was produced for at least 7 months. Morphometric, electroretinographic and optokinetic analyses highlighted that this continuous release of GDNF delayed photoreceptors (PRs) as well as retinal functions loss until at least 70 days of age in RCS rats. Unexpectedly, increasing the GDNF secretion level accelerated PR degeneration and the loss of electrophysiological responses. This is the first report: (i) demonstrating the efficacy of GDNF delivery through non-viral gene therapy in RP; (ii) establishing the efficacy of intravitreal administration of GDNF in RP associated with a mutation in the retinal pigment epithelium; and (iii) warning against potential toxic effects of GDNF within the eye/retina.