arxiv: v3 [math.na] 15 Dec 2009

Transcription

1 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING M.A. EBRAHIMI, MICHAEL HOLST, AND EVELYN LUNASIN arxiv: v3 [mat.na] 15 Dec 9 ABSTRACT. In tis paper we investigate te use of te D Navier-Stokes-Voigt (NSV) model for use in algoritms and explore its limits in te context of image inpainting. We begin by giving a brief review of te work of Bertalmio et. al. in 1 on exploiting an analogy between te image intensity function for te image inpainting problem and te stream function in D incompressible fluid. An approximate solution to te inpainting problem was ten obtained by numerically approximating te steady state solution of te D Navier-Stokes vorticity transport equation, and simultaneously solving te Poisson problem between te vorticity and stream function, in te region to be inpainted. Tis elegant approac allows one to produce an approximate solution to te image inpainting problem by using tecniques from computational fluid dynamics (CFD). Recently, te tree-dimensional (3D) Navier-Stokes-Voigt (NSV) model of viscoelastic fluid, was suggested by Cao, et. al. as an inviscid regularization to te 3D Navier-Stokes equations. We give some background on te NSV matematical model, describe wy it is a good candidate sub-grid-scale turbulence model, and propose tis model as an alternative for image inpainting. We describe an implementation of inpainting use te NSV model, and present numerical results comparing te resulting images wen using te NSE and NSV for inpainting. Our results sow tat te NSV model allows for a larger time step to converge to te steady state solution, yielding a more efficient numerical process wen automating te inpainting process. We compare quality of te resulting images using subjective measure (uman evaluation) and objected measure (by calculating te peak signal-to-noise ratio (PSNR), also known as peak signal-to-reconstructed measure). We also present some new teoretical results based on energy metods comparing te sufficient conditions for stability of te discretization sceme for te two model equations. Tese teoretical and numerical studies sed some ligt on wat can be expected from tis category of approac wen automating te inpainting problem. CONTENTS 1. Introduction. Image Inpainting Using Fluid Models 3.1. Inpainting Using te Navier-Stokes Equations 4.. Alpa Models and te Navier-Stokes-Voigt (NSV) Equations 5.3. Te Alpa Model as Sub-Grid-Scale Turbulence Model 7.4. Navier-Stokes-Voigt for Image Inpainting 7 3. Numerical Comparison NSE and NSV for Image Inpainting Stability Beavior of NSE and NSV for Inpainting Comparing te Efficiency of NSE and NSV for Inpainting 1 4. Some Supporting Teoretical Results Uniqueness of Steady State Solutions for Boundary Driven Flows Stability Analysis of te Sceme based on D NSE Stability Analysis of te Sceme based on D NSV Stability Analysis of a Semi-Implicit Sceme based on D NSE 1 5. Summary 3 6. Acknowledgments 4 References 4 Date: December 15, 9. Key words and prases. Navier-Stokes-Voigt, sub-grid scale turbulence model, fluid mecanics, image inpainting. MH was supported in part by NSF Awards and 915. ME and EL were supported in part by NSF Award

2 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN 1. INTRODUCTION Image inpainting (inpainting for sort) is te process of correcting a damaged image by filling in te missing or altered data of an image wit a better suited data for tat region. Te idea is to fill in te damaged part of an image using te information from surrounding areas. Te goal of researcers in tis area is to develop algoritms for automatic digital inpainting wic mimics te basic manual tecniques used by a professional restorer. In 1, Bertalmio, et. al. built a metod based on an analogy between image intensity and te stream function in a two-dimensional (D) incompressible fluid. Te solution to te inpainting problem is ten obtained by numerically approximating te steady state solution of te D Navier-Stokes vorticity transport equation, for some small viscosity, and simultaneously solving te Poisson problem between te vorticity and stream function in te region to be inpainted. Tis approac enables automation of te inpainting process by wic inpainting is done using tecniques from computational fluid dynamics (CFD). However, difficulties wic arise in CFD are also inerited. For inpainting, te viscosity in te fluid model sould be as small as possible to preserve edges. However, CFD simulation of ig Reynolds number flows (flows wit very small viscosity) requires very fine mes in order to resolve te wide range of scales of motion contributing to te dynamics of te flow. Sub-grid scale modeling is an alternative tat reduces te computational requirements wen simulating turbulent flows. Recently, te tree-dimensional (3D) Navier-Stokes-Voigt (NSV) model of viscoelastic fluid, α u t + u t ν u + (u )u + p = f, u =, (1.1) wit initial condition u(x, ) = u in (x) was suggested by Cao, et. al. as an inviscid regularization to te 3D Navier-Stokes equations (NSE), were te lengt-scale α is considered as te regularizing parameter, u is te velocity of te fluid wit viscosity ν >, p is te pressure and f is te body force. Note tat wen α =, we recover te Navier-Stokes equations of motion. Te system (1.1) was first introduced in 1973 by Oskolkov (see [3, 31]) as a model of a motion of linear, viscoelastic fluid. In tat setting, α is tougt of as a lengt-scale parameter caracterizing te elasticity of te fluid. Te 3D NSV model is sown to be globally well-posed [3, 31] and as a finite-dimensional global attractor [1], making it an attractive sub-grid scale turbulence model for numerical simulation. In te presence of pysical boundaries, te above regularization of te NSE is different in nature from oter alpa regularization models in [7, 8, 9, 1, 11, 13,, 4, 9] (see also Section.), because it does not require additional boundary conditions. In tis paper, we investigate a new approac wic arises by combining te recent tecnique in [4], wic uses ideas from classical fluid dynamics to propagate isopote lines, and recent results in [8,, 1], wic studies te NSV model of viscoelastic fluid as a candidate sub-grid scale turbulence model for purposes of numerical simulation. Oter PDEs ave been proposed recently for example in [6, 34] wic are possibly considered to be more efficient and peraps easier to implement. Of particular interest to us in tis area is to explore different sub-grid scale turbulent models applied to image inpainting. For te reasons outlined above, we ave cosen tis particular turbulence model to study in te context of inpainting. Using tis new model, we will study te effect of te regularization parameter α (wic one can tink of as te filter widt) on quality and efficiency of inpainting automation. As part of our investigation of D NSV for inpainting, we examine te differences between te D NSV model and te D NSE

