Uniformly best wavenumber approximations by spatial central difference operators: An initial investigation
|
|
- Louise Hill
- 5 years ago
- Views:
Transcription
1 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 by central difference schemes is resented. A central difference stencil is derived based on the theorem and is comared with disersion relation reserving schemes and with classical central differences for a relevant test roblem. 1 Introduction Modelling wave roagation over sizeable intervals using finite differences is a common roblem in fields ranging from aeroacoustics to seismology. For high frequency roblems the numerical error may over time be dominated by inaccurate aroximations of the disersion relation, leading to errors in hase and grou velocity, unless the satial increment, x is very small. A remedy, resented in [1] for central differences, is to erturb the classical schemes by an extra arameter but decreasing the formal accuracy. The new arameter is used to minimise the disersion error in the L 2 [,π/2] norm. Such schemes are nown as Disersion Relation Preserving (DRP). For other aroaches based on similar ideas, see e.g. [2, 3, 4]. The DRP aroach is disadvantageous in that it rovides no means of obtaining wavenumber-secific error bounds. For roblems involving a range of wavenumbers it is more convenient to minimise the disersion error in the L -norm. In this aer we resent a characterisation theorem for best uniform wavenumber aroxima- Vitor Linders Linöing University, Deartment of Mathematics, Comutational Mathematics, SE Linöing, Sweden, vitor.linders@liu.se Jan Nordström Linöing University, Deartment of Mathematics, Comutational Mathematics, SE Linöing, Sweden, jan.nordstrom@liu.se 1
2 2 Vitor Linders and Jan Nordström tions. New central difference schemes are derived and comared with their classical and DRP counterarts. 2 Central difference schemes We begin by demonstrating why classical central differences are subotimal for wavenumber aroximation. Consider a central difference stencil of order 2, (u x ) j = 1 x =1 (u j+ u j ) + O( x 2 ). The numerical wavenumber of this scheme is (see e.g. [1]) ξ c = 2 =1 sin(ξ ) (1) where ξ = xκ and κ is the exact wavenumber of the roagating solution. Here we let ξ [,ξ max ] [,π]. A Taylor exansion reveals that to obtain desired accuracy, must satisfy The following observation is useful. For the roof, see [5]. Lemma 1. Consider the function f (x) = 1 2 =1 2 =.. (2) T (x) where satisfies (2) and T (x) is the th order Chebyshev olynomial of the first ind, uniquely defined through the relation T (cos(φ)) = cos(φ). Then f (x) = d (1 x) for some d that deends exclusively on. Theorem 1. Classical central difference stencils underestimate the seed of roagating solutions. Proof. Let E c (ξ ) ξ ξ c = ξ 2 =1 sin(ξ )
3 Title Suressed Due to Excessive Length 3 be the error function associated with a classical central difference stencil of order O( x 2 ). Note from the definition of f and the Chebyshev olynomials that de c dξ = f (cos(ξ )). It follows that de c dξ = only when cos(ξ ) = 1. It is thus clear that E c(ξ ) has an inflection oint at ξ = and no other extrema in the domain of interest. Consequently E c (ξ ) is monotonic and since E c () = and E c (π) = π it is also increasing. We thus have ξ c ξ with equality only at ξ =. Therefore the classical central difference stencils underestimate the analytic wavenumber. It follows now that the relative error in the hase seed is v v v and so the hase seed is underestimated. = ξ c ξ 1 Let us now, lie for DRP schemes, erturb the stencil by adding an additional coefficient, a without increasing the accuracy. The coefficients of the new stencil must solve the system (2) for each a, though this system is now underdetermined. Calling the coefficients of the new system a, = 1,..., + 1, we must have a linear deendence of the first coefficients on the added arameter, a. We write a = + c a, = 1,...,. In view of (3) the numerical wavenumber of the erturbed stencil is ξ = 2 =1 Let us define the error function of this scheme as a sin(ξ ), ξ ξ max π. (3) E(ξ ) ξ ξ = ξ 2 a sin(ξ ). (4) Our goal is to choose a such as to minimise the magnitude of any extrema of E(ξ ). In order to do so we will extend Lemma 1 to the erturbed stencil. For a detailed roof, see [5]. Lemma 2. Let =1 g (x) = 1 2 a T (x) where a = + c a are the coefficients of the erturbed central difference stencil as defined reviously, and T (x) is the th order Chebushev olynomial. Then [ ( )] g (x) = (1 x) (1 x) d c a + d 1 a c =1
4 4 Vitor Linders and Jan Nordström where d and are defined as before. Corollary 1. E(ξ ) has at most one extremum in the oen interval (,ξ max ] and it is located at ξ = ξ r = arccos ( 1 d d ]) [1 c. a For good aroximations ξ r is a minimum. This occurs only when a and c have the same sign and a c, where equality holds only for classical stencils. Again, for the roof we refer to [5]. From the above corollary we conclude that we can have E(ξ ) = E only at two ossible oints, namely at ξ r or ξ max, i.e. E = max{ E(ξ max ), E(ξ r ) }. (5) Of course E(ξ r ) and E(ξ max ) deend on how we choose a and in view of (5) it is of interest to investigate this deendency. Our goal is to find the choice of a that minimises (5). In fact we have Theorem 2. Consider a oint central difference scheme of order O( x 2 ) with numerical wavenumber ξ defined as in (3), and a corresonding error function E(ξ ) = ξ ξ. The stencil that uniformly minimises the error function, i.e. min ξ E = min a R ξ ξ, is uniquely characterised by the roerty E(ξ r ) + E(ξ max ) =. (6) The roof is found in [5] where it is also demonstrated how this result generalises to an arbitrary number of free arameters, a,...,a +n. 3 A numerical examle For a given ξ max solving (6) is a simle matter. Even for general ξ max good estimates may be found by relacing E(ξ r ) by a suitable olynomial aroximation. For the case = 1 this is shown in Fig. 1 for ξ max [,π/2]. Here e = min ξ E. To illustrate the strength of Theorem 2 we consider a rofoundly olychromatic solution to the advection equation over a eriodic domain: u t + u x =, x 3, t u(x,) = ex( 32(x 1/2) 2 ). This ulse is narrow and its Fourier transform is wide resulting in a significant contribution from a broad range of wavenumbers. The dominating wavenumber is
5 Title Suressed Due to Excessive Length Ξ max Fig. 1: Error of best uniform wavenumber aroximation for given ξ max [,π/2] (orange dots) and general solution using third degree olynomial aroximation (blac line). κ =. Contributions from larger wavenumbers decay exonentially but slowly. It maes sense to use a scheme that accurately aroximates the disersion relation near ξ = κ x = and for some suitably chosen region of larger wavenumbers. From Fig. 1 we see that if we choose ξ max = π/4.785 we will have an error in the disersion relation of around 1 3. Solving (6) gives the scheme a 1 = , a 2 = In Fig. 2 the disersion error of the new scheme is lotted as a function of ξ. The aroximation starts to deviate from the exact result for ξ ξ max. The errors of a fourth order classical central difference scheme and of a five-oint DRP scheme [1] are also included. As exected from Theorem 1 the classical stencil underestimates the disersion relation, seen by the ositive sign of the error. For the shown range of ξ it seems that the DRP scheme overestimates the disersion relation whereas our new scheme stays within tight error bounds. We set x = 1/12 and integrate in time using the classical fourth order Runge- Kutta scheme with time ste t = 1 3 so the contribution from the temoral discretisation is small. The exact and numerical solutions are shown in Fig. 3 (to) together with the error as a function of time (bottom). All numerical solutions quicly diserse into a train of ulses of decaying amlitude trailing behind the main ea. As exected from Fig. 2, the DRP scheme overestimates the seed of some ulses. Our new scheme does this as well but to a much reduced extent. The smaller ulse train behind the DRP solution with resect to the new scheme may be attributed to a better aroximation for very high wavenumbers. However, since the contribution of these wavenumbers are comarably small, the resulting error remains larger for the DRP scheme as comared with the classical and the new stencil.
6 6 Vitor Linders and Jan Nordström.1.5 New Classical DRP E = e 4 E(ξ) ξ Fig. 2: Disersion error for the new scheme, the classical 4 th order stencil, and a five-oint DRP scheme. 4 Extension to multile dimensions Extending the new stencil to multile dimensions is in rincile straight forward. As an examle, for the advection roblem in 2D we may, after discretising in sace, write v t + (D x I y )v + (I x D y )v = where v is a grid vector aroximating the true solution, D x,y are eriodic oerators containing the central difference stencil and oerating on a given cartesian grid in the x and y direction resectively. Here I x,y are identity matrices of aroriate dimensions and denotes the Kronecer roduct. For this situation, the solution roagates at an angle θ with resect to the x- axis. It should be noted that the stencils resented here are otimal for the onedimensional roblem, that is for the cases when θ = nπ/4, n =,1,2,3. For any other angle the stencils will be subotimal since the numerical disersion relation deends on the direction of roagation. In other words, these stencils may be sensitive to numerical anisotroy. For a comrehensive overview of methods that handles this issue, see e.g. [6]. At resent we shall not consider this henomenon further.
7 Title Suressed Due to Excessive Length New DRP Classical Exact 3 u(x, t) x.35.3 New DRP Classical Error develoment in time.25 Error Time Fig. 3: (To) Exact solution and numerical aroximations after 2 time stes. (Bottom) Corresonding errors. 5 Conclusion We have roved a characterisation theorem for best uniform wavenumber aroximations by central difference stencils with one free arameter. The best aroximation is unique and may be easily obtained numerically for a given range of wavenumbers. This allows for accurate aroximations of roblems of high frequency waves, or multi-frequency solutions, with a relatively coarse satial mesh.
8 8 Vitor Linders and Jan Nordström References 1. C.K.W. Tam, J.C. Web, Disersion-Relation-Preserving Finite Difference Schemes for Comutational Acoustics, J. Comut. Phys. 17, (1993) 2. D.W. Zing, H. Lomax, H. Jurgens, High-accuracy finite-difference schemes for linear wave roagation, SIAM J. Sci. Comute. 17, (1996) 3. D.W. Zingg, H. Lomax, H. Jurgens, An otimized finite-difference scheme for wave roagation roblems, AIAA aer 93(459) (1993) 4. C. Bogey, C. Bailly, A family of low disersive and low dissiative exlicit schemes for flow and noise comutations, J. Comute. Phys. 194, (24) 5. V. Linders, J. Nordström, Uniformly best wavenumber aroximations by satial central difference oerators, LiTH-MAT-R, 214:17, 214, Deartment of Mathematics, Linöing University 6. A. Sescu, Numerical anisotroy in finite differencing, Advances in Difference Equations 215:9 (215), Mississii State University
Radial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More informationAll-fiber Optical Parametric Oscillator
All-fiber Otical Parametric Oscillator Chengao Wang Otical Science and Engineering, Deartment of Physics & Astronomy, University of New Mexico Albuquerque, NM 87131-0001, USA Abstract All-fiber otical
More informationApplications to stochastic PDE
15 Alications to stochastic PE In this final lecture we resent some alications of the theory develoed in this course to stochastic artial differential equations. We concentrate on two secific examles:
More informationVarious Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems
Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various
More informationEquivalence of Wilson actions
Prog. Theor. Ex. Phys. 05, 03B0 7 ages DOI: 0.093/te/tv30 Equivalence of Wilson actions Physics Deartment, Kobe University, Kobe 657-850, Jaan E-mail: hsonoda@kobe-u.ac.j Received June 6, 05; Revised August
More informationAn Estimate For Heilbronn s Exponential Sum
An Estimate For Heilbronn s Exonential Sum D.R. Heath-Brown Magdalen College, Oxford For Heini Halberstam, on his retirement Let be a rime, and set e(x) = ex(2πix). Heilbronn s exonential sum is defined
More informationParticipation Factors. However, it does not give the influence of each state on the mode.
Particiation Factors he mode shae, as indicated by the right eigenvector, gives the relative hase of each state in a articular mode. However, it does not give the influence of each state on the mode. We
More informationSums of independent random variables
3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where
More informationLORENZO BRANDOLESE AND MARIA E. SCHONBEK
LARGE TIME DECAY AND GROWTH FOR SOLUTIONS OF A VISCOUS BOUSSINESQ SYSTEM LORENZO BRANDOLESE AND MARIA E. SCHONBEK Abstract. In this aer we analyze the decay and the growth for large time of weak and strong
More informationThe Algebraic Structure of the p-order Cone
The Algebraic Structure of the -Order Cone Baha Alzalg Abstract We study and analyze the algebraic structure of the -order cones We show that, unlike oularly erceived, the -order cone is self-dual for
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
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 informationResearch Article An iterative Algorithm for Hemicontractive Mappings in Banach Spaces
Abstract and Alied Analysis Volume 2012, Article ID 264103, 11 ages doi:10.1155/2012/264103 Research Article An iterative Algorithm for Hemicontractive Maings in Banach Saces Youli Yu, 1 Zhitao Wu, 2 and
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationPulse Propagation in Optical Fibers using the Moment Method
Pulse Proagation in Otical Fibers using the Moment Method Bruno Miguel Viçoso Gonçalves das Mercês, Instituto Suerior Técnico Abstract The scoe of this aer is to use the semianalytic technique of the Moment
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More informationNumerical Methods for Particle Tracing in Vector Fields
On-Line Visualization Notes Numerical Methods for Particle Tracing in Vector Fields Kenneth I. Joy Visualization and Grahics Research Laboratory Deartment of Comuter Science University of California, Davis
More informationApproximation of the Euclidean Distance by Chamfer Distances
Acta Cybernetica 0 (0 399 47. Aroximation of the Euclidean Distance by Chamfer Distances András Hajdu, Lajos Hajdu, and Robert Tijdeman Abstract Chamfer distances lay an imortant role in the theory of
More informationCoding Along Hermite Polynomials for Gaussian Noise Channels
Coding Along Hermite olynomials for Gaussian Noise Channels Emmanuel A. Abbe IG, EFL Lausanne, 1015 CH Email: emmanuel.abbe@efl.ch Lizhong Zheng LIDS, MIT Cambridge, MA 0139 Email: lizhong@mit.edu Abstract
More informationImproved Bounds on Bell Numbers and on Moments of Sums of Random Variables
Imroved Bounds on Bell Numbers and on Moments of Sums of Random Variables Daniel Berend Tamir Tassa Abstract We rovide bounds for moments of sums of sequences of indeendent random variables. Concentrating
More informationLINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL
LINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL Mohammad Bozorg Deatment of Mechanical Engineering University of Yazd P. O. Box 89195-741 Yazd Iran Fax: +98-351-750110
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
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 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 informationGeneration of Linear Models using Simulation Results
4. IMACS-Symosium MATHMOD, Wien, 5..003,. 436-443 Generation of Linear Models using Simulation Results Georg Otte, Sven Reitz, Joachim Haase Fraunhofer Institute for Integrated Circuits, Branch Lab Design
More informationLecture 6. 2 Recurrence/transience, harmonic functions and martingales
Lecture 6 Classification of states We have shown that all states of an irreducible countable state Markov chain must of the same tye. This gives rise to the following classification. Definition. [Classification
More informationAdditive results for the generalized Drazin inverse in a Banach algebra
Additive results for the generalized Drazin inverse in a Banach algebra Dragana S. Cvetković-Ilić Dragan S. Djordjević and Yimin Wei* Abstract In this aer we investigate additive roerties of the generalized
More informationStochastic integration II: the Itô integral
13 Stochastic integration II: the Itô integral We have seen in Lecture 6 how to integrate functions Φ : (, ) L (H, E) with resect to an H-cylindrical Brownian motion W H. In this lecture we address the
More informationε(ω,k) =1 ω = ω'+kv (5) ω'= e2 n 2 < 0, where f is the particle distribution function and v p f v p = 0 then f v = 0. For a real f (v) v ω (kv T
High High Power Power Laser Laser Programme Programme Theory Theory and Comutation and Asects of electron acoustic wave hysics in laser backscatter N J Sircombe, T D Arber Deartment of Physics, University
More informationSTABILITY OF PERIODIC WAVES FOR THE GENERALIZED BBM EQUATION
Dedicated to Professor Philie G. Ciarlet on his 70th birthday STABILITY OF PERIODIC WAVES FOR THE GENERALIZED BBM EQUATION MARIANA H R GU We consider the generalized Benjamin-Bona-Mahony (gbbm) equation
More informationPositivity, local smoothing and Harnack inequalities for very fast diffusion equations
Positivity, local smoothing and Harnack inequalities for very fast diffusion equations Dedicated to Luis Caffarelli for his ucoming 60 th birthday Matteo Bonforte a, b and Juan Luis Vázquez a, c Abstract
More informationAsymptotically Optimal Simulation Allocation under Dependent Sampling
Asymtotically Otimal Simulation Allocation under Deendent Samling Xiaoing Xiong The Robert H. Smith School of Business, University of Maryland, College Park, MD 20742-1815, USA, xiaoingx@yahoo.com Sandee
More information3.4 Design Methods for Fractional Delay Allpass Filters
Chater 3. Fractional Delay Filters 15 3.4 Design Methods for Fractional Delay Allass Filters Above we have studied the design of FIR filters for fractional delay aroximation. ow we show how recursive or
More informationDomain Dynamics in a Ferroelastic Epilayer on a Paraelastic Substrate
Y. F. Gao Z. Suo Mechanical and Aerosace Engineering Deartment and Princeton Materials Institute, Princeton University, Princeton, NJ 08544 Domain Dynamics in a Ferroelastic Eilayer on a Paraelastic Substrate
More informationAN OPTIMAL CONTROL CHART FOR NON-NORMAL PROCESSES
AN OPTIMAL CONTROL CHART FOR NON-NORMAL PROCESSES Emmanuel Duclos, Maurice Pillet To cite this version: Emmanuel Duclos, Maurice Pillet. AN OPTIMAL CONTROL CHART FOR NON-NORMAL PRO- CESSES. st IFAC Worsho
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationON POLYNOMIAL SELECTION FOR THE GENERAL NUMBER FIELD SIEVE
MATHEMATICS OF COMPUTATIO Volume 75, umber 256, October 26, Pages 237 247 S 25-5718(6)187-9 Article electronically ublished on June 28, 26 O POLYOMIAL SELECTIO FOR THE GEERAL UMBER FIELD SIEVE THORSTE
More informationThe Poisson Regression Model
The Poisson Regression Model The Poisson regression model aims at modeling a counting variable Y, counting the number of times that a certain event occurs during a given time eriod. We observe a samle
More informationFor q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i
Comuting with Haar Functions Sami Khuri Deartment of Mathematics and Comuter Science San Jose State University One Washington Square San Jose, CA 9519-0103, USA khuri@juiter.sjsu.edu Fax: (40)94-500 Keywords:
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationLocal Discontinuous Galerkin Methods for the Khokhlov Zabolotskaya Kuznetzov Equation
J Sci Comut (7 73:593 66 DOI.7/s95-7-5-z Local Discontinuous Galerkin Methods for the Khokhlov Zabolotskaya Kuznetzov Equation Ching-Shan Chou Weizhou Sun Yulong Xing He Yang Received: 5 Setember 6 / Revised:
More informationarxiv: v2 [math.na] 6 Apr 2016
Existence and otimality of strong stability reserving linear multiste methods: a duality-based aroach arxiv:504.03930v [math.na] 6 Ar 06 Adrián Németh January 9, 08 Abstract David I. Ketcheson We rove
More informationSpectral Properties of Schrödinger-type Operators and Large-time Behavior of the Solutions to the Corresponding Wave Equation
Math. Model. Nat. Phenom. Vol. 8, No., 23,. 27 24 DOI:.5/mmn/2386 Sectral Proerties of Schrödinger-tye Oerators and Large-time Behavior of the Solutions to the Corresonding Wave Equation A.G. Ramm Deartment
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 informationMontgomery self-imaging effect using computer-generated diffractive optical elements
Otics Communications 225 (2003) 13 17 www.elsevier.com/locate/otcom Montgomery self-imaging effect using comuter-generated diffractive otical elements J urgen Jahns a, *, Hans Knuertz a, Adolf W. Lohmann
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 informationFocal Waveform of a Prolate-Spheroidal IRA
Sensor and Simulation Notes Note 59 February 6 Focal Waveform of a Prolate-Sheroidal IRA Carl E. Baum University of New Mexico Deartment of Electrical and Comuter Engineering Albuquerque New Mexico 873
More informationREFLECTION AND TRANSMISSION BAND STRUCTURES OF A ONE-DIMENSIONAL PERIODIC SYSTEM IN THE PRESENCE OF ABSORPTION
Armenian Journal of Physics, 0, vol. 4, issue,. 90-0 REFLECTIO AD TRASMISSIO BAD STRUCTURES OF A OE-DIMESIOAL PERIODIC SYSTEM I THE PRESECE OF ABSORPTIO A. Zh. Khachatrian State Engineering University
More informationSpectral Analysis by Stationary Time Series Modeling
Chater 6 Sectral Analysis by Stationary Time Series Modeling Choosing a arametric model among all the existing models is by itself a difficult roblem. Generally, this is a riori information about the signal
More informationarxiv: v1 [quant-ph] 22 Apr 2017
Quaternionic Quantum Particles SERGIO GIARDINO Institute of Science and Technology, Federal University of São Paulo (Unifes) Avenida Cesare G. M. Lattes 101, 147-014 São José dos Camos, SP, Brazil arxiv:1704.06848v1
More informationVariable Selection and Model Building
LINEAR REGRESSION ANALYSIS MODULE XIII Lecture - 38 Variable Selection and Model Building Dr. Shalabh Deartment of Mathematics and Statistics Indian Institute of Technology Kanur Evaluation of subset regression
More informationON MINKOWSKI MEASURABILITY
ON MINKOWSKI MEASURABILITY F. MENDIVIL AND J. C. SAUNDERS DEPARTMENT OF MATHEMATICS AND STATISTICS ACADIA UNIVERSITY WOLFVILLE, NS CANADA B4P 2R6 Abstract. Two athological roerties of Minkowski content
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
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 informationHENSEL S LEMMA KEITH CONRAD
HENSEL S LEMMA KEITH CONRAD 1. Introduction In the -adic integers, congruences are aroximations: for a and b in Z, a b mod n is the same as a b 1/ n. Turning information modulo one ower of into similar
More informationSolution sheet ξi ξ < ξ i+1 0 otherwise ξ ξ i N i,p 1 (ξ) + where 0 0
Advanced Finite Elements MA5337 - WS7/8 Solution sheet This exercise sheets deals with B-slines and NURBS, which are the basis of isogeometric analysis as they will later relace the olynomial ansatz-functions
More informationSUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS *
Journal of Comutational Mathematics Vol.8, No.,, 48 48. htt://www.global-sci.org/jcm doi:.48/jcm.9.-m6 SUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More informationProof: We follow thearoach develoed in [4]. We adot a useful but non-intuitive notion of time; a bin with z balls at time t receives its next ball at
A Scaling Result for Exlosive Processes M. Mitzenmacher Λ J. Sencer We consider the following balls and bins model, as described in [, 4]. Balls are sequentially thrown into bins so that the robability
More informationJournal of Computational and Applied Mathematics. Numerical modeling of unsteady flow in steam turbine stage
Journal of Comutational and Alied Mathematics 234 (2010) 2336 2341 Contents lists available at ScienceDirect Journal of Comutational and Alied Mathematics journal homeage: www.elsevier.com/locate/cam Numerical
More informationA. G. Falkowski Chrysler Corporation, Detroit, Michigan 48227
Three-ass mufflers with uniform erforations A. Selamet V. Easwaran The Ohio State University, Deartment of Mechanical Engineering Center for Automotive Research, 06 West 18th Avenue, Columbus, Ohio 4310
More information2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized
BIT 38 :2 (998), 283{292. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES II KARL ENTACHER y Deartment of Mathematics, University of Salzburg, Hellbrunnerstr. 34 A-52 Salzburg,
More informationJournal of Mathematical Analysis and Applications
J. Math. Anal. Al. 44 (3) 3 38 Contents lists available at SciVerse ScienceDirect Journal of Mathematical Analysis and Alications journal homeage: www.elsevier.com/locate/jmaa Maximal surface area of a
More informationA CONCRETE EXAMPLE OF PRIME BEHAVIOR IN QUADRATIC FIELDS. 1. Abstract
A CONCRETE EXAMPLE OF PRIME BEHAVIOR IN QUADRATIC FIELDS CASEY BRUCK 1. Abstract The goal of this aer is to rovide a concise way for undergraduate mathematics students to learn about how rime numbers behave
More informationSCHUR S LEMMA AND BEST CONSTANTS IN WEIGHTED NORM INEQUALITIES. Gord Sinnamon The University of Western Ontario. December 27, 2003
SCHUR S LEMMA AND BEST CONSTANTS IN WEIGHTED NORM INEQUALITIES Gord Sinnamon The University of Western Ontario December 27, 23 Abstract. Strong forms of Schur s Lemma and its converse are roved for mas
More informationAutomatic Generation and Integration of Equations of Motion for Linked Mechanical Systems
Automatic Generation and Integration of Equations of Motion for Linked Mechanical Systems D. Todd Griffith a, John L. Junkins a, and James D. Turner b a Deartment of Aerosace Engineering, Texas A&M University,
More informationOn Optimization of Power Coefficient of HAWT
Journal of Power and Energy Engineering, 14,, 198- Published Online Aril 14 in Scies htt://wwwscirorg/journal/jee htt://dxdoiorg/1436/jee1448 On Otimization of Power Coefficient of HAWT Marat Z Dosaev
More informationMobius Functions, Legendre Symbols, and Discriminants
Mobius Functions, Legendre Symbols, and Discriminants 1 Introduction Zev Chonoles, Erick Knight, Tim Kunisky Over the integers, there are two key number-theoretic functions that take on values of 1, 1,
More informationCOMMUNICATION BETWEEN SHAREHOLDERS 1
COMMUNICATION BTWN SHARHOLDRS 1 A B. O A : A D Lemma B.1. U to µ Z r 2 σ2 Z + σ2 X 2r ω 2 an additive constant that does not deend on a or θ, the agents ayoffs can be written as: 2r rθa ω2 + θ µ Y rcov
More informationA Note on the Positive Nonoscillatory Solutions of the Difference Equation
Int. Journal of Math. Analysis, Vol. 4, 1, no. 36, 1787-1798 A Note on the Positive Nonoscillatory Solutions of the Difference Equation x n+1 = α c ix n i + x n k c ix n i ) Vu Van Khuong 1 and Mai Nam
More informationBoundary regularity for elliptic problems with continuous coefficients
Boundary regularity for ellitic roblems with continuous coefficients Lisa Beck Abstract: We consider weak solutions of second order nonlinear ellitic systems in divergence form or of quasi-convex variational
More information#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS
#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS Ramy F. Taki ElDin Physics and Engineering Mathematics Deartment, Faculty of Engineering, Ain Shams University, Cairo, Egyt
More informationMODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS
MODEL-BASED MULIPLE FAUL DEECION AND ISOLAION FOR NONLINEAR SYSEMS Ivan Castillo, and homas F. Edgar he University of exas at Austin Austin, X 78712 David Hill Chemstations Houston, X 77009 Abstract A
More informationConstruction algorithms for good extensible lattice rules
Construction algorithms for good extensible lattice rules Harald Niederreiter and Friedrich Pillichshammer Dedicated to Ian H. Sloan on the occasion of his 70th birthday Abstract Extensible olynomial lattice
More informationarxiv: v1 [hep-th] 6 Oct 2017
Dressed infrared quantum information Daniel Carney, Laurent Chaurette, Domini Neuenfeld, and Gordon Walter Semenoff Deartment of Physics and Astronomy, University of British Columbia, BC, Canada We study
More informationThe ALaDyn Code Technical Report
The ALaDyn Code Technical Reort P. Londrillo, A.Sgattoni and F. Rossi INAF and INFN, Sezione Bologna Diartimento di Fisica Università di Bologna May 18, 2011 Contents 1 The Maxwell-Vlasov system using
More informationMultiplicity of weak solutions for a class of nonuniformly elliptic equations of p-laplacian type
Nonlinear Analysis 7 29 536 546 www.elsevier.com/locate/na Multilicity of weak solutions for a class of nonuniformly ellitic equations of -Lalacian tye Hoang Quoc Toan, Quô c-anh Ngô Deartment of Mathematics,
More informationQuantitative estimates of propagation of chaos for stochastic systems with W 1, kernels
oname manuscrit o. will be inserted by the editor) Quantitative estimates of roagation of chaos for stochastic systems with W, kernels Pierre-Emmanuel Jabin Zhenfu Wang Received: date / Acceted: date Abstract
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 informationA Closed-Form Solution to the Minimum V 2
Celestial Mechanics and Dynamical Astronomy manuscrit No. (will be inserted by the editor) Martín Avendaño Daniele Mortari A Closed-Form Solution to the Minimum V tot Lambert s Problem Received: Month
More informationEstimation of Separable Representations in Psychophysical Experiments
Estimation of Searable Reresentations in Psychohysical Exeriments Michele Bernasconi (mbernasconi@eco.uninsubria.it) Christine Choirat (cchoirat@eco.uninsubria.it) Raffaello Seri (rseri@eco.uninsubria.it)
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Jul 1998
arxiv:cond-mat/98773v1 [cond-mat.stat-mech] 5 Jul 1998 Floy modes and the free energy: Rigidity and connectivity ercolation on Bethe Lattices P.M. Duxbury, D.J. Jacobs, M.F. Thore Deartment of Physics
More informationIrregular reflection of spark-generated shock pulses from a rigid surface: Mach-Zehnder interferometry measurements in air
Irregular reflection of sark-generated shock ulses from a rigid surface: Mach-Zehnder interferometry measurements in air Maria M. Karzova, 1,a) Thomas Lechat, 2 Sebastien Ollivier, 3 Didier Dragna, 2 Petr
More informationThe Mathematics of Thermal Diffusion Shocks
The Mathematics of Thermal Diffusion Shocks Vitalyi Gusev,WalterCraig 2, Roberto LiVoti 3, Sorasak Danworahong 4, and Gerald J. Diebold 5 Université du Maine, av. Messiaen, 7285 LeMans, Cedex 9 France,
More informationINDIRECT PLANETARY CAPTURE VIA PERIODIC ORBITS ABOUT LIBRATION POINTS
INDIRECT PLANETARY CAPTURE VIA PERIODIC ORBITS ABOUT LIBRATION POINTS Li Xiangyu 1,2, Qiao Dong 1,2, Cui Pingyuan 1,2 (1. Institute of Dee Sace Exloration Technology, Beijing Institute of Technology, Beijing,
More informationOptimal Recognition Algorithm for Cameras of Lasers Evanescent
Otimal Recognition Algorithm for Cameras of Lasers Evanescent T. Gaudo * Abstract An algorithm based on the Bayesian aroach to detect and recognise off-axis ulse laser beams roagating in the atmoshere
More informationStudy on a Ship with 6 Degrees of Freedom Inertial Measurement System and Related Technologies
Oen Access Library Journal Study on a Shi with 6 Degrees of Freedom Inertial Measurement System and Related echnologies Jianhui Lu,, Yachao Li, Shaonan Chen, Yunxia Wu Shandong rovince Key Laboratory of
More informationEffective conductivity in a lattice model for binary disordered media with complex distributions of grain sizes
hys. stat. sol. b 36, 65-633 003 Effective conductivity in a lattice model for binary disordered media with comlex distributions of grain sizes R. PIASECKI Institute of Chemistry, University of Oole, Oleska
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationDetermination of the Best Apodization Function and Grating Length of Linearly Chirped Fiber Bragg Grating for Dispersion Compensation
84 JOURNAL OF COMMUNICATIONS, VOL. 7, NO., NOVEMBER Determination of the Best Aodization Function and Grating Length of Linearly Chired Fiber Bragg Grating for Disersion Comensation Sher Shermin A. Khan
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationThe non-stochastic multi-armed bandit problem
Submitted for journal ublication. The non-stochastic multi-armed bandit roblem Peter Auer Institute for Theoretical Comuter Science Graz University of Technology A-8010 Graz (Austria) auer@igi.tu-graz.ac.at
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationFactors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier
Australian Journal of Basic and Alied Sciences, 5(12): 2010-2020, 2011 ISSN 1991-8178 Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doed Fiber Amlifier
More informationAn Inverse Problem for Two Spectra of Complex Finite Jacobi Matrices
Coyright 202 Tech Science Press CMES, vol.86, no.4,.30-39, 202 An Inverse Problem for Two Sectra of Comlex Finite Jacobi Matrices Gusein Sh. Guseinov Abstract: This aer deals with the inverse sectral roblem
More information