ENSEEIHT, 2, rue Camichel, Toulouse, France. CERFACS, 42, avenue Gaspard Coriolis, Toulouse, France
|
|
- Avice Cole
- 5 years ago
- Views:
Transcription
1 Adaptive Observations And Multilevel Optimization In Data Assimilation by S. Gratton, M. M. Rincon-Camacho and Ph. L. Toint Report NAXYS May 203 ENSEEIHT, 2, rue Camichel, 3000 Toulouse, France CERFACS, 42, avenue Gaspard Coriolis, 3057 Toulouse, France University o Namur, 6, rue de Bruxelles, B5000 Namur (Belgium)
2 Adaptive Observations And Multilevel Optimization In Data Assimilation Serge Gratton, M. Monserrat Rincon-Camacho, and Philippe L. Toint May 2, 203 Abstract We propose to use a decomposition o large-scale incremental our dimensional (4D-Var) data assimilation problems in order to make their numerical solution more eicient. This decomposition is based on exploiting an adaptive hierarchy o the observations. Starting with a low-cardinality set and the solution o its corresponding optimization problem, observations are adaptively added based on a posteriori error estimates. The particular structure o the sequence o associated linear systems allows the use o a variant o the conjugate gradient algorithm which eectively exploits the act that the number o observations is smaller than the size o the vector state in the 4D-Var model. The method proposed is justiied by deriving the relevant error estimates at dierent levels o the hierarchy and a practical computational technique is then derived. The new algorithm is tested on a D-wave equation and on the Lorenz-96 system, the latter one being o special interest because o its similarity with Numerical Weather Prediction (NWP) systems. Keywords: Data assimilation, adaptive observations, numerical algorithms, multilevel optimization, a posteriori errors. Introduction In data assimilation, a substantial body o research has been conducted in the ield o adaptive observations. The main idea in this area is to adapt the set o observations used in the model calibration to improve computational eiciency while not sacriicing accuracy. Several techniques have been considered, depending on the practical context o the relevant application. In a Kalman ilter ramework, Leutbecher (2003b) introduces a Hessian reduced rank estimate whose aim is to predict expected change in the square norm o orecast error due to intermittent modiications o the observation network; its application to the Lorenz96 model (see Lorenz, 995) is presented in Leutbecher (2003a). Lorenz and Emanuel (998) present dierent strategies or weather orecasting which target the locations where the irst-guess errors are largest, and veriy i the including additional data improve the analyses and orecasts. Berliner, Lu and Snyder (999) introduce rigorous mathematical criteria to select uture data in a dynamical process, based on the uncertainty o the current analysis, the dynamics o error evolution, the orm and errors o observations and data assimilation. In Hansen and Smith (2000) the perormance o various observation targeting strategies is compared or both direct insertion and ensemble Kalman ilter assimilation. The authors conclude that the optimal adaptive observation strategies depend on the combination o model, assimilation scheme and observational network. Daescu and Navon (2004) describe a method based on a periodic update o the adjoint sensitivity ield which takes into account the ENSEEIHT, 2, rue Camichel, 3000 Toulouse, France and CERFACS, 42, avenue Gaspard Coriolis, 3057 Toulouse, France. (serge.gratton@enseeiht.r). CERFACS, 42, avenue Gaspard Coriolis, 3057 Toulouse, France. (monserratrc@ceracs.r). Namur Center or Complex Systems, University o Namur, 6, rue de Bruxelles, B-5000 Namur, Belgium. (philippe.toint@unamur.be).
3 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 2 interaction between time-distributed adaptive and routine observations. Cardinali, Pezzulli and Andersson (2004) introduce an inluence matrix whose purpose is to identiy inluential data in a 4D-Var ramework. They observe that low-inluence data occur in areas with a large number o observations while high-inluence data points occur in regions with sparse ones or in dynamically active areas. Finally, Liu, Kalnay, Miyoshi and Cardinali (2009) describe, in a Kalman ilter setting, the impact o the sel-sensitivity and the cross sensitivity (o-diagonal elements o inluence matrix). Our motivation or the present paper lies in several o the conclusions drawn in this literature. In particular, our work aims at a good numerical exploitation o the conclusion that the spatial distribution o observations may be inluential. More speciically, we introduce new techniques or adaptively selecting observations in a 4D-Var setting, based on mathematically sound a posteriori error estimates within a dynamic hierarchy. The technique used or constructing the adaptive set o observations is inspired by adaptive inite elements techniques studied in Section 2.2 (p. 00) o Rincon-Camacho (20). Related techniques or adaptive multigrid can be ound in Brandt (973) and McCormick (984). The paper is structured as ollows: Section 2 introduces the 4D-Var vocabulary and the considered hierarchy o observations. Section 3 then covers the associated error estimates and Section 4 presents the new adaptive algorithm. Section 5 discusses its application to a one-dimensional nonlinear wave equation and to the Lorenz96 model. Conclusions and perspectives are discussed in Section 6. 2 The problem and associated observation hierarchy Consider the nonlinear least-squares problems in 4D-Var data assimilation, whose objective is to ind an initial vector state at an initial time denoted as x = x(t 0 ) IR n. The structure o the problem is as ollows: min 2 x x b 2 x IR n B + 2 N t j=0 H j (x(t j )) y j 2 R j (2.) where the squared norm x 2 M is induced by the inner product xt Mx or a symmetric positive deinite matrix M IR l l and a vector x IR l. Here x b IR n is the background vector, which is an a priori estimate. The vector y j IR mj is the vector o observations at time t j and H j is the operator modeling the observation process at the same time. The state vector x(t j ) satisies the nonlinear model o evolution x(t j ) = M 0 j [x(t 0 )]. The matrix B IR n n is a symmetric positive deinite matrix representing the background-error covariance and the matrix R j IR mj mj is also a symmetric positive deinite matrix representing the observation-error covariance at time t j. The nonlinear least-squares problem (2.) is solved iteratively, and, at iteration k, the problem is simpliied by linearizing the nonlinear observation operator H j (x(t j )) at the iterate x k (which or the moment we denote only by x), leading to the optimization problem min 2 x x b + δx 2 δx IR n B + 2 H δx d 2 R (2.2) where R = diag(r 0, R,..., R Nt ) IR m m, d IR m denotes the concatenated misits over time y j H j (x(t j )) and H IR m n is the concatenated version o the linearized observation process d 0 H 0 d d =. H M 0 IRm, H =. IRn m. d Nt H Nt M 0 Nt Here H j and M 0 j are a (possibly approximate) linearization o the observation operator H j and the model M 0 j around x = x k (t 0 ) respectively.
4 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 3 Diverse techniques or solving problem (2.) have been proposed, see or instance Courtier, Thépaut and Hollingsworth (994) or the so called incremental method which is equivalent to applying a truncated Gauss-Newton iteration to problem (2.) (see Gratton, Lawless and Nichols, 2007). The general algorithm is the ollowing: Algorithm 2. Incremental 4D-Var. Initialize x 0 = x b IR n and set k = Compute x k (t j ) rom x k (t 0 ) = x k by running the nonlinear model M rom time t 0 to t Nt. 3. Calculate the vectors d k,j = y j H j (x k (t j )) or j =,..., N t. 4. Find an approximate solution δx k o the minimization problem (2.2), where x = x k, H = H k and d = d k. 5. Update the current solution as x k+ = x k + δx k. 6. Set k := k +. I convergence is not achieved return to Step 2. This algorithm is known as the outer loop and the minimization step (Step 4) as the inner loop. Our interest is to reduce the cost o this optimization problem. Termination criteria or the outer and inner loops o this algorithm are discussed in Gratton et al. (2007). It is most natural to solve the subproblem in Step 4 o Algorithm 2. directly in the space o dimension n, the size o x. This technique is reerred to as the primal approach. At variance, the quadratic optimization problem (2.2) may be rewritten in a space o dimension m, the number o observations. This is known as the dual approach and it is especially useul when m is signiicantly smaller than n, in which case the second term in (2.2) is o relatively low rank compared to the irst. A irst approach o this type is the Physical-space Statistical Analysis System (PSAS) method (see Courtier, 997), where the low rank observation term in (2.2) is handled by using a Sherman-Morrison ormula and applying the standard conjugate-gradient algorithm (see Hestenes and Stieel, 952, or Golub and Van Loan, 996, Section 0.2, p. 520) to the reduced system. This method has the drawback o not maintaining the inherent monotonicity o the conjugate-gradient algorithm in IR n, thereby making any stopping rule o the inner minimization diicult to deine (see El Akkroui, Gauthier, Pellerin and Buis, 2008, or Gratton, Gürol and Toint, 203). Fortunately, a better alternative is available, where the conjugate-gradient algorithm is reormulated in IR m using the inner product deined by the metric HBH T in order to guarantee the desired monotonicity properties without incurring additional cost. This technique, known as the Restricted Preconditioned Conjugate Gradient method (RPCG, see Gratton and Tshimanga, 2009 and Gratton et al., 203), provides an eicient numerical procedure where substantial computing gains are obtained when m n. We reer the reader to the cited publications or details. Our aim in this work is to make the best possible use o this dual technique and to propose an adaptive observations strategy or solving the data assimilation problem (2.). Suppose that we have a large set o m observations O which can be decomposed into a hierarchical collection o sets {O i } r i=0, each with cardinality m i, such that O i O i+ and m i < m i+ (i = 0,..., r ) where, by convention, O r = O and m r = m. To each set o observations O i we associate a misit vector y i IR mi (by selecting the relevant components in the vector y), and the corresponding observation-error covariance matrix R i. Given {O i } r i=0, we may thereore consider the hierarchical collection o minimization problems min 2 x x b 2 x IR n B + 2 H i (x) y i 2, i = 0,..., r, (2.3) R i where the vector H i (x) denotes the nonlinear observation operator concatenated over time. Thus our objective is to construct the collection {O i } r i=0 such that the sequential solution o the problems (2.3) is obtained at signiicantly lower cost compared to solving (2.2) directly, while at the same time maintaining equivalent accuracy requirements.
5 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 4 3 Mathematical analysis Our adaptive method is based on the exploitation o a posteriori bounds or the error between the solutions obtained using ew observations or many. In particular, we are interested in comparing the accuracy o the solutions o problems (2.3) or O i and O i+. For simplicity o notation, we denote them as O c a set with m c observations and O a set containing m observations such that m c < m and O c O, where the indices c and stand or coarse and ine, respectively. The ine optimization problem (2.3) is given, or the starting point x, by min 2 x + δx x b 2 δx IR n B + 2 H δx d 2 R (3.4) where d = y H (x) and H is the linearized version o the observation process at x involving the set o observations O. We may reormulate (3.4) as a convex quadratic problem with linear equality constraints given by min 2 x + δx x b 2 δx IR n, v IR m B + 2 v 2, subject to v R = H δx d, (3.5) see Gratton et al. (203). The Lagrange unction corresponding to this reormulated problem is given by L(δx, v, λ ) de = 2 x + δx x b 2 B + 2 v 2 λ T R (H δx v d ). Using this unction, the optimality conditions or problem (3.5) can then be expressed as which leads to the system δx L(δx, v, λ ) = B (x + δx x b ) H T λ = 0, v L(δx, v, λ ) = R (v ) + λ = 0, λ L(δx, v, λ ) = v H δx + d = 0, H T λ = B (x + δx x b ), λ = R (H δx d ). Gratton et al. (203) and Gratton and Tshimanga (2009) show that the solution o this system can be obtained by solving (3.6) (R H BH T + I m )λ = R (d H (x b x)), (3.7) or λ and then substituting in the second equation o (3.6) or δx and v. We now compare the solution (δx, λ ) to the solution o the coarse level minimization problem min 2 x + δx c x b 2 δx c IR n B + 2 Γ (H δx c d ) 2, (3.8) Rc where Γ is the restriction operator Γ : IR m IR mc rom the ine observation space to the coarse one. As above, we reormulate (3.8) as a convex quadratic problem with linear equality constraints, and obtain min 2 x + δx c x b 2 δx c IR n, v c IR mc B + 2 v c 2, subject to v Rc c = Γ (H δx c d ), (3.9) whose Lagrangian unction is given by L(δx c, v c, λ c ) de = 2 x + δx c x b 2 B + 2 v c 2 λ T Rc c (Γ H δx c v c Γ d ). The optimality conditions now become H T ΓT λ c = B (x + δx c x b ), λ c = R c Γ (H δx c d ) (3.0)
6 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 5 These conditions may again be solved by computing λ c such that (Rc Γ H BH T Γ T + I mc )λ c = Rc Γ (d H (x b x)). (3.) Since the dimension o the system (3.) is smaller than that o the system (3.7), solving problem (3.8) is (oten much) cheaper than solving problem (3.4). The RPCG algorithm mentioned in the previous section derives its eiciency by using the ormulation (3.) rather than (3.7), see Gratton et al. (203) and Gratton and Tshimanga (2009). Ater obtaining the solution (δx c, λ c ), our objective is now to compute the dierence between the (still unknown) λ and Π c λ c in order to identiy the set o observations where this dierence is large. Here Π c is the prolongation operator given by Π c de = σ Γ T (3.2) or some σ > 0. Our intention is then to construct the ine problem rom the coarse one by adding to the coarse the observations that are singled out by this comparison. We start by computing the dierence between λ and Π c λ c in an appropriate norm. Using (3.6) and (3.0), we deine and E 2 E de = λ Π c λ c 2 R = R (λ Π c λ c ), R (H δx d ) Π c λ c = λ Π c λ c, H δx + d R Π c λ c, de = λ Π c λ c 2 H BH T = H T λ H T Π cλ c 2 B By adding errors E and E 2 we obtain that = B(H T λ H T Π cλ c ), B (x + δx x b ) H T c Π c λ c = λ Π c λ c, H (x x b ) + H δx H BH T Π cλ c. E + E 2 = λ Π c λ c, H δx + d R Π c λ c As a consequence, we deduce that + λ Π c λ c, H (x x b ) + H δx H BH T Π cλ c + λ Π c λ c, H δx c H δx c = λ Π c λ c, d R Π c λ c H δx c + λ Π c λ c, H (x x b + δx c ) H BH T Π cλ c. E + E 2 = R /2 (λ Π c λ c ), R /2 (d H δx c R Π c λ c ) Thus, we obtain that + H BH T (λ Π c λ c ), (H BH T ) (H (x x b + δx c ) H BH T Π cλ c ) = R /2 (λ Π c λ c ), R /2 (d H δx c R Π c λ c ) + B /2 H T (λ Π c λ c ), B /2 H T (H BH T ) (H (x x b + δx c ) H BH T Π cλ c ). E + E 2 λ Π c λ c R d H δx c R Π c λ c R + H T λ H T Π cλ c B H T (H BH T ) (H (x x b + δx c ) H BH T Π cλ c ) B λ Π c λ c R d H δx c R Π c λ c R + λ Π c λ c H BH T H (x x b + δx c ) H BH T Π cλ c (H BH T ).
7 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 6 The Cauchy-Schwarz inequality then yields that [ ] /2 E + E 2 λ Π c λ c 2 R + λ Π c λ c 2 H BH T [ ] /2 d H δx c R Π c λ c 2 + H R (x x b + δx c ) H BH T Π cλ c 2 (H BH T ) = [ E + E 2 ] /2 [ ] /2 d H δx c R Π c λ c 2 + H R (x x b + δx c ) H BH T Π cλ c 2 (H. BH T ) This inequality leads to the a posteriori error bound E + E 2 d H δx c R Π c λ c 2 R + H (x x b + δx c ) H BH T Π c λ c 2 (H BH T ). (3.3) The question then arises o how to compute the second term in the right hand side o this inequality, i.e. H (x x b + δx c ) H BH T Π cλ c 2 (H BH T ) = H (x x b + δx c ) H BH T Π cλ c, (H BH T ) ( H (x x b + δx c ) H BH T Π cλ c ) = H (x x b + δx c ) H BH T Π cλ c, (H BH T ) H (x x b + δx c ) Π c λ c. But the irst equation o (3.0) multiplied by H B gives that (3.4) H BH T Γ T λ c = H (x x b + δx c ), (3.5) and hence that Thereore Γ T λ c = (H BH T ) H (x x b + δx c ). (3.6) H (x x b + δx c ) H BH T Π cλ c 2 (H BH T ) because o (3.2). We note also that = H (x x b + δx c ) H BH T Π cλ c, Γ T λ c Π c λ c = ( σ ) H (x x b + δx c ) H BH T Π cλ c, Π c λ c, (3.7) λ Π c λ c 2 R +H BH T = λ Π c λ c 2 R + λ Π c λ c 2 H BH T = E + E 2. (3.8) As a consequence, the let-hand side o this inequality (the sought error estimate) can be bounded above using inequality (3.3), giving the ollowing a posteriori error. Theorem 3. Let δx be the solution to the problem (3.5) and λ the corresponding Lagrange multiplier to the constraint in (3.5) such that the couple (δx, λ ) satisies (3.6). Analogously, let δx c be the solution to (3.9) and λ c the corresponding Lagrange multiplier such that (δx c, λ c ) satisies (3.0). Then the a posteriori error bound satisies the inequality λ Π c λ c 2 R +H BH T d H δx c R Π c λ c 2 R + ( σ ) H (x x b + δx c ) H BH T Π cλ c, Π c λ c. (3.9) Note that the bound (3.9) does not involve the computation o λ or δx. The use o an a posteriori error bound depending on a primal-dual problem in adaptive inite elements was studied in Rincon-Camacho (20). In the next section we describe how to make use o the bound (3.9) in order to algorithmically construct the ine problem rom the coarse one.
8 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 7 4 A new adaptive algorithm Starting rom a small set o observations O c, our goal is to add only signiicant observations to produce O so that the a posteriori error (3.9) is reduced. Our strategy is to deine an auxiliary set o potential ine observations Õ rom which the observations in O are selected. However, describing our strategy (and algorithm) requires additional assumptions on the hierarchy o (potential) observations. More speciically, we complete our assumptions as ollows. The observations correspond to localizations in some underlying continuous measurable observation space. The coarse observation set partitions the observation space in a inite number o cells {c j } pc j= o measures {w j } pc j=. The auxiliary set Õ is constructed by considering all observations in O c with the addition o a single additional potential observation point in the interior o each cell. The cell is said to be associated with this additional potential observation. There exists a prolongation operator Π c rom O c to Õ such that, or each potential observation o j in Õ \ O c, Π c deines the value o this observation only in terms o the observations o the associated cell c j. As expected, we deine Π c = σ ΓT. We illustrate these assumptions by an example: suppose that the observation space is the plane and the coarse observation set O c is the rectangular grid shown in Figure 4. (a): the cells are elementary rectangles in this grid, whose measure is given by their surace. We may then deine Õ as the grid shown in Figure 4. (b), which we obtained by locally adding a new potential observation in the center o each rectangle and our additional ones on its boundary (eectively doubling the mesh in every direction). The observations in O (as shown in Figure 4. (c)) can then be extracted rom Õ. (a) Coarse observation set O c (b) Auxiliary observation set Õ (c) Fine observation set O Figure 4.: Auxiliary observations In this example, the restriction operator Γ can be deined as the usual ull weighting operator associated with bilinear interpolation prolongations (which justiies the introduction o the
9 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 8 boundary points). This ull-weighting restriction operator is given, on every grid node, by the stencil (4.20) 6 2 In order to achieve our goal to select important observations rom Õ, we need to compute localized error indicators (3.9). We deine them as η j ν j de = w j ( d H δx c R Πc λ c ) j, R ( d H δx c R Πc λ c ) j, de = ( σ )w j ( H (x x b + δx c ) H B H T Π c λ c ) j, Π c λ c j, where d, H, R are constructed rom the set o observations Õ and where the symbol j denotes the restriction o the associated quantity to the cell c j. Note that the η s correspond to the irst term on the right-hand side o (3.9) while the ν s correspond the the second term. In order to decide which cells will be chosen to include a new interior potential observation, we use the bulk criterion (also known as Dörler marking, see, or instance, Dörler, 996, Morin, Nochetto and Siebert, 2000, Logg, Mardal and Wells, 202). For some constants θ, θ 2 (0, ), we construct the sets S η and S ν such that p θ η j p η k and θ 2 ν j ν k. (4.2) j= k S η j= k S ν In practice this construction is carried out by progressively constructing each set using a greedy heuristic which includes irst the non-included cell with maximal indicator value. Once S η and S ν are constructed, we decide that a cell k o O c is reined i k S η or k S ν, meaning that the observations associated with the corresponding cell in Õ are added to the set O c to construct the new set O. More ormally, de O = O c. (4.22) k S η S ν o k Thus, starting with a small set o observations O 0, we progressively add observations using the method just described, resulting in the ollowing algorithm. Algorithm 4. An algorithm using adaptive observations. Set i = 0, initialize x and the coarse observation set O Given the set o observations O i, construct the auxiliary set Õ i+ such that the conditions described at the beginning o this section hold. 3. Find the solution (δx i, λ i ) to the problem min 2 x i + δx i x b 2 δx i IR n B + 2 Γ i+ ( H i+ δx i d i+ ) 2, (4.23) R i by approximately solving the system (R i Γ i Hi+ B H i+ Γ T T i+ + I mi )λ i = R i Γ i+ ( d i+ H i+ (x b x i )) using RPCG and then setting δx i = x b x i + B H T i+ Γ T i λ i. 4. For each cell c j o observation set O i compute the error indicators η j = w j ( d i+ H i+ δx i R i+ Πi λ i ) j, R i+ ( d i+ H i+ δx i R i+ Πi λ i ) j, ν j = ( σ i+ )w j ( H i+ (x i x b + δx i ) H i+ B H T i+ Π i λ i ) j, Π i λ i j.
10 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 9 5. Build the sets S η and S ν such that p i θ p i η k, and θ 2 k S η using the bulk criterion. η j ν j j= j= k S ν η k. 6. Construct the set O i+ as. O i+ := O i k S η S ν o k 7. Update x x + δx i. 8. Increment i and return to Step 2. We note that the computation o (δx i, λ i ) in Step 3 corresponds to applying the RCPG algorithm to (4.23), thereby making this computation essentially dependent on m, the number o observations, which the algorithm maintains as small as necessary by design. We also note that the bulk o the remaining work is in the computation o the error indicators in Step 4, as this requires the product o H i+ with δx i. Observe that, because o (3.2), HT i+ Πi λ i = σ i+ HT i+ ΓT i+ λ i and that both this quantity and the product H i+ (x i x b ) are already available as by-products o Step 3. Similar methods or constructing adaptive grids are the multi-level adaptive technique (MLAT) studied in Brandt (973) and the ast adaptive composite grid (FAC) presented in McCormick (984). 5 Numerical experiments 5. Two test problems This section is devoted to showing the perormance o our new adaptive algorithm on two test cases. The irst one is a one-dimensional wave equation system, which we reer to as the D-Wave model rom now on. The dynamics on this model are governed by the ollowing nonlinear wave equations: u(z, t) 2 u(z, t) + (u) = 0, t z u(0, t) = u(, t) = 0, (5.24) u(z, 0) = u 0 (z), u(z, 0) = 0, t 0 t T, 0 z, where we have chosen (u) = µe ηu. In this case we look or the initial unction u 0 (z), which corresponds to x in the data assimilation problem (2.). We illustrate in Figure 5.2 (a) the initial state vector u 0 (x = u 0 ) and the evolution o the system in Figure 5.2 (b) (view rom the top) where the space domain corresponds to the horizontal axis and the time domain to the vertical axis. Our second example is the model reerred as Lorenz96 presented in Lorenz (995). The variable ū is a vector o N-equally spaced entries around a circle o constant latitude, i.e. ū(t) = (u (t), u 2 (t),..., u N (t)). The N-dimensional system is determined by the ollowing N equations du j dt = κ ( u j 2u j + u j u j+ u j + F ), j =,... N, (5.25)
11 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 0 (a) Initial x 0 := u 0 (z) (b) Dynamical system over space and time Figure 5.2: Nonlinear D Wave equation where F and κ are constants independent o j. To orm a cyclic chain, we set u N = u 0, u = u N and u N+ = u. This system is known to have a chaotic behaviour over time depending on the parameters N, F and κ (see or instance Karimi and Paul, 2009) For a given set o parameters N, F and κ or which a stable behavior is observed, i.e. where data assimilation can be perormed, we consider then the ollowing dynamical system: du j+θ dt = κ ( u j+θ 2u j+θ + u j+θ u j+θ+ u j+θ + F ), j =,... N, θ =,..., Θ, (5.26) where θ and Θ are integers. Thus, the new size o the vector ū(t) is N Θ, which may be speciied as large as needed in our numerical experiments. The dynamical system is plotted or N = 40, F = 8, κ = 20 and Θ = 0, using the coordinate graph described in Figure 5.3 (a). For an initial state x = ū(0) as that shown in Figure 5.3 (b) the system develops over time as described in Figure 5.3 (c) (view rom the top). As we observe that the system becomes chaotic ater a certain time, we consider a reduced window o assimilation plotted in Figure 5.3 (d). 5.2 Results We now provide numerical results using Algorithm 4., which uses the RPCG algorithm in Step 3, as was mentioned in Section 2. When tuning the accuracy parameter or stopping the inner iterations, we noticed that it was suitable to choose a value in the middle o the possible range. In the case o the DWave equation the range o the parameter is more or less deined by [0, 0 4 ], thus we choose 0 2, while the range is approximately [0, 0 6 ] in the case o the Lorenz96 model, and we chose the value 0 4 in our experiment. For the DWave problem, we depict the background vector in a black dashed line together with the true solution as a red line, in Figure 6.4 (a). The result given by our algorithm is plotted as a blue line in Figure 6.4 (b). The background vector and the true solution are plotted in Figure 6.5 (a) and the result given by our algorithm is presented in Figure 6.5 (b) or the Lorenz96 model. In order to observe the adaptive nature o the algorithm, or an intermediate iteration i, we display two consecutive sets o observations O i and O i+ in Figures 6.6 (a)-(b) and 6.7 (a)-(b) or the DWave and the Lorenz96 respectively. In order to appreciate the impact o the local error indicators deined in Step 4, we also illustrate, in Figures 6.6 (e) and 6.7 (e), the local behaviour o the error between the prolongation o the current λ i to the set O i+ and the true λ i+, as given by ɛ j = ( λ i+ Π i λ i ) j, [( R i+ + H i+ B H i+)( λ T i+ Π i λ i )] j, together with that o the local error indicators themselves (η j and ν j are displayed Figures 6.6 (c)-(d) and 6.7 (c)-(d), respectively).
12 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA (a) Coordinate system (b) Initial x 0 := (u (0), u 2 (0),..., u N (0)) (c) Dynamical system over space and time (d) Window o assimilation Figure 5.3: Lorenz96 problem In the case o the DWave equation we observe in Figure 5.2 (a) that the peak in the middle o the signal produces a dynamical reaction over space and time, shown in part (b) o the same igure. The errors η j and ν j in Figure 6.6 (c) and (d) are larger in the regions o strong dynamical activity, whose identiication is clear in Figure 6.6 (a)-(b) which shows the evolution o the observation sets. In Figure 6.6 (e), we also observe that the regions where the dierence ɛ j is large coincide with those where η j and ν j are also large in Figure 6.6 (c)-(d). In the case o the Lorenz96 model where the initial signal has also a peak in the middle (Figure 5.3 (b)), the high dynamics are on the let part o the space and time graph (Figure 5.3 (d)). In this case the quantities η j and ν j also provide correct indication o where more observations are needed, see Figure 6.7 (c)-(d). We also note that the distance ɛ j in Figure 6.7 (c) is consistent with η j and ν j. Indeed, the sets o observations in Figure 6.7 (a)-(b) show that the observations concentrate on the let part o the graph where the highly dynamic nature o the model is most evident. In Figures 6.8 and 6.9, we compare the perormance o the new algorithm with the simple use o uniorm observations and a benchmark method where a pre-established hierarchy o uniorm distributed and progressively denser observations sets is used in Algorithm 4. (skipping Steps 4-6). This last method bypasses the new adaptive eatures o the method and its associated computing cost. For this comparison, we deine, or a given iterate x resulting rom Step 7 o Algorithm 4., the cost unction as the value o the unction in (2.3) when i = r, i.e. when all the possible observations are used. We then plot the evolution o this cost unction (in logarithmic scale) against the number o observations used and the associated lop (loating-point operations) counts or the three algorithms. (The curves or the uniorm case do not correspond to an algorithm, but show the accuracy and computational costs associated with directly solving the problem on a uniorm grid or each speciic size.) We observe in both cases that the new method achieves a smaller error (plotted in red) than
13 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 2 that obtained by using uniormly distributed observations (plotted in blue) or the benchmark (plotted in black), which is explained by the act that most observations used by the adaptive algorithm are in regions o space where their contribution to accuracy is largest. Computational experience not reported here also indicates that the eect o the choice o starting value o x does not aect signiicantly the hierarchy o observation sets beyond the act that observations concentrate in regions o high dynamical activity. We also ound that [0.4, 0.7] appears to be an adequate range or the parameters θ and θ 2 in (4.2), yielding a satisactory rate o inclusion o new observations at each iteration. Both results on accuracy and computational costs are thereore highly encouraging. 6 Conclusions and perspectives A method which identiies the inluential data on a 4D-Var data assimilation problem is proposed, as well as an algorithm that exploits this identiication to improve on computing eiciency. The cpu-time gains are obtained or two reasons, the irst being that the available number o observations is used very eectively, and the second the act that the cost o the subproblem solution is signicantly reduced by the use o dual-space conjugate-gradient techniques like RPCG. Numerical experience has been presented on two nonlinear test problems, and the results are so ar encouraging. Further reinements o the algorithm could be considered, such as the use o adaptive preconditioners in the subproblem solver, and continued experience with the method is o course desirable to assess its true potential. Extensions o these ideas in other domains are also possible: we think in particular o data assimilation problems in rameworks where each state o the system is itsel an image on which adaptive reconstruction techniques could be applied. Acknowledgments. The authors would like to thank S. Gürol and X. Vasseur at CERFACS or helpul scientiic discussions.
14 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 3 (a) Background vector and true initial u(0) (b) True u(0) and algorithm solution Figure 6.4: Results: Nonlinear D Wave equation (a) Background vector and true initial u(0) (b) True u(0) and algorithm solution Figure 6.5: Results: Lorenz96
15 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 4 (a) O i (b) O i+ (c) η j (d) ν j (e) ɛ j Figure 6.6: Observations set and adaptive errors (a) O i (b) O i+ (c) η j (d) ν j (e) ɛ j Figure 6.7: Observations set and adaptive errors
16 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 5 (a) Cost unction versus number o observations (b) Cost unction versus lops Figure 6.8: Perormance o the algorithm on the nonlinear wave equation (a) Cost unction versus number o observations (b) Cost unction versus lops Figure 6.9: Perormance o the algorithm on the Lorenz96 problem
17 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 6 Reerences L.M. Berliner, Z.Q. Lu, and C. Snyder. Statistical design or adaptive weather observations. Journal o the Atmospheric Sciences, 56(5), , 999. A. Brandt. Multi-level adaptive technique (MLAT) or ast numerical solution to boundary value problems. in Proceedings o the Third International Conerence on Numerical Methods in Fluid Mechanics, pp , Heidelberg, Berlin, New York, 973. Springer Verlag. C. Cardinali, S. Pezzulli, and E. Andersson. Inluence-matrix diagnostic o a data assimilation system. Quarterly Journal o the Royal Meteorological Society, 30(603), , Ph. Courtier. Dual ormulation o our-dimensional variational assimilation. Quarterly Journal o the Royal Meteorological Society, 23, , 997. Ph. Courtier, J.-N. Thépaut, and A. Hollingsworth. A strategy or operational implementation o 4D-Var using an incremental approach. Quarterly Journal o the Royal Meteorological Society, 20, , 994. D.N. Daescu and I.M. Navon. Adaptive observations in the context o 4D-Var data assimilation. Meteorology and Atmospheric Physics, 85(4), , W. Dörler. A convergent adaptive algorithm or Poisson s equation. SIAM Journal on Numerical Analysis, 33(3), 06 24, 996. A. El Akkroui, P. Gauthier, S. Pellerin, and S. Buis. Intercomparison o the primal and dual ormulations o variational data assimilation. Quarterly Journal o the Royal Meteorological Society, 34(633), , G. H. Golub and C. F. Van Loan. Matrix Computations. Johns Hopkins University Press, Baltimore, third edn, 996. S. Gratton and J. Tshimanga. An observation-space ormulation o variational assimilation using a restricted preconditioned conjugate-gradient algorithm. Quarterly Journal o the Royal Meteorological Society, 35, , S. Gratton, S. Gürol, and Ph. L. Toint. Preconditioning and globalizing conjugate gradients in dual space or quadratically penalized nonlinear-least squares problems. Computational Optimization and Applications, 54(), 25, 203. S. Gratton, A. Lawless, and N. K. Nichols. Approximate Gauss-Newton methods or nonlinear least-squares problems. SIAM Journal on Optimization, 8, 06 32, J.A. Hansen and L.A. Smith. The role o operational constraints in selecting supplementary observations. Journal o the Atmospheric Sciences, 57(7), , M. R. Hestenes and E. Stieel. Methods o conjugate gradients or solving linear systems. Journal o the National Bureau o Standards, 49, , 952. A. Karimi and M.R. Paul. Extensive chaos in the Lorenz-96 model. Arxiv preprint arxiv: , M. Leutbecher. Adaptive observations, the Hessian metric and singular vectors. in Proc. ECMWF Seminar on Recent developments in data assimilation or atmosphere and ocean, Reading, UK, pp. 8 2, 2003a. M. Leutbecher. A reduced rank estimate o orecast error variance changes due to intermittent modiications o the observing network. Journal o the Atmospheric Sciences, 60(5), , 2003b.
18 Gratton, Rincon-Camacho, Toint: Adaptive observations and optimization in DA 7 J. Liu, E. Kalnay, T. Miyoshi, and C. Cardinali. Analysis sensitivity calculation in an ensemble Kalman ilter. Quarterly Journal o the Royal Meteorological Society, 35(644), , A. Logg, K.A. Mardal, and G.N. Wells. Automated solution o dierential equations by the inite element method, Vol. 84 o Lecture Notes in Computational Science and Engineering. Springer Verlag, Heidelberg, Berlin, New York, 202. E.N. Lorenz and K.A. Emanuel. Optimal sites or supplementary weather observations: Simulation with a small model. Journal o the Atmospheric Sciences, 55(3), , 998. E. N. Lorenz. Predictability: a problem partly solved. in Seminar on Predictability, Vol., pp. 8, Reading, UK, 995. ECMWF. S. T. McCormick. Fast adaptive composite grid (FAC) methods: Theory or the variational case. in Deect correction methods, pp Springer Verlag, Heidelberg, Berlin, New York, 984. P. Morin, R.H. Nochetto, and K.G. Siebert. Data oscillation and convergence o adaptive FEM. SIAM Journal on Numerical Analysis, 38(2), , M. M. Rincon-Camacho. Adaptive methods or Total Variation based image restoration. PhD thesis, University o Graz, Graz, Austria, 20.
Observation Thinning In Data Assimilation Computations
Observation Thinning In Data Assimilation Computations Serge Gratton, Monserrat Rincon-Camacho, Ehouarn Simon, and Philippe L. Toint June 23, 204 Abstract We propose to use an observation-thinning method
More informationData Assimilation for Weather Forecasting: Reducing the Curse of Dimensionality
Data Assimilation for Weather Forecasting: Reducing the Curse of Dimensionality Philippe Toint (with S. Gratton, S. Gürol, M. Rincon-Camacho, E. Simon and J. Tshimanga) University of Namur, Belgium Leverhulme
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationOn the Development of Implicit Solvers for Time-Dependent Systems
School o something FACULTY School OF OTHER o Computing On the Development o Implicit Solvers or Time-Dependent Systems Peter Jimack School o Computing, University o Leeds In collaboration with: P.H. Gaskell,
More informationAnalysis Scheme in the Ensemble Kalman Filter
JUNE 1998 BURGERS ET AL. 1719 Analysis Scheme in the Ensemble Kalman Filter GERRIT BURGERS Royal Netherlands Meteorological Institute, De Bilt, the Netherlands PETER JAN VAN LEEUWEN Institute or Marine
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationAn Ensemble Kalman Smoother for Nonlinear Dynamics
1852 MONTHLY WEATHER REVIEW VOLUME 128 An Ensemble Kalman Smoother or Nonlinear Dynamics GEIR EVENSEN Nansen Environmental and Remote Sensing Center, Bergen, Norway PETER JAN VAN LEEUWEN Institute or Marine
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More information(One Dimension) Problem: for a function f(x), find x 0 such that f(x 0 ) = 0. f(x)
Solving Nonlinear Equations & Optimization One Dimension Problem: or a unction, ind 0 such that 0 = 0. 0 One Root: The Bisection Method This one s guaranteed to converge at least to a singularity, i not
More informationLecture 8 Optimization
4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional
More informationRobust Residual Selection for Fault Detection
Robust Residual Selection or Fault Detection Hamed Khorasgani*, Daniel E Jung**, Gautam Biswas*, Erik Frisk**, and Mattias Krysander** Abstract A number o residual generation methods have been developed
More informationVariational data assimilation
Background and methods NCEO, Dept. of Meteorology, Univ. of Reading 710 March 2018, Univ. of Reading Bayes' Theorem Bayes' Theorem p(x y) = posterior distribution = p(x) p(y x) p(y) prior distribution
More informationParallel Algorithms for Four-Dimensional Variational Data Assimilation
Parallel Algorithms for Four-Dimensional Variational Data Assimilation Mie Fisher ECMWF October 24, 2011 Mie Fisher (ECMWF) Parallel 4D-Var October 24, 2011 1 / 37 Brief Introduction to 4D-Var Four-Dimensional
More informationDATA ASSIMILATION IN A COMBINED 1D-2D FLOOD MODEL
Proceedings o the International Conerence Innovation, Advances and Implementation o Flood Forecasting Technology, Tromsø, DATA ASSIMILATION IN A COMBINED 1D-2D FLOOD MODEL Johan Hartnac, Henri Madsen and
More informationOPTIMAL PLACEMENT AND UTILIZATION OF PHASOR MEASUREMENTS FOR STATE ESTIMATION
OPTIMAL PLACEMENT AND UTILIZATION OF PHASOR MEASUREMENTS FOR STATE ESTIMATION Xu Bei, Yeo Jun Yoon and Ali Abur Teas A&M University College Station, Teas, U.S.A. abur@ee.tamu.edu Abstract This paper presents
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More informationAdditional exercises in Stationary Stochastic Processes
Mathematical Statistics, Centre or Mathematical Sciences Lund University Additional exercises 8 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
More informationNONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION
NONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION Steven Ball Science Applications International Corporation Columbia, MD email: sball@nmtedu Steve Schaer Department o Mathematics
More informationIterative Methods for Stokes/Darcy Coupling
Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch
More informationSupplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values
Supplementary material or Continuous-action planning or discounted ininite-horizon nonlinear optimal control with Lipschitz values List o main notations x, X, u, U state, state space, action, action space,
More informationRevision of TR-09-25: A Hybrid Variational/Ensemble Filter Approach to Data Assimilation
Revision of TR-9-25: A Hybrid Variational/Ensemble ilter Approach to Data Assimilation Adrian Sandu 1 and Haiyan Cheng 1 Computational Science Laboratory Department of Computer Science Virginia Polytechnic
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More informationA Penalized 4-D Var data assimilation method for reducing forecast error
A Penalized 4-D Var data assimilation method for reducing forecast error M. J. Hossen Department of Mathematics and Natural Sciences BRAC University 66 Mohakhali, Dhaka-1212, Bangladesh I. M. Navon Department
More informationThe Ascent Trajectory Optimization of Two-Stage-To-Orbit Aerospace Plane Based on Pseudospectral Method
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (014) 000 000 www.elsevier.com/locate/procedia APISAT014, 014 Asia-Paciic International Symposium on Aerospace Technology,
More informationAPPROXIMATE NULLSPACE ITERATIONS FOR KKT SYSTEMS IN MODEL BASED OPTIMIZATION. min. s.t. c(x, p) = 0. c(x, p) = 0
Forschungsbericht Nr. 6-5, FB IV Mathematik/Inormatik, Universität Trier, 6 PPROXIMTE NULLSPCE ITERTIONS FOR KKT SYSTEMS IN MODEL BSED OPTIMIZTION KZUFUMI ITO, KRL KUNISCH, ILI GHERMN, ND VOLKER SCHULZ
More information2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.
.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics Spline interpolation was originally developed or image processing. In GIS, it is mainly used in visualization o spatial
More informationTelescoping Decomposition Method for Solving First Order Nonlinear Differential Equations
Telescoping Decomposition Method or Solving First Order Nonlinear Dierential Equations 1 Mohammed Al-Reai 2 Maysem Abu-Dalu 3 Ahmed Al-Rawashdeh Abstract The Telescoping Decomposition Method TDM is a new
More informationSensitivity analysis in variational data assimilation and applications
Sensitivity analysis in variational data assimilation and applications Dacian N. Daescu Portland State University, P.O. Box 751, Portland, Oregon 977-751, U.S.A. daescu@pdx.edu ABSTRACT Mathematical aspects
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. OPTIM. Vol. 18, No. 1, pp. 106 13 c 007 Society for Industrial and Applied Mathematics APPROXIMATE GAUSS NEWTON METHODS FOR NONLINEAR LEAST SQUARES PROBLEMS S. GRATTON, A. S. LAWLESS, AND N. K.
More informationMath 2412 Activity 1(Due by EOC Sep. 17)
Math 4 Activity (Due by EOC Sep. 7) Determine whether each relation is a unction.(indicate why or why not.) Find the domain and range o each relation.. 4,5, 6,7, 8,8. 5,6, 5,7, 6,6, 6,7 Determine whether
More informationObservation Impact Assessment for Dynamic. Data-Driven Coupled Chaotic System
Applied Mathematical Sciences, Vol. 10, 2016, no. 45, 2239-2248 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.65170 Observation Impact Assessment for Dynamic Data-Driven Coupled Chaotic
More informationROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS
ROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS I. Thirunavukkarasu 1, V. I. George 1, G. Saravana Kumar 1 and A. Ramakalyan 2 1 Department o Instrumentation
More informationDifferentiation. The main problem of differential calculus deals with finding the slope of the tangent line at a point on a curve.
Dierentiation The main problem o dierential calculus deals with inding the slope o the tangent line at a point on a curve. deinition() : The slope o a curve at a point p is the slope, i it eists, o the
More informationThe achievable limits of operational modal analysis. * Siu-Kui Au 1)
The achievable limits o operational modal analysis * Siu-Kui Au 1) 1) Center or Engineering Dynamics and Institute or Risk and Uncertainty, University o Liverpool, Liverpool L69 3GH, United Kingdom 1)
More informationA solid-fluid mixture theory of porous media
A solid-luid mixture theory o porous media I-Shih Liu Instituto de Matemática, Universidade Federal do Rio de Janeiro Abstract The theories o mixtures in the ramework o continuum mechanics have been developed
More informationRoberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points
Roberto s Notes on Dierential Calculus Chapter 8: Graphical analysis Section 1 Extreme points What you need to know already: How to solve basic algebraic and trigonometric equations. All basic techniques
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
More informationReview of Prerequisite Skills for Unit # 2 (Derivatives) U2L2: Sec.2.1 The Derivative Function
UL1: Review o Prerequisite Skills or Unit # (Derivatives) Working with the properties o exponents Simpliying radical expressions Finding the slopes o parallel and perpendicular lines Simpliying rational
More informationA general hybrid formulation of the background-error covariance matrix for ensemble-variational ocean data assimilation
.. A general hybrid formulation of the background-error covariance matrix for ensemble-variational ocean data assimilation Anthony Weaver Andrea Piacentini, Serge Gratton, Selime Gürol, Jean Tshimanga
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationLeast-Squares Spectral Analysis Theory Summary
Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,
More information(C) The rationals and the reals as linearly ordered sets. Contents. 1 The characterizing results
(C) The rationals and the reals as linearly ordered sets We know that both Q and R are something special. When we think about about either o these we usually view it as a ield, or at least some kind o
More informationExtreme Values of Functions
Extreme Values o Functions When we are using mathematics to model the physical world in which we live, we oten express observed physical quantities in terms o variables. Then, unctions are used to describe
More informationRobust feedback linearization
Robust eedback linearization Hervé Guillard Henri Bourlès Laboratoire d Automatique des Arts et Métiers CNAM/ENSAM 21 rue Pinel 75013 Paris France {herveguillardhenribourles}@parisensamr Keywords: Nonlinear
More informationSTAT 801: Mathematical Statistics. Hypothesis Testing
STAT 801: Mathematical Statistics Hypothesis Testing Hypothesis testing: a statistical problem where you must choose, on the basis o data X, between two alternatives. We ormalize this as the problem o
More informationAn Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator
An Alternative Poincaré Section or Steady-State Responses and Biurcations o a Duing-Van der Pol Oscillator Jang-Der Jeng, Yuan Kang *, Yeon-Pun Chang Department o Mechanical Engineering, National United
More informationIn the derivation of Optimal Interpolation, we found the optimal weight matrix W that minimizes the total analysis error variance.
hree-dimensional variational assimilation (3D-Var) In the derivation of Optimal Interpolation, we found the optimal weight matrix W that minimizes the total analysis error variance. Lorenc (1986) showed
More informationBANDELET IMAGE APPROXIMATION AND COMPRESSION
BANDELET IMAGE APPOXIMATION AND COMPESSION E. LE PENNEC AND S. MALLAT Abstract. Finding eicient geometric representations o images is a central issue to improve image compression and noise removal algorithms.
More informationANALYSIS OF p-norm REGULARIZED SUBPROBLEM MINIMIZATION FOR SPARSE PHOTON-LIMITED IMAGE RECOVERY
ANALYSIS OF p-norm REGULARIZED SUBPROBLEM MINIMIZATION FOR SPARSE PHOTON-LIMITED IMAGE RECOVERY Aramayis Orkusyan, Lasith Adhikari, Joanna Valenzuela, and Roummel F. Marcia Department o Mathematics, Caliornia
More informationRobin-Robin preconditioned Krylov methods for fluid-structure interaction problems
MOX Report No. 19/2008 Robin-Robin preconditioned Krylov methods or luid-structure interaction problems Santiago Badia, Fabio Nobile, Christian Vergara MOX, Dipartimento di Matematica F. Brioschi Politecnico
More informationBasic mathematics of economic models. 3. Maximization
John Riley 1 January 16 Basic mathematics o economic models 3 Maimization 31 Single variable maimization 1 3 Multi variable maimization 6 33 Concave unctions 9 34 Maimization with non-negativity constraints
More informationExponential and Logarithmic. Functions CHAPTER The Algebra of Functions; Composite
CHAPTER 9 Exponential and Logarithmic Functions 9. The Algebra o Functions; Composite Functions 9.2 Inverse Functions 9.3 Exponential Functions 9.4 Exponential Growth and Decay Functions 9.5 Logarithmic
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More informationOBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL
OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL Dionisio Bernal, Burcu Gunes Associate Proessor, Graduate Student Department o Civil and Environmental Engineering, 7 Snell
More informationFin efficiency of the newly developed Compartmented Coil of a Single Coil Twin Fan System
Fin eiciency o the newly developed Compartmented Coil o a Single Coil Twin Fan System ABSTRACT In predicting the perormance o any cooling coil, HVAC designers ace multiold challenges in designing the system
More informationThe University of Reading
The University of Reading Radial Velocity Assimilation and Experiments with a Simple Shallow Water Model S.J. Rennie 2 and S.L. Dance 1,2 NUMERICAL ANALYSIS REPORT 1/2008 1 Department of Mathematics 2
More information4. DATA ASSIMILATION FUNDAMENTALS
4. DATA ASSIMILATION FUNDAMENTALS... [the atmosphere] "is a chaotic system in which errors introduced into the system can grow with time... As a consequence, data assimilation is a struggle between chaotic
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationA Method for Assimilating Lagrangian Data into a Shallow-Water-Equation Ocean Model
APRIL 2006 S A L M A N E T A L. 1081 A Method or Assimilating Lagrangian Data into a Shallow-Water-Equation Ocean Model H. SALMAN, L.KUZNETSOV, AND C. K. R. T. JONES Department o Mathematics, University
More informationObjectives. By the time the student is finished with this section of the workbook, he/she should be able
FUNCTIONS Quadratic Functions......8 Absolute Value Functions.....48 Translations o Functions..57 Radical Functions...61 Eponential Functions...7 Logarithmic Functions......8 Cubic Functions......91 Piece-Wise
More informationOn the Efficiency of Maximum-Likelihood Estimators of Misspecified Models
217 25th European Signal Processing Conerence EUSIPCO On the Eiciency o Maximum-ikelihood Estimators o Misspeciied Models Mahamadou amine Diong Eric Chaumette and François Vincent University o Toulouse
More informationCONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a
Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions
More informationPreconditioning of Radial Basis Function Interpolation Systems via Accelerated Iterated Approximate Moving Least Squares Approximation
Preconditioning o Radial Basis Function Interpolation Systems via Accelerated Iterated Approximate Moving Least Squares Approximation Gregory E. Fasshauer and Jack G. Zhang Abstract The standard approach
More informationProbabilistic Model of Error in Fixed-Point Arithmetic Gaussian Pyramid
Probabilistic Model o Error in Fixed-Point Arithmetic Gaussian Pyramid Antoine Méler John A. Ruiz-Hernandez James L. Crowley INRIA Grenoble - Rhône-Alpes 655 avenue de l Europe 38 334 Saint Ismier Cedex
More informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
More informationA Local Ensemble Transform Kalman Filter Data Assimilation System for the Global FSU Atmospheric Model
Journal o Mechanics Engineering and Automation 5 (2015) 185-196 doi: 10.17265/2159-5275/2015.03.008 D DAVID PUBLISHING A Local Ensemble Transorm Kalman Filter Data Assimilation System or the Global FSU
More informationEstimation and detection of a periodic signal
Estimation and detection o a periodic signal Daniel Aronsson, Erik Björnemo, Mathias Johansson Signals and Systems Group, Uppsala University, Sweden, e-mail: Daniel.Aronsson,Erik.Bjornemo,Mathias.Johansson}@Angstrom.uu.se
More informationAH 2700A. Attenuator Pair Ratio for C vs Frequency. Option-E 50 Hz-20 khz Ultra-precision Capacitance/Loss Bridge
0 E ttenuator Pair Ratio or vs requency NEEN-ERLN 700 Option-E 0-0 k Ultra-precision apacitance/loss ridge ttenuator Ratio Pair Uncertainty o in ppm or ll Usable Pairs o Taps 0 0 0. 0. 0. 07/08/0 E E E
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationThe concept of limit
Roberto s Notes on Dierential Calculus Chapter 1: Limits and continuity Section 1 The concept o limit What you need to know already: All basic concepts about unctions. What you can learn here: What limits
More informationSimulation of electromagnetic devices using coupled models
Manuscript received Aug. 2, 2007; revised Nov. 30, 2007 Simulation o electromagnetic devices using coupled models ION CÂRSEA Department o Computer Engineering and Communication University o Craiova Str.
More informationOn High-Rate Cryptographic Compression Functions
On High-Rate Cryptographic Compression Functions Richard Ostertág and Martin Stanek Department o Computer Science Faculty o Mathematics, Physics and Inormatics Comenius University Mlynská dolina, 842 48
More informationCHAPTER 8 ANALYSIS OF AVERAGE SQUARED DIFFERENCE SURFACES
CAPTER 8 ANALYSS O AVERAGE SQUARED DERENCE SURACES n Chapters 5, 6, and 7, the Spectral it algorithm was used to estimate both scatterer size and total attenuation rom the backscattered waveorms by minimizing
More informationImprovement of Sparse Computation Application in Power System Short Circuit Study
Volume 44, Number 1, 2003 3 Improvement o Sparse Computation Application in Power System Short Circuit Study A. MEGA *, M. BELKACEMI * and J.M. KAUFFMANN ** * Research Laboratory LEB, L2ES Department o
More informationYURI LEVIN AND ADI BEN-ISRAEL
Pp. 1447-1457 in Progress in Analysis, Vol. Heinrich G W Begehr. Robert P Gilbert and Man Wah Wong, Editors, World Scientiic, Singapore, 003, ISBN 981-38-967-9 AN INVERSE-FREE DIRECTIONAL NEWTON METHOD
More informationFundamentals of Data Assimilation
National Center for Atmospheric Research, Boulder, CO USA GSI Data Assimilation Tutorial - June 28-30, 2010 Acknowledgments and References WRFDA Overview (WRF Tutorial Lectures, H. Huang and D. Barker)
More informationComputing proximal points of nonconvex functions
Mathematical Programming manuscript No. (will be inserted by the editor) Warren Hare Claudia Sagastizábal Computing proximal points o nonconvex unctions the date o receipt and acceptance should be inserted
More informationOptimal state estimation for numerical weather prediction using reduced order models
Optimal state estimation for numerical weather prediction using reduced order models A.S. Lawless 1, C. Boess 1, N.K. Nichols 1 and A. Bunse-Gerstner 2 1 School of Mathematical and Physical Sciences, University
More informationarxiv: v1 [math.na] 29 Oct 2018
A Multirate Approach or Fluid-Structure Interaction Computation with Decoupled Methods Lian Zhang 1, Mingchao Cai 2 Mo Mu 3, arxiv:1810.11933v1 [math.na] 29 Oct 2018 Abstract We investigate a multirate
More informationSymbolic-Numeric Methods for Improving Structural Analysis of DAEs
Symbolic-Numeric Methods or Improving Structural Analysis o DAEs Guangning Tan, Nedialko S. Nedialkov, and John D. Pryce Abstract Systems o dierential-algebraic equations (DAEs) are generated routinely
More informationA REPORT ON PERFORMANCE OF ANNULAR FINS HAVING VARYING THICKNESS
VOL., NO. 8, APRIL 6 ISSN 89-668 ARPN Journal o Engineering and Applied Sciences 6-6 Asian Research Publishing Networ (ARPN). All rights reserved. A REPORT ON PERFORMANCE OF ANNULAR FINS HAVING VARYING
More informationFeasibility of a Multi-Pass Thomson Scattering System with Confocal Spherical Mirrors
Plasma and Fusion Research: Letters Volume 5, 044 200) Feasibility o a Multi-Pass Thomson Scattering System with Conocal Spherical Mirrors Junichi HIRATSUKA, Akira EJIRI, Yuichi TAKASE and Takashi YAMAGUCHI
More informationAnalytical expressions for field astigmatism in decentered two mirror telescopes and application to the collimation of the ESO VLT
ASTRONOMY & ASTROPHYSICS MAY II 000, PAGE 57 SUPPLEMENT SERIES Astron. Astrophys. Suppl. Ser. 44, 57 67 000) Analytical expressions or ield astigmatism in decentered two mirror telescopes and application
More informationSparsity-Aware Pseudo Affine Projection Algorithm for Active Noise Control
Sparsity-Aware Pseudo Aine Projection Algorithm or Active Noise Control Felix Albu *, Amelia Gully and Rodrigo de Lamare * Department o Electronics, Valahia University o argoviste, argoviste, Romania E-mail:
More information2 Frequency-Domain Analysis
2 requency-domain Analysis Electrical engineers live in the two worlds, so to speak, o time and requency. requency-domain analysis is an extremely valuable tool to the communications engineer, more so
More informationControlling the Heat Flux Distribution by Changing the Thickness of Heated Wall
J. Basic. Appl. Sci. Res., 2(7)7270-7275, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal o Basic and Applied Scientiic Research www.textroad.com Controlling the Heat Flux Distribution by Changing
More informationLevenberg-Marquardt methods based on probabilistic gradient models and inexact subproblem solution, with application to data assimilation
Levenberg-Marquardt methods based on probabilistic gradient models and inexact subproblem solution, with application to data assimilation E. Bergou S. Gratton L. N. Vicente June 26, 204 Abstract The Levenberg-Marquardt
More informationToday. Introduction to optimization Definition and motivation 1-dimensional methods. Multi-dimensional methods. General strategies, value-only methods
Optimization Last time Root inding: deinition, motivation Algorithms: Bisection, alse position, secant, Newton-Raphson Convergence & tradeos Eample applications o Newton s method Root inding in > 1 dimension
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationCOMPARATIVE EVALUATION OF ENKF AND MLEF FOR ASSIMILATION OF STREAMFLOW DATA INTO NWS OPERATIONAL HYDROLOGIC MODELS
COMPARATIVE EVALUATION OF ENKF AND MLEF FOR ASSIMILATION OF STREAMFLOW DATA INTO NWS OPERATIONAL HYDROLOGIC MODELS Arezoo Raieei Nasab 1, Dong-Jun Seo 1, Haksu Lee 2,3, Sunghee Kim 1 1 Dept. o Civil Eng.,
More informationAPPENDIX 1 ERROR ESTIMATION
1 APPENDIX 1 ERROR ESTIMATION Measurements are always subject to some uncertainties no matter how modern and expensive equipment is used or how careully the measurements are perormed These uncertainties
More informationAnalysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits
Computer Applications in Electrical Engineering Vol. 4 Analysis o the regularity pointwise completeness pointwise generacy o descriptor linear electrical circuits Tadeusz Kaczorek Białystok University
More informationIonosphere Ray-Tracing of RF Signals and Solution Sensitivities to Model Parameters
Ionosphere Ray-Tracing o RF Signals and Solution Sensitivities to Model Parameters Mark L. Psiaki Virginia Tech, Blacksburg, VA 2461, USA BIOGRAPHY Mark L. Psiaki is Proessor and Kevin Croton Faculty Chair
More informationarxiv: v3 [math.oc] 8 May 2018
Noname manuscript No will be inserted by the editor Generalized Sel-Concordant Functions: A Recipe or Newton-Type Methods Tianxiao Sun Quoc Tran-Dinh arxiv:70304599v3 [mathoc] 8 May 08 Received: date /
More informationGaussian Process Regression Models for Predicting Stock Trends
Gaussian Process Regression Models or Predicting Stock Trends M. Todd Farrell Andrew Correa December 5, 7 Introduction Historical stock price data is a massive amount o time-series data with little-to-no
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationA recursive model-based trust-region method for derivative-free bound-constrained optimization.
A recursive model-based trust-region method for derivative-free bound-constrained optimization. ANKE TRÖLTZSCH [CERFACS, TOULOUSE, FRANCE] JOINT WORK WITH: SERGE GRATTON [ENSEEIHT, TOULOUSE, FRANCE] PHILIPPE
More informationA PROBABILISTIC POWER DOMAIN ALGORITHM FOR FRACTAL IMAGE DECODING
Stochastics and Dynamics, Vol. 2, o. 2 (2002) 6 73 c World Scientiic Publishing Company A PROBABILISTIC POWER DOMAI ALGORITHM FOR FRACTAL IMAGE DECODIG V. DRAKOPOULOS Department o Inormatics and Telecommunications,
More informationOptimal Control. with. Aerospace Applications. James M. Longuski. Jose J. Guzman. John E. Prussing
Optimal Control with Aerospace Applications by James M. Longuski Jose J. Guzman John E. Prussing Published jointly by Microcosm Press and Springer 2014 Copyright Springer Science+Business Media New York
More informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More information