3 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 3 for inpainting. We look at semi-implicit forward-time upwind metod for bot NSV and NSE and compare teir stability and efficiency. We will refer to tese numerical metods by teir corresponding matematical model: NSV and NSE numerical model. Our numerical results sow te NSV model, in comparison to NSE, yields a more stable solution to te inpainting process. Tat is, te NSV can be computed wit larger time-steps, reducing te computational expense in te automated inpainting procedure. However, we also study te efficiency of te calculation of te model equations as well as te quality of te resulting images. We ope tat tese teoretical and numerical studies sed some ligt on wat can be expected from tis category of approac wen automating te inpainting problem. Tis paper is organized as follows. In Section we review te work of Bertalmio et. al. in [4] on exploiting an analogy between image intensity and te stream function in D incompressible fluid. In Section. we give some background on te NSV matematical model, discuss wy it is a good candidate sub-grid-scale turbulence model, and argue tat it makes for a plausible alternative model for inpainting. In Section 3 we describe our implementation, and ten present our numerical results, comparing te resulting images wen using te NSE and NSV for inpainting. Te results sow te NSV model allows for larger time steps to converge to te steady state solution yielding a more efficient numerical process wen automating te inpainting process. We compare quality of te resulting images using subjective measure (uman evaluation) and objected measure (by calculating te peak signal-to-noise ratio (PSNR), also known as peak signal-to-reconstructed measure). We also give approximate operation counts needed by te two models to converge to steady state solution for some fixed tolerance. In Section 4 we present some supporting teoretical results. In Section 4.1, we summarize some existing results on te dependence of uniqueness of te solution to te inpainting problem on te image at te boundary, viscosity and on te size of te inpainting region. In Sections 4., 4.3, and 4.4, we establis some new results comparing te NSE and NSV in terms of stability conditions. We describe te discretization under consideration, and following [33], derive sufficient conditions for stability based on energy metods, establising results tat support our numerical results. In Section 5 we summarize our results and outline some future directions for tis work.. IMAGE INPAINTING USING FLUID MODELS In tis section we present a summary of te approac in [4] for automating te inpainting process motivated by ideas from classical fluid dynamics. Te idea is to propagate isopote lines from te exterior into te region to be inpainted, wic we denote as Ω, using te D NSE for fluid dynamics. Trougout te paper, we denote te digital grayscale image by I, an m n matrix wit gray-scale value to 55 at eac pixel, and D be te set of points (x, y) were I is defined. Te value represents black and te value 55 represents wite. Te values in between are different gray levels between black and wite. To solve te inpainting problem, a common tecnique for restorers (see e.g.,[5]) is to extend edges from te boundary of Ω, filling in te intra-region wit te correct color gradient. For example, in Figure 1, te direction of te isopotes (level lines of equal color gradients) determines were te smootness sould propagate. To see ow we can automate tis tecnique, we first introduce some basic matematical concepts. Matematically, te direction of isopotes (level curves of equal gray-levels) can be represented as I, wic indicates te direction of zero cange, and te smootness of te image,

4 4 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN FIGURE 1. Tis image sows te direction of te isopotes and were te smootness sould propagate. can be represented by I, were is te usual Laplacian operator (see [4] and references terein)..1. Inpainting Using te Navier-Stokes Equations. In [5], Bertalmio, et. al. proposed tat te solution I can be approximated wit te steady solution to te PDE of te form I t = I I + ν (g( I ) I), (.1) for small ν >, were te anisotropic diffusion is added to preserve te edges. We will list below a few diffusivity functions g tat are used in te literature. To be more precise, te main goal is to find te solution I suc tat te level curves of I almost parallel to I, tat is, I I. (.) Bertalmio et. al. [4] exploited an analogy between te image of intensity function I and te stream function in a D incompressible fluid. To see tis, we recall te vorticitystream formulation of te D NSE, ω + v ω = ν ω. (.3) t Here v = Ψ is te velocity, were Ψ is te stream function and ω = v is te vorticity. If te viscosity ν is zero, te steady state solution for te stream-function Ψ in D satisfies Ψ Ψ =. (.4) Te similarity between (.) and (.4) allows one to develop metods for te inpainting problem using tecniques from fluid dynamics. For te image inpainting problem, instead of simulating (.3), we compute te following numerically: ω + v ω = ν (g( ω ) ω), (.5) t were g accounts for te anisotropic diffusion (edge preserving diffusion) to sarpen te image. At eac time step, we solve te Poisson equation: I = ω, I Ω = I, (.6) were te boundary values I are derived from te values of I in D Ω. We recall tat oter PDEs ave been proposed recently (cf. [6, 34]) wic ave been sown to more efficient and peraps easier to implement. Here, our main interest is to explore different sub-grid scale turbulent models applied to image inpainting. Te idea is to test different sub-grid scale turbulence models instead of simulating te NSE (or te equation in (.5)) for purposes of finding te solution to te image inpainting problem

5 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 5 for reduced computational requirements. Modifying te function g appropriately in (.5) can also add to te accuracy of preserving te edges and ence give a more accurate image inpainting. In (.5) we noted tat for te image inpainting, te dissipation term in te NSE was modified to accommodate anisotropic diffusion. If g = 1, we recover te usual NSE wit isotropic diffusion. Tere as been extensive analysis of various types of anisotropic filtering, see for example [35]. In our numerical simulation of equation (.5) we ave used several different diffusivity functions, namely, g(s) = [ 1 + ( s k ) ] 1, g(s) = exp ( s k ), g(s) = [ 1 + ( s k )] 1, were k is a predefined diffusion parameter. For te particular test cases performed, we observed no significant difference in te numerical solution to te image inpainting problem. Original image wit noise Recovered image wit NSE: =, dt = y 15 y x x (a) (b) FIGURE. Tis is an example of image inpainting resulting from te steady state solution of te Navier-Stokes equations (wit anisotropic diffusion to preserve edges) wit viscosity ν =. Te time-step dt was set to.1 to guarantee convergence of te numerical solution. Just as an illustration of te idea, we give a simple numerical experiment. In Figure.1a, we took te poto of Alexander Ebraimi wit a noise mask as a test case. We run te Navier-Stokes equations (NSE) wit anisotropic diffusion until its steady state [4], and recovered te picture in Figure.1b. Te iterative inpainting process written in MATLAB took 4 iterations wit a tolerance of.1, in about 1.84 seconds on a laptop computer wit a Intel(R) 1.6 GHz CPU. Te grid for te inpainting region is pixels... Alpa Models and te Navier-Stokes-Voigt (NSV) Equations. Te NSV equations (1.1) were suggested as a regularization model for te 3D NSE, were te lengtscale α is considered as te regularizing parameter. Te system (1.1) was first introduced in 1973 by Oskolkov (see [3, 31]) as a model of a motion of linear, viscoelastic fluid. In tat setting, α is tougt of as a lengt-scale parameter caracterizing te elasticity of te fluid. As noted in [, 1], te addition of te term α u t regularizes te 3D NSE wic makes it globally well-posed [8, 3] and it canges te parabolic caracter of te original Navier-Stokes equations were in tis case one does not observe an immediate smooting of te solutions as usually expected in parabolic PDEs. Instead, te equations beave like a damped yperbolic system. Despite te damped yperbolic

