Parameter estimation for Euler equations with uncertain inputs
|
|
- Jerome Terry
- 5 years ago
- Views:
Transcription
1 Parameter estimation for Euler equations wit uncertain inputs Sergiy Zuk Abstract Te paper presents a new state estimation algoritm for 2D incompressible Euler equations wit periodic boundary conditions and uncertain but bounded inputs and initial conditions. Te algoritm converges (in L 2 -sense) to a least squares estimator given incomplete and noisy observations. It is also sown experimentally tat te proposed algoritm applies to several conservation laws developing sock discontinuities. Te results are illustrated by numerical examples. I. INTRODUCTION Data assimilation algoritms represent a backbone of modern cyber-pysical systems. Tese algoritms allow one to optimally combine a priori knowledge encoded in equations of matematical pysics wit a posteriori information in te form of sensor readings. Weater forecasting is one of many examples were data assimilation is applied to improve predictions generated by ydrodynamical models. Te divergence-free Euler equations provide te most basic model of incompressible flows of omogeneous inviscid fluids [7], yet tis model may provide insigts in studies of turbulence. We refer te reader to [2] for a detailed overview of various matematical questions related to Euler equations. In tis paper we design a new state estimator for 2D Euler equations subject to uncertain but bounded input and unknown initial conditions. We rely on a vorticity-stream formulation of te Euler equations. Te resulting vorticity equation describes te rotation of te vorticity of te fluid velocity field (see Section II). We stress tat te omogenous vorticity equation possesses a nice property: te L 2 -norm of te vorticity is conserved over time. Te latter, in fact, allows one to prove existence and uniqueness IBM Researc, Dublin, Ireland, {sergiy.zuk,tigran}@ie.ibm.com Tigran Tcrakian of te solution of 2D Euler equations globally (in time) [7, Cor.3.3, p.6]. Assuming periodic boundary conditions, we apply Fourier-Galerkin (FG) approximation: we project te vorticity equation onto a subspace generated by {e inx } n N and obtain an ODE for te projection coefficients, a FG model. Note tat Fourier-Galerkin approximation possesses a spectral convergence rate provided te solution of te Euler equation is smoot []. Altoug te main results of tis work may be derived for te case of bounded domains wit non-penetration boundary conditions or unbounded domains (R 2 ), we focus on smoot periodic flows to simplify te convergence analysis. We refer te reader to Section II, were oter types of boundary conditions and weak solutions of Euler equations are discussed. Given an ODE representing projection coefficients, we design a discrete-in-time FG model suc tat te L 2 -norm of te discrete vorticity is conserved over time (see Section III-B). Tis property allows us to prove te convergence of te approximations provided by te discrete FG model to a solution of te continuous FG model. We ten derive te state estimator (in te form of a minimax filter) for te discrete FG model and note tat te combination of te mentioned convergence proof wit results of [6] may be used to derive a continuous formulation of te minimax filter (see Section III- C). To te best of our knowledge, tis result is new. A similar approac as been used to design data assimilation algoritms for scalar macroscopic traffic flow models [9] and flood models (StVenant equations) []. We stress tat te proposed metod sows very good performance not only for smoot periodic flows but also for conservation laws wit sock discontinuities and non-periodic boundary conditions (see Section IV). Data assimilation for systems of yperbolic conservation laws based on te calculus of variations
2 was proposed in [3], were te autors adopt te strategy: optimize, ten discretize so tat te estimate of te initial density does not depend on a discretization metod. In contrast, te present paper solves te filtering problem; tat is, te state estimate at time instant, t, depends on te observation at te same time instant and te previous estimate only. A comparison of classical filtering algoritms (extended Kalman filter and ensemble Kalman filter) for scalar conservation laws wit quadratic non-linearity may be found in [4]. Adaptive parameter estimators for yperbolic equations were considered for instance in [5]. Tis paper is organized as follows. Te next section presents a brief overview of Euler equations. Section II-A presents discrete in time FG model for Euler equations; te state estimation algoritm is sown in section III-C. Numerical experiments are given in section IV. Section V contains concluding remarks. Notation. Let Ω denote a subset of R 2 wit boundary Ω, and n is a normal vector for Ω poining outside, Ω T := Ω [, T ]. C s (Ω) denotes a space of continuously differentiable functions on Ω (up to order s), L 2 (Ω) is te space of square-integrable functions on Ω, H (Ω) is a Sobolev space of functions wit weak derivatives of L 2 (Ω)-class. div( u) = x u + x2 u 2, curl( u) = x u 2 x2 u, u = ( x u, x2 u), u = ( x2 u, x u) u v denotes te canonical inner product of R 2, H 2 := H H denotes te cartesian product of H wit itself. L 2 (t, t, H) := {f : f(t) H and T f(t) 2 H < + }. We write u = v a.e. on Ω if u(x) = v(x) for almost all x Ω. II. EULER EQUATIONS Assume tat ω verifies te vorticity-stream formulation of te Euler equation: t ω + u ω = f, ψ = ω, u = ū + ψ, ω() = curl( u ), (x, t) Ω T := [ r, r] 2 [, T ], () were ū R 2 is a given vector, u C 2 (Ω) 2 is a smoot -periodic vector-function and f C (Ω T ) as zero mean, f(x, t)dx =. It is not ard to Ω prove (see [7, Prop 2.4,p.5]) tat u = (u, u 2 ) C 2 (Ω T ) 2 is a smoot -periodic solution of te incompressible Euler equation on Ω T : d u + ( u ) u + p = f, div( u) =, u() = u + ū, (2) were f C 2 (Ω T ) and curl( f) = f, and te pressure p is a function of u and f. We note tat Euler equation (2) as te unique smoot solution u for Ω = R 2 provided f = and u is a smoot function suc tat div( u ) = and curl( u ) L (R 2 ) (see [7, L.3.2, p.93 and Cor.3.3, p.6]). For te case of bounded domains wit smoot boundary one may also derive existence of smoot solutions for (2) provided u and f are smoot and u n = on Ω (see [, p.356]). Weak solutions (or Sobolev space solutions) for (2) were constructed in [2] provided curl( u ) is bounded, and f C(, T, L p (Ω)) is so tat curl( f) is bounded. Solutions of L 2 (, T, L 2 (R 2 ) 2 )-class corresponding to so called vortex seets were discussed in [7, p.33]. A. Fourier-Galerkin approximation Define b( u, w, v) := ( u w)vdx. Ω Assume tat div( u) = and u, w, v are smoot -periodic functions on Ω. We find integrating by parts: b( u, w, v) = b( u, v, w). (3) Now, we apply Fourier spectral metod to construct a finite dimensional approximation of (). Let s Z 2 and define φ s (x) := eiπr s x 2r. It is known tat {φ s } s Z 2 is a total ortonormal system in L 2 (Ω) and (φ s, φ k ) L2 (Ω) = δ s,k δ s2,k 2 were s i, k i are components of s, k. Now, we multiply () by φ s
3 and integrate over Ω. We obtain a weak vorticity formulation: d (ω, φ s) L2 (Ω) + b( u, ω, φ s ) = (f, φ s ) L2 (Ω). Let us now set ω N (x, t) := k ω i N k(t)φ k (x) wit ω k (t) := (ω(t), φ k ) L 2 (Ω) and plug ω N into te above formulation. We get te following ODE for te coefficients: ω s + b( u, φ k, φ s )ω k (t) = (f, φ s ) L 2 (Ω). k i N Define ω(t) := {ω s (t)} s,2 N, f(t) := {(f(t), φ s ) L2 (Ω)} s,2 N and set B( ω) := {b( u, φ k, φ s )} s,2, k,2 N. Now we use tis notation to rewrite te finite dimensional vorticity formulation in te vector form: d ω + B( ω) ω = f wit initial condition ω() = ω were ω := {ω s ()} s,2 N. We stress tat B( ω) = B ( ω) as b( u, φ k, φ s ) = b( u, φ s, φ k ) by (3) so tat B( w) is skew-symmetric. A. Problem statement Let ω solve III. MAIN RESULTS d ω + B( ω) ω = f, ω() = ω, (4) and assume tat a vector-function y is observed in te following form ( k,2 M): y k (t) = (H k, φ s ) L 2 (Ω)ω s (t) + η k (t), (5) s,2 N were H k L 2 (Ω) and η = {η k } k,2 M is a measurable vector-function modelling noise. Our aim is to construct a state estimate ˆω(T ) for ω(t ) given data y and assuming tat ω S ω + T f Q f + η R η (6) provided S, Q, R are symmetric positive definite matrices of appropriate dimensions. Let us stress tat if ω is te solution of (4), ω N corresponds to ω and ω solves () ten according to [, T.5.] we ave: ω N (, t) ω(, t) L2 (Ω) e g(t) N 2l ω(, ) 2 H l (Ω) + e g(t) N 2 l max τ t ω(, τ) H l (Ω), g(t). (7) Tis so called spectral convergence rate justifies our coice of te state equation. Namely, (4) is a standard finite dimensional Fourier-Galerkin model wit uncertain input f. Now, following te idea of [4], [5] we may incorporate te effect of te unresolved modes (te projection error) by simply adding anoter model error term e and introduce an additional algebraic equation to filter out inadmissible e (see [8], [3]) so tat te exact projection coefficients of ω will be among te solutions of (4). On te oter and, te above convergence rate estimate suggests tat te effect of e is negligible for reasonably large N and l (so for smoot ω) and it may be terefore absorbed by f (by increasing te size of te ellipsoid (6)). B. Discrete Fourier-Galerkin model Assume f =. If we multiply () by ω and integrate over Ω we get: 2 t ω 2 L 2 (Ω) = (f, ω) L 2 (Ω) as b( u, ω, ω) = by (3). Hence, L 2 (Ω)-norm of ω is conserved. We stress tat (4) as te same property, namely ω 2 R is conserved as B( ω) is skewsymmetric. In wat follows we propose a metod 2N+ wic approximates ω(j) by w j := w(j), j =, m, := T m, m N and te norm of w j is conserved. By Newton-Leibniz formula we get: ω((j + )) = ω(j) (j+) j B( ω(τ)) ω(τ)dτ Define B j := B( ω(j)). Approximating te integral by te trapezoidal rule one gets: (I+ 2 B j+) ω((j+)) = (I 2 B j) ω(j)+o( 3 ). ω N (x, t) := k i N ω k(t)φ k (x) were te coefficients ω k are components of ω
4 By noting tat B j+ = B j +B( d ω (j))+o(2 ) we can simplify te above equation compromising te order of te approximation: specifically, approximating B j+ by B j + B( d ω (j)) we reduce te order down to O( 2 ) (locally); if we simply take B j we get O()-approximation. In wat follows we stick to te latter and define w j as a solution of te linear system: (I + 2 B j) w j+ = (I 2 B j) w j, w = ω. (8) Note tat (I + 2 B j) is invertible as B j is skewsymmetric and K j := (I + 2 B j) (I 2 B j) is te Caley transform of te skew-symmetric matrix B j. Hence, K j is an ortogonal matrix and so w j 2 R 2N+ = w 2 R 2N+. For j =,..., m and j t < (j + ) we define te following functions: U (m) wj++ wj (t) := 2, V (m) (t) = w j and W (m) (t) = w j (t j)b( w j )U (m) (t). Since w j 2 R = w 2N+ 2 R it follows tat 2N+ V (m) L2 (,T ), U (m) L2 (,T ) C < + and so W (m) L2 (,T ) C 2 < +. Terefore, te sequences of piecewise constant functions, {V (m) } m N, {U (m) } m N contain weakly convergent subsequences in L 2 (, T ). Te same olds true for te sequence {W (m) } m N. We denote te convergent subsequences by {V (mi) }, {U (mi) } and {W (mi) } respectively and let V, U and W be teir limiting functions. We claim tat all te tree mentioned sequences converge strongly and V = U = W. Indeed, dw (mi) dw (mi) = B(V (mi) (t))u (mi) (t) and so is bounded in L 2 (, T ). Hence, {W (mi) } weakly converges in H (, T ) tat implies strong convergence in L 2 (, T ). Now, V (mi) U (mi) L2 (,T ) 3 2 m 2 w 2 C 2 R 2N+ 2 were C := 2N+ j= B(e j ) 2, e j is j-t canonical basis vector. Hence, by taking weak limit (m ) and using weak lowersemicontinuity of L 2 -norm we get: V = U. On te oter and W (mi) V (mi) 2 L 2 (,T ) 3 3 m j= B( w j)u (mi) (t) 2 R 2N+ 3 3 C w 4 R 2N+. By te same argument we get: W = V. Taking weak limits in dw (mi) = B(V (mi) (t))u (mi) (t) we find tat: dw = B(W )W (t). Now, recalling te spectral convergence rate estimate (7) and noting tat tere exists unique smoot solution of () we deduce tat any convergent subsequence of U (m), V (m) and W (m) as te same limit. Te latter proves tat te entire sequences U (m), V (m) and W (m) are weakly convergent and sare te same limiting function W wic is te unique solution of (4). Let us summarize te above results. Proposition : If f = ten (i) (4) as a unique solution ω suc tat ω() R 2N+ = ω(t) R 2N+, (ii) a sequence of piecewise constant functions V (m) (t) = w(j), j t < (j + ) converges to ω in L 2 (, T ) provided w(j) solves (8) and := T m, m N. We note tat it is not ard to generalize point (ii) of te latter proposition to te case f. We omit tis generalization ere for te sake of space. C. State estimator Following [6] we introduce te following system of linear Hamiltonian equations: ( I 2 B j 2 H RH )( Uj+ 2 Q I+ 2 Bj 2 H RH ) ( I+ )( V j+ = 2 B j IPj ) 2 Q I, 2 Bj (9) were P j := U j V j for j > and P = S, and B j := B(ˆω j ), and ˆω j solves te following system
5 of linear equations: (I + 2 B j + 2 P j+h RH)ˆω j+ = (I 2 B j 2 P jh RH)ˆω j () + 2 P j+ H R y j+ + P j H R y j 2 were H := {(H k, φ s ) L 2 (Ω)} k,2, s,2 N. By combining te idea of te proof of Proposition wit a well-known stabilization effect brougt by P j, te solution of approximated Riccati equation (see [6] for te details), one may prove te following proposition. Proposition 2: Let U j, V j, P j and ˆω j be defined as above and := T m, m N. A sequence of piecewise constant functions V (m) (t) = V j, U (m) (t) = U j, ˆω (m) (t) = ˆω j for j t < (j+) converges to U, V and ˆω in L 2 (, T ). Moreover, U, V and ˆω represent te unique solution of te following system of equations: dˆω = B(ˆω)ˆω + P H RH( y H ˆω), ˆω() =, du = B (ˆω)U + H RHV, U() = I, dv = Q U B(ˆω)V, V () = S, P = V U. () In fact, te latter proposition presents a continuous version of te minimax filter (equation ()) and equations (9)-() define an approximation for te filter. In te following section we illustrate te convergence properties of tis approximation on numerical examples. IV. NUMERICAL EXPERIMENTS A. Euler equations For te divergence-free Euler system, we look for solutions in te space, span{φ k (x) := φ k (x)φ k2 (y)} k=...n x,k 2=...N y, were ψ k (z) = sin(k πz/l), L = 2r is te lengt of eac side of te spatial domain, and N x and N y are te numbers of basis functions in eac of te spatial, (a) at t =.5 (b) at t =.5 (c) at t =. (d) at t =. (e) at t = 3. (f) at t = 3. Fig. : s and trut for different times dimensions. Te coice of basis satisfies te boundary conditions, ω = on Ω. Te corresponding Fourier-Galerkin model (4) is subject to te perturbation vector f wic as only two non-zero entries, meaning tat in effect, we perturb only two modes: one low, and one ig frequency mode. For te filter, we first generate observations y by running (4) te Euler system forward in timewit initial condition represented by a translation and scaling of Gaussian distributions. Te filter is initialized to zero, meaning we assume no knowledge of te initial condition. Observations are taken on an evenly spaced 5 5 grid away from te boundary, were te state is fixed at zero. For bot te generation of observations and te filtering, we coose N x = N y = 6. Figure sows te estimate and trut at different times. Te sparsity of te observations is apparent from Figure a,
6 u u.4.2 n=, m=.4.2 n=2, m=2.2 time =. Spectral coefficient Spectral coefficient Observations Perturbed observations.2 Time Time.6 (a) Projection coefficients for k = k 2 = (b) Projection coefficients for k = k 2 = n=3, m=3.4.2 Spectral coefficient Relative L 2 norm of estimation error x Fig. 3: True traffic state versus estimate at t= Time Time time = 3. (c) Projection coefficients for(d) Relative L 2 -error of te k = k 2 = 3 estimate.2 Observations Perturbed observations.8 wic is from an early point in te estimation. In Figure c, we see tat by t =., te estimate does still not fully mimic te trut, but tat te flow is being captured. In Figure e, te estimate appears to capture te trut quite well. Figures 2a - 2c sow te estimated and true projection coefficients for some of te low wave numbers. B. Ligtill-Witam-Ricards (LWR) model In tis section we apply te above metod to scalar conservation laws. Recall from [9] tat te standard equilibrium traffic-flow model, LWR model consists of a scalar conservation law, t u(x, t) + x f(u(x, t)) =, (2) wit initial data u (x) = u(x, ) were u : R R + R is te traffic density, x R and t R + are te independent variables, space and time respectively, and f : R R is te flux function. ) A typical flux function is f(u) = uv m ( u u m were te constants, V m and u m, are te maximum speed and te maximum density respectively. We impose periodic boundary conditions on te interval (, 2π). Unlike te Euler equation presented above, x Fig. 4: True traffic state versus estimate at t=3 tis model develops sock discontinuities even subject to periodic boundary conditions and smoot initial condition. Te filter as been applied to te Fourier-Galerkin model 2 wic as been used to generate observations. Te filter as been initialized to. Te estimation results are presented on te figures 3-5. We refer te reader to [9] for te furter details. C. StVenant equations Tis final example sows tat te proposed state estimator may andle systems of conservation laws 2 Te model ad artificial viscosity term activated on iger order modes to allow for correct sock tracking
7 u time = 6..2 Observations Perturbed observations x Fig. 7:, u and te estimate after 2 time steps Fig. 5: True traffic state versus estimate at t=6 Fig. 8:, u and te estimate after 2 time steps Fig. 6: Initial conditions for, u and te estimate wit non-periodic boundary conditions. Te standard equilibrium flood model consists of a system of scalar conservation laws: t + x (u) =, t (u) + x (u 2 + g2 2 ) = (3) wit boundary conditions u(, t) = u l (t) and (, t) = l (t) on (, ), were is te fluid dept, u is te averaged velocity and g is te gravitational constant. Te Discountinuos Galerkin metod as been applied to te above equation to generate observations. Tere was no perturbation, so f =. Te initial conditions for te model and filter are presented at Figure 6. Te estimation results are presented on te figures 7-8. We refer te reader to [] for te furter details. V. CONCLUSION Te paper presents te discrete and continuous versions of te minimax state estimator for Euler equations. Te estimator is derived for a Fourier- Galerkin model wic approximates smoot solutions of te Euler equation wit spectral convergence rate. A very curious topic for te future researc is to develop te idea of [9] for weak solutions of Euler equations, tat is to combine te presented state estimation approac and vanising viscosity metod [] to design a state estimator for L solutions of Euler equations in vorticitystream formulation. We stress tat te presented metod and convergence results apply to Navier- Stokes equations in dimension 2 witout major modifications.
8 REFERENCES [] C. Bardos and E. Tadmor. Stability and spectral convergence of fourier metod for nonlinear problems: on te sortcomings of te 2/3 de-aliasing metod. Numerisce Matematik, 29(4), 25. [2] C. Bardos and E. Titi. Euler equations of incompressible ideal fluids. Uspeki Matematiceskik Nauk, 62(3):5 46, 27. [3] Claude Bardos and Olivier Pironneau. Data assimilation for conservation laws. Metods and Applications of Analysis, 2(2):3 34, [4] S. Blandin, A. Couque, A. Bayen, and D. Work. On sequential data assimilation for scalar macroscopic traffic flow models. Pysica D, 22. [5] M. Demetriou and H. Banks. Adaptive parameter estimation of yperbolic distributed parameter systems: Nonsymmetric damping and slowly time varying systems. ESAIM: Control, Optimisation and Calculus of Variations, 3:33 62, 998. [6] J. Frank and S. Zuk. Symplectic möbius integrators for lq optimal control problems. In Proc. IEEE Conference on Decision and Control, 24. [7] A. Majda and A. Bertozzi. Vorticity and incompressible flow. Cambridge Univ. Press, 22. [8] V. Mallet and S. Zuk. Reduced minimax filtering by means of differential-algebraic equations. In 5t Int. Conf. on Pysics and Control, 2. lib.pyscon.ru. [9] T. Tcrakian and S. Zuk. A macroscopic traffic data assimilation framework based on Fourier-Galerkin metod and minimax estimation. IEEE Transactions on Intelligent Transportation Systems, (99): 3, 24. special issue. [] R. Temam. Navier-Stokes equations: Teory and Numerical Analysis. AMS Celsea Publising, 2. [] S. Tirupati, S. Zuk, T. Tcrakian, and S. McKenna. Data assimilation for d sallow water equations: a minimax approac. in preparation, 25. [2] V. Yudovic. Non-stationary flow of an incompressible liquid. Z. Vycisl. Mat. Mat. Fiz., pages 32 66, 963. [3] S. Zuk. Kalman duality principle for a class of illposed minimax control problems wit linear differentialalgebraic constraints. Applied Matematics and Optimisation, 68(2):289 39, 23. [4] S. Zuk. Minimax projection metod for linear evolution equations. In Proc. IEEE Conference on Decision and Control, ieeexplorer.ieee.org, 23. [5] S. Zuk, J. Frank, I. Herlin, and R. Sorten. Data assimilation for linear parabolic equations: minimax projection metod. SIAM J. Sci. Comp., 25. to appear.
LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationarxiv: v3 [math.na] 15 Dec 2009
THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING M.A. EBRAHIMI, MICHAEL HOLST, AND EVELYN LUNASIN arxiv:91.4548v3 [mat.na] 15 Dec 9 ABSTRACT. In tis paper we investigate te use of te D Navier-Stokes-Voigt
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationSimulations of the turbulent channel flow at Re τ = 180 with projection-based finite element variational multiscale methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Met. Fluids 7; 55:47 49 Publised online 4 Marc 7 in Wiley InterScience (www.interscience.wiley.com). DOI:./fld.46 Simulations of te
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationMODIFIED DIFFERENTIAL EQUATIONS. Dedicated to Prof. Michel Crouzeix
ESAIM: PROCEEDINGS, September 2007, Vol.21, 16-20 Gabriel Caloz & Monique Dauge, Editors MODIFIED DIFFERENTIAL EQUATIONS Pilippe Cartier 1, Ernst Hairer 2 and Gilles Vilmart 1,2 Dedicated to Prof. Micel
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More informationMinimal surfaces of revolution
5 April 013 Minimal surfaces of revolution Maggie Miller 1 Introduction In tis paper, we will prove tat all non-planar minimal surfaces of revolution can be generated by functions of te form f = 1 C cos(cx),
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More information, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim
Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationRecent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm
Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
More informationDecay of solutions of wave equations with memory
Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationUnconditional long-time stability of a velocity-vorticity method for the 2D Navier-Stokes equations
Numerical Analysis and Scientific Computing Preprint Seria Unconditional long-time stability of a velocity-vorticity metod for te D Navier-Stokes equations T. Heister M.A. Olsanskii L.G. Rebolz Preprint
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationAdvancements In Finite Element Methods For Newtonian And Non-Newtonian Flows
Clemson University TigerPrints All Dissertations Dissertations 8-3 Advancements In Finite Element Metods For Newtonian And Non-Newtonian Flows Keit Galvin Clemson University, kjgalvi@clemson.edu Follow
More informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationHigh-Order Energy and Linear Momentum Conserving Methods for the Klein-Gordon Equation
matematics Article Hig-Order Energy and Linear Momentum Conserving Metods for te Klein-Gordon Equation He Yang Department of Matematics, Augusta University, Augusta, GA 39, USA; yang@augusta.edu; Tel.:
More informationA Demonstration of the Advantage of Asymptotic Preserving Schemes over Standard Finite Volume Schemes
A Demonstration of te Advantage of Asymptotic Preserving Scemes over Standard Finite Volume Scemes Jocen Scütz Berict Nr. 366 Juni 213 Key words: conservation laws, asymptotic metods, finite volume metods,
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationSOLUTIONS OF FOURTH-ORDER PARTIAL DIFFERENTIAL EQUATIONS IN A NOISE REMOVAL MODEL
Electronic Journal of Differential Equations, Vol. 7(7, No., pp.. ISSN: 7-669. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu (login: ftp SOLUTIONS OF FOURTH-ORDER PARTIAL
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. MATH. ANAL. ol. 37, No. 3, pp. 737 751 c 2005 Society for Industrial and Applied Matematics DIFFEOMORPHISMS AND NONLINEAR HEAT FLOWS L. C. EANS, O. SAIN, AND W. GANGBO Abstract. We sow tat te gradient
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More informationA Polynomial Adaptive LCP Scheme for Viscous Compressible Flows
A Polynomial Adaptive LCP Sceme for Viscous Compressible Flows J.S. Cagnone, B.C. Vermeire, and S.. Nadaraja Department of Mecanical Engineering, McGill University, Montreal, Canada, H3A 2S6 Email: jean-sebastien.cagnone@mail.mcgill.ca
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationSmoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations
396 Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 3, 396 412 ttp://dx.doi.org/10.15388/na.2014.3.6 Smootness of solutions wit respect to multi-strip integral boundary conditions for nt
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume XX Number 0 ages 45 c 03 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND LORNTZ
More informationON A PROBLEM OF SHALLOW WATER TYPE
24-Fez conference on Differential Equations and Mecanics Electronic Journal of Differential Equations, Conference 11, 24, pp. 19 116. ISSN: 172-6691. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu
More information2009 Elsevier Science. Reprinted with permission from Elsevier.
P. Hansbo and M. Juntunen. 009. Weakly imposed Diriclet boundary conditions for te Brinkman model of porous media flow. Applied Numerical Matematics, volume 59, number 6, pages 174 189. doi:10.1016/j.apnum.008.07.003.
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationEXTENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS WITH APPLICATION TO AN AEROACOUSTIC PROBLEM
SIAM J. SCI. COMPUT. Vol. 26, No. 3, pp. 821 843 c 2005 Society for Industrial and Applied Matematics ETENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationResearch Article New Results on Multiple Solutions for Nth-Order Fuzzy Differential Equations under Generalized Differentiability
Hindawi Publising Corporation Boundary Value Problems Volume 009, Article ID 395714, 13 pages doi:10.1155/009/395714 Researc Article New Results on Multiple Solutions for Nt-Order Fuzzy Differential Equations
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationarxiv: v2 [math.na] 11 Dec 2016
Noname manuscript No. will be inserted by te editor Sallow water equations: Split-form, entropy stable, well-balanced, and positivity preserving numerical metods Hendrik Ranoca arxiv:609.009v [mat.na]
More informationAnalysis of time-dependent Navier-Stokes flow coupled with Darcy
Analysis of time-dependent Navier-Stokes flow coupled wit Darcy flow Ayçıl Çeşmelioğlu and Béatrice Rivière Abstract Tis paper formulates and analyzes a weak solution to te coupling of time-dependent Navier-Stokes
More informationEntropy and the numerical integration of conservation laws
Pysics Procedia Pysics Procedia 00 2011) 1 28 Entropy and te numerical integration of conservation laws Gabriella Puppo Dipartimento di Matematica, Politecnico di Torino Italy) Matteo Semplice Dipartimento
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationA Hybridizable Discontinuous Galerkin Method for the Compressible Euler and Navier-Stokes Equations
A Hybridizable Discontinuous Galerkin Metod for te Compressible Euler and Navier-Stokes Equations J. Peraire and N. C. Nguyen Massacusetts Institute of Tecnology, Cambridge, MA 02139, USA B. Cockburn University
More informationThe Convergence of a Central-Difference Discretization of Rudin-Osher-Fatemi Model for Image Denoising
Te Convergence of a Central-Difference Discretization of Rudin-Oser-Fatemi Model for Image Denoising Ming-Jun Lai 1, Bradley Lucier 2, and Jingyue Wang 3 1 University of Georgia, Atens GA 30602, USA mjlai@mat.uga.edu
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationGELFAND S PROOF OF WIENER S THEOREM
GELFAND S PROOF OF WIENER S THEOREM S. H. KULKARNI 1. Introduction Te following teorem was proved by te famous matematician Norbert Wiener. Wiener s proof can be found in is book [5]. Teorem 1.1. (Wiener
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationTHE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718XX0000-0 THE DISCRETE PLATEAU PROBLEM: CONVERGENCE RESULTS GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We solve te problem of
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More informationLarge eddy simulation of turbulent flow downstream of a backward-facing step
Available online at www.sciencedirect.com Procedia Engineering 31 (01) 16 International Conference on Advances in Computational Modeling and Simulation Large eddy simulation of turbulent flow downstream
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
More information