Numerical Methods: Structured vs. unstructured grids. General Introduction: Why numerical methods? Numerical methods and their fields of application
|
|
- Adelia Murphy
- 5 years ago
- Views:
Transcription
1 Numerical Methods: Structured vs. unstructured grids The goals o this course General : Why numerical methods? Numerical methods and their ields o alication Review o inite dierences Goals: Understanding the basics o all the inite s (dierences, elements, volumes The beauty: in the linear limit they are all really the same
2 Why numerical methods? Eamle: seismic wave roagation Seismometers homogeneous medium elosion In this case there are analytical solutions? Are they useul?
3 Why numerical methods? Eamle: seismic wave roagation Seismometers layered medium elosion... in this case quasi-analytical solutions eist, alicable or eamle or layered sediments...
4 Why numerical methods? Eamle: seismic wave roagation Seismometers long wavelength erturbations elosion in this case high-requency aroimations can be used (ray theory
5 Why numerical methods Eamle: seismic wave roagation Seismometers Generally heterogeneous medium elosion we need numerical solutions! we need grids! And big comuters
6 Partial Dierential Equations in Geohysics t P c s c ( s y z ressure acoustic wave seed sources The acoustic wave equation - seismology - acoustics - oceanograhy - meteorology t C k C v C RC C tracer concentration k diusivity v low velocity R reactivity sources Diusion, advection, Reaction - geodynamics - oceanograhy - meteorology - geochemistry - sedimentology - geohysical luid dynamics
7 Numerical methods: ields o alication Finite dierences - time-deendent PDEs - seismic wave roagation - geohysical luid dynamics - Mawell s equations - Ground enetrating radar -> robust, simle concet, easy to arallelize, regular grids, elicit method Finite elements - static and time-deendent PDEs - seismic wave roagation - geohysical luid dynamics - all roblems -> imlicit aroach, matri inversion, well ounded, irregular grids, more comle algorithms, engineering roblems Finite volumes - time-deendent PDEs - seismic wave roagation - mainly luid dynamics -> robust, simle concet, irregular grids, elicit method
8 Other Numerical methods: Particle-based methods - lattice gas methods - molecular dynamics - granular roblems - luid low - earthquake simulations -> very heterogeneous roblems, nonlinear roblems Boundary element methods - roblems with boundaries (ruture - based in analytical solutions - only discretization o lanes --> good or roblems with secial boundary conditions (ruture, cracks, etc Pseudosectral methods - orthogonal basis unctions - sectral accuracy o sace derivatives - wave roagation, GPR -> regular grids, elicit method, roblems with discontinuities
9 Numerical methods in all ields o Earth sciences Seismology Granular media - ruture Global dynamics Miing - Geochemistry Earthquake hysics Earth s magnetic ield Regional earthquakes
10 What is a inite dierence? Common deinitions o the derivative o (: lim 0 ( ( lim 0 ( ( lim 0 ( ( These are all correct deinitions in the limit ->0. But we want to remain FINITE
11 What is a inite dierence? The equivalent aroimations o the derivatives are: ( ( orward dierence ( ( backward dierence ( ( centered dierence
12 The big question: How good are the FD aroimations? This leads us to Taylor series...
13 Taylor Series... that leads to : 3 ( ( 1 ' '' ''' ( (! 3! (... ' ( O( The error o the irst derivative using the orward ormulation is o order. Is this the case or other ormulations o the derivative? Let s check!
14 Taylor Series... with the centered ormulation we get: 3 ( / ( / 1 ' ''' ( 3! (... ' ( O( The error o the irst derivative using the centered aroimation is o order. This is an imortant results: it DOES matter which ormulation we use. The centered scheme is more accurate!
15 Alternative Derivation o FD a a a ' d b b b ' d a b ( a b ( a b ' d Interolation Derivative a b 0 a b w1 0.5, w 0.5 Interolation weights w ' d 1, w 1 1 Derivative weights
16 Our irst FD algorithm! ( z y t s c P ressure c acoustic wave seed s sources Problem: Solve the 1D acoustic wave equation using the inite Dierence method. Solution: [ ] ( ( ( ( ( ( s t t c t
17 Problems: Stability ( t c [ ( ( ( ] ( t ( t s Stability: Careul analysis using harmonic unctions shows that a stable numerical calculation is subject to secial conditions (conditional stability. This holds or many numerical roblems. c ε 1
18 Problems: Disersion ( t c [ ( ( ( ] ( t ( t s True velocity Disersion: The numerical aroimation has artiicial disersion, in other words, the wave seed becomes requency deendent. You have to ind a requency bandwih where this eect is small. The solution is to use suicient grid oints er wavelength.
19 Our irst FD code! ( t c [ ( ( ( ] ( t ( t s % Time steing or i1:nt, end % FD dis(srint(' Time ste : %i',i; or j:n-1 d(j((j1-*(j(j-1/^; % sace derivative end new*-oldd*^; % time etraolation new(n/new(n/src(i*^; % add source term old; % time levels new; (10; % set boundaries ressure ree (n0; % Dislay lot(,,'b-' title(' FD ' drawnow
20 Our irst FD code! ( t c [ ( ( ( ] ( t ( t s % Time steing or i1:nt, % FD dis(srint(' Time ste : %i',i; Eercises (FD: 1. Increase the time ste manually and determine the stability limit numerically (c*/.. Make the medium heterogeneous. Put in a velocity contrast along the ais with a 30% erturbation. 3. Perturb the medium (e.g. 0% with random erturbations (dcrand([1 n]. What eect do you have on the roagating ulse? Do you think the result is accurate? 4. Comare the results with a simulation with time stes. or j:n-1 d(j((j1-*(j(j-1/^; % sace derivative end new*-oldd*^; % time etraolation new(n/new(n/src(i*^; % add source term old; % time levels new; (10; % set boundaries ressure ree (n0; end % Dislay lot(,,'b-' title(' FD ' drawnow
21 Finite Dierences - Summary Deending on the choice o the FD scheme (e.g. orward, backward, centered a numerical solution may be more or less accurate. Elicit inite dierence solutions to dierential equations are oten conditionally stable. The correct choice o the sace or time increment is crucial to enable accurate solutions. Sometimes it is useul to emloy so-called staggered grids where the ields are deined on searate grids which may imrove the overall accuracy o the scheme.
Numerical Methods in Geophysics. Introduction
: Why numerical methods? simple geometries analytical solutions complex geometries numerical solutions Applications in geophysics seismology geodynamics electromagnetism... in all domains History of computers
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationMagnetospheric Physics - Homework, 4/04/2014
Magnetosheric hysics - Homework, // 7. Fast wave, fast shock, erendicular shock, entroy, jum relation. Consider the fast erendicular shock. a) Determine the ositive root of the solution for the comression
More informationChapter 10. Classical Fourier Series
Math 344, Male ab Manual Chater : Classical Fourier Series Real and Comle Chater. Classical Fourier Series Fourier Series in PS K, Classical Fourier Series in PS K, are aroimations obtained using orthogonal
More informationNonlinear programming
08-04- htt://staff.chemeng.lth.se/~berntn/courses/otimps.htm Otimization of Process Systems Nonlinear rogramming PhD course 08 Bernt Nilsson, Det of Chemical Engineering, Lund University Content Unconstraint
More informationMHD Heat and Mass Transfer Forced Convection Flow Along a Stretching Sheet with Heat Generation, Radiation and Viscous Dissipation
Dhaka Univ. J. Sci. 59(): -8, (July) MHD Heat and Mass Transer Forced onvection Flo Along a Stretching Sheet ith Heat Generation, Radiation and Viscous Dissiation M. A. Samad, Sajid Ahmed and Md. Motaleb
More informationMeshless Methods for Scientific Computing Final Project
Meshless Methods for Scientific Comuting Final Project D0051008 洪啟耀 Introduction Floating island becomes an imortant study in recent years, because the lands we can use are limit, so eole start thinking
More informationChapter 2 Section 3. Partial Derivatives
Chapter Section 3 Partial Derivatives Deinition. Let be a unction o two variables and. The partial derivative o with respect to is the unction, denoted b D1 1 such that its value at an point (,) in the
More informationOn Closing the Constitutive Gap Between Forming and Crash Simulation
0 th International LS-DYNA Users Conerence etal Forming (3) On Closing the Constitutive Ga Between Forming and Crash Simulation F. Neukamm *,. Feucht *, A. Haue **, K. Roll * * Daimler AG, 7059 Sindelingen,
More informationPrimary Author: James D. McCalley, Iowa State University
Module T7: The ower Flow roblem 7 Module T7 The ower Flow roblem rimary Author: ames D. McCalley, Iowa State University Email Address: dm@iastate.edu Co-author: one rereuisite Cometencies:. Steady State
More information5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:
5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum
More informationResonances in high-contrast gratings with complex unit cell topology
Resonances in high-contrast gratings with comle unit cell toology Milan Maksimovic Focal-Vision & Otics, Oldenzaal, The Netherlands The XXI International Worksho on Otical Wave & Waveguide Theory and Numerical
More informationP043 Anisotropic 2.5D - 3C Finite-difference Modeling
P04 Anisotroic.5D - C Finite-difference Modeling A. Kostyukevych* (esseral echnologies Inc.), N. Marmalevskyi (Ukrainian State Geological Prosecting Institute), Y. Roganov (Ukrainian State Geological Prosecting
More informationOptical Fibres - Dispersion Part 1
ECE 455 Lecture 05 1 Otical Fibres - Disersion Part 1 Stavros Iezekiel Deartment of Electrical and Comuter Engineering University of Cyrus HMY 445 Lecture 05 Fall Semester 016 ECE 455 Lecture 05 Otical
More informationOn the relationship between sound intensity and wave impedance
Buenos Aires 5 to 9 Setember, 16 Acoustics for the 1 st Century PROCEEDINGS of the nd International Congress on Acoustics Sound Intensity and Inverse Methods in Acoustics: Paer ICA16-198 On the relationshi
More informationANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE BARCELONA PLASTIC DAMAGE MODEL
XII International Conerence on Comutational Plasticity. Fundamentals and Alications COMPLAS XII E. Oñate, D.R.J. Owen, D. Peric and B. Suárez (Eds) ANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE
More informationAccelerator School Transverse Beam Dynamics-2. V. S. Pandit
Accelerator School 8 Transverse Beam Dnamics- V. S. Pandit Equation of Motion Reference orbit is a single laner curve. Diole is used for bending and quadruole for focusing We use coordinates (r, θ, ) Diole
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More informationDifferentiation. The main problem of differential calculus deals with finding the slope of the tangent line at a point on a curve.
Dierentiation The main problem o dierential calculus deals with inding the slope o the tangent line at a point on a curve. deinition() : The slope o a curve at a point p is the slope, i it eists, o the
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationComputational Seismology: An Introduction
Computational Seismology: An Introduction Aim of lecture: Understand why we need numerical methods to understand our world Learn about various numerical methods (finite differences, pseudospectal methods,
More informationScaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling
Scaling Multile Point Statistics or Non-Stationary Geostatistical Modeling Julián M. Ortiz, Steven Lyster and Clayton V. Deutsch Centre or Comutational Geostatistics Deartment o Civil & Environmental Engineering
More informationTh P10 13 Alternative Misfit Functions for FWI Applied to Surface Waves
Th P0 3 Alternative Misfit Functions for FWI Alied to Surface Waves I. Masoni* Total E&P, Joseh Fourier University), R. Brossier Joseh Fourier University), J. Virieux Joseh Fourier University) & J.L. Boelle
More informationAn analytical approximation method for the stabilizing solution of the Hamilton-Jacobi equation based on stable manifold theory
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July -3, 27 ThA8.4 An analytical aroimation method for the stabilizing solution of the Hamilton-Jacobi
More informationComputer arithmetic. Intensive Computation. Annalisa Massini 2017/2018
Comuter arithmetic Intensive Comutation Annalisa Massini 7/8 Intensive Comutation - 7/8 References Comuter Architecture - A Quantitative Aroach Hennessy Patterson Aendix J Intensive Comutation - 7/8 3
More informationOptimal array pattern synthesis with desired magnitude response
Otimal array attern synthesis with desired magnitude resonse A.M. Pasqual a, J.R. Arruda a and P. erzog b a Universidade Estadual de Caminas, Rua Mendeleiev, 00, Cidade Universitária Zeferino Vaz, 13083-970
More informationNumerical Computation of the Eigenstructure of Cylindrical Acoustic Waveguides with Heated (or Cooled) Walls
Alied Mathematical Sciences, Vol. 3, 29, no. 7, 82-837 Numerical Comutation of the Eigenstructure of Clindrical Acoustic Waveguides with Heated (or Cooled) Walls Brian J. M c Cartin Alied Mathematics,
More informationApplied Fitting Theory VI. Formulas for Kinematic Fitting
Alied Fitting heory VI Paul Avery CBX 98 37 June 9, 1998 Ar. 17, 1999 (rev.) I Introduction Formulas for Kinematic Fitting I intend for this note and the one following it to serve as mathematical references,
More informationME 437/ME 537 PARTICLE TRANSPORT, DEPOSITION AND REMOVAL. Goodarz Ahmadi
ME 437/ME 537 PARTICLE TRANSPORT, DEPOSITION AND REMOVAL Goodarz Ahmadi Deartment o Mechanical and Aeronautical Engineering Clarkson University Potsdam, NY 13699-575 1 Air ollution and smog. Particle trajectories
More information3 More applications of derivatives
3 More alications of derivatives 3.1 Eact & ineact differentials in thermodynamics So far we have been discussing total or eact differentials ( ( u u du = d + dy, (1 but we could imagine a more general
More informationA BSS-BASED APPROACH FOR LOCALIZATION OF SIMULTANEOUS SPEAKERS IN REVERBERANT CONDITIONS
A BSS-BASED AROACH FOR LOCALIZATION OF SIMULTANEOUS SEAKERS IN REVERBERANT CONDITIONS Hamid Reza Abutalebi,2, Hedieh Heli, Danil Korchagin 2, and Hervé Bourlard 2 Seech rocessing Research Lab (SRL), Elec.
More informationDecoding First-Spike Latency: A Likelihood Approach. Rick L. Jenison
eurocomuting 38-40 (00) 39-48 Decoding First-Sike Latency: A Likelihood Aroach Rick L. Jenison Deartment of Psychology University of Wisconsin Madison WI 53706 enison@wavelet.sych.wisc.edu ABSTRACT First-sike
More information8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b).
8. THEOREM I the partial derivatives and eist near (a b) and are continuous at (a b) then is dierentiable at (a b). For a dierentiable unction o two variables z= ( ) we deine the dierentials d and d to
More informationA NEW GEOMETRIC ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD
Electronic Transactions on Numerical Analysis. Volume 45,. 330 341, 2016. Coyright c 2016,. ISSN 1068 9613. ETNA A NEW GEOMETRIC ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD WILLIAMS LÓPEZ
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More informationProblem set 6 for Quantum Field Theory course
Problem set 6 or Quantum Field Theory course 2018.03.13. Toics covered Scattering cross-section and decay rate Yukawa theory and Yukawa otential Scattering in external electromagnetic ield, Rutherord ormula
More informationThe achievable limits of operational modal analysis. * Siu-Kui Au 1)
The achievable limits o operational modal analysis * Siu-Kui Au 1) 1) Center or Engineering Dynamics and Institute or Risk and Uncertainty, University o Liverpool, Liverpool L69 3GH, United Kingdom 1)
More information. This is the Basic Chain Rule. x dt y dt z dt Chain Rule in this context.
Math 18.0A Gradients, Chain Rule, Implicit Dierentiation, igher Order Derivatives These notes ocus on our things: (a) the application o gradients to ind normal vectors to curves suraces; (b) the generaliation
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 80 29 JUNE 1998 NUMBER 26 Heavy Meson Decay Constants from Quenched Lattice QCD S. Aoki, 1, * M. Fukugita, 2 S. Hashimoto, 3, N. Ishizuka, 1,4 Y. Iwasaki, 1,4 K. Kanaya,
More informationOn the Fluid Dependence of Rock Compressibility: Biot-Gassmann Refined
Downloaded 0/9/3 to 99.86.4.8. Redistribution subject to SEG license or coyright; see Terms of Use at htt://library.seg.org/ On the luid Deendence of Rock Comressibility: Biot-Gassmann Refined Leon Thomsen,
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More informationBasic mathematics of economic models. 3. Maximization
John Riley 1 January 16 Basic mathematics o economic models 3 Maimization 31 Single variable maimization 1 3 Multi variable maimization 6 33 Concave unctions 9 34 Maimization with non-negativity constraints
More informationPMD monitoring in traffic-carrying optical systems and its statistical analysis
PMD monitoring in traic-carrying otical systems and its statistical analysis Juneng Jiang 1,, Sathyanarayanan Sundhararajan 1, Doug Richards, Steve Oliva and Rongqing Hui 1,* 1 Electrical Engineering &
More informationApplications of local fractional calculus to engineering in fractal time-space:
Applications o local ractional calculus to engineering in ractal time-space: Local ractional dierential equations with local ractional derivative Yang XiaoJun Department o Mathematics and Mechanics, China
More information( 1) ( 2) ( 1) nan integer, since the potential is no longer simple harmonic.
. Anharmonic Oscillators Michael Fowler Landau (para 8) considers a simple harmonic oscillator with added small potential energy terms mα + mβ. We ll simpliy slightly by dropping the term, to give an equation
More information9.3 Graphing Functions by Plotting Points, The Domain and Range of Functions
9. Graphing Functions by Plotting Points, The Domain and Range o Functions Now that we have a basic idea o what unctions are and how to deal with them, we would like to start talking about the graph o
More informationAcoustic forcing of flexural waves and acoustic fields for a thin plate in a fluid
Acoustic orcing o leural waves and acoustic ields or a thin plate in a luid Darryl MCMAHON Maritime Division, Deence Science and Technology Organisation, HMAS Stirling, WA Australia ABSTACT Consistency
More informationA Recursive Block Incomplete Factorization. Preconditioner for Adaptive Filtering Problem
Alied Mathematical Sciences, Vol. 7, 03, no. 63, 3-3 HIKARI Ltd, www.m-hiari.com A Recursive Bloc Incomlete Factorization Preconditioner for Adative Filtering Problem Shazia Javed School of Mathematical
More informationChapter 6. Thermodynamics and the Equations of Motion
Chater 6 hermodynamics and the Equations of Motion 6.1 he first law of thermodynamics for a fluid and the equation of state. We noted in chater 4 that the full formulation of the equations of motion required
More informationDirac s Hole Theory and the Pauli Principle: Clearing up the Confusion
Adv. Studies Theor. Phys., Vol. 3, 29, no. 9, 323-332 Dirac s Hole Theory and the Pauli Princile: Clearing u the Conusion Dan Solomon Rauland-Borg Cororation 82 W. Central Road Mount Prosect, IL 656, USA
More informationTime Frequency Aggregation Performance Optimization of Power Quality Disturbances Based on Generalized S Transform
Time Frequency Aggregation Perormance Otimization o Power Quality Disturbances Based on Generalized S Transorm Mengda Li Shanghai Dianji University, Shanghai 01306, China limd @ sdju.edu.cn Abstract In
More informationSemicontinuous filter limits of nets of lattice groupvalued
Semicontinuous ilter limits o nets o lattice grouvalued unctions THEMATIC UNIT: MATHEMATICS AND APPLICATIONS A Boccuto, Diartimento di Matematica e Inormatica, via Vanvitelli, I- 623 Perugia, Italy, E-mail:
More information+++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic +++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic
Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic Daniel Leitner, Johannes Krof, Günther Zauner, TU Vienna, Austria, dleitner@osiris.tuwien.ac.at Yuri Karov, Yuri Senichenkov, Yuri
More informationShort Solutions to Practice Material for Test #2 MATH 2421
Short Solutions to Practice Material for Test # MATH 4 Kawai (#) Describe recisely the D surfaces listed here (a) + y + z z = Shere ( ) + (y ) + (z ) = 4 = The center is located at C (; ; ) and the radius
More informationExtreme Values of Functions
Extreme Values o Functions When we are using mathematics to model the physical world in which we live, we oten express observed physical quantities in terms o variables. Then, unctions are used to describe
More informationFFTs in Graphics and Vision. Rotational and Reflective Symmetry Detection
FFTs in Grahics and Vision Rotational and Relectie Symmetry Detection Outline Reresentation Theory Symmetry Detection Rotational Symmetry Relectie Symmetry Reresentation Theory Recall: A rou is a set o
More informationL p Norms and the Sinc Function
L Norms and the Sinc Function D. Borwein, J. M. Borwein and I. E. Leonard March 6, 9 Introduction The sinc function is a real valued function defined on the real line R by the following eression: sin if
More informationChapter 6: Sound Wave Equation
Lecture notes on OPAC0- ntroduction to Acoustics Dr. Eser OLĞAR, 08 Chater 6: Sound Wave Equation. Sound Waves in a medium the wave equation Just like the eriodic motion of the simle harmonic oscillator,
More informationCaveats on the Implementation of the Generalized Material Point Method
Coyright c 2008 Tech Science Press CMES, vol.31, no.2,.85-106, 2008 Caveats on the Imlementation of the Generalized Material Point Method O. Buzzi 1, D. M. Pedroso 2 and A. Giacomini 1 Abstract: The material
More information( x) f = where P and Q are polynomials.
9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational
More informationToday. CS-184: Computer Graphics. Introduction. Some Examples. 2D Transformations
Toda CS-184: Comuter Grahics Lecture #3: 2D Transformations Prof. James O Brien Universit of California, Berkele V2005F-03-1.0 2D Transformations Primitive Oerations Scale, Rotate, Shear, Fli, Translate
More informationarxiv: v1 [physics.class-ph] 14 Mar 2008
Full-wave finite-difference time-domain simulation of electromagnetic cloaking structures arxiv:8.6v hysics.class-h Mar 8 Yan Zhao*, Christos Argyrooulos, and Yang Hao Queen Mary University of London,
More information1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum
Estimation of fluid mobility from frequency deendent azimuthal AVO a synthetic model study Yingrui Ren 1*, Xiaoyang Wu 2, Mark Chaman 1 and Xiangyang Li 2,3 1 University of Edinburgh, 2 British Geological
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. III Stability Theory - Peter C. Müller
STABILITY THEORY Peter C. Müller University of Wuertal, Germany Keywords: Asymtotic stability, Eonential stability, Linearization, Linear systems, Lyaunov equation, Lyaunov function, Lyaunov stability,
More informationLECTURE 3 BASIC QUANTUM THEORY
LECTURE 3 BASIC QUANTUM THEORY Matter waves and the wave function In 194 De Broglie roosed that all matter has a wavelength and exhibits wave like behavior. He roosed that the wavelength of a article of
More informationPlanar Transformations and Displacements
Chater Planar Transformations and Dislacements Kinematics is concerned with the roerties of the motion of oints. These oints are on objects in the environment or on a robot maniulator. Two features that
More informationSolved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.
Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the
More informationFUZZY CONTROL. Khurshid Ahmad, Professor of Computer Science, Department of Computer Science Trinity College, Dublin-2, IRELAND Nov 15-16th, 2011
: Mamdani & akagi-sugeno Controllers Khurshid Ahmad, Professor of Comuter Science, Deartment of Comuter Science rinit College, Dublin-, IRELAND Nov 5-6th, 0 htts://www.cs.tcd.ie/khurshid.ahmad/eaching/eaching.html
More informationUseful concepts associated with the Bernoulli equation. Dynamic
Useful concets associated with the Bernoulli equation - Static, Stagnation, and Dynamic Pressures Bernoulli eq. along a streamline + ρ v + γ z = constant (Unit of Pressure Static (Thermodynamic Dynamic
More informationSection 3.4: Concavity and the second Derivative Test. Find any points of inflection of the graph of a function.
Unit 3: Applications o Dierentiation Section 3.4: Concavity and the second Derivative Test Determine intervals on which a unction is concave upward or concave downward. Find any points o inlection o the
More informationCISE-301: Numerical Methods Topic 1:
CISE-3: Numerical Methods Topic : Introduction to Numerical Methods and Taylor Series Lectures -4: KFUPM Term 9 Section 8 CISE3_Topic KFUPM - T9 - Section 8 Lecture Introduction to Numerical Methods What
More informationStrong Lyapunov Functions for Systems Satisfying the Conditions of La Salle
06 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 6, JUNE 004 Strong Lyapunov Functions or Systems Satisying the Conditions o La Salle Frédéric Mazenc and Dragan Ne sić Abstract We present a construction
More information9 The Theory of Special Relativity
9 The Theory of Secial Relativity Assign: Read Chater 4 of Carrol and Ostlie (2006) Newtonian hysics is a quantitative descrition of Nature excet under three circumstances: 1. In the realm of the very
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationAvailable online at ScienceDirect. Procedia Engineering 102 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 102 (2015 ) 1366 1372 The 7th World Congress on Particle Technology (WCPT7) Numerical Modelling o ESP or Design Otimization
More informationNumerical Simulations of Methane Aromatization with and without a Ceramic Hydrogen Separation Membrane
Numerical Simulations o Methane Aromatization with and without a Ceramic ydrogen Searation Membrane Zuoan Li*, Christian ølseth, Selene ernandez Moreudo, and eidar augsrud University o Oslo, Deartment
More informationReview of Prerequisite Skills for Unit # 2 (Derivatives) U2L2: Sec.2.1 The Derivative Function
UL1: Review o Prerequisite Skills or Unit # (Derivatives) Working with the properties o exponents Simpliying radical expressions Finding the slopes o parallel and perpendicular lines Simpliying rational
More informationCan ISO MRS and ASTM HDB Rated Materials Be Harmonized
an ISO MRS and ASTM HDB Rated Materials Be Harmonized Dr. Z. Jimmy Zhou The Dow hemical omany, USA Dr. Eugene F. Palermo Plastics Pie Institute, Washington D, USA ABSTRAT ISO MRS and ASTM HDB rating systems
More informationUnsteady Flow of a Dusty Conducting Fluid through porous medium between Parallel Porous Plates with Temperature Dependent Viscosity and Heat Source
Volume Issue3 3- June www.ijsret.org ISSN 78-88 Unsteady Flow of a Dusty Conducting Fluid through orous medium between Parallel Porous Plates with Temerature Deendent Viscosity and Heat Source Shalini
More informationElectronic properties of Graphene and 2-D materials. Part 2
Electronic roerties of Grahene and -D materials z Part y The laws of Quantum Mechanics continued Time indeendent Schroedinger equation Eamles Particle in a well Harmonic oscillator Quantum Tunneling The
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
More informationBasic statistical models
Basic statistical models Valery Pokrovsky March 27, 2012 Part I Ising model 1 Definition and the basic roerties The Ising model (IM) was invented by Lenz. His student Ising has found the artition function
More informationPassive Identification is Non Stationary Objects With Closed Loop Control
IOP Conerence Series: Materials Science and Engineering PAPER OPEN ACCESS Passive Identiication is Non Stationary Obects With Closed Loo Control To cite this article: Valeriy F Dyadik et al 2016 IOP Con.
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationAleatory or epistemic? Does it matter?
Worksho on isk Accetance and isk Communication tanord University, March 6-7, 007 Aleatory or eistemic? Does it matter? Armen Der Kiureghian University o Caliornia, Berkeley Ove Ditlevsen Technical University
More informationf self = 1/T self (b) With revolution, rotaton period T rot in second and the frequency Ω rot are T yr T yr + T day T rot = T self > f self
Problem : Units : Q-a Mathematically exress the relationshi between the different units of the hysical variables: i) Temerature: ) Fahrenheit and Celsius; 2) Fahrenheit and Kelvin ii) Length: ) foot and
More informationCFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids
CFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids T. Toulorge a,, W. Desmet a a K.U. Leuven, Det. of Mechanical Engineering, Celestijnenlaan 3, B-31 Heverlee, Belgium Abstract
More informationAnalysis Scheme in the Ensemble Kalman Filter
JUNE 1998 BURGERS ET AL. 1719 Analysis Scheme in the Ensemble Kalman Filter GERRIT BURGERS Royal Netherlands Meteorological Institute, De Bilt, the Netherlands PETER JAN VAN LEEUWEN Institute or Marine
More informationA Critical Investigation of High-Order Flux Limiters In Multiphase Flow Problems
A Critical Investigation o High-Order Flux Limiters In Multiphase Flow Problems Chris Guenther Fluent In., 3647 Collins Ferry Rd., Morgantown, WV 26505, USA cpg@luent.com ABSTRACT. In recent years inite
More informationLecture 2 The First Law of Thermodynamics (Ch.1)
Lecture he First Law o hermodynamics (h.) Lecture - we introduced macroscopic parameters that describe the state o a thermodynamic system (including temperature), the equation o state (,,) 0, and linked
More informationDiffusion in generalized lattice-gas models
PHYSICAL REVIEW B VOLUME 58, NUMBER Diusion in generalized lattice-gas models O. M. Braun Institute o Physics, Ukrainian Academy o Sciences, UA-50 Kiev, Ukraine DECEMBER 998-II C. A. Sholl Department o
More informationInternational Journal of Mathematics Trends and Technology- Volume3 Issue4-2012
Effect of Hall current on Unsteady Flow of a Dusty Conducting Fluid through orous medium between Parallel Porous Plates with Temerature Deendent Viscosity and Thermal Radiation Harshbardhan Singh and Dr.
More informationDouble porosity modeling in elastic wave propagation for reservoir characterization
Stanord Exploration Project, Report 98, September 2, 1998, pages 1 243 Double porosity modeling in elastic wave propagation or reservoir characterization James G. Berryman 1 and Herbert F. Wang 2 keywords:
More information16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes
Infinite Series 6.2 Introduction We extend the concet of a finite series, met in Section 6., to the situation in which the number of terms increase without bound. We define what is meant by an infinite
More informationHigher order theory for analytic saddle point approximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface
Higher order theory for analytic saddle oint aroximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface P.F Daley ABSTACT The high frequency solution to the roblem of a Ρ P and Ρ S reflected
More informationNON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION
NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100
More informationUNDERSTANDING of unsteady combustion of solid
JOURNAL OF PROPULSION AND POWER Vol. 4, No. 5, September October 8 Analytical Solution or Pressure-Coupled Combustion Response Functions o Composite Solid Propellants Michael Shusser and Fred E. C. Culick
More informationComptes rendus de l Academie bulgare des Sciences, Tome 59, 4, 2006, p POSITIVE DEFINITE RANDOM MATRICES. Evelina Veleva
Comtes rendus de l Academie bulgare des ciences Tome 59 4 6 353 36 POITIVE DEFINITE RANDOM MATRICE Evelina Veleva Abstract: The aer begins with necessary and suicient conditions or ositive deiniteness
More information