6 6 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN beavior, it was sown in [] tat te solutions to te 3D NSV equations lying on te global attractor posses a dissipation range. Tis property supports te claim tat indeed, te NSV equations can be used as a sub-grid scale turbulence model. Tis type of inviscid regularization (tat is, a regularization tecnique witout introducing extra viscous or yperviscous terms) as been used for D surface quasi-geostropic (SQG) model ([3]), see also [3] for te Birkoff-Rott equation, induced by te D Euler-α equations for vortex seet dynamics and [5] for te 3D Magnetoydrodynamics (MHD) equations. Tere is an interesting literature beind te rediscovery of te NSV equation as a turbulence model. It can be traced back from te early study of 3D Navier-Stokes-α (NS-α) turbulence model in 1998 (also known as te viscous Camassa-Holm equations (VCHE) and Lagrangian averaged Navier-Stokes-α (LANS-α) model, wic can be written as [15, 16, 19, 8] t v ν v u v = p + f, u = v =, v = (I α )u, (.7) wit u(x, ) = u in (x). Te parameter α can be viewed as te lengt scale associated wit filter widt smooting v to obtain u, wit filter associated wit te Green s function (Bessel potential) of te Helmoltz operator (I α ). Te system is supplemented wit periodic boundary conditions in a box [, L] 3. Te inviscid and unforced version of te 3D NS-α was introduced in [19] based on te Hamilton variational principle subject to te incompressibility constraint div v =. By adding te viscous term ν v and te forcing f in an ad oc fasion, te autors in [9, 1, 11] and [16] obtain te NS-α system wic tey named, at te time, te viscous Camassa-Holm equations (VCHE), also known as te Lagrangian averaged Navier- Stokes-α model (LANS-α). In [9, 1, 11] it was found tat steady state solutions of te 3D NS-α model compared well wit averaged experimental data from turbulent flows in cannels and pipes for large Reynolds numbers. Tis led researcers [9, 1, 11] to suggest tat te NS-α model be used as a closure model for te Reynolds averaged equations. Since ten, an entire family of α - models as been found wic provide similar successful comparison wit empirical data among tese are te Clark-α model [7], te Leray-α model [13], te modified Leray-α model [] and te simplified Bardina model [8, 4] (see also [9] for a family of similar models). We focus our attention on te simplified Bardina model: t v ν v + (u )u = p + f, u = v =, v = u α u, (.8) wit initial condition u(x, ) = u in (x). Te equation above was introduced and studied in [4] supplemented wit periodic boundary conditions. Notice tat consistent wit all te oter alpa models, te above system is te Navier-Stokes system of equations wen α =, i.e. u = v. In [8], te viscous and inviscid simplified Bardina models were sown to be globally well-posed. It was also sown tat te viscous simplified Bardina model as a finite dimensional global attractor. In [8] it was observed tat te inviscid simplified Bardina model, is equivalent to te following modification of te 3D Euler

7 equations THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 7 α u t + u + (u )u + p = f, t u =, (.9) wit initial condition u(x, ) = u in. In particular, it is equal to te Euler equations wen α =. Following tis observation, Cao, et. al. [8] proposed te inviscid simplified Bardina model as regularization of te 3D Euler equations tat could be used for simulations of 3D inviscid flows. Inspired by te above model, Cao, et. al ten proposed te following regularization of te 3D Navier-Stokes equations α u t + u ν u + (u )u + p = f, t (.1) u =, wit initial condition u(x, ) = u in, and subject to eiter periodic boundary condition or te no-slip Diriclet boundary condition u Ω =. In te presence of pysical boundaries te above regularization (.1) of te Navier-Stokes equations, is different in nature from any of te oter alpa regularization models, because it does not require any additional boundary conditions. Tis newest addition to te family of α models was later on discovered to ave already existed in te literature known as te Navier-Stokes Voigt equations but as a model for viscoelastic fluids [3, 31]. As mentioned earlier, te α in tat setting represents te elasticity of te fluid..3. Te Alpa Model as Sub-Grid-Scale Turbulence Model. In addition to te remarkable matc of explicit analytical steady state solutions of te alpa models to te experimental data in cannels and pipes, te validity of te first alpa model, namely te NS-α model, as a sub-grid scale turbulence model was also tested numerically in [1, 8]. In te numerical simulation of te 3D NS-α model, te autors of [1, 17, 18, 8] sowed tat te large scale (to be more specific, tose scales of motion bigger tan te lengt scale α) features of a turbulent flow is captured. For scales of motion smaller tan te lengt scale α, te energy spectra decays faster in comparison to tat of NSE. Tis numerical observation as been justified analytically in [15]. In direct numerical simulation, te fast decay of te energy spectra for scales of motion smaller tan te supplied filter lengt represents reduced grid requirements in simulating a flow. Te numerical study of [1] gives te same results, wic also old for te Leray-α model in [13, 17]. In [6] it was observed tat te scaling exponent of te energy spectrum of te D Navier-Stokes-α model (NS-α), for wave numbers k suc tat kα 1, is k 7. A posteriori, it was observed tat tis scaling corresponding to tat predicted by assuming tat te dynamics for kα 1 was governed by te caracteristic time scale for flux of te conserved enstropy. Tis observation led researcers [6] to speculate tat (in general) te unknown scaling exponent for any α-model may be predicted by te dynamical time scales for te dominant conserved quantity for tat model in te regime kα 1. Tis speculation is verified to be correct in [7], were numerical simulations of te D Leray-α model were presented..4. Navier-Stokes-Voigt for Image Inpainting. Te main difference between NSV and oter alpa models is tat it needs no additional boundary conditions in te presence of pysical boundaries. Tis makes NSV a natural alternative to NSE for inpainting. In order to solve te inpainting problem we need to compute te solution to (.4). We approximate it wit te steady state solutions of (.5) wit viscosity ν >. As noted

8 8 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN in [4] te presence of viscosity is needed to ave a unique steady-state solution and its stability depends on ow big or small is te viscosity. It is well-known tat for very ig Reynolds number flows (ie very small viscosity ν), direct numerical simulation of te NSE, requires computational resources wic cannot be accessed easily. For purposes of direct numerical simulations, te NSV equations were proposed as a sub-grid scale turbulence model. Here we would like to investigate te effects of tis regularized equations to image inpainting. Notice tat te steady state solution to te NSV equations in (1.1) matces te steady state solution to te NSE. We expect tat by using te NSV, te convergence to te steady state solution will be at a muc faster rate reducing te computational cost. 3. NUMERICAL COMPARISON NSE AND NSV FOR IMAGE INPAINTING We simulate te inviscid D NSV wit anisotropic diffusion (artificial viscosity) α ω ω + + u ω = ν (g( ω ) ω), (3.1) t t were, ω = u. Notice tat wen α =, we recover (.5). Using a forward time up-wind finite-difference sceme [], we solve for ω (n+1) i,j in te discretized equation: (1 α )ω n+1 i,j νdt ( x (gωx n+1 ) + y (gωy n+1 )) = (1 α )ωi,j n [ + dt u n 1 sgn(u n 1) ωn i+ sgn(u n 1 ),j ωn i,j u n dx sgn(u n ) ωn i,j+ sgn(u n ) ωn i,j dy ]. (3.) Our numerical experiments can be summarized into tree categories: test for stability, test for efficiency and test for quality. In section 3.1, we present some test cases sowing tat for a given time stepping size, te gray-level blows up wen te inpainting is done using NSE wile te inpainting using NSV produced a stable solution. We give teoretical arguments in te appendix wic in some sense support tese numerical results. It is wortwile to mention tat using an explicit numerical sceme for bot te D NSE and D NSV equations, one can see a nice advantage of using NSV instead of te NSE in automating te inpainting problem in terms of CFL condition wen we for simplicity ignore te nonlinear term in te governing equations. In particular, te usual stiffness of te NSE by discretizing te linear part explicitly is no longer present for NSV. In te Appendix, we present rigorous argument for te stability condition wen te nonlinear term is taken into consideration. In section 3. we present some numerical studies using NSE and NSV for different values of α. We look at te effect of te parameter α on te quality of te resulting images and estimate te number of floating point operations (FLOPS) needed to converge to steady state. We perform tis numerical experiment for two different time-steps. Te quality of te resulting image is measured using PSNR wic will be defined in te next section Stability Beavior of NSE and NSV for Inpainting. In tis section we will present two test cases comparing te NSE and NSV bot wit anisotropic diffusion. Te results presented ere suggests tat for some fixed specified time-stepping size dt, te NSV produced a stable solution in te image inpainting iteration process wile te NSE wit anisotropic diffusion provided an unstable solution (gray-level blows up). A quick summary of results is provided in Table 1. For eac test cases, we fix te viscosity ν, tolerance tol, and te inpainting region Ω. In Figure 3a, te inpainting region was computed in a grid (in pixels) of 1 1. For dt =.1 and ν =, te NSE did not converge to te steady state solution due to large

9 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 9 TABLE 1. Comparison between NSE(α = ) and NSV. ν = viscosity coefficient, dt- size of time step, Ω (in pixels) = size of te inpainting region. Figure α ν dt Ω Recovered image 3a 1 1 Original image 3b.1 Does not converge 3c.5.1 Converges 3d.1 Converges 4a 64 1 Original image 4b.1 Does not converge 4c.9.1 Converges 4d.1 Converges time-stepping, as sown in Figure 3b. On te oter and, for te same dt =.1 and viscosity ν =, te steady state solution to NSV wit α =.5 gave a reasonable image inpainting as sown in Figure 3c. In Figure 3d, we sow te converged steady state solution of NSE for a smaller dt =.1. Using a subjective measure, Figures 3c and 3d are comparable in quality but Figure 3d required a larger time-step to converge to te steady state solution. y 4 6 Original image wit noise x 5 y No convergence wit NSE: =, =, dt= x 5 (a) (b) y Recovered image wit NSV: =.5, =, dt= x 5 y Recovered image wit NSE: =, =, dt= x 5 (c) (d) FIGURE 3. A comparison between te NSE and NSV. Figure a: Original image I wit te inpainting region. Figure 3b: For step size dt =.1 te NSE wit anisotropic diffusion does not converge to its steady state solution. Figure 3c: On te oter and te NSV wit anisotropic diffusion does converge for α =.5 and dt =.1. Figure d: For a smaller step size dt =.1, te NSE wit anisotropic diffusion does converge.

10 1 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN [t] In Figure 4a, we coose a sligtly different size of te inpainting region computed in a grid (in pixels) of For dt =.1 and ν =, te NSE did not converge to te steady state solution due to large time-stepping, as sown in Figure 4b. Similar to te first example, for te same dt =.1 and viscosity ν =, te steady state solution to te NSV model, wit α =.9, gave a reasonable result to te image inpainting problem as sown in Figure 4c. In Figure 4d, we sow te converged steady state solution of NSE for a smaller dt =.1. Figures 4c and 4d are comparable in quality but Figure 4d required a larger time-step to converge to te steady state solution. Original image wit noise 5 No convergence wit NSE: =, =, dt=.1 5 y 15 y x x (a) (b) y Recovered image wit NSV: =.9, =, dt= x y Recovered image wit NSE: =, =, dt= x (c) (d) FIGURE 4. A comparison between te NSE and NSV. Figure a: Original image I wit te inpainting region. Figure4b: For step size dt =.1 te NSE wit anisotropic diffusion does not converge to its steady state solution. Figure 4c: On te oter and te NSV wit anisotropic diffusion does converge for α =.9 and dt =.1. Figure 4d: For a smaller step size dt =.1, te NSE wit anisotropic diffusion does converge. 3.. Comparing te Efficiency of NSE and NSV for Inpainting. Te NSV converges to its steady state solution using a larger time-step or fewer number of iterations. However, it is not clear weter it is more efficient tan te NSE because of te extra computational cost, for te α term, for eac iteration. In addition, we also need to compare te quality of te resulting images. We do a FLOP counting for te two model equations until it reaced a steady state solution by using te Ligtspeed MATLAB toolbox. We approximated te FLOP by taking only in consideration te cost of matrix inversions. We use te FLOP inv function wic returns te number of FLOP necessary to invert an n n matrix. Note tat tis function returns te minimal number of FLOP for matrix inversion, as toug te best possible algoritm was used, no matter wic metod we

11 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 11 are actually using. Tis sould not be an issue for we use te same matrix inversion tecnique wen running te two model equations. Given tat Ω is N M, we take n = max(n, M). We ten normalize te FLOP count by dividing it by te number of pixels in Ω. Tat is, te FLOP count is in units of number of FLOP per pixel. In addition, we also compute te PSNR of te resulting images to give some measure of comparison for te quality of te resulting images in addition to subjective measure. Let P (i, j) be te original image tat contains N by M pixels and I(i, j) te reconstructed image. Te pixel values range between black () and wite (55). Te PSNR in decibels (db) are computed as follows: i ( ) 55 j (P (i, j) I(i, j)) P SNR = log 1, were RMSE =. RMSE N M Comparison between two PSNR values for different reconstructed image gives one mea- TABLE. Comparison between NSE(α = ) and NSV. ν = viscosity coefficient, dt- size of time step, FLOP count = FLOP count/(n*m), PSNR = Peak signalto-noise ratio. FLOPS α dt =.1 dt = e4 9.e3 1/ e e3 /3 4.15e3 3.75e e e3 4/ e e3 sure of quality. Reconstructed images wit iger PSNR are considered better. (However, as frequently mentioned in te literature, PSNR may not equate wit uman perception; we will view PSNR as a reasonable metric ere.) In any case, we now describe of our results based on observations in Table and Figure 5. In tis numerical experiment, we fix te value of ν =. (1) Comparing NSE (α = ) and NSV wit α = 1/3 for bot time steps dt =.1 and dt =.1 in Table 5, we notice te number of FLOPS needed by NSV to converge to steady state as been reduced to almost alf of te FLOPS required by te NSE to converge to its steady state solution. () In Figure 5a, for smaller dt =.1, te PSNR of NSV (α = 1/) is 45 db and is 1 point iger tan tat of NSE. In te case wen dt =.1, (Figure 5b) te PSNR of NSE is 46.5 db wile te NSV for α = 1/3 as a PSNR value of 46. (3) We also notice tat increasing te time step from dt =.1 to dt =.1 for te NSE canges te quality of te picture dramatically under subjective measure. In Figure 5, te PSNR for NSE is reduced from 46.5 db to 44 db wen dt is increased from.1 to.1. In te case of NSV for α = 1/3, te cange in quality of te picture under subjective measure is unnoticeable. In particular, te PSNR is canged from 46 db to 45 db as te dt is increased. In summary, one can sow tat for some value of α >, NSV gives an inpainting wic is comparable to NSE but requiring less resources. It is ten natural to ask if one can determine te value of α a priori wic will yield te smallest FLOP count wit maximum PSNR. Our preliminary studies suggest tat coice of α may be image dependent. For tis reason, one may prefer te NSE to avoid finding te optimal α. Tis suggests te development of scemes to adaptively determine optimized α.

12 1 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN PSNR dt =.1 a. NSE ( =) b. NSV ( =1/3, =) c. NSV ( =/3, =) d. NSV ( =1, =) e. NSV ( =4/3, =) f. NSV ( =5/3, =) g. NSV ( =6, =) no. of frames (iterations) (a) PSNR dt =.1 a. NSE ( =) b. NSV ( =1/3, =) c. NSV ( =/3, =) d. NSV ( =1, =) e. NSV ( =4/3, =) f. NSV ( =5/3, =) g. NSV ( =6, =) no. of frames (iterations) (b) FIGURE 5. PSNR of NSE vs NSV (for various values of α). Wen dt =.1, te PSNR of NSV (α=1/3) is iger tan te NSE and requires only about alf of te number of FLOPS needed for te NSE to converge to its steady state solution. Recovered Image = Recovered Image = 1/ y 15 1 y x x 4 (a) (b) FIGURE 6. An example of image recovered for (a) NSE and (b) NSV (α = 1/3) wit dt = SOME SUPPORTING THEORETICAL RESULTS In tis section we summarize some existing results, and ten establis some new results tat support te numerical results in te previous section. To tis end, let Ω = [, πl]. Te NSE of viscous incompressible flows, subject to periodic boundary condition on domain Ω, is written in te form: t u ν u + (u )u = p + f, (4.1) u =,

13 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 13 wit initial condition u(x, ) = u in (x), were u represents te unknown fluid velocity, p is te unknown pressure scalar, ν > is te constant kinematic viscosity, Te given forcing function f is assumed ere to be time independent and wit mean zero: Ω f(x)dx =. Te initial velocity uin is also assumed to ave zero mean, ence also u. Below we use will use standard notation in te context of te matematical teory of Navier-Stokes equations (NSE) (see, e.g., [14, 3, 33]). In particular, (1) We denote by L p and H m te usual Lebesgue and Sobolev spaces, respectively. And we denote by and (, ) te L norm and L inner product, respectively. () Let F be te set of all vector trigonometric polynomials wit periodic domain Ω. We ten set { } V = φ F : φ = and φ(x) dx =. Ω We set H and V to be te closures of V in L and H 1, respectively. We note by Rellic lemma (see, e.g., [1]) we ave te V is compactly embedded in H. (3) We denote by P σ : L H te Helmoltz-Leray ortogonal projection operator, and by A = P σ te Stokes operator subject to periodic boundary condition wit domain D(A) = (H (Ω)) V. We note tat in te space-periodic case, Au = P σ u = u, for all u D(A), wit A 1 self-adjoint positive definite and compact from H into H (cf. [14, 33]). Denote < L = λ 1 λ te eigenvalues of A, repeated according to teir multiplicities. (4) We recall te following two-dimensional Ladyzenskaya inequality: φ L 4 c φ 1/ L φ 1/ H 1, for every φ H 1 (Ω). (4.) Hereafter c will denote a generic dimensionless constant. (5) For w 1, w V, we define te bilinear form B(w 1, w ) = P σ ((w 1 )w ). (4.3) (6) By te Sobolev inequalities and compactness teorems, in two dimensions (or any dimensions less tan 4) we can define on V a trilinear continuous form b by setting b(u, v, w) = (B(u, v), w). (4.4) If u, v V ten b(u, v, v) =. (4.5) A special case we will need to establis for stability of D NSV is: b(u, u, Au) =, for all u D(A). (4.6) We finis by mentioning some useful inequalities for semi-discretizations. Let n be te spatial dimension. Te following are easily establised (cf. [3, 33]): { n ( u 1 = u i x + ) ( j u i x ) } j dx i,j=1 j Ω ( n ) (4.7) 1 4 u = S() u j=1 j

14 14 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN u = 1 n { 1 ( u 4 i x + ) ( j + u i x ) } j u i (x) dx i,j=1 j Ω ( n ) (4.8) 1 4 u = S() u j=1 j 4.1. Uniqueness of Steady State Solutions for Boundary Driven Flows. Te results in [4] is obtained by requiring te user to tune te parameters and related data in te inpainting region to suit te particular problem at and. In [] wile inpainting various sample images, wit similar inpainting region, viscosity ν and δt, it was found tat certain set of parameters produced stable solutions in some images and unstable solutions (gray-level blows up) in oters. Tey ave noted tat certain caracteristics of I near Ω as some effect on te maximum allowable stable coice of δt. In tis section we will present analytical arguments on te dependence of te solution on te image at te boundary, size of te inpainting region and viscosity. We give some ypotesis on ow te related data affects te convergence of te numerical solution. Note tat te steady state solution for te NSV is exactly equal to te steady state solution of NSE. Te main goal for tis study is to determine te relationsip between te viscosity, te image at te boundary and te size of te inpainting region. For example, we would like to determine for wic size of te inpainting region and norm of I on te boundary, do we get a unique steady state solution. In [4], for te Navier-Stokes based inpainting, a discussion on te uniqueness of steady state solution and its relevance to inpainting was presented. It was expected tat Navier- Stokes based inpainting may inerit some of te stability and uniqueness issues known for incompressible fluids, altoug te effect of anisotropic diffusion is still unclear. Te dependence of uniqueness in te viscosity of te fluid is discussed in [4]. In tis section, we present some rigorous arguments following te work in [3, 33]), modified sligtly to interpret it in te context of image inpainting. In tis section we sow te dependence of te uniqueness of steady solution on te viscosity, on te image at te boundary, and on te size of te inpainting region. We start ere by recalling some notation. Te notation used ere are similar to tose used in Section.. Let us denote by Ω a bounded domain of R of class C wic is filled wit an incompressible viscous fluid. For te particular application of NSE and NSV in image inpainting, we supplement it wit Diriclet boundary condition u = φ on Ω, in wic φ is independent of time. We consider ere te non-omogeneous steady state Navier-Stokes problem wic coincides wit te steady state Navier-Stokes-Voigt problem. Find u and p suc tat ν u + (u )u + p =, in Ω, (4.9) div u =, in Ω, (4.1) u = φ on Ω. (4.11) We assume tat φ is given as te trace on Ω of a function Φ, Φ H (Ω), div Φ =, Φ n ds =, (4.1) were n is te unit outward normal in Ω. Te idea is te following: Given te pysical data φ defined on Ω, we find an extension Φ of φ inside Ω satisfying (4.1). Under te above ypotesis tere exist at least one u H 1 and distribution p on Ω satisfying (4.9) (see [3, 33]) provided we coose te extension Φ L a (Ω) of φ suc tat Φ a is Ω

