Numerical issues in the implementation of high order polynomial multidomain penalty spectral Galerkin methods for hyperbolic conservation laws
|
|
- Merilyn Horn
- 6 years ago
- Views:
Transcription
1 Nuerical issues in the ipleentation of high order polynoial ultidoain penalty spectral Galerkin ethods for hyperbolic conservation laws Sigal Gottlieb 1 and Jae-Hun Jung 1, 1 Departent of Matheatics, University of Massachusetts at Dartouth North Dartouth, MA 747-3, U.S.A. Abstract. In this paper, we consider high order ulti-doain penalty spectral Galerkin ethods for the approxiation of hyperbolic conservation laws. This forulation has a penalty paraeter which can vary in space and tie, allowing for flexibility in the penalty forulation. This flexibility is particularly adventageous for probles with an inhoogeneous esh. We show that the discontinuous Galerkin ethod is equivalent to the ultidoain spectral penalty Galerkin ethod with a particular penalty ter. The penalty paraeter has an effect on both the accuracy and stability of the ethod. We exaine the nuerical issues which arise in the ipleentation of high order ulti-doain penalty spectral Galerkin ethods. The coefficient truncation ethod is proposed to prevent the rapid error growth due to round-off errors when high order polynoials are used. Finally, we show that an inconsistent evaluation of the integrals in the penalty ethod ay lead to growth of errors. Nuerical exaples for linear and nonlinear probles are presented. Key words: High Order Polynoial Galerkin Methods, Penalty Boundary Conditions, Discontinuous Galerkin Methods, Hyperbolic Conservation Laws, Round-off Errors, Truncation Methods. 1 Introduction Polynoial Galerkin ethods are widely used for the nuerical solution of hyperbolic conservation laws [1,, 9, 1]. These ethods seek a polynoial approxiation of the solution for which the projected residual of the differential equation to the polynoial space vanishes. Two such classes of ethods are the spectral Galerkin ethods (sgm) and the discontinuous Galerkin ethods (dgm). Traditionally, sgm have used high order polynoials on one eleent, while dgm use lower order polynoials on any eleents. However, ulti-doain sgm exist and are known as spectral eleent ethods. In order to increase the accuracy of the approxiation, these ethods can use ore saller eleents Corresponding author. Gottlieb) Eail addresses: jjung@uassd.edu (J.-H. Jung), sgottlieb@uassd.edu (S. Global Science Preprint
2 (h-refineent) or raise the degree of the polynoial in each eleent (p refineent). High order polynoials have nuerical issues such as sensitivity to roundoff errors, so it is iportant to carefully study their effects on accuracy and stability when used in ultidoain penalty spectral Galerkin ethods. The penalty forulation penalizes the boundary or interface conditions at each eleent by introducing a penalty ter which includes a penalty paraeter. This penalty paraeter confers a great deal of flexibility on the proble, as it can change over space and tie. We deonstrate the advantages of the flexibility in the choice of penalty paraeter, especially in the case where an inhoogeneous grid syste is used. An inhoogeneous grid can be due to a difference in eleent size or in polynoial order at each eleent, but this type of grid is subject to non-physical reflecting or dispersive odes which ay appear in the solution. By odifying the penalty conditions near the grid discontinuity, we show that the sgm can reduce the non-physical odes while coputing the other eleents efficiently and accurately. We further consider the effects of the penalty ethod on the stability and accuracy. In this context, we show that the dg forulation is a special case of the penalty ulti-doain sgm. Next we discuss the effect of round-off errors for high order sgms. These round-off error effects can arise fro the the ill-conditioned ass atrix for high order polynoials, and the nuerically inconsistent evaluations of the ass atrix and the load vector. The coefficient truncation ethod is introduced to reduce round-off errors. This ethod truncates high order coefficients which do not show rapid decay in the solution of the linear syste, and prevents error growth.the second round-off error effect we explore is the error resulting fro inconsistent evaluations of the two sides of the equation. While the right-hand-side of the linear syste is evaluated by quadrature, the left-hand side can usually be coputed exactly. In the case of orthogonal polynoials, this atrix is a diagonal atrix which siplifies the process of solving the syste. However, coputing one side of the equation exactly and the other by quadrature results in inconsistency errors which rise when the polynoial order is raised. We show these nuerical inconsistency errors, and their dependence on the penalty paraeter, in nuerical coputations. The paper is structured as follows. In Section, we forulate the ultidoain penalty sgm and deonstrate the effect of the penalty ters on the stability and accuracy. The penalty sgm with non-hoogeneous grid and the flexibility of the penalty ethods are discussed. The equivalence of the dgm to the sgm is shown as well. In Section 3, the effect of round-off errors on the high order sgm is discussed. The coefficient truncation ethod is introduced and nuerical consistency is studied for the reduction of round-off errors. In Section 4, we suarize our results and discuss future research directions. Penalty Spectral Galerkin ethod In this paper, we consider the one-diensional hyperbolic conservation law u t +f(u) x =, x Ω=[ 1,1], u:[ 1,1] R + R,t>, (.1)
3 3 with the initial condition u(x,)=g(x), and boundary conditions B ± u(x,t)=h ± (t), x Ω, where B ± are the boundary operators at the doain boundary Ω. For siplicity, we will consider ainly the case f (u), for which the boundary condition is u( 1,t)=h + (t). In practice, this case can be easily generalized using the flux splitting and treating each flux with the appropriate boundary condition. In 1988, Gottlieb and Funaro introduced a penalty boundary condition for the spectral collocation approxiation of Eq. (.1) [4] and there have has been uch research on the penalty collocation ethods [3,5 8]. The ain otivation of penalty boundary conditions for collocation ethods is that the differential equation is satisfied exactly at a given set of collocation points, while at the boundaries we add a ter which penalizes for the approxiation s distance fro the prescribed boundary conditions to ake the given equations satisfied asyptotically at any points. The nuerical approxiation is the polynoial U(x,t)= k= b k(t)p k (x), where P k (x) are the basis polynoials of degree k in x and b k (t) are the unknown expansion coefficients, which will be deterined. The spectral collocation penalty ethod leads to the requireent that U t +f(u) x =τq + (x) ( f (U( 1,t)) f ( h + (t) )), at each of the collocation points x=x j. The penalty paraeter τ can depend on x and t. Q + (x) is a polynoial of degree which vanishes at all the collocation points x j, except at the boundary point x = 1, so that the penalty ter is only applied to the boundary point. To forulate the penalty Galerkin spectral ethod, we assue a solution of the for U(x,t)= k= b k(t)p k (x) and find the expansion coefficients by requiring that the projection of the residual onto the solution subspace vanish. To satisfy the boundary conditions, we can ipose a penalty ter on the strong forulation (as in the collocation case), and then require the projection of the residual of the penalized equation to vanish. Alternatively, we can ipose the penalty ter after the projection. If we choose P k (x) tobea set of orthogonal polynoials, such that Ω P k(x)p j (x)dx = γ k δ jk, the penalty Galerkin forulation yields the syste γ j b j(t)+ x P j (x)dx=τ Ωf(U) j ( f(u( 1,t)) f(h + (t)) ) Q + (x)p j (x)dx Ω for j =,,, where the superscript denotes the derivative with respect to tie and γ j are the noralization factors. We have flexibility in the choice of the penalty paraeter τ j and Q +, provided only that the polynoial Q + ( 1) = 1. One way to accoplish this is to set Q+ (x) tobe polynoials of degree, that is, Q + (x)=1=p (x). In this case, the penalty ters appear
4 4 only in the equation for b (t), that is, γ b (t) = f(u) x dx+τ ( f(u( 1,t)) f(h + (t)) ) Ω = f(u( 1,t)) f(u(1,t))+τ ( f(u( 1,t)) f(h + (t)) ), γ j b j (t) = f(u) x P j (x)dx, j =1,,. Ω We refer to this forulation as the strong forulation. However, this choice of Q + (x)=1 does not lead to good stability properties for a linear hyperbolic wave equation, as we will discuss below. A different choice of Q + (x) ay resolve this proble. For exaple, if we let Q + 1 = P k (x), γ k then the integral in the penalty ter can be evaluated exactly 1 1 Q + 1 (x)p j (x)dx= P k (x)p j (x)dx=p j ( 1). γ l 1 k= 1 k= For the orthogonal polynoials such as Chebyshev or Legendre polynoials, P j ( 1) = ( 1) j. In this case the forulation becoes siilar to the collocation case above: γ j b j(t)= f(u) x P j (x)dx+τ ( f(u( 1,t)) f(h + (t)) ) ( 1) j, Ω for j =,,. We refer to this forulation as the weak forulation. For the linear hyperbolic wave equation this penalty ter yields a better stability profile. To see the difference between these two forulations, consider the siple advection equation u t +u x = with both the strong and weak penalty forulations. Figure 1 shows that the strong forulation (represented by ) has positive real eigenvalues which will lead to instability, while the eigenvalues for the weak forulation (represented by ) are all in the left half plane. This is a clear indication that the choice of the penalty polynoial, Q + (x) is critical for stability. Fro now on, we will consider only the weak forulation polynoial. The penalty paraeter τ also plays an iportant role in the stability of the ethod. Figure shows the eigenvalues with the dgm and the penalty sgm for two different values of τ. The top figures show the eigenvalues in the coplex plane for τ = 1(left) and τ = (right) for polynoial order = 8. In each figure, the sybols and represent the eigenvalues with the dgm and the sgm respectively. Fro the left top figure it is clear that the dgm and sgm are identical for τ = 1, but when τ is increased, as in the top right figure, the real part of one eigenvalue increases. This one eigenvalue alone is responsible for the growth in the spectral radius ρ. This eigenvalue will later be iplicated in penalty sgm s increased sensitivity to roundoff errors for large τ.
5 I(λ) I(λ) Re(λ) Re(λ) Figure 1: Eigenvalues in coplex plane with the strong penalty forulation with =8(left) and = (right) for τ = 1. The sybols and represent the eigenvalues with the weak and strong penalty ethods respectively. The left botto figure in Figure shows the spectral radius ρ increasing as a function of the polynoial order for τ = 1. The right botto figure shows the spectral radius ρ as a function of the penalty paraeter τ, for various polynoial orders, on a logarithic scale. The figure shows that ρ increases fast around τ = 1. The figures show that the stability of the forulation depends on the penalty paraeters τ j, as well as Q + (x). Later we will see that the choice of penalty paraeter will affect the accuracy as well as stability of the ethod..1 Multidoain spectral penalty Galerkin ethods To increase the order of accuracy in the penalty sgms, we can allow the polynoial order to rise or we can divide the doain into any saller subdoains, and apply the sgm in each doain. To forulate the ulti-doain ethod, we divide the doain Ω=[ 1,1] into N subdoains, or eleents, I l,l=1,...,n. For siplicity, we assue here that each eleent has the sae polynoial order,, but this is not necessary, or in fact advisable, in general. The solution in each doain is given by U l (x,t)= b l k (t)p k(ξ(x)), x I l = [ ] x l 1/,x l+1/, k= where I l is the lth eleent with the doain interval x l and x l 1/ = x l x l, x l+1/ = x l + x l where x l the cell center. Finally, ξ(x) is the linear ap ξ : x [ 1,1]. For each eleent, we have for each j =,,N k= db l k (t) dt I l P j (ξ(x))p k (ξ(x))dx+ I l P j (ξ(x))f(u l ) x dx=p j I l, (.)
6 I(λ) I(λ) Re(λ) Re(λ) = 4 = 8 = 1 = 16 Log 1 ρ Log 1 ρ τ Figure : Top: Eigenvalues in coplex plane with τ = 1(left) and τ = (right) for =8. The sybols and represent the eigenvalues with the dgm and the sgm respectively. Botto: The spectral radius ρ versus τ with the dgm(left) and the sgm(right) in logarithic scale. where P j I l is the penalty ter for the jth coefficient. Assuing, as above, that f (U), the penalty ter can take the for P j I = τ j ( f(u ( 1,t)) f(h + (t)) ) ( P j ( 1) ) P j I l = τ j f(u l (x l 1,t)) f(u l 1 (x l 1,t)) P j ( 1) l =1,...,N where the τ j s ay be different in each subdoain. This forulation easily generalizes to the case where f (U) is allowed to be negative. In that case, we split the flux f =f + +f into its positive ( df + df du ) and negative ( du ) parts. For scalar hyperbolic equations, f =f + when f U > and f =f when f U <. For the syste of Eq. (.1), f ± are defined as f ± = A ± du, where A is the Jacobian atrix, i.e., A= f U. The Jacobian A is then given by A ± =T Λ ± T 1,
7 7 where T is the siilarity transforation of A and Λ + and Λ are the atrices coposed of the positive and negative eigenvalues respectively with Λ=Λ + +Λ so that f =f + +f. In this case we can extend the penalty ters to include ters not seen in a characteristic decoposition, so that the ultidoain spectral Galerkin penalty ethod becoes for each interior doain, ( ) P j I l = τ j 1 f + (U l (x l 1/ )),t)) f + (U l 1 xk (x l 1/ )),t)) Q + (ξ)p j (ξ)dξ+ Ω ( ) τ j f (U l (x l 1/ )),t)) f (U l 1 xk (x l 1/ )),t)) Q + (ξ)p j (ξ)dξ+ Ω ( ) τ j 3 f + (U l (x l+1/ )),t)) f + (U l+1 xk (x l+1/ )),t)) Q (ξ)p j (ξ)dξ+ Ω ( ) τ j 4 f (U l (x l+1/ )),t)) f (U l+1 xk (x l+1/ )),t)) Q (ξ)p j (ξ)dξ, (.3) where τ j are the penalty paraeters, and the polynoials Q can be chosen as before. If τ j ==τ j 3, the above penalty forulation is basically the characteristic decoposition. The case where τ j 1 =τ j =τ j 3 =τ j 4 is equivalent to no flux-splitting. Adjusting the coefficients allows for flexibility in applying this ethod. The ajor advantage of the flexibility of this forulation is in the case of the non-hoogeneous grid, as we see in the following exaple. Exaple.1: Consider the following equation q t +f x =, (.4) where q =(u,v) T, and f =(v,u). The initial conditions are u(x,)=sin(ωπx) and v(x,)= with ω = 5. With these initial conditions, the exact solutions are given by u(x,t)= 1 ( sin(ωπ(x t)) sin(ωπ(x+t))) and v(x,t)= 1 ( sin(ωπ(x t))+sin(ωπ(x+t))). This equation was previously investigated by Hu and Atkins for the dgm in [11]. The positive and negative fluxes are given by f + =(u+v,u+v) T and f =(u v,u v) T. We can easily show that the sgm is equivalent to the dgm with the characteristic decoposition if the following penalty paraeters are taken, that is, for q 1 =u, τ 1 = 1 =τ 4, and τ ==τ 3, Ω and for q =v, τ 1 = 1 = τ 4, and τ ==τ 3. On the other hand, the case τ 1 =τ = 1 and τ 3=τ 4 = 1 are taken for q 1=u and τ 1 =τ = 1 and τ 3 =τ 4 = 1 for q =v, there is no flux-splitting. In our nuerical exaple, we consider the sgm approxiation with a non-hoogeneous grid. A total of 5 eleents are used. In the interval x =[ 1,], there are 47 eleents
8 8 u(x,t) x u(x,t) x Figure 3: Solution u(x,t) to Eq. (.4) at t =.15 with =5 with totally 5 eleents. Left: Inhoogeneous grid with no splitting penalty ethod. Right: Inhoogeneous grid with no splitting penalty ethod except 3 interface eleents near x=. The interface of two inhoogeneous edia is at x=. each of size x i = 47 1, i =1,,47. In the interval x =[,1], there are 3 eleents with of size x i = 1 3, i=48,49,5. In each eleent, polynoial order =5 is used. In each interval x=[ 1,] and x=[,1], the solutions are sooth and there is no need to use flux splitting at all. In that case, we use the penalty ethod such that there is no flux-splitting, i.e. τ 1 = τ = τ 3 = τ 4 = 1 for u and τ 1 = τ = 1 = τ 3 = τ 4 for v. If there is no flux splitting, the ethod is easily ipleented without coputing the local flux conditions for the characteristic decoposition. However, the inhoogeneity of the grid raises the issue of possible non-physical reflection waves at the grid discontinuity at x =. The penalty ethod can be easily ipleented to reduce the artificial reflections at the grid discontinuity x =. By adopting the characteristic flux splitting only for the doains near the grid discontinuity, one can reduce the agnitude of the non-physical reflections significantly. Furtherore, one can take an advantage of the no flux splitting ethod inside each hoogeneous doain. The left figure of Figure 3 shows the non-splitting penalty ethod at t =.15. As expected, there are sall fluctuations in the solution in the region x. These fluctuations are the non-physical reflecting solutions reflected at x= and propagating to the left. The right figure shows the result with the penalty ethod where only three eleents near x= and two boundary eleents, use the characteristic penalty ethod while the other eleents use the non-splitting ethod. For the nuerical experient, for exaple, I 46,I 47 and I 48 are ipleented based on the characteristic penalty ethod. Here note that the grid discontinuity x= exists between I 47 and I 47. We also note that we use the characteristic penalty ethod for I 1 and I 5 to iniize the boundary effects. As shown in the figures, the penalty ethod is efficiently flexible to adopt the grid inhoogeneity and obtains an accurate result without any considerable non-physical reflecting odes. In [3] the sae issue, but in the collocation ethod, was also briefly discussed.
9 9. The dgm as a special case of the the ulti-doain penalty sgm In this section, we will present the dgm and then show that it can be seen as an exaple of the ulti-doain penalty sgm with a particular penalty ter τ = 1. To forulate the dgm we begin with the Galerkin for for every eleent I l, db l k (t) P j (ξ(x))p k (ξ(x))dx+ P j (ξ(x))f(u l ) x dx=, j=,,. dt I l I l k= and integrate by parts to obtain, k= db l k (t) P j (ξ(x))p k (ξ(x))dx dt I l dp j (ξ(x)) I l dx f(u l )dx= P j (ξ(x))f(u l ) x l+ 1 x l 1. (.5) for j =,,. As usual, we assue that f (U) without loss of generality. Under this assuption, the boundary ter in the above equation can be evaluated by replacing the left flux with the incoing flux based on the characteristic direction for l>, B j I l := P j (ξ(x))f(u l (x,t)) Il =P j ( 1)f(U l 1 (x l 1,t)) P j (1)f(U l (x l+ 1,t)). In the case l = where the left-ost interval in considered, we replace the incoing value by the given boundary condition, B j I = P j ( 1)f(h + (t)) P j (1)f(U l (x l+ 1,t)). (.6) To see the relationship between the dgm and sgm, we define the auxiliary integrals F and G F = and we note that F +G = I l dp j (ξ(x))) I l dx d ( ) P j (ξ(x))f(u l (x,t)) dx f(u l )dx G = P j (ξ(x))f(u l ) x dx, I l dx=p j (1)f(U l (x l+ 1,t)) P j ( 1)f(U l (x l 1,t)). (.7) With these two integrals and the penalty and boundary ters, the dgm and sgm can be rewritten, for each doain l: (dg) M b F =B Il, (sg) M b +G =P Il, (.8) where the ass atrix M and the coefficient vector b are defined by ( T M jk = P j (ξ(x))p k (ξ(x))dx, b= b l (t),,bn(t)) l. I l
10 1 x 1 7 x Figure 4: Exaple.: On the horizonal axis is the value of τ while on the vertical axis the L error of the nuerical solution at tie t =.1, with N =1. Left: The errors as a function of τ for =3. Right: The errors as a function of τ for =6. To show that the dgm and sgm forulations are equivalent, we can rearrange the dg forulation Eq. (.8) for l> M b =F +B Il =P j (1)f(UN l (x l+ 1,t)) P j ( 1)f(UN l (x l 1,t)) G+B Il, M b +G = P j (1)f(U l (x l+ 1,t)) P j ( 1)f(U l (x l 1,t))+B Il = P j (1)f(U l (x l+ 1,t)) P j ( 1)f(U l (x l 1,t)) +P j ( 1)f(U l 1 (x l 1,t)) P j (1)f(U l (x l+ 1,t)) ( = P j ( 1) f(u l 1 (x l 1,t)) f(u l (x l 1,t)) which is equal to the penalty ter P Il with τ = 1. For the case l=, the left boundary ter U l 1 (x l 1,t) is replaced by the boundary condition h + (t). Thus, the dgm forulation is just a special case of the penalty ultidoain sgm. We saw in Figures 1 and that different values of the penalty paraeter yield different stability properties. Additionally, changing the value of τ will also change the size of the errors. To see this, consider the following exaple. Exaple.: Consider the sgm for u t +u x = with the initial condition u(x,)= sin(πx) and periodic boundary conditions. Using the ultidoain sgm with N =1 and the weak penalty forulation, we exaine the the effect of the penalty paraeter τ on the errors. Figure 4 shows the L errors for different τ and =3(left) and =6(right) at t=.1. The figures show that the optial value of τ is not τ = 1, i.e. the dgm is not the optial ethod within the class of sgms. Furtherore, we observe that the optial value of τ also depends on the polynoial order. Note that L errors are 1 and 1 9 for =3 and =6, respectively. )
11 11.3 The flexibility of the penalty paraeter One of the ain advantages of the penalty sgm is that the penalty paraeters can vary depending on the proble. To see how the penalty paraeter affects the peforance of the evaluation of hyperbolic conservation laws, consider the following nuerical exaple. Exaple.3: Consider Burgers equation u t +(u /) x =, x [ 1,1], t>, (.9) with the initial condition u(x,)=x. Then the exact solution u e (x,t) is given by u e (x,t)= x 1+t. With the given initial condition, there is no incoing boundary condition for the boundaries both at x= 1 and x=1. Also, if x>, the characteristic incoing boundary conditions are applied at the left boundary of each eleent, I l.ifx<, the characteristic incoing boundary conditions are applied at the right boundary of each eleent, I + l. Since the solution is a polynoial of degree only 1, the approxiation in each eleent can be sufficiently resolved with =1, i.e. U = 1 k= b k(t)p k (x). However, the expansion coefficients b k, for k>1 are not necessarily zero in the Galerkin approxiation. The weak sgm iposes the interface conditions for every ode for which the overall approxiation is affected by round-off errors when the high order approxiation is sought. The penalty sg forulation has ore flexibility than the dg forulation to deal with such issue by exploiting the penalty paraeters. Since the solution is only a polynoial of degree 1, we eploy the following penalty paraeters { τ j τ if j =,1 = otherwise. (.1) This procedure can be autoated by exaining the regularity of the coefficients, and then designing penalty paraeters that take this into account. Figure 5 shows the sg forulation with τ j = 1, j =,, which is equivalent to the dgm and the odified penalty forulation with the condition of Eq. (.1). The figures show the L and L errors at t =.15. We use the Legendre polynoials with the total nuber of eleents being 3, with =, τ = 5 and CFL=.1. As shown in the figure, the odified penalty ethod represented by yields the best results. We also copute the errors between the exact expansion coefficients and the approxiated coefficients. The exact expansion coefficients for each doain are given in Table 1 at any tie t. Figure 6 shows the errors in the expansion coefficients for both the dgm approxiations and the odified penalty approxiations at t=t f =.15. The figures show that the errors of each expansion coefficient are larger than 1 14 for the dgm approxiation for b k,k>1 while they are close or below 1 14 for the odified penalty sgm.
12 1 SGM dgm SGM dgm Log 1 L Log 1 L Figure 5: Exaple.3: The L (left) and L (right) errors versus for Burgers equation. N =3, =, τ = 5, CFL=.1, and t f =.15. The sybols and represent the dgm and the odified sgm. Eleent 1 Eleent Eleent 3 Eleent 1 Eleent Eleent 3 Log 1 b b exact Log 1 b b exact Figure 6: The errors between the exact and approxiated coefficients in Exaple.3. Left: dgm approxiation. Right: Modified penalty approxiation. The eleent 1 is in x [ 1, 1 3 ], the eleent in x [ 1 3, 1 3 ] and the eleent 3 in x [ 1,1]. 3 Table 1: The exact expansion coefficients for Exaple.3, for each of three doains at tie t. x [ 1, 1 3 ] x [ 1 3, 1 3 ] x [1 3,1] b (t) t 3 1+t b 1 (t) 3 1+t 3 1+t 3 1+t b i (t),i>1 3 Roundoff errors of high order sgm 3.1 The coefficient truncation ethod Although a larger polynoial order should produce a ore accurate solution, in practice we see that Galerkin approxiations with higher order polynoials are ore sensitive to
13 13 roundoff errors. This sensitivity to roundoff errors can destroy the accuracy of a solution for large polynoial values. In this section, we present nuerical exaples of this proble, and suggest the coefficient truncation ethod to resolve it. Exaple 3.1: Consider, once again, the linear advection equation U t +U x =, x [ 1,1], t>. The solution of this equation is approxiated using a dg forulation with the onoial basis function ψ i (x)=x i. Figure 7 shows the decays of the transfored vector b (left) and the obtained expansion coefficient vector x (right) for the last eleent which contains the the right boundary x = 1. The paraeters used for the nuerical approxiation are total nuber of doain = 5, =37, t f =.11, and CFL=.5. The left figure in Figure 7 shows that the eleent of b decays with until around 33 with the agnitude 1. Beyond >33 it is observed not to decrease. This is due to the large condition nuber, κ, of the stiffness atrix M, i.e. κ The right figure in Figure 7 shows the decay or growth of the evaluated expansion coefficient vector x. The figure clearly shows that the expansion coefficients grow rapidly with. Figure 8 deonstrates that this results in large errors. This siple exaple shows that for high order polynoials, the calculation of the coefficients is very sensitive to round-off error, and this affects the accuracy of the solution. To resolve this proble, we odify the penalty ethod with the truncation ethod introduced in [13] for use in the inverse polynoial reconstruction ethod (IPRM) [1, 14] to reduce deterioration of the error when the polynoial error was large. To apply the coefficient truncation ethod, we look at the coefficients resulting fro the interediate step in the Gaussian eliination. Suppose that the set of expansion coefficients is to be found fro the linear equation Mb =h, then the syste is solved by Gaussian eliination to yield an upper triangular atrix syste U b =c. The coefficients of this syste should be rapidly decaying, so to reduce the proble of round-off errors we ipose this requireent by setting c i = { ci if c i >ɛ t otherwise. (3.1) Here ɛ t is the tolerance level which is to be deterined depending on the decay rate of c. The otivation behind this ethod is that, as shown in Figure 1, the spectral radius ρ increases with for the Galerkin ethod and this akes the ethod sensitive to roundoff errors. The truncation ethod tries to reduce the agnitude of ρ by truncating the RHS of the linear syste of the Galerkin ethod by reducing the rank of the atrix. For exaple, if the truncation order is n, then the rank of the truncated atrix becoes +1 n. With the rank reduced, ρ is also reduced. Exaple 3.: To deonstrate the success of the truncation ethod in reducing the effects of round-off errors, we return to Exaple 3.1 above. We apply the truncation ethod with various tolerence levels, i.e. ɛ t =1 14 (+), ɛ =1 1 ( ), ɛ =1 1 ( ), and ɛ=1 ( ). The left figure in Figure 7 shows that the the eleents of b decay faster when the truncation error is applied for different values of ɛ. The right figure in Figure 7 shows the evaluated expansion coefficient vector x. The figure clearly shows that the expansion
14 14 Log 1 b No Truncation = 1 14 = 1 1 = 1 1 = Log 1 b No Truncation = 1 14 = 1 1 = 1 1 = Figure 7: Exaples 3.1 and 3.: Left: The decay of the eleents of the load vector c with and without the truncation ethod. Right: The decay of the obtained expansion coefficient vector x with and without the truncation ethod. The dgm is used for the advection equation with onoial basis functions. Total nuber of doain is 5, =37, t f =.11, and CFL=.5. Various tolerance levels ɛ t are used, i.e. ɛ t = 1 14,1 1,1 1,1. Log 1 L No Truncation = 1 14 = 1 1 = 1 1 = Log 1 Error 8 1 = 1 1 No Truncation = 1 14 = 1 1 = x Figure 8: Exaples 3.1 and 3.: Left: The L errors versus. Right: The pointwise errors when =37 is used with and without the truncation ethod. The dgm is used for the advection equation with onoial basis functions. Total nuber of doain is 5, t f =.11, and CFL=.5. Various tolerance levels ɛ t are used, i.e. ɛ t =1 14,1 1,1 1,1. coefficients grow with exponentially without the truncation ethod while they decay if the truncation ethod is applied. The figures also shows that the expansion coefficients becoe truncated to alost achine zeros or very close to zeros if the truncation ethod is applied. In the figure, these truncated values are not displayed as the plot uses the logarithic scale. The larger ɛ t is used, the faster the truncation of x occurs. The L errors with (left) and the pointwise errors with =37 (right) are given in Figure 8 with the sae sybols as in Figure 7. The left figure in Figure 8 shows the L decay or growth with with and without the truncation ethod. The figure shows that the L errors decay up to 7 and they start to grow beyond if the truncation ethod is not applied or the truncation ethod is applied with ɛ t =1 14. The truncation ethod with ɛ =1 1 or 1 1 yields the best results up to 4. If ɛ t =1 is used,
15 15 the L reains around 1 for large. It is observed that for >4, the L errors reains around 1 when ɛ t =1 while the other cases result in the growth of L errors. The right figure in Figure 8 shows the pointwise errors with = 37. The figure shows that the best result is obtained when ɛ t =1 1 is used. Figure 8 suggests that the tolerance level used with the truncation ethod can be chosen by observing the decay rate of c for a given. Also, one can use the different tolerance level ɛ t for different eleent since each eleent has different decay rate of c. 3. Nuerical consistency Another factor which increases the sensitivity to round-off errors is the inconsistent evaluation of the integrals in the forulation of the syste of ODEs for the approxiation coefficients. On the Galerkin forulationwe need to copute two sets of integrals, I 1 = Ω P k(x)p j (x)dx and I = Ω f(u) xp j (x)dx. While I 1 can be evaluated exactly, because the for of the polynoials is known explicitly in general, I can not generally be evaluated, and so quadrature ust be used. We define the integral operators E and Q, where E is the exact integration operator E(g)= Ωg(x), and Q soe quadrature rule such as Gauss quadrature rule used to evaluate the given integral Q(g)= M i= g iω i Ωg, where ω i are the weights. The Galerkin forulation can now be perfored in one of four ways: E(P l P j )=E(f(U) x P j )+P l, E(P l P j )=Q(f(U) x P j )+P l, (3.) Q(P l P j )=E(f(U) x P j )+P l, Q(P l P j )=Q(f(U) x P j )+P l. We say that the Galerkin forulation is nuerically consistent if the integrals are evaluated in the sae way on both sides of the equations (both by E or both by Q), and it is nuerically inconsistent if the two integrals are evaluated differently on the two sides of the equation. The phenoenon we are considering is essentially a nuerical one. For high polynoial order, the nuerically consistent forulation yields ore accurate results, and lower sensitivity to roundoff errors, than the inconsistent forulation. Consequently the above 4 different forulations, Eq. (3.), show different errors. Exaple 3.3: We first consider the following steady-state proble u t +u x =sin(πx), x [ 1,1],t>, with initial condition u(x,) = u s (x)+ɛ M δ(x), where ɛ M is achine accuracy, ɛ M =1 16 and δ(x) is the rando function, δ(x) 1 with the noral distribution. The boundary condition is u( 1,t)= 1 π cos( π), t>. The steady-state solution, us =li t u(x,t)= u s (x) for t is then given by u s (x,t)= π 1 cos(π(x t)). For the nuerical experient, the final tie is t f =1, the total nuber of doain N =3, the penalty paraeter τ = 5, and the CFL nuber CFL=.1 for dt=cfl dx= are used where each eleent has the sae eleent size and polynoial order.
16 16 SGM QQ SGM EQ dgm QQ dgm EQ SGM QQ SGM EQ dgm QQ dgm EQ Log 1 L 8 1 Log 1 L Figure 9: Left: L errors vs for the steady-state proble u t+u x = sin(πx). The sybols and denote the quadrature-quadrature (QQ) and exact-quadrature (EQ) forulations for the penalty sgm and and the QQ and EQ forulations for the dgm. The final tie t f =1 and the total nuber of doain is 3. The penalty parater is τ = 5 and the CFL nuber is CFL=.1, i.e. dt= Figure 9 shows the convergence of L and L errors with for the the weak penalty forulation and the dg forulation. The LHS or RHS are evaluated either using the quadrature rules labeled by Q or the exact forula denotedlabeled by E for the penalty forulation with τ = 5. As the figures show, both the L and L errors decay exponentially until around 7. For >7, these errors grow slightly due to round-off errors. The figures show that the results with the consistent evaluation of the stiffness atrix of the sgm (QQ) show better perforance for the weak penalty sgm when round-off error becoe doinant. We note that although the dgm does not suffer fro an inconsistent forulation, the consistent forulation for the penalty sgm yields the best results when is large. 4 Suary and conclusion We describe the forulation of the ulti-doain penalty sgm and deonstrated that the flexibility of the penalty forulation can be advantageous, because it allows us to tailor the penalty paraeter to atch the proble. This is especially relevant in the case where the sub-doains have different esh size or polynoial order, as the flexibility in the penalty paraeters can allow us to avoid costly flux splitting except near the grid discontinuity. However, different values of the penalty paraeter effect the stability, accuracy, and sensitivity to round-off errors. For exaple, the dgm is siply a sgm with a particular choice of penalty paraeter, which for a linear wave equation has nice stability properties, though it is not optial in ters of accuracy. We presented two nuerical issues which arise in the solution of hyperbolic conservation laws using the ulti-doain penalty sgm with high order polynoials. The first is the sensitivity of high order sgms to roundoff errors, which can potentially ruin the accu-
17 17 racy of the solution. To resolve this issue we introduce the coefficient truncation ethod which prevents the rapid growth of the errors with, and has a stabilizing effect on the ethod. The second nuerical issue is that consistent evaluation of the stiffness atrix and the load vector yields a better result when the high order polynoials are used with the penalty forulation. This sensitivity, too, depends on the penalty paraeter, and we note that the case of τ = 1, (the dgm case) does not see affected by this nuerical consistency issue. Future studies will center around ethods for choosing the penalty paraeter to optiize for accuracy and stability, as well as further developent of the coefficient truncation ethod for ulti-diensional probles. 5 Acknowledgeents The authors gratefully acknowledges the advice of David Gottlieb and Jennifer Ryan. This work has been supported by the NSF under Grant No. DMS References [1] B. Cockburn, C. Johnson, C.-W. Shu, E. Tador, Advanced nuerical approxiation of nonlinear hyperbolic equations, Lecture Notes in Matheatics, Springer, Berlin, [] B. Cockburn, G. Karniadakis, C.-W. Shu, The developent of discontinuous Galerkin ethods, in: Discontinuous Galerkin ethods (Newport, RI, 1999), Lecture Notes in Coputer Science and Engineering, vol. 11, Springer, Berlin, pp [3] W.-S. Don, D. Gottlieb, J.-H. Jung, A ultidoain spectral ethod for supersonic reactive flows, J. Coput. Phys. 19 (3) [4] D. Funaro, D.Gottlieb, A New Method of Iposing Boundary Conditions for Hyperbolic Equations, Math. Coput. 51 (1988) [5] D. Funaro, D. Gottlieb, Convergence Results for Pseudospectral Approxiations of Hyperbolic Systes by a Penalty-type Boundary Treatent, Math. Coput. 57 (1991) [6] J. S. Hesthaven, D. Gottlieb, A Stable Penalty Method for the Copressible Navier Stokes Equations: I. Open Boundary Conditions, SIAM Journal on Scientific Coputing 17 (1996) [7] J. S. Hesthaven, A Stable Penalty Method for the Copressible Navier Stokes Equations: II. One-Diensional Doain Decoposition Schees, SIAM Journal on Scientific Coputing, 18 (1997) [8] J. S. Hesthaven, S. Gottlieb, D. Gottlieb, Spectral ethods for tie-dependent probles (Cabridge, 7), Cabridge UP. [9] J. S. Hesthaven and T. Warburton,, High-Order Nodal Methods on Unstructured Grids. I. Tie-Doain Solution of Maxwell s Equations, J. Coput. Phys. 181(1), [1] F. X. Giraldo, J. S. Hesthaven, and T. Warburton,, Nodal High-Order Discontinuous Galerkin Method for the Spherical Shallow Water Equations, J. Coput. Phys. 181(), [11] F. Q. Hu, H. L. Atkins, Eigensolution analysis of the discontinuous Galerkin ethod with nonunifor grids: I. one space diension, J. Coput. Phy. 18 ()
18 18 [1] J.-H. Jung, B. D. Shizgal, Generalization of the inverse polynoial reconstruction ethod in the resolution of the Gibbs phenoena, J. Coput. Appl. Math. 17 (4) [13] J.-H. Jung, B. D. Shizgal, On the nuerical convergence with the inverse polynoial reconstruction ethod for the resolution of the Gibbs phenoenon, J. Coput. Phys. 4 (7) [14] B. D. Shizgal, J.-H Jung, Towards the resolution of the Gibbs phenoena, J. Coput. Appl. Math. 161 (3)
lecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 38: Linear Multistep Methods: Absolute Stability, Part II
lecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 3: Linear Multistep Methods: Absolute Stability, Part II 5.7 Linear ultistep ethods: absolute stability At this point, it ay well
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN
More informationSupporting Information for Supression of Auger Processes in Confined Structures
Supporting Inforation for Supression of Auger Processes in Confined Structures George E. Cragg and Alexander. Efros Naval Research aboratory, Washington, DC 20375, USA 1 Solution of the Coupled, Two-band
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationComparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations
International Journal of Applied Science and Technology Vol. 7, No. 3, Septeber 217 Coparison of Stability of Selected Nuerical Methods for Solving Stiff Sei- Linear Differential Equations Kwaku Darkwah
More informationEE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng
EE59 Spring Parallel LSI AD Algoriths Lecture I interconnect odeling ethods Zhuo Feng. Z. Feng MTU EE59 So far we ve considered only tie doain analyses We ll soon see that it is soeties preferable to odel
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationNumerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationStability Analysis of the Matrix-Free Linearly Implicit 2 Euler Method 3 UNCORRECTED PROOF
1 Stability Analysis of the Matrix-Free Linearly Iplicit 2 Euler Method 3 Adrian Sandu 1 andaikst-cyr 2 4 1 Coputational Science Laboratory, Departent of Coputer Science, Virginia 5 Polytechnic Institute,
More informationOn the characterization of non-linear diffusion equations. An application in soil mechanics
On the characterization of non-linear diffusion equations. An application in soil echanics GARCÍA-ROS, G., ALHAMA, I., CÁNOVAS, M *. Civil Engineering Departent Universidad Politécnica de Cartagena Paseo
More informationSOLUTIONS. PROBLEM 1. The Hamiltonian of the particle in the gravitational field can be written as, x 0, + U(x), U(x) =
SOLUTIONS PROBLEM 1. The Hailtonian of the particle in the gravitational field can be written as { Ĥ = ˆp2, x 0, + U(x), U(x) = (1) 2 gx, x > 0. The siplest estiate coes fro the uncertainty relation. If
More information. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe
PROPERTIES OF MULTIVARIATE HOMOGENEOUS ORTHOGONAL POLYNOMIALS Brahi Benouahane y Annie Cuyt? Keywords Abstract It is well-known that the denoinators of Pade approxiants can be considered as orthogonal
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationKernel Methods and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley ENSIAG 2 / osig 1 Second Seester 2012/2013 Lesson 20 2 ay 2013 Kernel ethods and Support Vector achines Contents Kernel Functions...2 Quadratic
More informationPattern Recognition and Machine Learning. Artificial Neural networks
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lessons 7 20 Dec 2017 Outline Artificial Neural networks Notation...2 Introduction...3 Key Equations... 3 Artificial
More informationOn the approximation of Feynman-Kac path integrals
On the approxiation of Feynan-Kac path integrals Stephen D. Bond, Brian B. Laird, and Benedict J. Leikuhler University of California, San Diego, Departents of Matheatics and Cheistry, La Jolla, CA 993,
More informationRECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE
Proceedings of ICIPE rd International Conference on Inverse Probles in Engineering: Theory and Practice June -8, 999, Port Ludlow, Washington, USA : RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS
More informationIntelligent Systems: Reasoning and Recognition. Artificial Neural Networks
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley MOSIG M1 Winter Seester 2018 Lesson 7 1 March 2018 Outline Artificial Neural Networks Notation...2 Introduction...3 Key Equations... 3 Artificial
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationREDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION
ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas
More informationSD Numerical Simulation Technique for Hydrodynamic Flow Gas- Solids Mixing Mantilla Núñez, Irla* 1 S. De Vicente C. 2
SD Nuerical Siulation Technique for Hydrodynaic Flow Gas- Solids Mixing Mantilla Núñez, Irla* S. De Vicente C. National University of Engineering, Lia, Perú, Polytechnic University of Madrid. *Lia 5, irlan@uni.edu.pe.
More informationFast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials
Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter
More information1 Bounding the Margin
COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #12 Scribe: Jian Min Si March 14, 2013 1 Bounding the Margin We are continuing the proof of a bound on the generalization error of AdaBoost
More informationA model reduction approach to numerical inversion for a parabolic partial differential equation
Inverse Probles Inverse Probles 30 (204) 250 (33pp) doi:0.088/0266-56/30/2/250 A odel reduction approach to nuerical inversion for a parabolic partial differential equation Liliana Borcea, Vladiir Drusin
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationQuantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search
Quantu algoriths (CO 781, Winter 2008) Prof Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search ow we begin to discuss applications of quantu walks to search algoriths
More informationCelal S. Konor Release 1.1 (identical to 1.0) 3/21/08. 1-Hybrid isentropic-sigma vertical coordinate and governing equations in the free atmosphere
Celal S. Konor Release. (identical to.0) 3/2/08 -Hybrid isentropic-siga vertical coordinate governing equations in the free atosphere This section describes the equations in the free atosphere of the odel.
More informationThe linear sampling method and the MUSIC algorithm
INSTITUTE OF PHYSICS PUBLISHING INVERSE PROBLEMS Inverse Probles 17 (2001) 591 595 www.iop.org/journals/ip PII: S0266-5611(01)16989-3 The linear sapling ethod and the MUSIC algorith Margaret Cheney Departent
More informationIntroduction to Robotics (CS223A) (Winter 2006/2007) Homework #5 solutions
Introduction to Robotics (CS3A) Handout (Winter 6/7) Hoework #5 solutions. (a) Derive a forula that transfors an inertia tensor given in soe frae {C} into a new frae {A}. The frae {A} can differ fro frae
More informationSharp Time Data Tradeoffs for Linear Inverse Problems
Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used
More informationarxiv: v1 [physics.comp-ph] 20 Oct 2017 Abstract
Convergence Properties & Extension to Non-Conservative Multidiensional ystes Haran Jackson Cavendish Laboratory, JJ Thoson Ave, Cabridge, UK, CB3 0HE arxiv:70.07603v [physics.cop-ph] 20 Oct 207 Abstract
More informationThis model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t.
CS 493: Algoriths for Massive Data Sets Feb 2, 2002 Local Models, Bloo Filter Scribe: Qin Lv Local Models In global odels, every inverted file entry is copressed with the sae odel. This work wells when
More informationOptimal Jamming Over Additive Noise: Vector Source-Channel Case
Fifty-first Annual Allerton Conference Allerton House, UIUC, Illinois, USA October 2-3, 2013 Optial Jaing Over Additive Noise: Vector Source-Channel Case Erah Akyol and Kenneth Rose Abstract This paper
More information3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme
P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University
More informationNumerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term
Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2
More informationComput. Methods Appl. Mech. Engrg.
Coput. Methods Appl. Mech. Engrg. 200 (2011) 2768 2781 Contents lists available at ScienceDirect Coput. Methods Appl. Mech. Engrg. journal hoepage: www.elsevier.co/locate/ca A posteriori error analysis
More informationPHY307F/407F - Computational Physics Background Material for Expt. 3 - Heat Equation David Harrison
INTRODUCTION PHY37F/47F - Coputational Physics Background Material for Expt 3 - Heat Equation David Harrison In the Pendulu Experient, we studied the Runge-Kutta algorith for solving ordinary differential
More informationOn the hyperbolicity of two- and three-layer shallow water equations. 1 Introduction
On the hyperbolicity of two- and three-layer shallow water equations M. Castro, J.T. Frings, S. Noelle, C. Parés, G. Puppo RWTH Aachen, Universidad de Málaga, Politecnico di Torino E-ail: frings@igp.rwth-aachen.de
More informationModeling Chemical Reactions with Single Reactant Specie
Modeling Cheical Reactions with Single Reactant Specie Abhyudai Singh and João edro Hespanha Abstract A procedure for constructing approxiate stochastic odels for cheical reactions involving a single reactant
More informationPhysics 139B Solutions to Homework Set 3 Fall 2009
Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about
More informationANALYSIS OF A NUMERICAL SOLVER FOR RADIATIVE TRANSPORT EQUATION
ANALYSIS OF A NUMERICAL SOLVER FOR RADIATIVE TRANSPORT EQUATION HAO GAO AND HONGKAI ZHAO Abstract. We analyze a nuerical algorith for solving radiative transport equation with vacuu or reflection boundary
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More informationESTIMATE OF THE TRUNCATION ERROR OF FINITE VOLUME DISCRETISATION OF THE NAVIER-STOKES EQUATIONS ON COLOCATED GRIDS.
ESTIMATE OF THE TRUNCATION ERROR OF FINITE VOLUME DISCRETISATION OF THE NAVIER-STOKES EQUATIONS ON COLOCATED GRIDS. Alexandros Syrakos, Apostolos Goulas Departent of Mechanical Engineering, Aristotle University
More informationHESSIAN MATRICES OF PENALTY FUNCTIONS FOR SOLVING CONSTRAINED-OPTIMIZATION PROBLEMS
R 702 Philips Res. Repts 24, 322-330, 1969 HESSIAN MATRICES OF PENALTY FUNCTIONS FOR SOLVING CONSTRAINED-OPTIMIZATION PROBLEMS by F. A. LOOTSMA Abstract This paper deals with the Hessian atrices of penalty
More informationInspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information
Cite as: Straub D. (2014). Value of inforation analysis with structural reliability ethods. Structural Safety, 49: 75-86. Value of Inforation Analysis with Structural Reliability Methods Daniel Straub
More informationGeneralized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili,
Australian Journal of Basic and Applied Sciences, 5(3): 35-358, 20 ISSN 99-878 Generalized AOR Method for Solving Syste of Linear Equations Davod Khojasteh Salkuyeh Departent of Matheatics, University
More informationA Simple Regression Problem
A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where
More informationA BLOCK MONOTONE DOMAIN DECOMPOSITION ALGORITHM FOR A NONLINEAR SINGULARLY PERTURBED PARABOLIC PROBLEM
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volue 3, Nuber 2, Pages 211 231 c 2006 Institute for Scientific Coputing and Inforation A BLOCK MONOTONE DOMAIN DECOMPOSITION ALGORITHM FOR A NONLINEAR
More information12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015
18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.
More informationNUMERICAL SOLUTION OF BACKWARD HEAT CONDUCTION PROBLEMS BY A HIGH ORDER LATTICE-FREE FINITE DIFFERENCE METHOD
Journal of the Chinese Institute of Engineers, Vol. 7, o. 4, pp. 6-60 (004) 6 UMERICAL SOLUTIO OF BACKWARD HEAT CODUCTIO PROBLEMS BY A HIGH ORDER LATTICE-FREE FIITE DIFFERECE METHOD Kentaro Iijia ABSTRACT
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More information2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all
Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either
More informationGeneralized eigenfunctions and a Borel Theorem on the Sierpinski Gasket.
Generalized eigenfunctions and a Borel Theore on the Sierpinski Gasket. Kasso A. Okoudjou, Luke G. Rogers, and Robert S. Strichartz May 26, 2006 1 Introduction There is a well developed theory (see [5,
More informationMSEC MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL SOLUTION FOR MAINTENANCE AND PERFORMANCE
Proceeding of the ASME 9 International Manufacturing Science and Engineering Conference MSEC9 October 4-7, 9, West Lafayette, Indiana, USA MSEC9-8466 MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL
More informationProjectile Motion with Air Resistance (Numerical Modeling, Euler s Method)
Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,
More informationA Generalized Permanent Estimator and its Application in Computing Multi- Homogeneous Bézout Number
Research Journal of Applied Sciences, Engineering and Technology 4(23): 5206-52, 202 ISSN: 2040-7467 Maxwell Scientific Organization, 202 Subitted: April 25, 202 Accepted: May 3, 202 Published: Deceber
More informationScattering and bound states
Chapter Scattering and bound states In this chapter we give a review of quantu-echanical scattering theory. We focus on the relation between the scattering aplitude of a potential and its bound states
More informationAnalysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning
Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics
More informationCompressive Distilled Sensing: Sparse Recovery Using Adaptivity in Compressive Measurements
1 Copressive Distilled Sensing: Sparse Recovery Using Adaptivity in Copressive Measureents Jarvis D. Haupt 1 Richard G. Baraniuk 1 Rui M. Castro 2 and Robert D. Nowak 3 1 Dept. of Electrical and Coputer
More informationRecovery of high order accuracy in radial basis function approximation for discontinuous problems
Recovery of high order accuracy in radial basis function approximation for discontinuous problems Chris L. Bresten, Sigal Gottlieb 1, Daniel Higgs, Jae-Hun Jung* 2 Abstract Radial basis function(rbf) methods
More informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationRESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS
BIT Nuerical Matheatics 43: 459 466, 2003. 2003 Kluwer Acadeic Publishers. Printed in The Netherlands 459 RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS V. SIMONCINI Dipartiento di
More informationOn Conditions for Linearity of Optimal Estimation
On Conditions for Linearity of Optial Estiation Erah Akyol, Kuar Viswanatha and Kenneth Rose {eakyol, kuar, rose}@ece.ucsb.edu Departent of Electrical and Coputer Engineering University of California at
More informationModelling diabatic atmospheric boundary layer using a RANS-CFD code with a k-ε turbulence closure F. VENDEL
Modelling diabatic atospheric boundary layer using a RANS-CFD code with a k-ε turbulence closure F. VENDEL Florian Vendel 1, Guillevic Laaison 1, Lionel Soulhac 1, Ludovic Donnat 2, Olivier Duclaux 2,
More informationProbability Distributions
Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationHyperbolic Horn Helical Mass Spectrometer (3HMS) James G. Hagerman Hagerman Technology LLC & Pacific Environmental Technologies April 2005
Hyperbolic Horn Helical Mass Spectroeter (3HMS) Jaes G Hageran Hageran Technology LLC & Pacific Environental Technologies April 5 ABSTRACT This paper describes a new type of ass filter based on the REFIMS
More informationDESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *
IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS
More informationSupplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators
Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical
More informationThe Weierstrass Approximation Theorem
36 The Weierstrass Approxiation Theore Recall that the fundaental idea underlying the construction of the real nubers is approxiation by the sipler rational nubers. Firstly, nubers are often deterined
More informationVulnerability of MRD-Code-Based Universal Secure Error-Correcting Network Codes under Time-Varying Jamming Links
Vulnerability of MRD-Code-Based Universal Secure Error-Correcting Network Codes under Tie-Varying Jaing Links Jun Kurihara KDDI R&D Laboratories, Inc 2 5 Ohara, Fujiino, Saitaa, 356 8502 Japan Eail: kurihara@kddilabsjp
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Notes for EE227C (Spring 2018): Convex Optiization and Approxiation Instructor: Moritz Hardt Eail: hardt+ee227c@berkeley.edu Graduate Instructor: Max Sichowitz Eail: sichow+ee227c@berkeley.edu October
More informationHybrid System Identification: An SDP Approach
49th IEEE Conference on Decision and Control Deceber 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Hybrid Syste Identification: An SDP Approach C Feng, C M Lagoa, N Ozay and M Sznaier Abstract The
More informationlecture 36: Linear Multistep Mehods: Zero Stability
95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,
More informationma x = -bv x + F rod.
Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous
More informationUCSD Spring School lecture notes: Continuous-time quantum computing
UCSD Spring School lecture notes: Continuous-tie quantu coputing David Gosset 1 Efficient siulation of quantu dynaics Quantu echanics is described atheatically using linear algebra, so at soe level is
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationMechanics Physics 151
Mechanics Physics 5 Lecture Oscillations (Chapter 6) What We Did Last Tie Analyzed the otion of a heavy top Reduced into -diensional proble of θ Qualitative behavior Precession + nutation Initial condition
More informationAnisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS D: APPLIED PHYSICS J. Phys. D: Appl. Phys. 5 00 007 04 PII: S00-770867-6 Anisotropic reference edia and the possible linearized approxiations for phase
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More informationExplicit solution of the polynomial least-squares approximation problem on Chebyshev extrema nodes
Explicit solution of the polynoial least-squares approxiation proble on Chebyshev extrea nodes Alfredo Eisinberg, Giuseppe Fedele Dipartiento di Elettronica Inforatica e Sisteistica, Università degli Studi
More informationEfficient Filter Banks And Interpolators
Efficient Filter Banks And Interpolators A. G. DEMPSTER AND N. P. MURPHY Departent of Electronic Systes University of Westinster 115 New Cavendish St, London W1M 8JS United Kingdo Abstract: - Graphical
More informationDETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION
DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION Masaki WAKUI 1 and Jun IYAMA and Tsuyoshi KOYAMA 3 ABSTRACT This paper shows a criteria to detect
More informationMulti-Scale/Multi-Resolution: Wavelet Transform
Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationAn Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles
An Inverse Interpolation Method Utilizing In-Flight Strain Measureents for Deterining Loads and Structural Response of Aerospace Vehicles S. Shkarayev and R. Krashantisa University of Arizona, Tucson,
More informationA Note on the Applied Use of MDL Approximations
A Note on the Applied Use of MDL Approxiations Daniel J. Navarro Departent of Psychology Ohio State University Abstract An applied proble is discussed in which two nested psychological odels of retention
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin Heilbronner Str. 10, D Berlin - Wilmersdorf
Konrad-Zuse-Zentru für Inforationstechnik Berlin Heilbronner Str. 10, D-10711 Berlin - Wilersdorf Folkar A. Borneann On the Convergence of Cascadic Iterations for Elliptic Probles SC 94-8 (Marz 1994) 1
More information13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices
CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay
More informationComputable Shell Decomposition Bounds
Coputable Shell Decoposition Bounds John Langford TTI-Chicago jcl@cs.cu.edu David McAllester TTI-Chicago dac@autoreason.co Editor: Leslie Pack Kaelbling and David Cohn Abstract Haussler, Kearns, Seung
More informationAN APPLICATION OF CUBIC B-SPLINE FINITE ELEMENT METHOD FOR THE BURGERS EQUATION
Aksan, E..: An Applıcatıon of Cubıc B-Splıne Fınıte Eleent Method for... THERMAL SCIECE: Year 8, Vol., Suppl., pp. S95-S S95 A APPLICATIO OF CBIC B-SPLIE FIITE ELEMET METHOD FOR THE BRGERS EQATIO by Eine
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More informationProc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES
Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co
More informationNecessity of low effective dimension
Necessity of low effective diension Art B. Owen Stanford University October 2002, Orig: July 2002 Abstract Practitioners have long noticed that quasi-monte Carlo ethods work very well on functions that
More informationOn the Use of A Priori Information for Sparse Signal Approximations
ITS TECHNICAL REPORT NO. 3/4 On the Use of A Priori Inforation for Sparse Signal Approxiations Oscar Divorra Escoda, Lorenzo Granai and Pierre Vandergheynst Signal Processing Institute ITS) Ecole Polytechnique
More information