Optimal Reconstruction of Material Properties in Complex Multiphysics Phenomena
|
|
- Myles Patterson
- 5 years ago
- Views:
Transcription
1 Annual Meeting of the Canadian Society of Applied and Industrial Mathematics CAIMS 23 Optimal Reconstruction of Material Properties in Complex Multiphysics Phenomena Vladislav Bukshtynov Department of Energy Resources Engineering Stanford University CAIMS SCMAI Quebec City June 7, 23
2 he Cecil Graham Doctoral Dissertation Award: Acknowledgments Prof. Bartosz Protas (McMaster University) Words of appreciation to: Dr. Oleg Volkov (postdoc at McMaster, now at Stanford) Prof. Dmitry Pelinovsky (McMaster) Profs. David Novog, Nikolas Provatas and James Wadsley (McMaster) Profs. Khalid Aziz and Lou Durlofsky (Stanford)
3 Generous Funding Provided by: School of Computational Science & Engineering and Department of Mathematics & Statistics, McMaster University Natural Sciences and Engineering Research Council of Canada (Collaborative Research and Development Program) Ontario Centres of Excellence Centre for Materials and Manufacturing General Motors of Canada Energy Resources Engineering Department and Smart Fields Consortium, Stanford University
4 Acknowledgments Motivation Simple Model Multiphysics Model Reduced-Order Models Special hanks to My Family My wife Katerina and daughter Ksenia Open Problem & Outlook
5 Agenda Motivation Generalized Framework (Simple Model Problem) Multiphysics Problem Reduced-Order Models Open Problem & Outlook
6 Motivation Generalized Framework (Simple Model Problem) Multiphysics Problem Reduced-Order Models Open Problem & Outlook
7 Reconstruction of Constitutive Relations: Motivation Application: Welding in automotive industry (with GM of Canada) Problem: Unknown material properties of alloys in liquid phase viscosity µ( ) thermal conductivity k( ) surface tension σ( ) all functions of temperature. Goal: How to reconstruct µ( ), k( ), σ( ), etc., based only on available measurements?
8 wo distinct formulations: Parameter Estimation as Inverse Problem a) material properties dependent on the space variable (i.e., the independent variable), and b) material properties dependent on the state variable (i.e., the dependent variable), less attention in the literature
9 Generalized Problem Setup o propose and validate a suite of computational techniques to reconstruct certain material properties based on partial (pointwise) measurements available for a particular process and a particular material Strategy: Formulate an inverse problem of the second type (b) as partial differential equation (PDE) constrained optimization problem which can be solved using an adjoint based gradient method. Propose an optimize then discretize approach where the optimality conditions and the optimization algorithm are both formulated based on the continuous (PDE) form of the problem.
10 Problem Extensions he approach is easily extendable to many problems arising in nonequilibrium thermodynamics and continuum mechanics and is independent of the domain geometry and discretization. simple model multiphysics problem: complications of mathematical character only, e.g. in derivation of the cost functional gradient adding time dependence and/or changing the model geometry affects only the computational part of the solution Examples: D model problem (heat conductivity) 2D multiphysics problem (viscosity in coupled Navier-Stokes and energy equations) reduced-order models of hydrodynamic instabilities (inertial manifolds) open problem in petroleum engineering
11 Motivation Generalized Framework (Simple Model Problem) Multiphysics Problem Reduced-Order Models Open Problem & Outlook
12 Statement of the Problem Assuming the constitutive relation in thermodynamics in a general form [ ] [ ] thermodynamic thermodynamic = k (state variables), flux force in heat conduction problem this relation takes the form of the BVP for steady state heat equation: [k( ) ] = g(x) in Ω, = on Ω spatial domain x Ω, where Ω R n, n =, 2, 3 state (temperature) (x) : Ω R distribution of heat sources g(x) : Ω R control (thermal diffusivity) k( ) : R R well-posedness of the problem: g(x) >, k( ) > measurements { i } M i= at a number of points {x i } M i= in the domain Ω
13 Parameter Estimation as Optimization Problem Reformulation as minimization problem: defining the cost functional J : X R as J (k) 2 M [ (x i ; k) ] 2 i i= optimal reconstruction ˆk to be obtained as an unconstrained minimizer of cost functional J ˆk = argmin J (k) k X minimizer ˆk can be computed with the gradient descent algorithm as ˆk = lim n k (n), where { k (n+) = k (n) τ (n) k J (k (n) ), n =,... k () = k minimizer ˆk can be found by using a set of adjoint based techniques
14 Extraction of Cost Functional Gradients Gâteaux (directional) differential: J (k; k ) = M [ (x i ; k) i ] (x i ; k, k ), i= where perturbation variable (x i ; k, k ) satisfies the perturbation system obtained from the governing equation. Riesz representation theorem (for some Hilbert space X ) J (k; k ) = k J, k X, e.g. setting X = L 2 gives J (k; k ) = k J k dx. Inconsistency: k is hidden in the system defining (x i ; k, k ) Solution: J (k; k ) requires a form consistent with a Riesz representation. Ω
15 Gradient via Adjoint Analysis Adjoint variable (Lagrange multiplier) satisfies the adjoint system k( ) = = M [ i (x)]δ(x x i ) in Ω, i= where δ( ) is the Dirac delta function. on Ω, Gâteaux differential (after integration by parts and Kirchhoff transformation) β [ ] J (k; k ) = χ [α, (x)](s) dx k (s) ds α Ω β [ χ [α, (x)](s) ] α Ω n dσ k (s) ds = k J, k X, {, s [α, β ], χ [α,β ](s) =, s / [ α, β ]. [Proved in Bukshtynov, Volkov & Protas (2)]
16 Cost Functional Gradient: Solution Commonly used structure adapted to steady-state problem after Chavent & Lemonnier (974, in French) L 2 k J (s) = δ( (x) s) dx, Ω where Dirac delta function δ( (x) s) is used to define a level set of the temperature field (x) Γ s {x Ω, (x) = s} Alternative (computationally more efficient) structure of gradient L 2 k J (s) = M i= χ [α, (x i )](s) i (x i ) k( (x i )) χ [α, (x)](s) Ω n dσ [Equivalence proved in Bukshtynov, Volkov & Protas (2)]
17 Differentiability of Map k Solving parameter estimation problem for [k( ) ] = g in Ω, = on Ω or in a variational formulation in the space = { H (Ω); = on Ω} k( ) u dx = g u dx,, u (2) Ω is equivalent to finding a solution to the operator equation F(k) =, where F : K L 2(Ω) and K is the set of constitutive relations Ω K = {k( ) piecewise C on L; < m k < k( ) < M k, L} () heorem [Bukshtynov, Volkov & Protas (2)] Assume that m k is sufficiently large and solutions of () satisfy L (Ω) <. hen the map k ( ; k) from K to L 2(Ω) defined by (2) is Fréchet differentiable in the norm H (I).
18 Main Complications. L 2 gradient L 2 k J (s) is in general a discontinuous function necessary regularity (smoothness) achieved by preconditioning
19 Main Complications. L 2 gradient L 2 k J (s) is in general a discontinuous function necessary regularity (smoothness) achieved by preconditioning 2. Min problem is in general nonconvex condition J (k; k ) = characterizes only a local, rather than global, minimizer choice of initial guess k
20 Main Complications. L 2 gradient L 2 k J (s) is in general a discontinuous function necessary regularity (smoothness) achieved by preconditioning 2. Min problem is in general nonconvex condition J (k; k ) = characterizes only a local, rather than global, minimizer choice of initial guess k 3. Gradient L 2 k J (s) is nontrivial only for s I, and sensitivity of the functional J (k) to perturbations of k can only be determined on the identifiability region I algorithm for shifting (extending) identifiability region I
21 Main Complications. L 2 gradient L 2 k J (s) is in general a discontinuous function necessary regularity (smoothness) achieved by preconditioning 2. Min problem is in general nonconvex condition J (k; k ) = characterizes only a local, rather than global, minimizer choice of initial guess k 3. Gradient L 2 k J (s) is nontrivial only for s I, and sensitivity of the functional J (k) to perturbations of k can only be determined on the identifiability region I algorithm for shifting (extending) identifiability region I 4. Positivity condition k( ) > (well-posedness) is required for all L augmenting functional J (k) with the barrier function straightforward approach using slack variables
22 Regularization via Preconditioning Riesz representation theorem J (k; k ) = L 2 k J, k L 2 = H k J, k H and Sobolev space H is endowed with the inner product β [ ] z, z 2 H (L) = z z 2 + l 2 z z 2 ds, z, z 2 H (L), s s α where l R is a length scale (smoothening) parameter. Sobolev gradient H k J can be determined as a solution of the elliptic BVP H k J l 2 d 2 ds 2 H k J = L 2 d ds H k J on (a, b), k J = for s = a, b.
23 Single Identifiability Interval Reconstruction k() L 2 J vs. H J (l =.5 and l =.2) ˆk( ) using H J
24 Single Identifiability Interval Reconstruction k() L 2 J vs. H J (l =.5 and l =.2) ˆk( ) using H J
25 Single Identifiability Interval Reconstruction k() L 2 J vs. H J (l =.5 and l =.2) ˆk( ) using H J
26 Single Identifiability Interval Reconstruction k() L 2 J vs. H J (l =.5 and l =.2) ˆk( ) using H J
27 Reconstruction with Extended Gradients k() H k J (s) = H k J ( α) for s [ a, α], H k J (s) = H k J ( β ) for s [ β, b ]
28 Reconstruction with Shifting Reconstruction of k( ) over the entire interval of interest L by determining new heat source distribution boundary conditions g (j+) (x) g (j) (x) h [k (j) ( (j) (x))],(j+) =,(j) + h to perform shifting to right (h > ) or/and left (h < ) with step h << k() 2 k()
29 Reconstruction with Shifting Reconstruction of k( ) over the entire interval of interest L by determining new heat source distribution boundary conditions g (j+) (x) g (j) (x) h [k (j) ( (j) (x))],(j+) =,(j) + h to perform shifting to right (h > ) or/and left (h < ) with step h << k() 2 k()
30 Reconstruction with Shifting Reconstruction of k( ) over the entire interval of interest L by determining new heat source distribution boundary conditions g (j+) (x) g (j) (x) h [k (j) ( (j) (x))],(j+) =,(j) + h to perform shifting to right (h > ) or/and left (h < ) with step h << k() 2 k()
31 Motivation Generalized Framework (Simple Model Problem) Multiphysics Problem Reduced-Order Models Open Problem & Outlook
32 General Statement Goal: o reconstruct the viscosity coefficient µ( ) in the momentum equation (Navier-Stokes equation) coupled with energy equation (in 2D or 3D): [ ] tu + u u + p µ( )[ u + ( u) ] = in Ω, subject to boundary and initial conditions u = in Ω, t + u [k ] = in Ω, u = u b = b on Ω, on Ω, u(, ) = u, (, ) = in Ω. spatial domain x Ω, where Ω R d, d = 2, 3 states: u(x) : Ω R n, (x) : Ω R control (viscosity) µ( ) : R R + measurements { i } M i= at a number of points {x i } M i= in the domain Ω
33 Computations for Iterative Reconstruction. Direct problem (forward in time) 2. Adjoint problem (backward in time) tu u u σ + u ( u) + = in Ω, t u [k ] + dµ d ( )[ u + ( u) ] : u = u = in Ω, M [ (x i ; µ) i ]δ(x x i ) in Ω, where σ p I + (θ 2 ( ) + m µ) [ u + ( u ) ] and subject to 3. Cost functional gradient tf L 2 µ J (s) = i= u =, = on Ω, u ( ; t f ) =, ( ; t f ) = in Ω. [Derived in Bukshtynov & Protas (23)] Ω δ( (x) s)[ u + ( u) ] : u dx dτ
34 Positivity Constraint µ( ) > Positivity condition µ( ) > required for all L (well-posedness) ˆµ = argmin µ X, J (µ) µ( )>, L Augmenting cost functional by adding penalization term (barrier function).
35 Positivity Constraint µ( ) > Positivity condition µ( ) > required for all L (well-posedness) ˆµ = argmin µ X, J (µ) µ( )>, L Augmenting cost functional by adding penalization term (barrier function). Straightforward approach: slack variable θ( ) : R R µ( ) = θ 2 ( ) + m µ, m µ R +, m µ.
36 Positivity Constraint µ( ) > Positivity condition µ( ) > required for all L (well-posedness) ˆµ = argmin µ X, J (µ) µ( )>, L Augmenting cost functional by adding penalization term (barrier function). Straightforward approach: slack variable θ( ) : R R µ( ) = θ 2 ( ) + m µ, m µ R +, m µ. New unconstrained (except governing PDE) optimization problem where µ( ) > is satisfied automatically with L 2 gradient ˆθ = argmin J (θ), J (θ) tf θ X 2 L 2 θ J (s) = 2 tf Ω M i= [ (x i ; θ) i ] 2 dτ δ( (x) s) θ(s) [ u + ( u) ] : u dx dτ
37 Setup for Numerical Simulations algebraic expression (Andrade Law) to represent viscosity µ( ) µ( ) = C e C 2/ geometry: simple 2D lid driven (shear driven) cavity flow thermal conductivity k =.2 = const local Re = 25 BCs: top boundary (moving lid) u x = u top = cos(αt), u y =, = 5 all other: u x = u y =, = 3 synthetic temperature measurements: 9 thermal sensors capture the temperature correspondent to µ( ) μ()
38 Setup for Numerical Simulations algebraic expression (Andrade Law) to represent viscosity µ( ) µ( ) = C e C 2/ geometry: simple 2D lid driven (shear driven) cavity flow thermal conductivity k =.2 = const local Re = 25 BCs: top boundary (moving lid) u x = u top = cos(αt), u y =, = 5 all other: u x = u y =, = 3 synthetic temperature measurements: 9 thermal sensors capture the temperature correspondent to µ( )
39 Vorticity ω Simulations for Direct and Adjoint Fields: Vorticity Adjoint vorticity ω x t =.5
40 Vorticity ω Simulations for Direct and Adjoint Fields: Vorticity Adjoint vorticity ω x 4 x t =.5 t =.5
41 Vorticity ω Simulations for Direct and Adjoint Fields: Vorticity Adjoint vorticity ω x 4 x 4 x t =.5 t =.5 t =.95
42 Simulations for Direct and Adjoint Fields: emperature emperature Adjoint temperature t =
43 Simulations for Direct and Adjoint Fields: emperature emperature Adjoint temperature t =.5 t =.5
44 Simulations for Direct and Adjoint Fields: emperature emperature Adjoint temperature t =.5 t =.5 t =.95
45 Integration over Level Sets (emperature Isolines) Cost functional gradient L 2 µ J (s) = tf Ω δ( (x) s)[ u + ( u ) ] : u dx dτ Integration over codimension manifolds defined by a level set function φ(x) Γ: φ(x)= h(x)dσ = δ(φ(x)) φ(x) h(x) dx Ω
46 Integration over Level Sets (emperature Isolines) Cost functional gradient L 2 µ J (s) = tf Ω δ( (x) s)[ u + ( u ) ] : u dx dτ Integration over codimension manifolds defined by a level set function φ(x) Γ: φ(x)= h(x)dσ = Ω δ(φ(x)) φ(x) h(x) dx Points to consider: regularity of isolines φ(x, s) = (x) s discrete definition of field (x) computational efforts
47 Integration Over Codimension Manifolds Approaches to evaluate h(x)dσ Γ: φ(x)= fall into two main groups: or Ω δ(φ(x)) φ(x) h(x) dx. reduction to a line (contour) integral x ( ) dy 2 h(x, y(x)) + dx x dx t t ( ) dx 2 ( ) dy 2 h(x(t), y(t)) + dt dt dt geometric integration: Min & Gibou (27, 28) 2. evaluation of an area integral regularization δ ɛ of the Dirac delta measure with a suitable choice ɛ: Zahedi & ornberg (2) approximations δ of the Dirac delta function: Smereka (26), owers (27, 29)
48 Approximation of Contour Integrals with Area Integrals
49 Approximation of Contour Integrals with Area Integrals heorem [Bukshtynov & Protas (23)] Method (*) using 2-D compound midpoint rule is 2 nd order accurate with respect to the discretization of the space domain and st order accurate with respect to the discretization of the state domain f (s) = Ω δ( (x) s)g(x) dx = h N s,h j= g(x j )hx hy + O(h2 ) + O(h )
50 Approximation of Contour Integrals: Performance 3 2 Ne /2 ξ 2 h.5 time, s N e 3 2 h h Relative error vs. discretization step h = x = y CPU time (sec) vs. number of computational elements N e Ne Line integration ( ) vs. Delta Function Approximation ( ) vs. Area Integration( )
51 Validation of Gradients κ(ε)..99 κ ε ε Directional differential J (µ; µ ) obtained in two ways: finite difference approximation and using adjoint field (different perturbations µ and µ 2) κ(ɛ) ɛ [J (µ + ɛµ ) J (µ)] + µj (s) µ (s) ds Ratio κ as t (- - - and )
52 Choice of Initial Guess μ().2.5 μ() μ().2.5 μ()
53 Choice of Initial Guess μ().2.5 μ() μ().2.5 μ()
54 Choice of Initial Guess μ().2.5 μ() μ().2.5 μ()
55 Choice of Initial Guess μ().2.5 μ() μ().2.5 μ()
56 Choice of Initial Guess μ().2.5 μ() μ().2.5 μ()
57 Choice of Initial Guess μ().2.5 μ() μ().2.5 μ()
58 Choice of Initial Guess μ().2.5 μ() μ().2.5 μ()
59 Choice of Initial Guess μ().2.5 μ() μ().2.5 μ()
60 Identifiability Region: Initial vs. Extended Cases μ() Reconstruction on I with initial guess µ = µ(α) + µ( β) 2
61 Identifiability Region: Initial vs. Extended Cases μ() Reconstruction on I with initial guess µ = µ(α) + µ( β) 2
62 Identifiability Region: Initial vs. Extended Cases.35 μ() μ() Reconstruction on I with initial guess µ = µ(α) + µ( β) 2 I = L with BCs for b assigned b top = b b else = a
63 Identifiability Region: Initial vs. Extended Cases.35 μ() μ() Reconstruction on I with initial guess µ = µ(α) + µ( β) 2 I = L with BCs for b assigned b top = b b else = a
64 Identifiability Region: Initial vs. Extended Cases.35 μ() μ() Reconstruction on I with initial guess µ = µ(α) + µ( β) 2 I = L with BCs for b assigned b top = b b else = a
65 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.5.4 μ()
66 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.5% α.5.4 μ()
67 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.5.4 μ()
68 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.3% α.5.4 μ()
69 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.5% α.5.4 μ()
70 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.5.4 μ()
71 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.6 λ = μ().3 μ()
72 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.6 λ = μ().3 μ()
73 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.6 λ = μ().3 μ()
74 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.6 λ = μ().3 μ()
75 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.6 λ = μ().3 μ()
76 Noise in Measurements Cost functional J λ (µ) J (µ) + λ 2 tf µ µ 2 Y(I) dτ.6 augmented with ikhonov regularization term, e.g. of H form, where λ 2 µ µ 2 Ḣ (I) = λ β ( dµ 2 ds d µ ) 2 ds ds noise =.% α.6 λ = μ().3 μ()
77 Motivation Generalized Framework (Simple Model Problem) Multiphysics Problem Reduced-Order Models Open Problem & Outlook
78 Optimal Reconstruction of Inertial Manifolds in Reduced-Order Models by Protas, Noack & Morzyński (22) admor et al. (2) Approximation within a low-dimensional space R N by evolution equation d dt a = f(a), a RN, f : R N R N Goal: parameter-free propagator f identification in a system with oscillatory dynamics
79 From -dim to dim = 3 Using POD modes as basis consider model reduction } tu + u u + p ν u = in Ω, = u = in Ω, ṙ = g (r) r, θ = g 2(r), a 3 = g 3(r), where r a 2 + a2 2 and θ arctan(a2/a). Measurements: ã(t) = [ã (t), ã 2(t), ã (t)] [ r(t), θ(t), ã (t)], for t [, ] Cost functionals for each optimization problem: J (g ) 2 J 2(g 2) 2 J 3(g 3) 2 [r(t) r(t)] 2 dt [ e ıθ(t) e ı θ(t) ] 2 dt [g 3(r(t)) ã (t)] 2 dt
80 Optimal Reconstructions Model: 2D flow around a circular cylinder g (r) in ṙ = g (r) r.5. g 2(r) in θ = g 2(r) g (r) g 2 (r) r Protas et al. (22) initial guess (mean-field model) [ ] ṙ(t) = σ β α M r 2 (t) r(t), r θ(t) = ω + γ α M r 2 (t) optimal reconstruction measurement data dr (r(t r dt i)), dθ dt (r(t i)), i =,..., N
81 Comparison of Output (state magnitude r = a 2 + a2 2 ) Protas et al. (22) r(t) t initial guess (mean-field model) optimal reconstruction measurement data
82 Comparison of Output (state variable a ) Protas et al. (22) a (t) t initial guess (mean-field model) optimal reconstruction measurement data
83 Motivation Generalized Framework (Simple Model Problem) Multiphysics Problem Reduced-Order Models Open Problem & Outlook
84 Parameter Reconstruction for Large Real Field Models Norne field in the Norwegian Sea operated by Statoil, Statoil (2) Norne field model ( 45, active cells, 3 wells, 9-year production history available)
85 Closed Loop Reservoir Modelling
86 History Matching as Optimization Problem Goal: to determine geological model parameters m = [m m 2... m Nm ] s.t. simulated results G sim (m) match observed (production) data d obs = [d d 2... d Nd ] and honor prior information about geological structure of the model. Optimal solution ˆm is sought as Sarma et al. (SPE 676) min m N d i= (G i sim(m) d i obs) 2 + R geology s.t. g(x; m) = m feasibly bounded set where R geology is geology pertaining regularization
87 Open Problem: Reconstruction of Relative Permeability K r (S) Mass balance equation for (O)il, (G)as and (W)ater phases, p = o, g, w (simplified black oil model) where Darcy velocity of phase p (φρpsp) + ρpup =, t u p = k Krp(S) µ p [ p p ρ pg] relative permeability, K r oil water water saturation, S w K r,water and K r,oil as functions of S water setup in tables approximations, e.g. Corey-type (black oil) [ K rw (S w ) = Krw S w S wi S wi S orw ] Nw Enhanced Oil Recovery (EOR) models (multicomponent) require more complicated dependencies K rp(s w ) and K rp(s g )
88 Conclusions Novel computational approach to reconstruct solution dependent constitutive relations based on available measurement data using a gradient based optimization technique Systematic experimental design procedure: reconstruction over a broader range of the state variable required degree of smoothness of the reconstructed material properties superior accuracy and efficiency of the novel numerical approach for evaluation of cost functional gradients
89 References. V. Bukshtynov, O. Volkov, B. Protas, On Optimal Reconstruction of Constitutive Relations, Physica D 24, , (2). 2. V. Bukshtynov, B. Protas, Optimal Reconstruction of Material Properties in Complex Multiphysics Phenomena, Journal of Computational Physics 242, , (23). 3. B. Protas, B. R. Noack and M. Morzyński, An Optimal Model Identification For Oscillatory Dynamics With A Stable Limit Cycle, (submitted), (22).
On Optimal Reconstruction of Constitutive Relations
On Optimal Reconstruction of Constitutive Relations Vladislav Bukshtynov a, Oleg Volkov b, and Bartosz Protas b, a School of Computational Engineering and Science, McMaster University, Hamilton, Ontario,
More informationProbing Fundamental Bounds in Hydrodynamics Using Variational Optimization Methods
Probing Fundamental Bounds in Hydrodynamics Using Variational Optimization Methods Bartosz Protas and Diego Ayala Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL:
More informationOn Some Variational Optimization Problems in Classical Fluids and Superfluids
On Some Variational Optimization Problems in Classical Fluids and Superfluids Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas
More informationNumerical Optimization of Partial Differential Equations
Numerical Optimization of Partial Differential Equations Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas Rencontres
More informationReconstruction of Concentration-Dependent. Material Properties in Electrochemical Systems
Reconstruction of Concentration-Dependent Material Properties in Electrochemical Systems RECONSTRUCTION OF CONCENTRATION-DEPENDENT MATERIAL PROPERTIES IN ELECTROCHEMICAL SYSTEMS BY ATHINTHRA KRISHNASWAMY
More informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationSimulation of Naturally Fractured Reservoirs with Dual Porosity Models
Simulation of Naturally Fractured Reservoirs ith Dual Porosity Models Pallav Sarma Prof. Khalid Aziz Stanford University SUPRI-HW Motivations! NFRs represent at least 0%of orld reserves, but difficult
More informationPOD for Parametric PDEs and for Optimality Systems
POD for Parametric PDEs and for Optimality Systems M. Kahlbacher, K. Kunisch, H. Müller and S. Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria DMV-Jahrestagung 26,
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More informationSolving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Solving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics Fikriye Yılmaz 1, Bülent Karasözen
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationNumerical Optimization of Partial Differential Equations
Numerical Optimization of Partial Differential Equations Part I: basic optimization concepts in R n Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada
More informationContinuum Modeling of Transportation Networks with Differential Equations
with Differential Equations King Abdullah University of Science and Technology Thuwal, KSA Examples of transportation networks The Silk Road Examples of transportation networks Painting by Latifa Echakhch
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationOn the well-posedness of the Prandtl boundary layer equation
On the well-posedness of the Prandtl boundary layer equation Vlad Vicol Department of Mathematics, The University of Chicago Incompressible Fluids, Turbulence and Mixing In honor of Peter Constantin s
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
More informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationElectromagnetically Induced Flows in Water
Electromagnetically Induced Flows in Water Michiel de Reus 8 maart 213 () Electromagnetically Induced Flows 1 / 56 Outline 1 Introduction 2 Maxwell equations Complex Maxwell equations 3 Gaussian sources
More informationA NUMERICAL ALGORITHM FOR SET-POINT REGULATION OF NON-LINEAR PARABOLIC CONTROL SYSTEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 11, Number 1, Pages 54 85 c 2014 Institute for Scientific Computing and Information A NUMERICAL ALGORITHM FOR SET-POINT REGULATION OF NON-LINEAR
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 16.323 Principles of Optimal Control Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 16.323 Lecture
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationMaximum Rate of Growth of Enstrophy in Solutions of the Fractional Burgers Equation
Maximum Rate of Growth of Enstrophy in Solutions of the Fractional Burgers Equation arxiv:161.9578v [physics.flu-dyn] 15 Jan 18 Dongfang Yun and Bartosz Protas Department of Mathematics and Statistics,
More informationDecay profiles of a linear artificial viscosity system
Decay profiles of a linear artificial viscosity system Gene Wayne, Ryan Goh and Roland Welter Boston University rwelter@bu.edu July 2, 2018 This research was supported by funding from the NSF. Roland Welter
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationUniversity of North Carolina at Charlotte Charlotte, USA Norwegian Institute of Science and Technology Trondheim, Norway
1 A globally convergent numerical method for some coefficient inverse problems Michael V. Klibanov, Larisa Beilina University of North Carolina at Charlotte Charlotte, USA Norwegian Institute of Science
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationDRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS
DRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS ADILBEK KAIRZHAN, DMITRY E. PELINOVSKY, AND ROY H. GOODMAN Abstract. When the coefficients of the cubic terms match the coefficients in the boundary
More informationIntroduction to Fourier Analysis
Lecture Introduction to Fourier Analysis Jan 7, 2005 Lecturer: Nati Linial Notes: Atri Rudra & Ashish Sabharwal. ext he main text for the first part of this course would be. W. Körner, Fourier Analysis
More informationHorizontal buoyancy-driven flow along a differentially cooled underlying surface
Horizontal buoyancy-driven flow along a differentially cooled underlying surface By Alan Shapiro and Evgeni Fedorovich School of Meteorology, University of Oklahoma, Norman, OK, USA 6th Baltic Heat Transfer
More informationInverse Gravimetry Problem
Inverse Gravimetry Problem Victor Isakov September 21, 2010 Department of M athematics and Statistics W ichita State U niversity W ichita, KS 67260 0033, U.S.A. e mail : victor.isakov@wichita.edu 1 Formulation.
More informationFluid dynamics of electrons in graphene. Andrew Lucas
Fluid dynamics of electrons in graphene Andrew Lucas Stanford Physics Condensed Matter Seminar, Princeton October 17, 2016 Collaborators 2 Subir Sachdev Harvard Physics & Perimeter Institute Philip Kim
More informationChapter 2 Distributed Parameter Systems: Controllability, Observability, and Identification
Chapter 2 Distributed Parameter Systems: Controllability, Observability, and Identification 2.1 Mathematical Description We introduce the class of systems to be considered in the framework of this monograph
More informationControllability of linear PDEs (I): The wave equation
Controllability of linear PDEs (I): The wave equation M. González-Burgos IMUS, Universidad de Sevilla Doc Course, Course 2, Sevilla, 2018 Contents 1 Introduction. Statement of the problem 2 Distributed
More informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationHydrodynamic Modes of Incoherent Black Holes
Hydrodynamic Modes of Incoherent Black Holes Vaios Ziogas Durham University Based on work in collaboration with A. Donos, J. Gauntlett [arxiv: 1707.xxxxx, 170x.xxxxx] 9th Crete Regional Meeting on String
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationPart 5: Penalty and augmented Lagrangian methods for equality constrained optimization. Nick Gould (RAL)
Part 5: Penalty and augmented Lagrangian methods for equality constrained optimization Nick Gould (RAL) x IR n f(x) subject to c(x) = Part C course on continuoue optimization CONSTRAINED MINIMIZATION x
More informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
More information5 Handling Constraints
5 Handling Constraints Engineering design optimization problems are very rarely unconstrained. Moreover, the constraints that appear in these problems are typically nonlinear. This motivates our interest
More information9. Boundary layers. Flow around an arbitrarily-shaped bluff body. Inner flow (strong viscous effects produce vorticity) BL separates
9. Boundary layers Flow around an arbitrarily-shaped bluff body Inner flow (strong viscous effects produce vorticity) BL separates Wake region (vorticity, small viscosity) Boundary layer (BL) Outer flow
More informationMumford-Shah Level-Set Ideas for Problems in Medical Imaging Paris Dauphine 2005
1/26 Mumford-Shah Level-Set Ideas for Problems in Medical Imaging Paris Dauphine 2005 Wolfgang Ring, Marc Droske, Ronny Ramlau Institut für Mathematik, Universität Graz General Mumford-Shah like 2/26 functionals
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More information7 EQUATIONS OF MOTION FOR AN INVISCID FLUID
7 EQUATIONS OF MOTION FOR AN INISCID FLUID iscosity is a measure of the thickness of a fluid, and its resistance to shearing motions. Honey is difficult to stir because of its high viscosity, whereas water
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationTwo-Scale Wave Equation Modeling for Seismic Inversion
Two-Scale Wave Equation Modeling for Seismic Inversion Susan E. Minkoff Department of Mathematics and Statistics University of Maryland Baltimore County Baltimore, MD 21250, USA RICAM Workshop 3: Wave
More informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
More informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
More informationSHAPE OPTIMIZATION OF ELLIPTIC PDE PROBLEMS ON COMPLEX DOMAINS
SHAPE OPTIMIZATION OF ELLIPTIC PDE PROBLEMS ON COMPLEX DOMAINS SHAPE OPTIMIZATION OF ELLIPTIC PDE PROBLEMS ON COMPLEX DOMAINS By KATSIARYNA NIAKHAI, M.Sc. A Thesis Submitted to the School of Graduate Studies
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationA collocation method for solving some integral equations in distributions
A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm
More informationFully implicit higher-order schemes applied to polymer flooding. K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF)
Fully implicit higher-order schemes applied to polymer flooding K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF) Motivation Challenging and complex to simulate water-based EOR methods, unresolved
More informationAdjoint Based Optimization of PDE Systems with Alternative Gradients
Adjoint Based Optimization of PDE Systems with Alternative Gradients Bartosz Protas a a Department of Mathematics & Statistics, McMaster University, Hamilton, Ontario, Canada Abstract In this work we investigate
More informationPDE-Constrained and Nonsmooth Optimization
Frank E. Curtis October 1, 2009 Outline PDE-Constrained Optimization Introduction Newton s method Inexactness Results Summary and future work Nonsmooth Optimization Sequential quadratic programming (SQP)
More informationSUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description of many particle
More informationComputation of operators in wavelet coordinates
Computation of operators in wavelet coordinates Tsogtgerel Gantumur and Rob Stevenson Department of Mathematics Utrecht University Tsogtgerel Gantumur - Computation of operators in wavelet coordinates
More informationPreliminary Exam 2018 Solutions to Morning Exam
Preliminary Exam 28 Solutions to Morning Exam Part I. Solve four of the following five problems. Problem. Consider the series n 2 (n log n) and n 2 (n(log n)2 ). Show that one converges and one diverges
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationSome Semi-Markov Processes
Some Semi-Markov Processes and the Navier-Stokes equations NW Probability Seminar, October 22, 2006 Mina Ossiander Department of Mathematics Oregon State University 1 Abstract Semi-Markov processes were
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationWaveform inversion and time-reversal imaging in attenuative TI media
Waveform inversion and time-reversal imaging in attenuative TI media Tong Bai 1, Tieyuan Zhu 2, Ilya Tsvankin 1, Xinming Wu 3 1. Colorado School of Mines 2. Penn State University 3. University of Texas
More informationPDE Methods for Mean Field Games with Non-Separable Hamiltonian: Data in Sobolev Spaces (Continued) David Ambrose June 29, 2018
PDE Methods for Mean Field Games with Non-Separable Hamiltonian: Data in Sobolev Spaces Continued David Ambrose June 29, 218 Steps of the energy method Introduce an approximate problem. Prove existence
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationFundamentals of Fluid Dynamics: Waves in Fluids
Fundamentals of Fluid Dynamics: Waves in Fluids Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/ tzielins/ Institute
More informationPDE Constrained Optimization selected Proofs
PDE Constrained Optimization selected Proofs Jeff Snider jeff@snider.com George Mason University Department of Mathematical Sciences April, 214 Outline 1 Prelim 2 Thms 3.9 3.11 3 Thm 3.12 4 Thm 3.13 5
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationParameter Estimation in Reservoir Engineering Models via Data Assimilation Techniques
Parameter Estimation in Reservoir Engineering Models via Data Assimilation Techniques Mariya V. Krymskaya TU Delft July 6, 2007 Ensemble Kalman Filter (EnKF) Iterative Ensemble Kalman Filter (IEnKF) State
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationComparison Results for Semilinear Elliptic Equations in Equimeasurable Domains
Comparison Results for Semilinear Elliptic Equations in Equimeasurable Domains François HAMEL Aix-Marseille University & Institut Universitaire de France In collaboration with Emmanuel RUSS Workshop on
More information1 Constant Pressure Inner Boundary Condition
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 8 Notes These notes are based on Rosalind Archer s PE281 lecture notes, with some revisions by Jim Lambers. 1 Constant Pressure Inner Boundary Condition
More informationA Sobolev trust-region method for numerical solution of the Ginz
A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012
More informationA globally convergent numerical method and adaptivity for an inverse problem via Carleman estimates
A globally convergent numerical method and adaptivity for an inverse problem via Carleman estimates Larisa Beilina, Michael V. Klibanov Chalmers University of Technology and Gothenburg University, Gothenburg,
More informationThermal Detection of Inaccessible Corrosion
Thermal Detection of Inaccessible Corrosion Matt Charnley and Andrew Rzeznik October 18, 212 Abstract In this paper, we explore the mathematical inverse problem of detecting corroded material on the reverse
More information5. A step beyond linearization: velocity analysis
5. A step beyond linearization: velocity analysis 1 Partially linearized seismic inverse problem ( velocity analysis ): given observed seismic data S obs, find smooth velocity v E(X), X R 3 oscillatory
More informationApplied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces
Applied Numerical Mathematics 111 (2017) 64 91 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum High-order numerical schemes based on difference potentials
More informationPenalty and Barrier Methods General classical constrained minimization problem minimize f(x) subject to g(x) 0 h(x) =0 Penalty methods are motivated by the desire to use unconstrained optimization techniques
More informationI. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010
: an : an (Joint work with A. Fasano) Dipartimento di Matematica U. Dini, Università di Firenze (Italy) borsi@math.unifi.it http://web.math.unifi.it/users/borsi porous EMS SCHOOL ON INDUSTRIAL MATHEMATICS
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationDissipative quasi-geostrophic equations with L p data
Electronic Journal of Differential Equations, Vol. (), No. 56, pp. 3. ISSN: 7-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Dissipative quasi-geostrophic
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
More information