15 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 15 sufficiently small for a > so tat b(v, Φ, v) ν v, for all v V. (4.13) Te construction of Φ satisfying (4.1) and (4.13) is presented in [33]. Knowing Φ and knowing u is a solution of (4.9) ten letting û = u Φ, equation (4.9) is equivalent to ν û + (û )û + (û )Φ + (Φ )û + p = ν Φ (Φ )Φ, in Ω, (4.14) We now state a uniqueness result. div û =, in Ω, (4.15) û = Teorem 4.1. Suppose tat te norm of Φ in L a (Ω) is sufficiently small so tat and ν is sufficiently large so tat on Ω.(4.16) b(v, Φ, v) ν v, for all v V, (4.17) ν > Cλ 1/ 1 f, (4.18) were f = ν Φ B(Φ, Φ), λ 1 is te smallest eigenvalue of te Stokes operator, and C is constant, ten (4.9) as a unique solution. Proof. See [3, 33]. We can recast te teorem above, in te context of image inpainting, as te dependence of te uniqueness of steady solution on te viscosity ν, on te image at te boundary, wic is related to Φ (since te stream function is related to te velocity field), and on te size of te inpainting region wic is related to λ 1/. 4.. Stability Analysis of te Sceme based on D NSE. In tis section we are concerned wit a discussion on te discretization of te Navier Stokes equations in two dimensions, subject to periodic boundary conditions, wit basic domain Ω = [, πl]. We study ere a full discretization of te equations, bot in space and time. We start wit a description of te approximating sceme and ten proceed to te study of te stability of tis sceme. Our study ere is based on te energy metods similar as in [33] wic leads to sufficient conditions for stability. We now begin by defining a Galerkin approximation of te separable normed space V. For reference we direct te reader to [33]. Let V, N, be an increasing sequence of finite-dimensional subspaces of V wose union is dense in V. For simplicity, assume tat V L (Ω), for all N. (4.19) Te space V is terefore equipped wit two norms: te norm induced by L (Ω) and its own norm. Since V is finite dimensional ten tese two norms must be equivalent. To be more precise we ave wit d independent of, u d u, for all u V, (4.) u S() u, for all u V. (4.1) Similarly, let D(A), N, be an increasing sequence of finite-dimensional subspaces of D(A) wose union is dense in D(A). For simplicity, assume tat D(A) V L (Ω), for all N. (4.) Te space D(A) is terefore equipped wit two norms: te norm induced by V and its own norm A. Since D(A) is finite dimensional normed space ten tese two

