Hyperbolic Systems of Equations Posed on Erroneous Curved Domains
|
|
- Frederick Garrison
- 5 years ago
- Views:
Transcription
1 Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE Linköping, Sween ( jan.norstrom@liu.se). b Department of Mathematics, Computational Mathematics, Linköping University, SE Linköping, Sween (samira.nikkar@liu.se). Abstract The effect of an inaccurate geometry escription on the solution accuracy of a hyperbolic problem is iscusse. The inaccurate geometry can for eample come from an imperfect CAD system, a faulty mesh generator, ba measurements or simply a misconception. We show that inaccurate geometry escriptions might lea to the wrong wave spees, a misplacement of the bounary conitions, to the wrong bounary operator an a mismatch of bounary ata. The errors cause by an inaccurate geometry escription may affect the solution more than the accuracy of the specific iscretization techniques use. In etreme cases, the orer of accuracy goes to zero. Numerical eperiments corroborate the theoretical results.. Erroneous computational omain Consier the following hyperbolic system of equations, in two space imensions, W t + ÂW + ˆBW y =, (, y) Ω, t (, T ], LW = g(, y, t), (, y) δω, t (, T ], () W = f (, y), (, y) Ω, t =, in which the solution is represente by the vector W = W(, y, t). Â an ˆB are constant symmetric M M matrices, Ω is the spatial omain with the bounary δω. The bounary operator L is efine on δω, f (, y) R M an g(, y, t) R M are the ata in the problem. Preprint submitte to Journal of Computational Physics January 3, 6
2 Equation () is transforme to curvilinear coorinates (ξ, η) by (V ξ, V η, V t ) = [J](V, V y, V t ) T, where [J] is the Jacobian matri of the transformation. The transforme problem is JW t + AW ξ + BW η =, (ξ, η) Φ, t (, T ], LW = g(ξ, η, t), (ξ, η), t (, T ], W = f (ξ, η), (ξ, η) Φ, t =, () where A = Jξ Â + Jξ y ˆB, B = Jη Â + Jη y ˆB an J = ξ y η η y ξ > is the eterminant of [J]. The energy metho together with the metric ientities an the use of the Green-Gauss theorem yiels t W(ξ, η, t) J = W T C W s, (3) where the norm is efine by V J = Φ V T J V ξ η. In (3), C = (A, B) n = αa + βb = XΛX T, n = (α, β) is the unit normal vector pointing outwar from Φ, X contains the eigenvectors as columns an Λ contains the eigenvalues on the main iagonal. For more etails see []. The problem () is well-pose if we impose characteristic bounary conitions an the following energy rate is obtaine t W(ξ, η, t) J + W T C + W s = g T Λ g s, (4) where Λ contains the negative eigenvalues of C an C + = XΛ + C T... Errors ue to the wrong position of bounary Consier two hyperbolic problems with solutions U an V pose on two nearby omains. The two omains, epicte in Figure, are both mappe to the unit square. We consier Ω c to be the correct omain an Ω e to be the erroneous one (subscripts c an e enote correct an erroneous, respectively). The transforme equations become J c U t + A c U ξ + B c U η =, J e V t + A e V ξ + B e V η =, (5) where the matrices are A i = [Jξ Â + Jξ y ˆB] i an B i = [Jη Â + Jη y ˆB] i, for i = c,e. Our first result follows immeiately. Since A c A e, B c B e we realize that: ) The wave-spees given by the eigenvalues of the matrices will iffer, an that V will have the wrong wave spees.
3 y η Ω c Φ Ω e ξ)ξ) Figure : A schematic of two ifferent geometry escriptions both mappe to the unit square. ξ Net, we consier the results of the energy metho given by (4) where W {U,V } an C i = α i A i + β i B i for i = c,e. A closer look at the bounary term C i reveals that the normals as well as the eigenvalues of C i are moifie by the erroneous geometry. This paves the way for the secon, thir an fourth conclusions: ) Error in the normal: This error is cause by imposing the wrong bounary operator ue to an erroneous normal. Consier for eample the soli wall no penetration conition (u,v) n = αu + βv = g = for the Euler equations. An erroneous normal will lea to the wrong bounary operator. This means that the bounary operator L in () is wrong while the ata g is correct. See case in Figure. 3) Error in the position: Here the error is ue to the misplacement of the bounary conition. Consier again the soli wall no penetration conition for the Euler equations, now with a correct bounary operator an ata, impose at the wrong position in space. In this case, the bounary operator L an the ata g in () are both correct. See case in Figure. 4) Error in the ata: Bounary conitions with ata from Ω c but impose at Ω e will also lea to inaccurate results. In this case, the ata g in () is wrong while the bounary operator L can be either correct or wrong. See case an 3 in Figure. We can also have a combination of all the errors. These errors might be, an often are [4], more important than the orer of accuracy of the specific iscretization techniques use (which will be shown later in the numerical eperiments section below). To further iscuss the effect of the wrong position of bounary conitions, consier () an two types of bounary ata, correct an erroneous enote by g c an g e, respectively. The correct an erroneous solutions are enote by U an V, respectively, an for simplicity we assume that they have the same initial 3
4 (c) = correct geometry (e) = erroneous geometry ne nc ne nc ne 3 nc Possible error sources: normal wrong at right position normal right at wrong position 3 normal wrong at wrong position Figure : A schematic of the erroneous an correct (n c = correct normal, n e = erroneous normal) bounary efinitions. ata f which vanishes close to the bounaries. By subtraction we obtain an error equation with zero initial ata an non-zero bounary ata. Let E = U V enote the error, an δg = g e g c = ( g) e + O( ) = O( ). (6) In (6), = (, y), = c e, y = y c y e, see Figure 3. (c) = correct geometry (e) = erroneous geometry (!, y! )!!!! (!, y! ) Figure 3: A blow-up of two nearby bounary points which result in two types of ata, g c = g( c,y c,t) an g e = g( e,y e,t). The energy metho applie to the transforme version of the error equation gives the corresponing energy rate to (4), as t E J e + E T C e + E s = δg T Λ e δg s, (7) 4
5 y y in which Λ e contain the negative eigenvalues of C e. From (7) we can conclue that the error is boune for long time calculations, see [3] for etails.. Numerical eperiments To quantify the error stemming from an erroneous geometry we consier the geometries in Figures 4 an 5. The geometry escription is given by.. a a' b.6.4 ' b'.. c c' Figure 4: The correct omain, Ω c Figure 5: The erroneous omain, Ω e. y a () = a sin(ωπ) +, Ω c : b (y) = a sin(ωπy) +, y c () = a sin(ωπ), (y) = a sin(ωπy), y a () = y a () â sin( ˆωπ), Ω e : b (y) = b (y) + â sin( ˆωπy), y c () = y c () â sin( ˆωπ), (y) = (y) + â sin( ˆωπy), (8) where a =., ω = an ˆω =. Moreover, â =. for q = an â = q for q {,, 3}. Note that Ω e in Figure 5 is eaggerate for a better visualization, an Ω c is inclue for clarity. It is not use uring the numerical eperiments. We consier the two-imensional constant coefficient symmetrize Euler equations in which V =[ρ, u, v, T ] T is the perturbe solution. The components ρ, u, v, an T are respectively the ensity, the an y velocity components an the temperature perturbation. The matrices  an ˆB are given in []. To verify the accuracy of our computational metho, we use the manufacture solution V = [sin(θ), cos(θ), sin(φ), cos(φ)] T, where θ(,t) = 5t, φ(y,t) = 5y t. The manufacture solution is injecte in the problem through a forcing function, bounary an initial ata. We use characteristic bounary conitions as in the erivations above. The convergence rates using a finite ifference operator on summation-by-parts form (SBP84) in space an time with weak 5
6 Convergence Rates, P Energy Norm of the Errors bounary an initial conitions, are shown in Table. SBP84 in space an time converges with 5th orer an the convergence rates are clearly correct []. Table : Convergence rates for the problem () with correct ata Number of gri points The ρ component The u component The v component The T component Equippe with a 5th orer convergent computational scheme (we use SBP84 in space an time, see [] for etails on this technique), we procee to investigate the problem with erroneous ata. We compute the erroneous solution by computing on the erroneous omain Ω e, taking the bounary ata from the nearby omain Ω c an finally subtracting the eact manufacture solution. The convergence rates for the erroneous solution are shown in Figure 6. We have consiere ifferent magnitues for the eviation of the erroneous omain from the correct omain in terms of = O(h q ) where h is the gri spacing an q {,,, 3}. As seen in Figure 6, imposing the bounary ata at an erroneous position affects the results an reuces the global orer accuracy of the scheme to q. Finally, in Figure 7 we show how the error behaves in long time calculations. It reaches an error boun as preicte above, since we use characteristic bounary conitions as shown in [3]. Convergence Rates for Erroneous Solutions Long Time Calculations q= q= q= q=3..9 Eact q= q= q= N Figure 6: Convergence rates for the erroneous solutions in Ω e with q {,,,3}, T= an sufficiently small time steps Time Figure 7: Error versus time for eact an erroneous solutions with q {,,3}, 3 an 5 gri points in space an time, respectively. 6
7 3. Summary an conclusions The effects of an erroneous geometry escription were iscusse by using the energy metho as the analytical tool. We erive an error estimate which escribes the effects of imposing correct bounary ata at erroneous bounary positions. We use a 5th orer provable stable an convergent high orer finite ifference metho as our computational tool. The results obtaine are vali for all iscretization techniques, that are sufficiently high orer accurate an stable. We conclue that an erroneous geometry escription might lea to the wrong wave spees, a misplacement of the bounary conitions, to the wrong bounary operator an a mismatch of bounary ata. The errors cause by an inaccurate geometry escription may affect the solution more than the accuracy of the specific iscretization techniques use. In etreme cases, the orer of accuracy goes to zero. The numerical eperiments corroborate this conclusion. 4. References [] S. Nikkar an J. Norström, Fully iscrete energy stable high orer finite ifference methos for hyperbolic problems in eforming omains, Journal of Computational Physics, 9():8-98 (5). [] M. Svär an J. Norström, Review of summation-by-parts schemes for initial-bounary-value problems, Journal of Computational Physics, Volume 68, pp. 738 (4). [3] J. Norström, Error boune schemes for time-epenent hyperbolic problems, SIAM Journal of Scinetific Computing, 3:46-59, (7). [4] I. Babuska an J. Chleboun, Effects of uncertanities in the omain on the solution of Neumann bounary value problems in two spatial imensions, Math. Comput., 7:339-37, (). 7
COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationStable finite difference schemes for the magnetic induction equation with Hall effect
Eigenössische Technische Hochschule Zürich Ecole polytechnique féérale e Zurich Politecnico feerale i Zurigo Swiss Feeral Institute of Technology Zurich Stable finite ifference schemes for the magnetic
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationAbstract A nonlinear partial differential equation of the following form is considered:
M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationLower bounds on Locality Sensitive Hashing
Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationLecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.
b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationThe Sokhotski-Plemelj Formula
hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationDistribution Theory for Discontinuous Test Functions and Differential Operators with Generalized Coefficients
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 1, 9733 1996 ARTICLE NO 56 Distribution Theory for Discontinuous Test Functions an Differential Operators with Generalize Coefficients P Kurasov* Department
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationDepartment of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ
Department of Physics University of arylan College Park, arylan PHYSICS 4 Prof. S. J. Gates Fall 5 Final Exam Dec. 6, 5 This is a OPEN book examination. Rea the entire examination before you begin to work.
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationMULTISCALE FRICTION MODELING FOR SHEET METAL FORMING
MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationBOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION CHEMOTAXIS SYSTEM WITH NONLINEAR DIFFUSION AND LOGISTIC SOURCE
Electronic Journal of Differential Equations, Vol. 016 (016, No. 176, pp. 1 1. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu BOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More information1. At time t = 0, the wave function of a free particle moving in a one-dimension is given by, ψ(x,0) = N
Physics 15 Solution Set Winter 018 1. At time t = 0, the wave function of a free particle moving in a one-imension is given by, ψ(x,0) = N where N an k 0 are real positive constants. + e k /k 0 e ikx k,
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationSpectral properties of a near-periodic row-stochastic Leslie matrix
Linear Algebra an its Applications 409 2005) 66 86 wwwelseviercom/locate/laa Spectral properties of a near-perioic row-stochastic Leslie matrix Mei-Qin Chen a Xiezhang Li b a Department of Mathematics
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationS10.G.1. Fluid Flow Around the Brownian Particle
Rea Reichl s introuction. Tables & proofs for vector calculus formulas can be foun in the stanar textbooks G.Arfken s Mathematical Methos for Physicists an J.D.Jackson s Classical Electroynamics. S0.G..
More informationMathematical Basics. Chapter Introduction and Definitions
Chapter 2 Mathematical Basics 2.1 Introuction an Definitions Flui mechanics eals with transport processes, especially with the flow- an molecule-epenent momentum transports in fluis. Their thermoynamic
More informationRank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col
Review of Linear Algebra { E18 Hanout Vectors an Their Inner Proucts Let X an Y be two vectors: an Their inner prouct is ene as X =[x1; ;x n ] T Y =[y1; ;y n ] T (X; Y ) = X T Y = x k y k k=1 where T an
More informationText S1: Simulation models and detailed method for early warning signal calculation
1 Text S1: Simulation moels an etaile metho for early warning signal calculation Steven J. Lae, Thilo Gross Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, 01187 Dresen, Germany
More informationDomain Decomposition and Model Reduction of Systems with Local Nonlinearities
Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen Abstract The goal of this paper is to combine balance truncation
More informationME338A CONTINUUM MECHANICS
global vs local balance equations ME338A CONTINUUM MECHANICS lecture notes 11 tuesay, may 06, 2008 The balance equations of continuum mechanics serve as a basic set of equations require to solve an initial
More informationMany problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx
Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus
More informationDesigning Information Devices and Systems II Spring 2018 J. Roychowdhury and M. Maharbiz Discussion 2A
EECS 6B Designing Information Devices an Systems II Spring 208 J. Roychowhury an M. Maharbiz Discussion 2A Secon-Orer Differential Equations Secon-orer ifferential equations are ifferential equations of
More informationCHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold
CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the
More informationdiv(λ(x) u) = f in Ω, (1) u = 0 on Ω, (2) where we denote by Ω = Ω \ Ω the boundary of the domain Ω, under the following assumptions:
Discretisation of heterogeneous an anisotropic iffusion problems on general nonconforming meshes SUSHI: a scheme using stabilisation an hybri interfaces 1 R. Eymar 2, T. Gallouët 3 an R. Herbin 4 Abstract:
More informationMath 1272 Solutions for Spring 2005 Final Exam. asked to find the limit of the sequence. This is equivalent to evaluating lim. lim.
Math 7 Solutions for Spring 5 Final Exam ) We are gien an infinite sequence for which the general term is a n 3 + 5n n + n an are 3 + 5n aske to fin the limit of the sequence. This is equialent to ealuating
More informationA Spectral Method for the Biharmonic Equation
A Spectral Metho for the Biharmonic Equation Kenall Atkinson, Davi Chien, an Olaf Hansen Abstract Let Ω be an open, simply connecte, an boune region in Ê,, with a smooth bounary Ω that is homeomorphic
More informationMagnetic helicity evolution in a periodic domain with imposed field
PHYSICAL REVIEW E 69, 056407 (2004) Magnetic helicity evolution in a perioic omain with impose fiel Axel Branenburg* Norita, Blegamsvej 17, DK-2100 Copenhagen Ø, Denmark William H. Matthaeus University
More informationKinematics of Rotations: A Summary
A Kinematics of Rotations: A Summary The purpose of this appenix is to outline proofs of some results in the realm of kinematics of rotations that were invoke in the preceing chapters. Further etails are
More informationThe Sokhotski-Plemelj Formula
hysics 24 Winter 207 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationarxiv: v4 [cs.ds] 7 Mar 2014
Analysis of Agglomerative Clustering Marcel R. Ackermann Johannes Blömer Daniel Kuntze Christian Sohler arxiv:101.697v [cs.ds] 7 Mar 01 Abstract The iameter k-clustering problem is the problem of partitioning
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationPure Further Mathematics 1. Revision Notes
Pure Further Mathematics Revision Notes June 20 2 FP JUNE 20 SDB Further Pure Complex Numbers... 3 Definitions an arithmetical operations... 3 Complex conjugate... 3 Properties... 3 Complex number plane,
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationPHYS 414 Problem Set 2: Turtles all the way down
PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
More informationChapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies
hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation
More informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationLeChatelier Dynamics
LeChatelier Dynamics Robert Gilmore Physics Department, Drexel University, Philaelphia, Pennsylvania 1914, USA (Date: June 12, 28, Levine Birthay Party: To be submitte.) Dynamics of the relaxation of a
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationTechnion - Computer Science Department - M.Sc. Thesis MSC Constrained Codes for Two-Dimensional Channels.
Technion - Computer Science Department - M.Sc. Thesis MSC-2006- - 2006 Constraine Coes for Two-Dimensional Channels Keren Censor Technion - Computer Science Department - M.Sc. Thesis MSC-2006- - 2006 Technion
More informationAn Eulerian level set method for partial differential equations on evolving surfaces
Computing an Visualization in Science manuscript No. (will be inserte by the eitor) G. Dziuk C. M. Elliott An Eulerian level set metho for partial ifferential equations on evolving surfaces Receive: ate
More informationSOLUTIONS for Homework #3
SOLUTIONS for Hoework #3 1. In the potential of given for there is no unboun states. Boun states have positive energies E n labele by an integer n. For each energy level E, two syetrically locate classical
More informationOn some parabolic systems arising from a nuclear reactor model
On some parabolic systems arising from a nuclear reactor moel Kosuke Kita Grauate School of Avance Science an Engineering, Wasea University Introuction NR We stuy the following initial-bounary value problem
More informationcosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) 2t 1 + t 2 cos x = 1 t2 sin x =
6.4 Integration using tan/ We will revisit the ouble angle ientities: sin = sin/ cos/ = tan/ sec / = tan/ + tan / cos = cos / sin / tan = = tan / sec / tan/ tan /. = tan / + tan / So writing t = tan/ we
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationA. Exclusive KL View of the MLE
A. Exclusive KL View of the MLE Lets assume a change-of-variable moel p Z z on the ranom variable Z R m, such as the one use in Dinh et al. 2017: z 0 p 0 z 0 an z = ψz 0, where ψ is an invertible function
More information