A Multiphysics Framework Using hp-finite Elements for Electromagnetics Applications
|
|
- Dulcie Montgomery
- 5 years ago
- Views:
Transcription
1 Univ. Carlos III of Madrid A Multiphysics Framework Using hp-finite Elements for Electromagnetics Applications D. Pardo, L. E. García-Castillo, M. J. Nam, C. Torres-Verdín Team: D. Pardo, M. J. Nam, V. Calo, L.E. García-Castillo, I. Garay, M. Paszynski, P. Matuszyk, L. Demkowicz, C. Torres-Verdín March 13th, 2009
2 overview 1. Motivation and Objectives: Joint Multiphysics Inversion 2. Simulation of Forward Problems Parallel Self-Adaptive Goal-Oriented hp-finite Element Method Electromagnetic Applications Sonic Applications 3. Inversion Library (Work in Progress) h-adaptive Newton s Method Implementation 4. Conclusions and Future Work
3 motivation and objectives Seismic Measurements Figure from the USGS Science Center for Coastal and Marine Geology 2
4 motivation and objectives Marine Controlled-Source Electromagnetics (CSEM) Figure from the UCSD Institute of Oceanography 3
5 motivation and objectives Multiphysics Logging Measurements OBJECTIVES: To determine payzones (porosity), amount of oil/gas (saturation), and ability to extract oil/gas (permeability). 4
6 motivation and objectives Main Objective: To Solve a Multiphysics Inverse Problem Software to solve the DIRECT problem is essential in order to solve the INVERSE problem. 5
7 simulation of forward problems (hp-fem) The h-finite Element Method 1. Convergence limited by the polynomial degree, and large material contrasts. 2. Optimal h-grids do NOT converge exponentially in real applications. 3. They may lock (100% error). The p-finite Element Method 1. Exponential convergence feasible for analytical ( nice ) solutions. 2. Optimal p-grids do NOT converge exponentially in real applications. 3. If initial h-grid is not adequate, the p-method will fail miserably. The hp-finite Element Method 1. Exponential convergence feasible for ALL solutions. 2. Optimal hp-grids DO converge exponentially in real applications. 3. If initial hp-grid is not adequate, results will still be great. 6
8 simulation of forward problems (hp-fem) Energy norm based fully automatic hp-adaptive strategy Coarse grids Fine grids (hp) (h/2, p + 1) global hp-refinement global hp-refinement SOL. METHOD ON FINE GRIDS: A TWO GRID SOLVER 7
9 simulation of forward problems (hp-fem) Motivation (Goal-Oriented Adaptivity) Test Problem Infinite Domain Resistivity = 1 Ohm m Frequency = 2 Mhz 50 meters Transmitter (T) Receiver (R) 8
10 simulation of forward problems (hp-fem) Motivation (Goal-Oriented Adaptivity) Infinite Domain Resistivity = 1 Ohm m Frequency = 2 Mhz Solution decays exponentially. E(T) E(R) 1060 Results using energy-norm adaptivity: Transmitter (T) 50 meters Receiver (R) Energy-norm error: 0.001% Relative error in the quantity of interest > %. Goal-oriented adaptivity is needed. Becker-Rannacher (1995,1996), Rannacher-Stuttmeier (1997), Cirak-Ramm (1998), Paraschivoiu-Patera (1998), Peraire-Patera (1998), Prudhomme-Oden (1999, 2001), Heuveline-Rannacher (2003), Solin-Demkowicz (2004). 9
11 simulation of forward problems (hp-fem) Motivation (Goal-Oriented Adaptivity) Infinite Domain Resistivity = 1 Ohm m Frequency = 2 Mhz Solution decays exponentially. E(T) E(R) 1060 Results using energy-norm adaptivity: Transmitter (T) 50 meters Receiver (R) Energy-norm error: 0.001% Relative error in the quantity of interest > %. Goal-oriented adaptivity is needed. Becker-Rannacher (1995,1996), Rannacher-Stuttmeier (1997), Cirak-Ramm (1998), Paraschivoiu-Patera (1998), Peraire-Patera (1998), Prudhomme-Oden (1999, 2001), Heuveline-Rannacher (2003), Solin-Demkowicz (2004). 10
12 simulation of forward problems (hp-fem) Motivation (Goal-Oriented Adaptivity) Transmitter (T) Infinite Domain Resistivity = 1 Ohm m Frequency = 2 Mhz 50 meters Receiver (R) Distance in z axis from transmitter to receiver (in m) Analytical solution Numerical solution Distance in z axis from transmitter to receiver (in m) Analytical solution Numerical solution Amplitude (V/m) Goal-oriented adaptivity is needed Phase (degrees) 11
13 simulation of forward problems (hp-fem) Mathematical Formulation (Goal-Oriented Adaptivity) We consider the following problem (in variational form): { Find L(Ψ), where Ψ V such that : b(ψ, ξ) = f(ξ) ξ V. We define residual r e (ξ) = b(e, ξ). We seek for solution G of: { Find G V V such that : G(r e ) = L(e). This is necessarily solved if we find the solution of the dual problem: { Find G V such that : Notice that L(e) = b(e, G). b(ψ, G) = L(Ψ) Ψ V. 12
14 simulation of forward problems (hp-fem) Algorithm for Goal-Oriented Adaptivity Solve DIRECT and DUAL problems on Grid hp. Solve DIRECT and DUAL problems on Grid h/2, p + 1. Compute e = Ψ h/2,p+1 Ψ hp, and ǫ = G h/2,p+1 G hp. Represent the error as: L(e) = b(e, ǫ) K b K(e, ǫ). Apply the fully automatic hp-adaptive algorithm. Solve DIRECT and DUAL problems on Grid hp. Solve DIRECT and DUAL problems on Grid h/2, p
15 simulation of forward problems (hp-fem) Axisymmetric Logging-While-Drilling (LWD) Simulation ENERGY-NORM HP-ADAPTIVITY Dhp90: A Fully automatic hp-adaptive Finite Element code p=8 p=7 p=6 p=5 p=4 p=3 p=2 p=
16 simulation of forward problems (hp-fem) Axisymmetric Logging-While-Drilling (LWD) Simulation GOAL-ORIENTED HP-ADAPTIVITY (Quadrilateral Elements) Dhp90: A Fully automatic hp-adaptive Finite Element code p=8 p=7 p=6 p=5 p=4 p=3 p=2 p=
17 simulation of forward problems (hp-fem) Axisymmetric Logging-While-Drilling (LWD) Simulation GOAL-ORIENTED HP-ADAPTIVITY (Zoom towards first receiver) Dhp90: A Fully automatic hp-adaptive Finite Element code p=8 p=7 p=6 p=5 p=4 p=3 p=2 p= E
18 simulation of forward problems (hp-fem) Non-Orthogonal System of Coordinates ζ 3 ζ 2 7 ζ 1 Fourier Series Expansion in ζ 2 DC Problems: σ u = f u(ζ 1, ζ 2, ζ 3 ) = σ(ζ 1, ζ 2, ζ 3 ) = f(ζ 1, ζ 2, ζ 3 ) = l= l= m= m= n= n= u l (ζ 1, ζ 3 )e jlζ 2 σ m (ζ 1, ζ 3 )e jmζ 2 f n (ζ 1, ζ 3 )e jnζ 2 Fourier modes e jlζ 2 are orthogonal high-order basis functions that are (almost) invariant with respect to the gradient operator. 17
19 simulation of forward problems (hp-fem) De Rham diagram De Rham diagram is critical to the theory of FE discretizations of multi-physics problems. IR W Q V id Π L 2 0 Π curl Π div P IR W p Q p V p W p 1 0. This diagram relates two exact sequences of spaces, on both continuous and discrete levels, and corresponding interpolation operators. 18
20 simulation of forward problems (hp-fem) Tri-Axial Induction Tool 19
21 simulation of forward problems (hp-fem) Tri-Axial Induction Tools in Deviated Wells (0, 30, and 60 degrees) 6 5 Real part of Hzz at 20 khz vertical 30 degrees 60 degrees 100 ohm m 6 5 Real part of Hxx at 20 khz vertical 30 degrees 60 degrees 100 ohm m 6 5 Real part of Hyy at 20 khz vertical 30 degrees 60 degrees 100 ohm m ohm m ohm m ohm m Depth (m) ohm m Depth (m) ohm m Depth (m) ohm m ohm m 2 1 ohm m 2 1 ohm m ohm m 4 20 ohm m 4 20 ohm m Re(Hzz) field (A/m) x 10 6 Re(Hxx) field (A/m) Re(Hyy) field (A/m) Triaxial tools are more sensitive to dip angle effects 20
22 simulation of forward problems (hp-fem) Tri-Axial Induction Tools in Deviated Wells (0, 30, and 60 degrees) 6 5 Imaginary part of Hzz at 20 khz 100 ohm m vertical 30 degrees 60 degrees 6 5 Imaginary part of Hxx at 20 khz 100 ohm m vertical 30 degrees 60 degrees 6 5 Imaginary part of Hyy at 20 khz 100 ohm m vertical 30 degrees 60 degrees ohm m ohm m ohm m Depth (m) ohm m Depth (m) ohm m Depth (m) ohm m ohm m 2 1 ohm m 2 1 ohm m ohm m 4 20 ohm m 4 20 ohm m Im(Hzz) field (A/m) x 10 5 Im(Hxx) field (A/m) Im(Hyy) field (A/m) Triaxial tools are more sensitive to dip angle effects 21
23 simulation of forward problems (hp-fem) DLL Tool 22
24 simulation of forward problems (hp-fem) Groningen Effect 30 Comparison of Groningen Effects on LLd at 100 Hz 30 Comparison of Groningen Effects on LLd at 100 Hz 15 Casing 15 Casing Relative Depth z (m) 0 The earth surface is at z = 2 km or z = 200 m B: at ρ = 1 m, 1 m below the surface Relative Depth z (m) 0 The earth surface is at z = 2 km B: 1 m below the surface, at ρ = 1 m or ρ = 200 m Surface: at z = 2 km Surface: at z = 200 m DC Apparent Resistivity (Ω m) 60 B at ρ = 1 m B at ρ = 200 m DC Apparent Resistivity (Ω m) As we place the current return electrode B farther from the logging instrument, the Groningen effect diminishes 23
25 simulation of forward problems (hp-fem) DC DLL in Deviated Wells 4 2 LLd: Iso LLd: Ani LLs: Iso LLs: Ani Anisotropic Formation (Vertical Well) 4 2 Anisotropic Formation (60 degree Deviated Well) LLd: 60, Iso LLd: 60, Ani LLs: 60, Iso LLs: 60, Ani 0 Resistivity of Formation 0 Resistivity of Formation Relative Depth (m) Vertical Resistivity Vertical Resistivity Relative Depth (m) Vertical Resistivity Vertical Resistivity Apparent Resistivity (Ω m) Apparent Resistivity (Ω m) Anisotropy is better identified when using deviated wells 24
26 simulation of forward problems (hp-fem) Final hp-grid and solution 8 KHz, acoustics, open borehole setting (no logging instrument). 25
27 new library for inverse problems Notation: Variational Formulation (DC) B(u, v; σ) =< v, σ u > L 2 (Ω) (bilinear u, v) F i (v) =< v, f i > L 2 (Ω) + < v, g i > L 2 ( Ω) (linear v) L i (u) =< l i, u > L 2 (Ω) + < h i, u > L 2 ( Ω) (linear u) Direct Problem (homogeneous Dirichlet BC s): { Find û i V such that : Dual (Adjoint) Problem: { B(û i, v; σ) = F i (v) v V Find ˆv i V such that : B(u, ˆv i ; σ) = L i (u) u V 26
28 new library for inverse problems Notation: Variational Formulation (AC) B(E, F; σ) =< F, µ 1 E > L 2 (Ω) < F, (ω 2 ǫ jωσ)e > L 2 (Ω) F i (F) = jω < F, J imp i > L 2 (Ω) +jω < F, J imp S,i > L 2 ( Ω) L i (E) =< J adj i, E > L 2 (Ω) + < J adj S,i, E > L 2 ( Ω) Direct Problem (homogeneous Dirichlet BC s): { Find Ê i W such that : Dual (Adjoint) Problem: { B(Ê i, F; σ) = F i (F) F W Find ˆF i W such that : B(E, ˆF i ; σ) = L i (E) E W 27
29 new library for inverse problems Constrained Nonlinear Optimization Problem Cost Functional: { Find σ > 0 such that it minimizes C β (σ), where : C β (σ) = W m (L(ûσ) M) 2 l 2 + β R(σ σ 0 ) 2 L 2, where M i denotes the i-th measurement, M = (M 1,..., M n ) L i is the i-th quantity of interest, L = (L 1,..., L n ) n M 2 l 2 = M 2 i ; R(σ σ 0 ) 2 L 2 = (R(σ σ 0 )) 2 i=1 β is the relaxation parameter, σ 0 is given, W m are weights Main objective (inversion problem): Find ˆσ = min σ>0 C β(σ) 28
30 new library for inverse problems Solving a Constrained Nonlinear Optimization Problem We select the following deterministic iterative method: σ (n+1) = σ (n) + α (n) δσ (n) How to find a search direction δσ (n)? We will employ a change of coordinates and a truncated Taylor s series expansion. How to determine the step size α (n)? Either with a fixed size or using an approximation for computing L(σ (n) + α (n) δσ (n) ). How to guarantee that the nonlinear constraints will be satisfied? Imposing the Karush-Kuhn-Tucker (KKT) conditions or with a penalization method, or via a change of variables. 29
31 new library for inverse problems Change of coordinates: Search Direction Method h(s) = σ => Find ŝ = min h(s)>0 C β(s) Taylor s series expansion: A) C β (s + δs) C β (s) + δs C β (s) + 0.5δs 2 H Cβ (s) B) L(s + δs) L(s) + δs L(s), R(s + δs) = R(s) + δs R(s) Expansion A) leads to the Newton-Raphson method. Expansion B) leads to the Gauss-Newton method. Expansion A) with H Cβ = I leads to the steepest descent method. Higher-order expansions require from higher-order derivatives. 30
32 new library for inverse problems Computation of Jacobian Matrix Using the Fréchet Derivative: ( ) L i (û i ) ûi = B, ˆv i, h(s) +B s j s j ( ) ( ) ûi ûi L i = B, ˆv i, h(s) s j s j ( ) ˆvi F i = B s j ( û i, ˆv ) i, h(s) s j ( û i, ˆv ) i, h(s) s j +B ( û i, ˆv i, h(s) ) s j Therefore, we conclude: Jacobian Matrix = L i(û i ) s j = B ( û i, ˆv i, h(s) ) s j 31
33 new library for inverse problems Jacobian Function: One TX, one RX 32
34 new library for inverse problems Jacobian Function: One TX, one RX 33
35 new library for inverse problems Jacobian Function: One TX, one RX 34
36 new library for inverse problems Computation of Hessian Matrix Following a similar argument as for the Jacobian matrix, we obtain: 2 ( L i (û i ) ûi = B, ˆv i, h(s) ) ( B û i, ˆv i, h(s) ) ( ) B û i, ˆv i, 2 h(s) s j s k s j s k s j s k s j s k How do we compute û i and ˆv i? s j s j Find û ( ) i ûi such that : B, v i, h(s) s j s j Find ˆv i s j such that : B ( ˆvi s j, u i, h(s) ) = B = B ( û i, v i, h(s) ) s j ( ˆv i, u i, h(s) ) s j We can compute the Hessian matrix EXACTLY by just solving our original problem for different right-hand-sides, and performing additional integrations. v i u i 35
37 new library for inverse problems Hessian Function: One TX, one RX 36
38 new library for inverse problems Hessian Function: One TX, two RXs 37
39 new library for inverse problems Hessian Function: One TX, three RXs 38
40 new library for inverse problems Algorithms implemented within the inverse library JOINT MULTI PHYSICS INVERSION LIBRARY CONSTRAINED OPTIMIZATION SEARCH DIRECTIONS STEP SIZE 1. KKT Conditions 1. Steepest Descent 2. Penalization Method 2. Gauss Newton 3. Change of coordinates 3. Newton Raphson 1. Uniform step size 2. Variable step size (multiple algorithms) The inverse library is composed of multiple algorithms for imposing constraints, and finding search directions and corresponding step sizes. Jacobian and Hessian matrices are computed exactly by simply solving the dual (adjoint) formulation and performing additional integrations. The inverse library is compatible with multi-physics problems. 39
41 conclusions We have described an efficient numerical method for solving PDE s based on a self-adaptive goal-oriented hp refinement strategy. We are developing a multiphysics version of the code. We are building a new module for the resolution of inverse problems. Our main objective is to create a software infrastructure enabling solution of challenging multiphysics inverse problems with applications to geophysics (hydrocarbon detection and monitoring, etc.), aeronautics and medicine. To achieve this objective, we need Ph.D. students, post-doctoral fellows, experienced researchers, and collaborators in different areas (inversion, solvers, etc). 40
Multiphysics, Inversion, and Petroleum
CIC biogune Research Line II: Multiphysics, Inversion, and Petroleum Research Professor at BCAM Team: D. Pardo, I. Garay, I. Andonegui, J. Álvarez Collaborators: P. de la Hoz, M. Paszynski, L.E. García-Castillo,
More informationSIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS
SIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS hp-fem team: D. Pardo, M. J. Nam, L. Demkowicz, C. Torres-Verdín, V. M. Calo, M. Paszynski, and P. J. Matuszyk
More informationFast Inversion of Logging-While-Drilling (LWD) Resistivity Measurements
Fast Inversion of Logging-While-Drilling (LWD) Resistivity Measurements David Pardo 1 Carlos Torres-Verdín 2 1 University of the Basque Country (UPV/EHU) and Ikerbasque, Bilbao, Spain. 2 The University
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationA Nested Dissection Parallel Direct Solver. for Simulations of 3D DC/AC Resistivity. Measurements. Maciej Paszyński (1,2)
A Nested Dissection Parallel Direct Solver for Simulations of 3D DC/AC Resistivity Measurements Maciej Paszyński (1,2) David Pardo (2), Carlos Torres-Verdín (2) (1) Department of Computer Science, AGH
More informationSimulation of 3D DC Borehole Resistivity Measurements with a Goal- Oriented hp Finite-Element Method. Part I: Laterolog and LWD
Journal of the Serbian Society for Computational Mechanics / Vol. 1 / No. 1, 2007 / pp. 62-73 Simulation of 3D DC Borehole Resistivity Measurements with a Goal- Oriented hp Finite-Element Method. Part
More informationSIMULATION OF RESISTIVITY LOGGING-WHILE-DRILLING (LWD) MEASUREMENTS USING A SELF-ADAPTIVE GOAL-ORIENTED HP FINITE ELEMENT METHOD.
SIMULATION OF RESISTIVITY LOGGING-WHILE-DRILLING (LWD) MEASUREMENTS USING A SELF-ADAPTIVE GOAL-ORIENTED HP FINITE ELEMENT METHOD. D. PARDO,, L. DEMKOWICZ, C. TORRES-VERDÍN, AND M. PASZYNSKI, Abstract.
More informationSimulation of Borehole Resistivity Tools Through Metal Casing at Different Frequencies Using a Goal Oriented hp Finite Element Method (FEM).
Simulation of Borehole Resistivity Tools Through Metal Casing at Different Frequencies Using a Goal Oriented hp Finite Element Method (FEM). D. Pardo, C.Torres-Verdín, L. Demkowicz Abstract We simulate
More informationSimulation of marine controlled source electromagnetic measurements using a parallel fourier hp-finite element method
Comput Geosci (2011 15:53 67 DOI.07/s596-0-9195-1 ORIGINAL PAPER Simulation of marine controlled source electromagnetic measurements using a parallel fourier hp-finite element method David Pardo Myung
More informationSimulation of DC dual-laterolog measurements in complex formations: a Fourier series
Simulation of DC dual-laterolog measurements in complex formations: a Fourier series approach with non-orthogonal coordinates and self-adapting finite elements Myung Jin Nam 1, David Pardo *, and Carlos
More informationFourier Series Expansion in a Non-Orthogonal System of Coordinates for the Simulation of 3D Alternating Current Borehole Resistivity Measurements
Fourier Series Expansion in a Non-Orthogonal System of Coordinates for the Simulation of 3D Alternating Current Borehole Resistivity Measurements D. Pardo[a], C. Torres-Verdín[a], M. J. Nam[a], M. Paszynski[b],
More informationFourier Series Expansion in a Non-Orthogonal System of Coordinates for the Simulation of 3D AC Borehole Resistivity Measurements
Fourier Series Expansion in a Non-Orthogonal System of Coordinates for the Simulation of 3D AC Borehole Resistivity Measurements D. Pardo[a], C. Torres-Verdín[a], M. J. Nam[a], M. Paszynski[b], and V.
More informationPML ENHANCED WITH A SELF-ADAPTIVE GOAL-ORIENTED hp-finite ELEMENT METHOD: SIMULATION OF THROUGH-CASING BOREHOLE RESISTIVITY MEASUREMENTS
SIAM J. SCI. COMPUT. Vol. 0, No. 0, pp. 000 000 c XXXX Society for Industrial and Applied Mathematics PML ENHANCED WITH A SELF-ADAPTIVE GOAL-ORIENTED hp-finite ELEMENT METHOD: SIMULATION OF THROUGH-CASING
More informationMultiscale Analysis of Vibrations of Streamers
Multiscale Analysis of Vibrations of Streamers Leszek Demkowicz Joint work with S. Prudhomme, W. Rachowicz, W. Qiu and L. Chamoin Institute for Computational Engineering and Sciences (ICES) The University
More informationA goal-oriented hp-adaptive finite element method with electromagnetic applications. Part I: Electrostatics
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2006; 65:1269 1309 Published online 19 September 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.1488
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 66, No. 6, pp. 2085 206 c 2006 Society for Industrial and Applied Mathematics TWO-DIMENSIONAL HIGH-ACCURACY SIMULATION OF RESISTIVITY LOGGING-WHILE-DRILLING (LWD) MEASUREMENTS
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More informationThe Normal Equations. For A R m n with m > n, A T A is singular if and only if A is rank-deficient. 1 Proof:
Applied Math 205 Homework 1 now posted. Due 5 PM on September 26. Last time: piecewise polynomial interpolation, least-squares fitting Today: least-squares, nonlinear least-squares The Normal Equations
More informationSimulation based optimization
SimBOpt p.1/52 Simulation based optimization Feb 2005 Eldad Haber haber@mathcs.emory.edu Emory University SimBOpt p.2/52 Outline Introduction A few words about discretization The unconstrained framework
More informationNonlinear Optimization: What s important?
Nonlinear Optimization: What s important? Julian Hall 10th May 2012 Convexity: convex problems A local minimizer is a global minimizer A solution of f (x) = 0 (stationary point) is a minimizer A global
More informationFast inversion of 3D Borehole Resistivity Measurements using Model Reduction Techniques based on 1D Semi-Analytical Solutions.
Fast inversion of 3D Borehole Resistivity Measurements using Model Reduction Techniques based on 1D Semi-Analytical Solutions. Aralar Erdozain Supervised by: Hélène Barucq, David Pardo and Victor Péron
More informationJoint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research
Joint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research SUMMARY New-generation electromagnetic (EM) and sonic logging tools are
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationICES REPORT Solution of Coupled Acoustic-Elastic Wave Propagation Problems with Anelastic Attenuation Using Automatic
ICES REPORT 11-03 January 2011 Solution of Coupled Acoustic-Elastic Wave Propagation Problems with Anelastic Attenuation Using Automatic by P. J. Matuszyk, L. F. Demkowicz, C. Torres-Verdin The Institute
More informationApproximate- vs. full-hessian in FWI: 1D analytical and numerical experiments
Approximate- vs. full-hessian in FWI: 1D analytical and numerical experiments Raul Cova and Kris Innanen ABSTRACT Feasibility of using Full Waveform Inversion (FWI) to build velocity models has been increasing
More informationMyung Jin Nam 1, David Pardo 2, and Carlos Torres-Verdín 1 ABSTRACT INTRODUCTION
GEOPHYSICS, VOL. 7, NO. 1 JANUARY-FEBRUARY 9 ; P. E31 E3, 13 FIGS., 1 TABLE. 1.119/1.3681 Simulation of DC dual-laterolog measurements in complex formations: A Fourier-series approach with nonorthogonal
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 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 informationOptimization Methods
Optimization Methods Decision making Examples: determining which ingredients and in what quantities to add to a mixture being made so that it will meet specifications on its composition allocating available
More informationNon-linear least squares
Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined
More informationContents. Preface. 1 Introduction Optimization view on mathematical models NLP models, black-box versus explicit expression 3
Contents Preface ix 1 Introduction 1 1.1 Optimization view on mathematical models 1 1.2 NLP models, black-box versus explicit expression 3 2 Mathematical modeling, cases 7 2.1 Introduction 7 2.2 Enclosing
More informationA new inversion method for dissipating electromagnetic problem
A new inversion method for dissipating electromagnetic problem Elena Cherkaev Abstract The paper suggests a new technique for solution of the inverse problem for Maxwell s equations in a dissipating medium.
More informationAchieving depth resolution with gradient array survey data through transient electromagnetic inversion
Achieving depth resolution with gradient array survey data through transient electromagnetic inversion Downloaded /1/17 to 128.189.118.. Redistribution subject to SEG license or copyright; see Terms of
More informationG008 Advancing Marine Controlled Source Electromagnetics in the Santos Basin, Brazil
G008 Advancing Marine Controlled Source Electromagnetics in the Santos Basin, Brazil M.P. Buonora* (Petrobras), A. Zerilli (WesternGeco Electromagnetics), T. Labruzzo (WesternGeco Electromagnetics) & L.F.
More informationConstrained Optimization
1 / 22 Constrained Optimization ME598/494 Lecture Max Yi Ren Department of Mechanical Engineering, Arizona State University March 30, 2015 2 / 22 1. Equality constraints only 1.1 Reduced gradient 1.2 Lagrange
More informationNon-polynomial Least-squares fitting
Applied Math 205 Last time: piecewise polynomial interpolation, least-squares fitting Today: underdetermined least squares, nonlinear least squares Homework 1 (and subsequent homeworks) have several parts
More informationAnisotropic meshes for PDEs: a posteriori error analysis and mesh adaptivity
ICAM 2013, Heraklion, 17 September 2013 Anisotropic meshes for PDEs: a posteriori error analysis and mesh adaptivity Simona Perotto MOX-Modeling and Scientific Computing Department of Mathematics, Politecnico
More informationNumerical optimization
Numerical optimization Lecture 4 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 2 Longest Slowest Shortest Minimal Maximal
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon s method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationNumerical optimization. Numerical optimization. Longest Shortest where Maximal Minimal. Fastest. Largest. Optimization problems
1 Numerical optimization Alexander & Michael Bronstein, 2006-2009 Michael Bronstein, 2010 tosca.cs.technion.ac.il/book Numerical optimization 048921 Advanced topics in vision Processing and Analysis of
More informationSpring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014)
Spring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014) The first problem in this assignment is a paper-and-pencil exercise
More informationMulti-order Vector Finite Element Modeling of 3D Magnetotelluric Data Including Complex Geometry and Anisotropy. Aixa M.
Multi-order Vector Finite Element Modeling of 3D Magnetotelluric Data Including Complex Geometry and Anisotropy Aixa M. Rivera Ríos Thesis submitted for the degree of Doctor of Philosophy in The University
More informationEstimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS
Estimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS Summary Towed streamer EM data was acquired in October
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationFinite Element Simulations of Logging-While-Drilling and Extra-Deep Azimuthal Resistivity Measurements using Non-Fitting Grids
Finite Element Simulations of Logging-While-Drilling and Extra-Deep Azimuthal Resistivity Measurements using Non- Grids Théophile Chaumont-Frelet, D Pardo, Angel Rodriguez-Rozas To cite this version: Théophile
More informationCalculation of Sound Fields in Flowing Media Using CAPA and Diffpack
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack H. Landes 1, M. Kaltenbacher 2, W. Rathmann 3, F. Vogel 3 1 WisSoft, 2 Univ. Erlangen 3 inutech GmbH Outline Introduction Sound in Flowing
More informationOutput-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics: Overview and Recent Results
Output-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics: Overview and Recent Results Krzysztof J. Fidkowski University of Michigan David L. Darmofal Massachusetts Institute of
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 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 informationISM206 Lecture Optimization of Nonlinear Objective with Linear Constraints
ISM206 Lecture Optimization of Nonlinear Objective with Linear Constraints Instructor: Prof. Kevin Ross Scribe: Nitish John October 18, 2011 1 The Basic Goal The main idea is to transform a given constrained
More informationUncertainty Quantification in Discrete Fracture Network Models
Uncertainty Quantification in Discrete Fracture Network Models Claudio Canuto Dipartimento di Scienze Matematiche, Politecnico di Torino claudio.canuto@polito.it Joint work with Stefano Berrone, Sandra
More informationCS 542G: Robustifying Newton, Constraints, Nonlinear Least Squares
CS 542G: Robustifying Newton, Constraints, Nonlinear Least Squares Robert Bridson October 29, 2008 1 Hessian Problems in Newton Last time we fixed one of plain Newton s problems by introducing line search
More information2.098/6.255/ Optimization Methods Practice True/False Questions
2.098/6.255/15.093 Optimization Methods Practice True/False Questions December 11, 2009 Part I For each one of the statements below, state whether it is true or false. Include a 1-3 line supporting sentence
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 informationInterior-Point Methods for Linear Optimization
Interior-Point Methods for Linear Optimization Robert M. Freund and Jorge Vera March, 204 c 204 Robert M. Freund and Jorge Vera. All rights reserved. Linear Optimization with a Logarithmic Barrier Function
More informationApplied Mathematics 205. Unit I: Data Fitting. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit I: Data Fitting Lecturer: Dr. David Knezevic Unit I: Data Fitting Chapter I.4: Nonlinear Least Squares 2 / 25 Nonlinear Least Squares So far we have looked at finding a best
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 informationAppendix A Taylor Approximations and Definite Matrices
Appendix A Taylor Approximations and Definite Matrices Taylor approximations provide an easy way to approximate a function as a polynomial, using the derivatives of the function. We know, from elementary
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationReducing Uncertainty through Multi-Measurement Integration: from Regional to Reservoir scale
Reducing Uncertainty through Multi-Measurement Integration: from Regional to Reservoir scale Efthimios Tartaras Data Processing & Modeling Manager, Integrated Electromagnetics CoE, Schlumberger Geosolutions
More informationThe Finite Element Method
Chapter 1 The Finite Element Method 1.1 Introduction This document provides a brief introduction to the finite element method and illustrates how the method is implemented in oomph-lib. The first few sections
More informationENGINEERING MATHEMATICS I. CODE: 10 MAT 11 IA Marks: 25 Hrs/Week: 04 Exam Hrs: 03 PART-A
ENGINEERING MATHEMATICS I CODE: 10 MAT 11 IA Marks: 25 Hrs/Week: 04 Exam Hrs: 03 Total Hrs: 52 Exam Marks:100 PART-A Unit-I: DIFFERENTIAL CALCULUS - 1 Determination of n th derivative of standard functions-illustrative
More informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Advisor: Howard Elman Department of Computer Science Sept. 29, 2015 Tengfei
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
More informationA multigrid integral equation method for large-scale models with inhomogeneous backgrounds
IOP PUBLISHING JOURNAL OF GEOPHYSICS AND ENGINEERING J. Geophys. Eng. (28) 438 447 doi:1.188/1742-2132//4/7 A integral equation method for large-scale models with inhomogeneous backgrounds Masashi Endo,
More informationAnisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements
Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements K. Key (Scripps Institution of Oceanography), Z. Du* (PGS), J. Mattsson (PGS), A.
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationnonrobust estimation The n measurement vectors taken together give the vector X R N. The unknown parameter vector is P R M.
Introduction to nonlinear LS estimation R. I. Hartley and A. Zisserman: Multiple View Geometry in Computer Vision. Cambridge University Press, 2ed., 2004. After Chapter 5 and Appendix 6. We will use x
More informationIntroduction into Implementation of Optimization problems with PDEs: Sheet 3
Technische Universität München Center for Mathematical Sciences, M17 Lucas Bonifacius, Korbinian Singhammer www-m17.ma.tum.de/lehrstuhl/lehresose16labcourseoptpde Summer 216 Introduction into Implementation
More informationFinite Elements for Magnetohydrodynamics and its Optimal Control
Finite Elements for Magnetohydrodynamics and its Karl Kunisch Marco Discacciati (RICAM) FEM Symposium Chemnitz September 25 27, 2006 Overview 1 2 3 What is Magnetohydrodynamics? Magnetohydrodynamics (MHD)
More informationSeismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics
Seismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics Rune Mittet, EMGS SUMMARY are very useful for the interpretation of seismic data. Moreover, these
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2
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 informationComputational Modeling of Geoelectrical Soundings using PML-FDTD
120 Computational Modeling of Geoelectrical Soundings using PML-FDTD Luísa F. Ribeiro 1 and Marcela S. Novo 2 1,2 Federal University of Bahia UFBA, Electrical Engineering Department, Salvador BA, Brazil
More informationLINEAR AND NONLINEAR PROGRAMMING
LINEAR AND NONLINEAR PROGRAMMING Stephen G. Nash and Ariela Sofer George Mason University The McGraw-Hill Companies, Inc. New York St. Louis San Francisco Auckland Bogota Caracas Lisbon London Madrid Mexico
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Taylor s Theorem Can often approximate a function by a polynomial The error in the approximation
More informationManuscript received by the Editor 1 October 2009; revised manuscript received 2 January 2010; published online 16 June 2010.
GEOPHYSICS, VOL. 75, NO. MAY-JUNE ; P. F8F95, FIGS..9/. Simulation of triaxial induction measurements in dipping, invaded, and anisotropic formations using a Fourier series expansion in a nonorthogonal
More informationKeywords: Electric Machines, Rotating Machinery, Stator faults, Fault tolerant control, Field Weakening, Anisotropy, Dual rotor, 3D modeling
Analysis of Electromagnetic Behavior of Permanent Magnetized Electrical Machines in Fault Modes M. U. Hassan 1, R. Nilssen 1, A. Røkke 2 1. Department of Electrical Power Engineering, Norwegian University
More information1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that
Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a
More informationINTERIOR-POINT METHODS FOR NONCONVEX NONLINEAR PROGRAMMING: CONVERGENCE ANALYSIS AND COMPUTATIONAL PERFORMANCE
INTERIOR-POINT METHODS FOR NONCONVEX NONLINEAR PROGRAMMING: CONVERGENCE ANALYSIS AND COMPUTATIONAL PERFORMANCE HANDE Y. BENSON, ARUN SEN, AND DAVID F. SHANNO Abstract. In this paper, we present global
More informationA MATLAB Implementation of the. Minimum Relative Entropy Method. for Linear Inverse Problems
A MATLAB Implementation of the Minimum Relative Entropy Method for Linear Inverse Problems Roseanna M. Neupauer 1 and Brian Borchers 2 1 Department of Earth and Environmental Science 2 Department of Mathematics
More informationA multigrid method for large scale inverse problems
A multigrid method for large scale inverse problems Eldad Haber Dept. of Computer Science, Dept. of Earth and Ocean Science University of British Columbia haber@cs.ubc.ca July 4, 2003 E.Haber: Multigrid
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationSite Characterization & Hydrogeophysics
Site Characterization & Hydrogeophysics (Source: Matthew Becker, California State University) Site Characterization Definition: quantitative description of the hydraulic, geologic, and chemical properties
More informationOptimization and Root Finding. Kurt Hornik
Optimization and Root Finding Kurt Hornik Basics Root finding and unconstrained smooth optimization are closely related: Solving ƒ () = 0 can be accomplished via minimizing ƒ () 2 Slide 2 Basics Root finding
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 informationMorphing ensemble Kalman filter
Morphing ensemble Kalman filter and applications Center for Computational Mathematics Department of Mathematical and Statistical Sciences University of Colorado Denver Supported by NSF grants CNS-0623983
More informationPrimal-Dual Interior-Point Methods for Linear Programming based on Newton s Method
Primal-Dual Interior-Point Methods for Linear Programming based on Newton s Method Robert M. Freund March, 2004 2004 Massachusetts Institute of Technology. The Problem The logarithmic barrier approach
More informationAM 205: lecture 19. Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods
AM 205: lecture 19 Last time: Conditions for optimality Today: Newton s method for optimization, survey of optimization methods Optimality Conditions: Equality Constrained Case As another example of equality
More informationGradient Descent. Dr. Xiaowei Huang
Gradient Descent Dr. Xiaowei Huang https://cgi.csc.liv.ac.uk/~xiaowei/ Up to now, Three machine learning algorithms: decision tree learning k-nn linear regression only optimization objectives are discussed,
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationChapter 2. Optimization. Gradients, convexity, and ALS
Chapter 2 Optimization Gradients, convexity, and ALS Contents Background Gradient descent Stochastic gradient descent Newton s method Alternating least squares KKT conditions 2 Motivation We can solve
More informationFEM: Domain Decomposition and Homogenization for Maxwell s Equations of Large Scale Problems
FEM: and Homogenization for Maxwell s Equations of Large Scale Problems Karl Hollaus Vienna University of Technology, Austria Institute for Analysis and Scientific Computing February 13, 2012 Outline 1
More informationReduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems. Laura Iapichino Stefan Trenz Stefan Volkwein
Universität Konstanz Reduced-Order Multiobjective Optimal Control of Semilinear Parabolic Problems Laura Iapichino Stefan Trenz Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 347, Dezember 2015
More informationA Line search Multigrid Method for Large-Scale Nonlinear Optimization
A Line search Multigrid Method for Large-Scale Nonlinear Optimization Zaiwen Wen Donald Goldfarb Department of Industrial Engineering and Operations Research Columbia University 2008 Siam Conference on
More informationA self-adaptive goal-oriented hp finite element method with electromagnetic applications. Part II: Electrodynamics
Comput. Methods Appl. Mech. Engrg. xxx (007) xxx xxx www.elsevier.com/locate/cma A self-adaptive goal-oriented finite element method with electromagnetic applications. Part II: Electrodynamics D. Pardo
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More informationInexact Algorithms in Function Space and Adaptivity
Inexact Algorithms in Function Space and Adaptivity Anton Schiela adaptive discretization nonlinear solver Zuse Institute Berlin linear solver Matheon project A1 ZIB: P. Deuflhard, M. Weiser TUB: F. Tröltzsch,
More information