16 16 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN norms must be equivalent. To be more precise we ave wit d independent of, u d A u, for all u D(A), (4.3) A u cs() u, for all u D(A). (4.4) Te constant S(), wic depends on is sometimes called te stability constant since it plays a major important role in obtaining necessary conditions on te stability of numerical approximations. Usually S(), as. Let tere be given a trilinear continuous form on V, say b (u, v, w ) wic satisfies te following: (1) For all u, v V, bot of te following old: were at least, b (u, v, v ) =, (4.5) b (u, u, v ) d 1 u v d 1 S () u u v () For all u, v, w V, S 1 () u u v, (4.6) S 1 () d 1 S (). (4.7) b (u, v, w ) d 1 u v w. (4.8) We divide te interval [, T ] into N intervals of equal lengt k = T/N. We associate wit k and te function f, te elements f 1,..., f N : f m = 1 k mk (m 1)k f(t) dt, m = 1,..., N; f m L (Ω). (4.9) We denote by u te ortogonal projection of te initial condition u onto D(A) in L (Ω). Wen u,..., um 1, are known, u m is te solution in D(A) of 1 k (um u m 1, v ) + ν((u m 1, v )) + b (u m 1, u m 1, v ) = (f m, v ) (4.3) for all v D(A). Lemma 4.. We assume tat k and satisfy (1) ks () 1 δ for some δ, < δ < 1, 4ν () ks () 1, (3) ks1()s () νδ, 8d d 5 were d is as in (4.) and ( ) d d 5 = u + d + 1 δ T f(s) ds, (4.31) ten, te u m νd given by (4.3) remain bounded in te following sense u m d 5 m = 1,..., N, (4.3) k A u m 1 d 5 νδ, (4.33) N ( ) δ T u m u m 1 d f(s) ds. (4.34) δ

