A posteriori error analysis of multiscale operator decomposition methods for multiphysics models
|
|
- MargaretMargaret Todd
- 5 years ago
- Views:
Transcription
1 A posteriori error analysis of multiscale operator decomposition metods for multipysics models D. Estep, V. Carey 2, V. Ginting 3, S. Tavener 2, T. Wildey 4 5 Department of Matematics and Department of Statistics, Colorado State University, Fort Collins, CO Department of Matematics, Colorado State University, Fort Collins, CO Department of Matematics, University of Wyoming, Laramie, WY Te Institute for Computational Engineering and Sciences, Te University of Texas, Austin, TX estep@mat.colostate.edu Abstract. Multipysics, multiscale models present significant callenges in terms of computing accurate solutions and for estimating te error in information computed from numerical solutions. In tis paper, we describe recent advances in extending te tecniques of a posteriori error analysis to multiscale operator decomposition solution metods. Wile te particulars of te analysis vary considerably wit te problem, tere are several key ideas underlaying a general approac to treat operator decomposition multiscale metods tat is underdevelopment. We explain tese ideas in te context of tree specific examples.. Multiscale, multipysics problems Multipysics, multiscale models tat couple different pysical processes interacting across a wide range of scales abound in te application areas of te SCIDAC program. Suc models present significant callenges for computing accurate numerical solutions, wic makes it very important to obtain accurate estimates of error in computed information. But, multipysics, multiscale problems also raise significant callenges for error estimation. In tis paper, we describe recent advances in computing a posteriori error estimates for multiscale operator decomposition, wic is a powerful solution metod for multipysics, multiscale problems tat is in widespread use in ig performance computing... Tree examples of multipysics models We briefly describe tree examples of multipysics penomena tat arise in SCIDAC application domains. Tese examples illustrate te myriad ways in wic pysical processes are coupled. Example.. A termal actuator. Tis is a MEMS (microelectronic mecanical switc) device modeled by a system of tree coupled equations, eac representing a distinct pysical 5 Tis work was supported in part by DOE (DE-FG2-4ER2562, DE-FG2-5ER25699, DE-FC2-7ER5499), NASA (NNG4GH63G), NSF (DMS-7832, DMS-7535, DGE , MSPA-CSE , ECCS- 7559), INL (69249), LLNL (B57339) and SNL (PO299784)
2 process acting on its own scale, (σ u ) =, x Ω, (κ(u 2 ) u 2 ) = σ( u u ), x Ω, (λ tr(e)i + 2µE β(u 2 u 2,ref )I ) =, E = ( () u 3 + u ) 3 /2, x Ω Te first is an electrostatic current equation governing potential u (current J = σ u ), te second is a steady-state energy equation for te governing temperature u 2, and te tird equation is linear elasticity describing te steady-state displacement u 3. Tis is an example of parameter passing, in wic te solution of one component is used to compute te parameters and/or data for anoter component. Note tat u can be calculated independently, wile computing u 2 requires u. Computing u 3 requires u 2 and terefore u. Example.2. Te Brusselator problem. Tis is a model of cemical dynamics [], u t k 2 u = α (β + )u x 2 + u 2 u 2, x (, ), t >, u 2 t k 2 u 2 2 = βu x 2 u 2 u 2, x (, ), t >, u (, t) = u (, t) = α, u 2 (, t) = u 2 (, t) = β/α, t >, u (x, ) = u, (x), u 2 (x, ) = u 2, (x), x (, ), (2) were u and u 2 are te concentrations of species and 2, respectively. Tis problem combines different pysics - in tis case, reaction and diffusion - in one equation. Te generic picture for a reaction-diffusion equation is a relatively fast, destabilizing reaction component interacting wit a relatively slow, stabilizing diffusion component. Tus, te pysical components ave bot different scales and different stability properties. Example.3. Conjugate eat transfer between a fluid and solid object. We consider te flow of a eat-conducting Newtonian fluid past a solid cylinder. Te model consists of te eat equation in te solid and te equations governing te conservation of momentum, mass and energy in te fluid, were we apply te Boussinesq approximation to te momentum equations in te fluid. Te temperature field is advected by te fluid and couples back to te momentum equations troug te buoyancy term. Let Ω F and Ω S be polygonal domains in R 2 wit boundaries Ω F and Ω S intersecting along an interface Γ I = Ω S Ω F. Te complete coupled problem is µ u + ρ (u ) u + p + ρ βt F g = ρ ( + βt ) g, x Ω F, u =, x Ω F, k F T F + ρ c p (u T F ) = Q F, x Ω F, (3) T S = T F, interface x Γ I, k F (n T F ) = k S (n T S ), k S T S = Q S, x Ω S, were we supplement tese equations by imposing Diriclet and Neumann conditions for te velocity field respectively on boundaries Γ u,d and Γ u,n and Diriclet and Neumann conditions for te temperature fields in te fluid and te solid respectively on boundaries Γ TF,D, Γ TF,N, Γ TS,D, and Γ TS,N. Here, ρ and T are reference values for te density and temperature respectively, µ is te molecular viscosity, β is te coefficient of termal expansion, c p is te specific eat, k F and k S are te termal conductivities of te fluid and solid respectively, Q F and Q S are source terms and n is te unit normal vector directed into te fluid. Note tat u is a vector. Tis presents a class of problems were pysics in different pysical domains are coupled across a common boundary.
3 .2. Callenges and goals for multiscale, multipysics models Callenges arising in te numerical solution of multiscale, multipysics problems include producing accurate numerical solutions efficiently, dealing wit complex stability properties, and transforming information across scales. Anoter complication is te range of applications of multipysics models, including model prediction, sensitivity analysis, and parameter optimization, tat require computation of solutions corresponding to wide range of data and parameters. We expect te solution beavior to vary significantly and te ability to obtain accurate numerical solutions terefore to vary as well. On te oter and, predictive science and engineering ave come to rely on ig performance simulation of complex pysical penomena. In tis context, it is critically important to accurately quantify te numerical error in any computed information. Anoter important goal is determining efficient ways to compute specific information accurately. 2. Multiscale operator decomposition Multiscale operator decomposition is a widely used tecnique for solving multipysics, multiscale models in SCIDAC simulations. Te general approac is to decompose te multipysics and/or multiscale problem into components involving simpler pysics over a relatively limited range of scales, and ten to seek te solution of te entire system troug some sort of iterative procedure involving numerical solutions of te individual components. We illustrate in Fig.. In general, Fully Coupled Multiscale Operator Decomposition Solution Pysics Pysics Pysics 2 Pysics 3 Pys upscale downscale upscale downscale upscale downscale upscale downscale Pysics 2 Pysics 2 Pysics 2 2 Pysics 2 3 Figure. Left: Illustration of a multiscale, multipysics model. Rigt: Illustration of a multiscale operator decomposition solution. different components are solved wit different numerical metods as well as wit different scale discretizations. Reasons tat tis approac is appealing include; Tere is generally a good understanding of ow to solve a broad spectrum of single pysics problems accurately and efficiently. It provides an alternative to accommodating multiple scales in one discretization. It allows for efficient ig performance simulations. However, multiscale operator decomposition generally as significant impact on te accuracy and stability of numerical solutions. 3. Analysis of multiscale operator decomposition In general, multiscale, multipysics problems ave complex stability properties tat defy accurate description via classic a priori analysis. We use duality and adjoint operators to quantify stability a posteriori. We focus on computing a particular quantity of interest determined by a linear functional. We combine tese tools wit variational analysis to produce accurate estimates of te effects of perturbation and error. By numerically solving te adjoint problem, we can compute very accurate error estimates.
4 Duality and adjoint operators ave a long istory of application in model sensitivity analysis and optimization, dating back to Lagrange, see [2, 3, 4, 5] for example. Te application of tese tools to a posteriori error estimation as a more recent istory, see [6, 7, 8, 9,,, 2, 3]. Our main purpose is to describe recent advances in extending te tecniques of a posteriori error analysis to multiscale operator decomposition solutions of multipysics, multiscale problems. Wile te particulars of te analysis vary considerably wit te problem, tere are several key ideas underlaying a general approac to treat operator decomposition multiscale metods tat is under development, including: We identify auxiliary quantities of interest associated wit information passed between pysical components and solve auxiliary adjoint problems to estimate te error in tose quantities. We deal wit scale differences by introducing projections between discrete spaces used for component solutions and estimate te effects of tose projections. Te standard linearization argument used to define an adjoint operator associated wit error analysis for a nonlinear problem may fail, requiring a new approac. Te adjoint operator associated wit a multiscale operator decomposition solution metod is often significantly different tan te adjoint associated wit te original problem, wic affects te stability of te metod. In practice, solving te adjoint associated wit te original fully-coupled problem may present te same kinds of computational callenges posed by te original problem, so attention must be paid to te solution of te adjoint problem. We explain tese ideas in te context of tree examples. 3.. Multiscale decomposition of triangular systems of elliptic problems We capture te essential features of te termal actuator model in Example. using a oneway coupled system a u + b u + c u = f (x), x Ω, a 2 u 2 + b 2 u 2 + c 2 u 2 = f 2 (x, u, Du ), x Ω, (4) u = u 2 =, x Ω, were a i, b i, c i, f i are smoot functions, wit a, a 2 α > on a bounded domain Ω in R N wit boundary Ω, and α is a constant. We introduce te finite element space S, (Ω) H (Ω), corresponding to a discretization T, of Ω for te first component, and anoter finite element space S,2 (Ω), on a different mes T,2, for te second component. To evaluate integrals involving functions defined on different meses, we introduce projections Π i j from S,i to S,j, e.g. interpolants or an L 2 ortogonal projection. We also define te bilinear forms, A i (u i, v i ) (a i u i, v i )+(b i (x) u i, v i )+(c i u i, v i ), i =, 2. We present te multiscale operator decomposition algoritm in Alg.. Algoritm Multiscale Operator Decomposition for Triangular Systems of Elliptic Equations Construct discretizations T,, T,2 and corresponding spaces S,, S,2 Compute U S, (Ω) satisfying A (U, v ) = (f, v ), for all v S, (Ω). Compute U 2 S,2 (Ω) satisfying A 2 (U 2, v 2 ) = (f 2 (x, Π 2 U, Π 2 DU ), v 2 ), for all v 2 S,2 (Ω). We observe tat any errors made in te solution of te first component affect te solution of te second component. Tis turns out to be a crucial fact for a posteriori error analysis.
5 Example 3.. We solve a system u = sin(4πx) sin(πy), x Ω u 2 = b u =, x Ω, u = u 2 =, x Ω, b = 2 π ( ) 25 sin(4πx), Ω = [, ] [, ], (5) sin(πx) using a standard piecewise linear, continuous finite element metod to compute te quantity of interest u 2 (.25,.25). We solve for u first and ten solve for u 2 using independent meses. Using uniform meses, evaluating te standard a posteriori error estimate for te second component problem wile ignoring any effect arising from error in te solution of te first component, we estimate te error in te quantity of interest to be.42. However, te true error is.48, so tere is discrepancy of.6 ( 3%) in te estimate. Tis is a consequence of ignoring te error transferred from te solution of te first component. If we adapt te mes for te solution of te second component based on te a posteriori error estimate of te error in tat component wile neglecting te effects of te decomposition, te discrepancy becomes alarmingly worse. For example, we can refine te mes until te estimate of te error in te second component is.. But, we find tat te true error is.2244! See Fig Figure 2. Multiscale operator decomposition solutions (u on te left, u 2 in te middle) of te component problems of (5). We solve for u 2 using a mes adapted to control te pointwise error of te numerical solution of u 2 wile ignoring te effect of errors in te approximation of u. On te rigt, we plot te pointwise error of te computed u Description of te a posteriori analysis We seek te error in a quantity of interest given by a functional of e 2 = U 2 u 2. Since we introduce some additional, auxiliary quantities of interest, we denote te primary quantity of interest by (ψ, e) = ( ψ () 2, e 2). We use te adjoint operators A i (φ i, v i ) = (a i φ i, v i ) (div(b i φ i ), v i ) + (c i φ i, v i ), i =, 2. We also use te linearization Lf 2 (u )(u U ) = f 2 u (u s + U ( s)) ds. Noting tat te solution of te first adjoint component is not needed to compute te quantity of interest, we define te primary adjoint problem to be A 2(φ () 2, v 2) = (ψ () 2, v 2), for all v 2 W 2 (Ω). Te error representation formula for te transfer error is ( Df2 (U ) e, Π 2 φ () 2 ) ( + Df2 (U ) e, (I Π 2 )φ () ) 2, (6)
6 were Π i j is a projection from discretization i to discretization j. Te first term in (6) is te error contribution arising from te transfer wile te second term is significant wen te approximation spaces are different. Te transfer error is a (nominally nonlinear) functional of te error in u, defining an auxiliary quantity of interest. We approximate it by a linear functional and define te auxiliary quantity of interest, ( f2 (u ) f 2 (U ), φ () ) ( 2 Df2 (U ) e, Π 2 φ () ) ( (2) 2 = ψ, e ). We define te corresponding transfer error adjoint problem A (φ (2), v ) = ( ψ (2), v ) for all v W 2 (Ω). (7) Te additional term ( Df 2 (U ) e, (I Π 2 )φ () ) 2 is a linear functional, so we define anoter auxiliary quantity of interest ( (3) ψ, e ) ( = Df2 (U ) e, (I Π 2 )φ () ) 2 and te corresponding adjoint problem Finally, we define te weak residuals A (φ (3), v ) = ( ψ (3), v ) for all v W 2 (Ω). (8) R (U, χ; ν) = (f (ν), χ) A (U, χ), R 2 (U 2, χ; ν) = (f 2 (Π 2 ν), χ) A 2 (U 2, χ). Te final error representation is terefore [4] Teorem 3.2 (ψ, e) = R 2 (U 2, (I Π 2 )φ () 2 ; U ) + R (U, (I Π )(φ (2) + φ (3) )) + ( Π 2 f 2 (U ) f 2 (Π 2 U ), φ () ) ( 2 + (I Π 2 )f 2 (U ), φ () ) (9) 2. We empasize tat evaluating te integrals in (9) is far from trivial. We ave used Monte- Carlo tecniques wit good results, see [4]. Example 3.3. In Example 3., we estimate te contributions to te error reported in tat computation using te relevant portions of (9). To produce te adaptive mes results sown in Fig. 2, we construct te adapted mes using equidistribution based on a bound derived from te first term in (9), i.e. neglecting te terms tat estimate te transfer error. Instead, we consider te system (5) for te quantity of interest equal to te average value of U 2. We begin wit te same initial coarse meses as in Fig. 2, but add te transfer error expression to te mes refinement criterion. Adapting te mes so tat te total error in te quantity of interest for U 2 as error estimates less tan 4 yields te meses sown in Fig. 3. We see tat te first component solve requires significantly more refinement tan te second component Multiscale decomposition of reaction-diffusion problems Operator splitting for reaction-diffusion problems is te classic example of multiscale operator decomposition. We consider u = ɛ u + f(u), x Ω, < t, t suitable boundary conditions, x Ω, < t, () u(, ) = u ( )
7 Figure 3. Te adapted meses resulting from te full estimate tat accounts for primary and transfer errors. Te transfer error dominates and te mes for te first component is refined wile te mes for te second component is not. Noneteless, te simulation yields an accurate value for te quantity of interest. It is interesting to compare tis result wit tat sown in Fig. 2. Since te quantity of interest is computed from te second component, intuition migt suggest tat te mes for te second component needs refinement. But, it is te mes for te first component tat is critical for determining accuracy. were Ω R d is a spatial domain and f is a smoot function. Employing te metod of lines, we discretize in space using a continuous, piecewise linear finite element metod wit M elements and obtain te initial value problem: find y R M suc tat ẏ = Ay(t) + F (y(t)), < t T, () y() = y, were A is an l l constant matrix representing a diffusion component and F (y) = (F (y), F 2 (y),, F l (y)) is a vector of nonlinear functions representing a reaction component. We next discretize in time using a discontinuous (dg) Galerkin metod [5, 8]. Wit appropriate coices of quadrature, we can treat standard integration metods. We first discretize [, T ] into = t < t < t 2 < < t N = T wit diffusion time steps t n N n=, t n = t n t n. For eac diffusion step, we coose a (small) time step s n = t n /M n wit s = max n N s n, and te nodes t n = s,n < s,n < < s Mn,n = t n. We define te time intervals I n = [t n, t n ] and I m,n = [s m,n, s m,n ]. See Fig. 4. Te finite element approximate solutions are sougt in Diffusion Integration: Reaction Integration: t t t 2 t 3 t 4 t t t 5 2 t t t s... s, s M,... s,2 s M2,2 s 2 s 3... s,3 s M3,3... s,4 s M4,4 Figure 4. Discretization of time used for multiscale operator splitting piecewise polynomial spaces for te diffusion and reaction components respectively, V (q d) = U : U In P (q d) (I n ), V (q r) (I n ) = U : U Im,n P (q r) (I m,n ), for n =,, N and, for eac n, m M n. P (q) (I) denotes te space of polynomials in R l of degree q on interval I. We let U n +, denote te left- and rigt-and limits of U at t n and [U] n = U n + Un te jump value of U at t n.
8 wit Let Ỹ (t) be te piecewise continuous finite element approximation of te operator splitting Ỹ (t) = t n t t n Ỹ n + t t n t n Ỹ n, t n t t n. Te nodal values Ỹn are obtained from te following procedure: Algoritm 2 Multiscale Operator Splitting for Reaction-Diffusion Equations Set Ỹ = y for n =,, N do Set Y r,n = Ỹn for m =,, M n do Compute Y r Im,n P (q r) (I m,n ) satisfying I m,n ( ) Y, W dt + ( [Y r ] m,n, W m + ) = (F (Y r ), W ) dt W P (qr) (I m,n ) (2) I m,n end for Set Y d n = Y r M n,n, compute Y d In P (q d) (I n ) satisfying Set ỹ n = y d (t n ) end for ( ( ) Y d, V ) dt + [Y d ] n, V n + = (AY d, V ) dt V P (qd) (I n ) (3) I n I n Example 3.4. We illustrate te instability of operator splitting applied to te Brusselator problem (2). We apply a standard first order splitting sceme to a space discretization of te Operator Split Solution Spatial Location L 2 norm of error t = 6.4 t = 6 t = 32-4 t = 64 t = slope Time Step Size Figure 5. Left: Typical instability arising from multiscale operator splitting applied to Brusselator problem. Solution is sown at time 8. Rigt: Plots of te error in te L 2 norm versus time step size at different times. Brusselator model wit 5 discrete points wit α =.6, β = 2, k = k 2 =.25 consisting of te trapezoidal sceme for te diffusion wit time step of.2 and backward Euler sceme for te reaction wit time step of.4. On te left of Fig. 5, we sow a numerical solution tat exibits nonpysical oscillations tat developed after some time. On te rigt, we sow plots of te error versus time steps at different times. Tere is a critical time step above wic te instability develops. Moreover, canging te space discretization does not improve te accuracy. In [6], it is demonstrated tat a finer spatial discretization for a constant time step size leads to significantly more error in te long time solution.
9 3.2.. Description of te ybrid a posteriori-a priori error analysis An accurate a posteriori estimate of te error must account for te stability effects arising from operator splitting. Te standard approac used to define an adjoint operator for analyzing nonlinear problems based on linearization of a perturbation equation fails. It turns out tat te adjoint operators associated wit te original problem and te multiscale operator decomposition metod are fundamentally different. In te estimate below, tis difference takes te form of residuals between certain adjoint operators. A practical difficulty wit suc a result is tat solving te adjoint for te fully coupled problem poses te same multipysics callenges as solving te original forward problem. We terefore develop a new ybrid a priori - a posteriori estimate tat combines a computable leading order expression obtained using a posteriori arguments wit a provably iger order bound obtained using a priori convergence result. Note tat Ỹ is a consistent numerical solution for an analytic version of te operator splitting and te expression for its error can be estimated using te standard a posteriori error analysis. We let ϑ d define te adjoint solution associated wit te diffusion component (3) and ϑ r define te adjoint solution associated wit te reaction component (2) satisfying respectively ϑ d = A ϑ d (t), t n > t t n, ϑ d (t n ) = ψ n, n = N,,, ϑ r = ( ˆF (y r, Y r )) ϑ r (t), s m,n > t s m,n, ϑ r (s m,n ) = ψm,n, r n = N,,, m = M n,,, wit ψ r M n,n = ϑ d+ n and ψ r m,n = ϑ r m,n for m < M n. Tus ϑ r is continuous across te internal reaction time nodes s m,n, m =,, M n. Here, ˆF (y r, Y r ) = F (sy r + ( s)y r ) ds. Te nonlinearity complicates te analysis because te coice of trajectory around wic to linearize is not clear. We cannot use te standard approac of linearizing te error representation because of te operator splitting. Instead, we assume tat bot te original problem and te operator split version ave a common solution and we linearize eac problem in a neigborood of tis common solution. For example, we assume tat y = is a steady state solution of bot problems, wic can be acieved by assuming tat Homogeneity Assumption: F () =, and we linearize in a region around. In terms of applications to reaction-diffusion problems, tere are matematical reasons for tis assumption and it is satisfied in a great many applications. We can modify te analysis to allow for linearization around any known common solution, see [7]. For n =,, N, we define te tree adjoint problems. Te diffusion problem is simplest because it is linear, ϕ d = A ϕ d (t), t n > t t n, ϕ d (t n ) = ψn. d (4) It is convenient to let Φ d n denote te solution operator, so ϕ d (t n ) = Φ d nψ d n. We require two adjoint problems to treat te reaction component. Te difference between te problems is te function around wic tey are linearized, ϕ r = F (Ỹ ) ϕ r (t), t n > t t n, ϕ r (t n ) = ψ r n. ϕ r 2 = F (Y r ) ϕ r 2 (t), t n > t t n, ϕ r 2 (t n ) = ψ r n. If Φ r n(z) denotes te solution operator for te problem linearized around a function z, ten we ave ϕ r (t n ) = Φ r n(ỹ )ψr n and ϕ r 2 (t n ) = Φ r n(y r )ψ r n. We can now prove [7]. (5)
10 Teorem 3.5 A ybrid a posteriori - a priori error estimate (ỸN y N, ψ N ) = were + E = 2 t n N n= M n n= m= + N ( ) ( Y d AY d, ϑ d Πϑ d ) dt + ([Y d ] n, ϑ d+ n Πϑ d+ n ) I n ( ) ( Y r F (Y r ), ϑ r Πϑ r ) dt + ([Y r ] m,n, ϑ r+ m,n Πϑr+ m,n ) I m,n N (Ỹn, (E + E 2 )ψ n ) + O( t qd+2 ) + O( t s qr+ ), n= (A F(Ỹ ) F(Ỹ )A ), F(Ỹ ) = I n F (Ỹ ) dt, E 2 = ( ) Φ r n(ỹ ) Φr n(y r ) Φ d n. Te first expression on te rigt is te error introduced by te numerical solution of te diffusion component wile te second expression is te error introduced by te numerical solution of te reaction component. Te tird expression on te rigt measures te effects of operator splitting, were E is a leading order estimate for te effects of operator splitting wile E 2 accounts for issues arising from te differences in linearizing around te global computed solution as opposed to te solution of te reaction component. Bot of tese quantities are scaled by te solution itself, so tat tese effects become negligible wen te solution approaces zero. Finally, te remaining terms represent bounds on terms tat are not computable but are iger order. In practice, we neglect tose terms wen computing an estimate. Example 3.6. We consider te Brusselator problem (2) wit α = 2, β = 5.45, k =.8, k 2 =.4 and initial conditions u (x, ) = α +. sin(πx) and u 2 (x, ) = β/α +. sin(πx), wic yields an oscillatory solution. In tis case, te reaction is very mildly unstable, wit at most polynomial rate accumulation of perturbations as time passes. We use a 32 node spatial finite element discretization, resulting in an differential equation system wit dimension 62. We note tat in original form, te reaction terms do not satisfy te requirement F () = so we linearize around te steady state solution c wit wit c i = α for i =,, N e and c i = β/α for i = N e,, 2N e 2, so tat F (c) =. Fig. 6 compares te errors computed using t =. and M = reaction time steps to te ybrid a posteriori error estimates. We sow results for [, 2], wen te solution is still in a transient stage, and at T = 4 wen te solution as become periodic. All te results sow tat te exact and estimated errors are in remarkable agreement Multiscale decomposition of a fluid-solid conjugate eat transfer problem We next consider te multiscale decomposition solution of te eat transfer problem described in Example.. Te weak formulation of (3) consists of computing u V F, p L 2 (Ω F ), T F W F and T S W S suc tat a (u, v) + c (u, u, v) + b(v, p) + d(t F, v) = (f, v), b(u, q) =, (6) a 2 (T F, w F ) + c 2 (u, T F, w F ) + a 3 (T S, w S ) = (Q F, w F ) + (Q S, w S ), for all v V F,, q L 2 (Ω F ), w F W F, and w S W S,, were
11 Error and Estimate Species Species Component Error and Estimate Species Species Time Error and Estimate Species Species Component Figure 6. Brusselator results. Left: Comparison of errors against te spatial location at T = 2. Middle: Time istory of errors at te midpoint location on [, 2]. Rigt: Comparison of errors against te spatial location at T = 4. Te dotted line is te exact error and te (+) is te estimated error and a (u, v) = Ω F µ( u : v) dx, a 2 (T F, w F ) = Ω F k F ( T F w F ) dx, a 3 (T S, w S ) = Ω S k S ( T S w S ) dx, b(v, q) = Ω F ( v)q dx, c (u, v, z) = Ω F ρ (u ) v z dx, c 2 (u, T, w) = Ω F ρ c p (u T ) w dx, d(t, v) = Ω F ρ βt g v dx, f = ρ ( + βt ) g, V F = v H (Ω F ) v = g u,d on Γ u,d, VF, = v V F v = on Γ u,d, W F = w H (Ω F ) w = g TF,D on Γ TF,D, WF, = w W F w = on Γ TF,D, W S = w H (Ω S ) w = g TS,D on Γ TS,D, WS, = w W S w = on Γ TS,D. To discretize, we construct independent locally-quasi-uniform triangulations T F, and T S, of Ω F and Ω S respectively. We use te piecewise polynomial spaces V F = v V F v continuous on Ω F, v i P 2 (K) for all K τ F,, Z = z Z z continuous on Ω F, z P (K) for all K τ F,, W F = w W F w continuous on Ω F, w P 2 (K) for all K τ F,, W S = w W S w continuous on Ω S, w P 2 (K) for all K τ S,, and te associated subspaces VF, = v V v = on Γ u,d, WF, = w WF w = on Γ TF,D, WS, = w WS w = on Γ TS,D and w = on Γ I, were P q (K) denotes te space of polynomials of degree q on an element K. Note tat WS, is different tan WF, in an important way since Γ T S,D does not include Γ I. We let π V, π WF, π WS, and π Z be projections into VF, W F, W S and Z respectively. We also use π WF and π WS to denote projections into WF and W S respectively along te interface Γ I. To compute a stable solution of te fluid equations, we coose VF and Z to be te Taylor- Hood finite element pair satisfying te discrete inf-sup condition. In general, te convergence of te iteration defined by tis Algoritm depends on te values of k S and k F along te interface and te geometry of eac region. Often, a relaxation is used to elp improve convergence properties, wic affects te analysis [8, 9].
12 Algoritm 3 Multiscale Decomposition Metod for Conjugate Heat Transfer k = wile ( T k S π S T k F Γ I > T OL) do k = k+ Compute T k S, W S Compute u k V F, pk k suc tat T S, = π W S T k F, along te interface Γ I and a 3 (T k S,, w) = (Q S, w), w W S,, (7) Z and T k F, W F suc tat a (u k, v) + c (u k, uk k, v) + b(v, pk ) + d(t F,, v) = (f, v), v V F,, b(u k, q) =, q Z, (8) a 2 (T k F,, w) + c 2(u k, T k F,, w) = (Q F, w) (k S (n T k S, ), w) Γ I, w WF,. end wile Description of an a posteriori error analysis We define te errors e u = u u k e p = p p k, e T F = T F T k F,, and e T S = T S T k S,. Te adjoint problem for te quantity of interest (ψ, e) = (ψ u, e u ) + (ψ p, e p ) + (ψ TF, e TF ) + (ψ TS, e TS ) is µ φ + c (φ) + z + c 2u (θ F ) = ψ u, x Ω F, φ = ψ p, x Ω F, k F θ F + c 2T (θ F ) + ρ β(g φ) = ψ TF, x Ω F, θ F = θ S, x Γ I, k F (n θ F ) = k S (n θ S ), k S θ S = ψ TS, x Ω S, Here, we ave used te linearizations φ =, x Γ u,d, µ φ n =, x Γ u,n, θ F =, x Γ TF,D, k F (n θ F ) =, x Γ TF,N, θ S =, x Γ TS,D, k S (n θ S ) =, x Γ TS,N. c (φ) = 2 ρ (u + u ) φ 2 ρ (u + u ) φ 2 ρ ( (u + u )) φ, c 2u(θ) = 2 ρ c p (T + T ) θ, c 2T (θ) = 2 ρ c p (u + u ) θ 2 ρ c p ( (u + u )) θ. We solve (9) numerically using an iterative operator decomposition approac as for te forward problem independently of te forward iterations. In [8, 9], we derive estimates tat only require adjoint solutions of te two component problems. To write out te a posteriori error representation, we introduce an additional projection π W S : H 2 W S, defined suc tat for any node x i. π W S θ S (x i ) = π WS θ S (x i ) for x i / Γ I, and oterwise. We can now prove [9],, (9)
13 Teorem 3.7 (ψ, e) = (f, φ π V φ) a (u k, φ π V φ ) c (u k, uk, φ π V φ) b(φ π V φ, p ) d(t k F,, φ π V φ) b(u k, z π Zz) (2) + (Q F, θ F π WF θ F ) a 2 (T k F,, θ F π WF θ F ) c 2 (u k, T k F,, θ F π WF θ F ) + (Q S, θ S π WS θ S ) a 3 (T k S,, θ S π WS θ S ) (2) + ( T k S, π ST k F,, k S(n θ S ) ) + ( π Γ I S T k F, T k F,, k S(n θ S ) ) (22) Γ I + ( k S (n T k S, ), π ) W F θ F + ( ) ( k Q Γ I S, π WS θ S a3 T S,, π ) W S θ S. (23) Te contributions to te error are Equations (2)-(2) represents te contribution of te discretization error arising from eac component solve. Equation (22) represents te contribution from te iteration. Te first term in (23) represents contribution of te transfer error wile te remaining terms represent te contribution arising from projections between two different discretizations. Example 3.8. For te flow past a cylinder, we solve te steady non-dimensionalized Boussinesq equations in te fluid domain and te non-dimensional eat equation in te solid domain. To simulate te flow of water past a cylinder made from stainless steel, we set te dimensionless constants P r = 6.6 and k R = 3, and coose te inflow velocity and te temperature gradient so tat Re = 75, P e = 495, F r =., Ra = 5. Te temperature gradient is imposed by setting different temperatures along te top and bottom boundaries, wit a linear temperature gradient on te inflow boundary, and an adiabatic condition on te outflow boundary. We sow results for two quantities of interest. Te first is te temperature in a small region in te wake, located approximately one cannel widt downstream of te center of te cylinder and /4 of a cannel widt below te upper wall. Te second is temperature at a small region in te center of te cylinder. In eac case, we base adaptivity on an element tolerance of 8. We sow te final adaptive meses for te flow and solid in Fig. 7. For te first quantity of interest, te flow mes is most refined near te region of interest and upstream of te region of interest, locating more elements between te cylinder and te top wall tan te cylinder and te bottom wall since te flow advecting eat to te region of interest passes above rater tan below te cylinder. For te solid, te mes is igly refined along te top in order to increase te accuracy of te normal derivative tat is computed in te solid and used as a boundary condition in te fluid computation. Evidently, te normal derivatives elsewere on te interface ave less of an influence on te first quantity of interest. For te second quantity of interest, te mes is igly refined upstream of te cylinder. We note tat te refinement downstream of te cylinder corresponds closely to te recirculation region, and te mes refinement is sligtly asymmetric about te midplane of te cannel due to te asymmetric initial mes. Te mes in te solid is refined uniformly near te boundary, reflecting te fact tat te error in te finite element flux makes a significant contribution to te error in te quantity of interest Loss of order and flux correction Te meses sown in Fig. 7 are igly refined near te interface. Tis reflects te fact tat tere is significant error in te numerical flux passed between te components. It turns out tat tis pollutes te entire computation, so tat overall te metod loses an entire order of accuracy.
14 Figure 7. Upper: Final adaptive meses in te fluid and solid wen te quantity of interest is te temperature in a small region in te wake above te cylinder. Lower: Final adaptive meses in te fluid and solid wen te quantity of interest is te temperature in a small region in te center of te solid. One way to compensate for te loss of order is by refining te mes locally near te interface. Anoter way is to compute te particular information, in tis case te flux on te interface, more accurately. It turns out tat we can adapt a post-processing tecnique called flux correction developed originally by Weeler [2] and Carey [2, 22] to recover boundary flux values wit increased accuracy. We denote te set of elements in τ S, tat intersect te interface boundary by τ Γ I S, = K τ S, K Γ I, and consider te corresponding finite element space Σ = v P 2 (K) wit K τ Γ I S,, v(η i) = if η i / Γ I, were η i denotes te nodes of element K. Te degrees of freedom correspond to te nodes on te boundary. We compute σ k Σ satisfying ( σ k, v ) Γ I = (Q S, v) a 3 (T k S,, v), for all v Σ, were T k S, solves (7). Since te dimension of te problem scales wit te number of nodes on a boundary, it is relatively inexpensive to solve. In te Algoritm, (8) is replaced by: Compute u k VF, pk Z and T k F, W F suc tat a (u k, v) + c (u k, uk k, v) + b(v, pk ) + d(t F,, v) = (f, v), b(u k, q) =, a 2 (T k F,, w) + c 2(u k, T k F,, w) = (Q F, w) ( σ k, w ), v V Γ I F,, q Z, w WF,. (24) It turns out tat using te recovered boundary flux leads to a cancelation of te transfer error term in te error representation formula (23), wic is te source of te loss of order. We can prove tat using te recovered flux recovers te expected cubic order of convergence, see [9]. Example 3.9. Te recovered accuracy is easily demonstrated by considering te adapted meses produced by te modified algoritm. We sow te final adaptive meses for te solid in Fig. 8. Tere is no mes refinement near te boundaries, indicating tat te flux error is no longer dominant. 4. Conclusion Multipysics, multiscale models present significant callenges in terms of computing accurate solutions and for estimating te error in information computed from numerical solutions. In tis
15 Figure 8. Left: Final adaptive mes in te solid wen te quantity of interest is te temperature in a small region in te wake above te cylinder. Rigt: Final adaptive mes in te solid wen te quantity of interest is te temperature in a small region in te center of te solid. paper, we describe recent advances in extending te tecniques of a posteriori error analysis to multiscale operator decomposition solution metods. Wile te particulars of te analysis vary considerably wit te problem, tere are several key ideas underlaying a general approac to treat operator decomposition multiscale metods tat is under development. We explain tese ideas in te context of tree specific examples. In tis presentation, we ave minimized te effects arising from te solution of nonlinear and/or fully coupled systems by carefully coosing te models and results tat are discussed. Referring back to Fig., we generally expect multiscale operator decomposition to require a number of iterations between te pysical components. Tis raises additional issues tat are discussed in [8, 9, 23, 7, 24]. References [] Prigogine I and Lefever R 968 J. Cem. Pys [2] Marcuk G I, Agoskov V I and Sutyaev V P 996 Adjoint equations and perturbation algoritms in nonlinear problems (Boca Raton, FL: CRC Press) ISBN [3] Lanczos C 997 Linear Differential Operators (Dover Publications) [4] Cacuci D 997 Sensitivity and Uncertainty Analysis: Teory vol I (Capman & Hall/CRC) [5] Marcuk G I 995 Adjoint equations and analysis of complex systems (Kluwer) [6] Eriksson K, Estep D, Hansbo P and Jonson C 996 Computational differential equations (Cambridge: Cambridge University Press) ISBN ; [7] Eriksson K, Estep D, Hansbo P and Jonson C 995 Acta numerica, 995 Acta Numer. (Cambridge: Cambridge Univ. Press) pp 5 58 [8] Estep D, Larson M G and Williams R D 2 Mem. Amer. Mat. Soc. 46 viii+9 ISSN [9] Becker R and Rannacer R 2 Acta Numerica 2 [] Giles M and Süli E 22 Acta Numerica [] Bangert W and Rannacer R 23 Adaptive Finite Element Metods for Differential Equations (Birkauser Verlag) [2] Parascivoiu M, Peraire J and Patera A 997 Comput. Metods Appl. Mec. Engrg [3] Bart T J 24 A-Posteriori Error Estimation and Mes Adaptivity for Finite Volume and Finite Element Metods (Lecture Notes in Computational Science and Engineering vol 4) (New York: Springer) [4] Carey V, Estep D and Tavener S 26 SINUM To appear [5] Estep D and Stuart A M 22 Mat. Comp (electronic) ISSN [6] Ropp D L, Sadid J N and Ober C C 24 J. Comput. Pys ISSN [7] Ginting V, Estep D, Sadid J and Tavener S 26 SINUM To appear [8] Estep D, Tavener S and Wildey T 26 SINUM To appear [9] Estep D, Tavener S and Wildey T 28 JCP Submitted [2] Weeler M 974 Numer. Mat [2] GF Carey SS Cow M S 985 Comp. Met. in Applied Mec. and Engr [22] Carey G 982 Comp. Met. in Applied Mec. and Engr [23] Carey V, Estep D and Tavener S 28 In preparation [24] Ginting V, Estep D and Tavener S 28 IMA J. Numer. Analysis Submitted
Adjoint-Fueled Advances in Error Estimation for Multiscale, Multiphysics Systems
Adjoint-Fueled Advances in Error Estimation for Multiscale, Multiphysics Systems Donald Estep University Interdisciplinary Research Scholar Department of Statistics Department of Mathematics Center for
More informationAdjoint-Fueled Advances in Error Estimation for Multiscale, Multiphysics Systems
Adjoint-Fueled Advances in Error Estimation for Multiscale, Multiphysics Systems Donald Estep University Interdisciplinary Research Scholar Department of Statistics Center for Interdisciplinary Mathematics
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 informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
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 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 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 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 informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
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 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 informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
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 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 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 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 informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
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 informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
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 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 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 informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
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 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 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 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 informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
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 informationLEAST-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 informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
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 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 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 informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
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 informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
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 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 informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
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 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 informationarxiv: v1 [math.na] 27 Mar 2018
Numerical study of goal-oriented error control for stabilized finite element metods M. P. Brucäuser, K. Scwegler, M. Bause Helmut Scmidt University, Faculty of Mecanical Engineering, Holstenofweg 85, 2243
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationarxiv: v1 [physics.flu-dyn] 3 Jun 2015
A Convective-like Energy-Stable Open Boundary Condition for Simulations of Incompressible Flows arxiv:156.132v1 [pysics.flu-dyn] 3 Jun 215 S. Dong Center for Computational & Applied Matematics Department
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationEffect of Numerical Integration on Meshless Methods
Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions
More informationChapters 19 & 20 Heat and the First Law of Thermodynamics
Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,
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 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 informationERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS. 1. Introduction. We consider the parabolic PDE of the form,
ERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS K. CHRYSAFINOS AND NOEL. J. WALKINGTON Abstract. We analyze te classical discontinuous Galerkin metod for a general parabolic
More informationA MESH ADAPTATION METHOD FOR 1D-BOUNDARY LAYER PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 3, Number 4, Pages 408 428 c 2012 Institute for Scientific Computing and Information A MESH ADAPTATION METHOD FOR 1D-BOUNDARY LAYER
More informationChapter 1. Density Estimation
Capter 1 Density Estimation Let X 1, X,..., X n be observations from a density f X x. Te aim is to use only tis data to obtain an estimate ˆf X x of f X x. Properties of f f X x x, Parametric metods f
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 informationA Modified Distributed Lagrange Multiplier/Fictitious Domain Method for Particulate Flows with Collisions
A Modified Distributed Lagrange Multiplier/Fictitious Domain Metod for Particulate Flows wit Collisions P. Sing Department of Mecanical Engineering New Jersey Institute of Tecnology University Heigts Newark,
More informationThe entransy dissipation minimization principle under given heat duty and heat transfer area conditions
Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area
More informationA posteriori error estimation of approximate boundary fluxes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2008; 24:421 434 Published online 24 May 2007 in Wiley InterScience (www.interscience.wiley.com)..1014 A posteriori error estimation
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
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 informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
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 informationCLEMSON U N I V E R S I T Y
A Fractional Step θ-metod for Convection-Diffusion Equations Jon Crispell December, 006 Advisors: Dr. Lea Jenkins and Dr. Vincent Ervin Fractional Step θ-metod Outline Crispell,Ervin,Jenkins Motivation
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 informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
More informationThird order Approximation on Icosahedral Great Circle Grids on the Sphere. J. Steppeler, P. Ripodas DWD Langen 2006
Tird order Approximation on Icosaedral Great Circle Grids on te Spere J. Steppeler, P. Ripodas DWD Langen 2006 Deasirable features of discretisation metods on te spere Great circle divisions of te spere:
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 informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationExplicit Hyperbolic Reconstructed Discontinuous Galerkin Methods for Time-Dependent Problems
AIAA AVIATION Forum June 25-29 218 Atlanta Georgia 218 Fluid Dynamics Conference 1.2514/6.218-427 Explicit Hyperbolic Reconstructed Discontinuous Galerkin Metods for Time-Dependent Problems Jialin Lou
More informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationAPPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction
Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an
More informationEfficient, unconditionally stable, and optimally accurate FE algorithms for approximate deconvolution models of fluid flow
Efficient, unconditionally stable, and optimally accurate FE algoritms for approximate deconvolution models of fluid flow Leo G. Rebolz Abstract Tis paper addresses an open question of ow to devise numerical
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 informationLecture 2: Symplectic integrators
Geometric Numerical Integration TU Müncen Ernst Hairer January February 010 Lecture : Symplectic integrators Table of contents 1 Basic symplectic integration scemes 1 Symplectic Runge Kutta metods 4 3
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 informationarxiv: v1 [math.na] 9 Mar 2018
A simple embedded discrete fracture-matrix model for a coupled flow and transport problem in porous media Lars H. Odsæter a,, Trond Kvamsdal a, Mats G. Larson b a Department of Matematical Sciences, NTNU
More informationEffects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates with Ramped Wall Temperature
Volume 39 No. February 01 Effects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates wit Ramped Wall Temperature S. Das Department of Matematics University of Gour Banga Malda 73
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 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 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 informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
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 informationAn approximation method using approximate approximations
Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More informationAPPROXIMATION BY QUADRILATERAL FINITE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More information