17 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 17 Proof. We replace by v by A u m 1 in (4.3); due to te identity (a b, b) = a b a b, (4.35) we find u m u m 1 u m u m 1 + kν A u m 1 = k(f m, A u m 1 ) kd f m A u m 1 (4.36) kν A u m 1 + k d ν f m. We would like to majorize u m um 1 in (4.36). Let v = A u m A u m 1 in (4.3). Tis gives u m u m 1 = kν ( A u m 1, A(u m u m 1 ) ) kb (u m 1, u m 1, A u m u m 1 ) (4.37) + k(f m, A (u m u m 1 )) =: I 1 + I + I 3. We successively majorize I 1, I and I 3 using (4.1), (4.4), (4.6), and Caucy-Scwarz inequality I 1 kν A u m 1 A (u m u m 1 ) (4.38) kνs() A u m 1 u m u m 1 (4.39) 1 4 um u m 1 + 4k ν S () A u m 1, (4.4) I ks 1 () u m 1 u m 1 ks()s 1 () u m um u m 1 A (u m u m 1 ) (4.41) u m u m 1 (4.4) u m 1 + 4k d S ()S1() u m 1 u m 1 (4.43) 1 4 um u m 1 + 4k d S ()S1() u m 1 A u m 1, (4.44) I 3 k f m A (u m u m 1 ) ks() f m u m u m 1 (4.45) 1 4 um u m 1 + 4k S () f m. (4.46) Denoting as Θ = 4k d S 1()S () u m 1 4k ν S ( ) A u m 1 u m u m 1 Going back to (4.36) we ave u m u m 1 A u m 1, we ave tat (4.37) becomes + Θ kν(1 δ) A u m 1 + Θ. (4.47) + k(νδ 4kd S1()S () u m 1 ( ) d k ν + 4kS () f m ( 1 kd ν + 1 δ ) f m νd ( ) d kd + 1 δ f m. νd ) A u m 1 (4.48)

18 18 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN Summing up (4.48) for m = 1 to r we get ( u r + k νδ 4kd S1()S () u m 1 ) A u m 1 µ r, (4.49) were ( ) d µ r = u + kd + 1 δ f m. Now using condition (iii) in Lemma 4., we will prove by te metod of induction tat u r + kνδ A u m 1 µ r, r = 1 to N. (4.5) First observe tat ( ) d µ r µ N = u + kd + 1 δ N f m νd m 1 ( ) d u + d (4.51) + 1 δ T f(s) ds =: d 5 To establis te basis for induction (r = 1) we write (4.48) for m = 1 and use (iii) in Lemma 4. to get νd νd u 1 + kνδ A u u + 4k d S1()S () u A u ( ) d + kd + 1 δ f 1 µ 1 + kνδ A u. Tus equation (4.5) for r = 1 is satisfied. By induction on r, assume tat (4.5) olds up to te order r 1. Note tat by te recurrence ypotesis νd (4.5) u r 1 µ r 1 µ N d 5. (4.53) Tus by (4.49), we ave u r + kνδ A u m 1 µ r + 4k d S1()S () u m 1 Hence, u r + kνδ µ r + 4k d S 1()S ()d 5 µ + kνδ A u m 1. A u m 1 (4.54) A u m 1 µ r. (4.55)

19 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 19 Tis gives (4.33). It remains to prove (4.34). From (4.47), applying (ii) and (iii) from Lemma 4., we get u m u m 1 k ν S ( ) A u m 1 + 4k d S1()S () u m 1 A u m 1 + 4k S () f m kν(1 δ) A u m 1 + 4k νδ u m 1 A u m 1 + 4k f m 8d 5 kν(1 δ) A u m 1 + kνδ A u m 1 + 4k f m kν( δ) A u m 1 + 4k f m. (4.56) Summing it up over all m = 1 to N and using (4.33) we find (4.34) Stability Analysis of te Sceme based on D NSV. Te preliminary setup is te same as tose of D NSE in te previous section. Wen u,..., um 1, are known, u m is te solution in D(A) of for all v V. 1 k (um u m 1, v ) + α k ( A u m A u m 1 Lemma 4.3. We assume tat k and satisfy, v ) + ν((u m 1, v )) +b (u m 1, u m 1, v ) = (f m, v ), (4.57) (1) ks () 1 δ for some δ, < δ < 1, 4ν () ks1() νδ, 8d 6 ( ) d d 6 = u + α u T + ν + 4T f(s) ds, (4.58) and d is as in (4.), ten, te u m given by (4.57) remain bounded in te following sense u m + α u m d 6 m = 1,..., N, (4.59) k u m 1 d 6 νδ, (4.6) N ( ) ( ) u m δ T u m 1 + α u m u m 1 d 6 + 4T f(s) ds. (4.61) δ Proof. We replace by v by u m 1 Scwarz, we get u m u m 1 in (4.57); due to (4.), (4.1),(4.35) and Caucy- u m u m 1 + kν u m 1 + α ( u m u m 1 u m u m 1 = k(f m, A u m 1 ) kd f m u m 1 kν u m 1 + k d ν f m, ) (4.6)

20 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN tat is, ( u m u m 1 ) + α ( u m u m 1 ( u m u m 1 ) + α u m u m 1 ) (4.63) k d ν f m. We would like to majorize te term ( u m um 1 + α u m um 1 ) in (4.63). Let v = u m um 1 in (4.57). Tis gives u m u m 1 + α u m u m 1 = kν (( )) u m 1, u m u m 1 (4.64) kb (u m 1, u m 1, u m u m 1 ) + k(f m, u m u m 1 ) =: I 1 + I + I 3. We successively majorize I 1, I and I 3 using repeatedly (4.1), (4.6), and Caucy- Scwarz inequality and get exactly te estimates as in page 34 of [33]. Terefore (4.64) becomes u m u m 1 + α u m u m 1 4k ν S ( ) u m 1 + 4k S1() u m 1 u m 1 + 4k f m (4.65) kν(1 δ) u m 1 + 4k S1() u m 1 u m 1 + 4k f m. Going back to (4.63) we ave u m u m 1 + α ( ) u m u m 1 ( d k k + k(νδ 4kS1() u m 1 ) ν + 4k f m ( d ν + 4T ) u m 1 ) f m (since k T ). (4.66) Summing up (4.65) for m = 1 to r we get ( u r + α u r + k νδ 4kS 1 () u m 1 ) u m 1 µ r, (4.67) were ( ) d µ r = u + α u + k ν + 4T f m. Now using condition (ii) in Lemma 4.3, we will prove by te metod of induction tat u r + α u r + kνδ u m 1 µ r, r = 1 to N. (4.68) First observe tat µ r µ N = u + α u + k ( ) d N ν + 4T f m =: d 6. (4.69) m 1

21 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 1 To establis te basis for induction (r = 1) we write (4.66) for m = 1 and use (ii) in Lemma 4.3 to get u 1 + α u 1 + kνδ u u + α u, ( ) d + 4k S1() u u + k ν + 4T f 1 µ 1 + kνδ (4.7) u. Tus equation (4.68) for r = 1 is satisfied. By induction on r, assume tat (4.68) olds up to te order r 1. Note tat by te recurrence ypotesis u r 1 + α u r 1 µ r 1 µ N d 6. (4.71) Tus by (4.67), we ave Hence, u r 1 + α u r + kνδ u m 1 u r + α u r + kνδ µ r + 4k S1() u m 1 u m 1 µ r + 4k S 1()d 6 µ + kνδ u m 1. u m 1 (4.7) u m 1 µ r. (4.73) Tis gives us (4.6). It remains to prove (4.61). From (4.65) and applying condition (i) and (ii) from Lemma 4.3, we get u m u m 1 + α u m u m 1 kν(1 δ) u m 1 + 4k S1() u m 1 u m 1 + 4k f m ( kν 1 δ ) u m 1 + 4kT f m. (4.74) Summing it up over all m = 1 to N and using (4.6) we find (4.61) Stability Analysis of a Semi-Implicit Sceme based on D NSE. We present ere general results for stability analysis for te full NSE equations vorticity formulation (implicit only on te linear part of te operator). Te setup is as follows: Wen u,..., um 1, are known, u m is te solution in D(A) of 1 k (um u m 1, v ) + ν((u m, v )) + b (u m 1, u m 1, v ) = (f m, v ), (4.75) for all v V. Lemma 4.4. We assume tat k and satisfy (1) ks 4 () d, () ks 1()S () d, were d = and d =, ten, te u m given by (4.75) remain bounded in te following sense u m d 7, m =,..., N, N N (4.76) u m u m 1 d 7, k A u m d 7,

22 M.A. EBRAHIMI, M. HOLST, AND E. LUNASIN were d 7 is some constant depending only on te data d and d. Proof. We let v = A u m in (4.75). Using te identity (4.35), we get u m u m 1 + u m u m 1 + kν A u m ( ) = kb u m 1, u m 1, A u m + k(f m, A u m ) (4.77) := I 1 + I. ( ) We would like bound te terms I 1 and I. Te identity b u m 1, u m 1, Au m 1 = togeter wit (4.) and (4.6) gives and Hence, I 1 ks 1 () u m 1 u m 1 ks 1 ()S() u m 1 k S 1()S () u m 1 A(u m u m 1 ) u m 1 u m 1 u m u m um u m 1, (4.78) ( ) d I kd f m A u m k ν f m + ν Au m. (4.79) u m u m um u m 1 + kν A u m k S1()S () u m 1 u m 1 kd ν f m. We sum up tese inequalities for m = 1 to r, to get were, u r + 1 λ r = u = kd ν u m u m 1 + kν k S 1()S () m= A u m u m 1 u m 1 λ r, (4.8) (4.81) f m + k S1()S () u u. (4.8) We assume tat kd d S1()S ()λ N ν δ, (4.83) for some fixed δ, δ ν. If tis olds ten one can sow recursively tat u r + 1 u m u m 1 + kδ A u m λ r, r = 1,..., N. (4.84) m 1 Clearly, letting m = 1 in (4.8) sows tat (4.84) is true for r = 1. Let us assume tat (4.84) is valid up to te order r 1, we want to sow tat (4.84) is valid for r. Observe

23 THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING 3 tat by te inductive ypotesis u m λ m λ N. Hence, by te condition (4.83) k S1()S () u m 1 u m 1 d k S1()S () u m 1 A u m 1 m= m= d d k S 1()S ()λ N k(ν δ) A u m 1. m 1 A u m 1 (4.85) We apply tis upper bound into (4.81), we get (4.84) for te integer r. To complete te proof it suffices to sow tat conditions (i) and (ii) in Lemma 4.4 ensure te condition (4.83). We recall tat since u u for all and N T k f m f(s) ds, ten, T λ N u + d f(s) ds + k S ν 1()S () u u d 1 + k S1()S () u u. Hence, if ks 1S () d, ten (4.86) kd d S 1()S ()λ N d (d 1 + d d u u ) ν δ, (4.87) provided d, d are sufficiently small. 5. SUMMARY Te NSV model of viscoelastic incompressible fluid as been proposed as a regularization of te 3D NSE for purposes of direct numerical simulation. In tis work, we ave sown one of te benefits of using te D NSV turbulence model for small regularization parameter α, instead of te D NSE to reduce computational expense wen automating te inpainting process. To be more precise, one can find a parameter α > in wic te D NSV gives a solution to te image inpainting problem comparable (bot using subjective and objective measure) to tat of te solution of te D NSE but only requires a time step muc larger in comparison to tat of D NSE. Tat is, te D NSV converge to te steady state solution wit a muc larger time step and ence solves te image inpainting problem using less computational resources. In te numerical experiments, we found tat after accounting for te relative costs of te two metods, te D NSV gives a solution to te inpainting problem wic matces te quality of te image produced by wen NSE is used but using less resources. In future work we would like to investigate oter PDEs wic can be used instead of te D NSE and D NSV wen solving te inpainting process. In particular, we would like to use a PDE to solve te inpainting problem witout te addition of anisotropic diffusion. Te dependence of te stability of solutions on certain caracteristics of te image near te boundary is also of major interest in tis topic. We would like to investigate te dependence of te uniqueness of steady solution on te viscosity, on te image at te boundary, and on te size of te inpainting region. We ave presented some basic results on te uniqueness of steady state solution of te D steady state applied in