Higher-Order Accurate Schemes for a Constant Coefficient Singularly Perturbed Reaction-Diffusion Problem of Boundary Layer Type
|
|
- Arnold James
- 5 years ago
- Views:
Transcription
1 Applied Mathematical Sciences, Vol., 6, no. 4, HIKARI Ltd, Higher-Order Accurate Schemes for a Constant Coefficient Singularly Perturbed Reaction-Diffusion Problem of Boundary Layer Type Brian J. M c Cartin Applied Mathematics, Kettering University 7 University Avenue, Flint, MI USA Copyright c 6 Brian J. M c Cartin. This article is distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract This paper addresses the numerical approximation of linear, steady, one-dimensional, constant coefficient singular perturbation problems of reaction-diffusion type which exhibit boundary layers. Since the diffusion and reaction terms are of opposite sign, this is a special case of the linear source problem []. Herein, an original derivation of the hyperbolically fitted discretization [8], second-order accurate uniformly in the small parameter, due to Miller [7] is presented. It is then shown how to apply the Mehrstellenverfahren technique of L. Collatz [] to produce formal fourth-order accuracy. This is then the springboard to schemes of uniform fourth-order and formal sixth-order accuracy. Mathematics Subject Classification: 65L, 65L, 65L Keywords: singular perturbation, reaction-diffusion, boundary layer Introduction A singular perturbation problem is said to be of reaction-diffusion type if the order of the problem is reduced by two when the small parameter is set Deceased January 9, 6
2 64 Brian J. McCartin to zero []. This is to be distinguished from a convection-diffusion problem [9] where the corresponding reduction in order is only one. Depending on the relative signs of the diffusion and reaction terms, the solution to such a reaction-diffusion equation may exhibit either boundary layer behavior (the subject of the present paper) or become highly oscillatory [6]. While the most general of such problems are nonlinear, transient and multidimensional, this paper focuses on the linear, steady-state, one-dimensional version with constant coefficients. The variable coefficient problem is the subject of [5]. Within these confines, we present new methods of numerical approximation that are uniformly fourth-order accurate and formally sixth-order accurate. This will be achieved by coupling the method of exponential fitting [3] with the Mehrstellenverfahren technique of L. Collatz []. In the development that follows, the focus will be on the boundary layer case. We begin by examining the local truncation error of central differencing as applied to this problem as well as that of Miller s scheme [7]. It is then shown how to apply the Mehrstellenverfahren technique (Hermitian method) of L. Collatz [] to produce formal fourth-order accuracy. These considerations reveal that the respective local truncation errors predict the known asymptotic behaviors of these varied methods. Motivated by this observation, we next provide a new derivation of Miller s scheme for the case of constant coefficients. The approach taken herein amounts to the annihilation of all objectionable terms (i.e. terms that misbehave in the boundary layer) in the local truncation error. Moreover, this approach indicates how to achieve uniform fourth-order accuracy in the small parameter ɛ for constant coefficients by combining exponential fitting with Mehrstellenverfahren. This then provides the springboard to a scheme of formal sixth-order accuracy. Finally, we turn to a comparison of central differencing, Miller s scheme, Mehrstellenverfahren as well as our new schemes on a constant coefficient boundary layer problem with known exact solution. These results are indeed very encouraging. Reaction-Diffusion Equation The subject of our investigation is the steady, one-dimensional, reactiondiffusion equation with constant coefficients. The corresponding self-adjoint two-point boundary value problem may be cast in the form ɛ u (x) + c u(x) = f(x); u() = u L, u() = u R, () where we assume that < ɛ and c >. Physically, ɛ is the diffusion coefficient while c is the reaction coefficient. Since ɛ is a small parameter, we
3 Constant coefficient reaction-diffusion problem of boundary layer type 65 have here a singular perturbation problem. When (as herein) ɛ >, boundary layers may develop while, for ɛ <, the solution may be highly oscillatory [6]. In what follows, all functions appearing (especially, the source term f(x)) are assumed to possess sufficient smoothness to justify any operations to which we may subject them. Also, we define the parameter z := h c/ɛ, which measures the strength of reaction relative to diffusion when scaled by the mesh width h := /N of Figure. For the problems of primary interest in this study, we have z. Lastly, in the ensuing derivations, we will restrict ourselves to ɛ > (boundary layer case) where singular perturbation theory [8] reveals that the boundary layers may be as thin as O( ɛ) and that in such a boundary layer u = O(), u = O(/ ɛ), and u = O(/ɛ). Figure : Computational Grid 3 The Central Difference Scheme If we discretize the two-point boundary value problem, Equation (), using central differences, we obtain ɛ ui+ u i + u i h + c u i = f i. () The local truncation error of this scheme, which is easily established using Taylor series, is given by L.T.E. = h [ f i }{{} O() c i u i }{{} O(/ɛ) ] + O(h 4 ). (3)
4 66 Brian J. McCartin In order to numerically resolve any boundary layers that arise, we require that h ɛ (boundary layer width). Perusal of the above local truncation error then reveals that it is effectively O() due to the presence of the term containing u i. Thus, in spite of the fact that this scheme is formally second-order accurate, it is not uniformly accurate to any degree in the small parameter ɛ. 4 The Miller Scheme Miller [7] has developed a scheme for Equation () which is second-order uniformly accurate for constant c. This scheme is defined as follows: ɛ z [cosh z ] ui+ u i + u i h + c u i = f i. (4) The appearance of hyperbolic functions in this context is referred to as exponential fitting [8]. The local truncation error of this scheme is given by L.T.E. = h [ f i ] + O(h 4 ). (5) }{{} Notice that, in comparison to Equation (3), the undesirable term involving u i has been eliminated. Thus, the local truncation error is O(h ), independently of ɛ. The Miller scheme is thereby seen to be second-order accurate, uniformly in ɛ. O() 5 Mehrstellenverfahren By incorporating the principal term of the local truncation error given by Equation (5) into the Miller scheme, Equation (4), we obtain the formally fourth-order accurate scheme: ɛ z [cosh z ] ui+ u i + u i + c u h i = f i + h f i. (6) The inclusion of higher-order information concerning the source term f(x) is the hallmark of the Mehrstellenverfahren technique (Hermitian method) of L. Collatz []. The local truncation error of this scheme is given by L.T.E. = h c [ 4 ɛ u i }{{} O(/ɛ) + c 6 uiv i ɛ }{{} O(/ɛ ) 5 uvi i ] + O(h 6 ). (7) }{{} O(/ɛ 3 )
5 Constant coefficient reaction-diffusion problem of boundary layer type 67 Observe that, for fixed ɛ, the principal term of the local truncation error is now O(h 4 ) so that this Mehrstellenverfahren scheme is formally fourth-order accurate. However, with h ɛ in order to resolve the boundary layer, this leading term is actually O(h 4 /ɛ ) = O(h ). Thus, this Mehrstellenverfahren scheme is not uniformly fourth-order accurate in the small parameter ɛ. If evaluation of f (x) is either difficult or impossible, then the further approximation f i f i+ f i + f i h (8) may be made without loss of formal order of accuracy. 6 An Original Derivation of the Miller Scheme Motivated by the above observed correlation between the behavior of the local truncation error in the boundary layer and the order of uniform accuracy of the corresponding scheme, we next present a new derivation of the Miller scheme, Equation (4), for the case of constant c. We begin by considering the second central difference ɛ ui+ u i + u i h = ɛ n= hn (n)! u(n) i. (9) Keeping only the first term in the above sum would produce the central difference approximation. Although this would be formally second-order accurate, the neglected terms in the series have components that display excessive growth in the boundary layer. Specifically, { [ ɛ h n u (n) = ɛ h n ( ) ( ) ] } c c n ɛ f (n ) + f (n 4) + + f + ( cɛ ɛ ɛ )n u ( ) c n ( ) c n = [h n f (n ) + + h n f ] + h n (ɛu ), () ɛ ɛ }{{}}{{} O(h ) O() where we have used Equation () recursively. With h ɛ, detailed examination of each of the terms on the right hand side of Equation () reveals that only the last one is O(), the remaining terms being O(h ) or higher. Hence, if we want to achieve uniform second-order accuracy then we need to include all such terms when deriving an approximation from Equation (9). Thus, we make the boundary layer approximation ( ) c n u (n) u ɛ ()
6 68 Brian J. McCartin which, when inserted into Equation (9), yields the approximation Solving for u i u i+ u i + u i h n= and summing the infinite series yields u i ( ) hn c n (n)! u i. () ɛ z [cosh z ] ui+ u i + u i h, (3) which is precisely Miller s approximation. In [5], this markedly new derivation of Miller s scheme will provide the insight needed to generalize it, in a straightforward fashion, to the general case of variable c. However, we now show how to capitalize on this new perspective to achieve several enhancements to Miller s scheme for the case of constant c. 7 Fitted Mehrstellenverfahren For the case of constant c, the above alternative derivation of Miller s scheme hints at how to achieve uniform fourth-order accuracy without expanding the difference stencil. Simply include the penultimate term of Equation () when making the boundary layer approximation, Equation (). Any attempt to eke additional accuracy from Equation () in this fashion would indeed require an expansion of the difference stencil beyond three points due to the need to approximate fi iv or higher derivatives of the source term. The resulting modified Miller scheme for constant c is given by z u i + u i ɛ +c u [cosh z ] ui+ h i = f i + c ɛ [ The local truncation error of this scheme is given by z [cosh z ] ] f i. (4) L.T.E. = h4 36 [ f i iv ] + O(h 6 ). (5) }{{} Thus, the local truncation error is O(h 4 ), independently of ɛ. The Fitted Mehrstellenverfahren scheme is thereby seen to be fourth-order accurate, uniformly in ɛ! 8 The New Scheme in Equation (4), the further approxi- In order to obviate the need for f i mation f i O() f i+ f i + f i h (6)
7 Constant coefficient reaction-diffusion problem of boundary layer type 69 is invoked without any loss of either formal or uniform accuracy. The resulting new scheme z u i + u i ɛ [cosh z ] ui+ h +c u i = f i + c ɛ [ z ] f i + f i fi+ [cosh z ] h. (7) will be referred to as the McCartin scheme in what follows. The local truncation error of this scheme is given by L.T.E. = h4 4 [ f i iv ] + O(h 6 ). (8) }{{} Thus, the local truncation error is once again O(h 4 ), independently of ɛ. The McCartin scheme is likewise also seen to be fourth-order accurate, uniformly in ɛ! O() 9 The Blended Scheme Comparison of the principal parts of the local truncation errors of the Fitted Mehrstellenverfahren scheme, Equation (5), and the McCartin scheme, Equation (8), leads to the surprising conclusion that an appropriate linear combination of these two schemes will yield a blended scheme which is formally sixth-order accurate. Specifically, 3 5 L.T.E. F itted Mehrstellenverfahren + 5 L.T.E. McCartin = O(h 6 ). (9) The resulting new scheme z ɛ [cosh z ] u i+ u i + u i h +cu i = f i + ɛ c [ z ] 3h f i + (f i+ f i + f i ) [cosh z ] 5h. () will be referred to as the Blended McCartin scheme in what follows. Boundary Layer Test Problem We now present the results of some numerical experiments aimed at demonstrating the efficacy of the new schemes developed above. First, we present a model problem for which we know the exact solution. Then, this test problem is used to compare our schemes with those commonly used in the literature for the boundary layer case. These new schemes provide greatly improved results in all instances examined.
8 7 Brian J. McCartin ε =.3963 ε =.953 ε = ε =.4888 ε =.444 ε =.7 Figure : Birth of a Boundary Layer (solid=u(x; ɛ), dashed=u (x)) The constant coefficient two-point boundary value problem, ɛ u + c u = e αx sin (βx); u() = u L, u() = u R, () possesses the analytical solution determined by the method of undetermined coefficients c/ɛ x u(x; ɛ) = c e + c e c/ɛ x + A e αx sin (βx) + B e αx cos (βx), () where and A = c ɛ(α + β ) [c ɛ(α + β )] + [ɛαβ] ; B = ɛαβ [c ɛ(α + β )] + [ɛαβ], c = u R + (B u L ) e c/ɛ e α [A sin β + B cos β] sinh ( c/ɛ ) ; c = u L B c,
9 Constant coefficient reaction-diffusion problem of boundary layer type 7 while the reduced problem, corresponding to ɛ = and ignoring the boundary conditions, has solution u (x) = c eαx sin (βx). (3) Figure displays the development of the boundary layers as ɛ + for the case of c =, u L =, u R =, α =, β = π. The dashed curve denotes the solution to the reduced problem, Equation (3), while the solid curves correspond to the solution of the boundary value problem, Equation (), for ɛ = m (m = 8,..., 3). This test problem will now be used as the basis for comparison of the methods described above. In the maximum norm, we have [4, p. 57] u(x; ɛ) u h (x; ɛ) Ch p, (4) where h := /N is the mesh parameter, C is the asymptotic error constant and p is the asymptotic rate of convergence which is estimated as [4, p. 59] ( ) un/ u p p N := log. (5) u N u Figures 3 through 8 display, for each method, the exact (solid line) and numerical (circles) solutions on the left with the corresponding error (dots) on the right for ɛ = 3 and N = 3 where h = /N. Tables through 6 contain the corresponding maximum error and convergence rate estimated from Equation (5) for ɛ = m (m = 8,..., 3) and N = n (n = 5,..., 9). Table clearly shows the formal second-order accuracy of the Central Difference scheme as h for fixed ɛ. However, the nonuniform nature of this convergence is evidenced by the growth of the errors as ɛ for fixed h. This contrasts sharply with Table where the uniform second-order accuracy of the Miller scheme for h with fixed ɛ as well as ɛ with fixed h is prominently displayed. Likewise, the results of Table 3 clearly indicate the formal fourth-order accuracy of the Mehrstellenverfahren scheme while those of Table 4 testify to the uniform fourth-order accuracy of the Fitted Mehrstellenverfahren scheme. Table 5 reinforces the fact that the uniformity of the fourth-order convergence of the McCartin scheme is not sacrificed by making the approximation of Equation (6). Finally, the nonuniform sixth-order accuracy of the Blended McCartin scheme is in evidence in Table 6. The erratic results for N = 5 are almost certainly a numerical artifact due to finite-precision arithmetic which is unsurprising since they involve errors on the order of %!
10 7 Brian J. McCartin Conclusion In the preceding sections, new higher-order accurate schemes for the steady, one-dimensional, constant coefficient reaction-diffusion equation of boundary layer type have been presented. These included schemes of uniform fourthorder and formal sixth-order accuracy. Detailed derivations as well as extensive comparisons to existing schemes have been provided. The computational efficacy of the McCartin scheme is underscored by the fact that the maximum error for ɛ = 3 with N = is 3 6 which is less than half of the corresponding error of the Miller scheme with N = 5 thereby indicating that only N grid points need be used for comparable accuracy. These boundary layer test case results provide strong numerical evidence for the validity of the following: Conjecture (Equivalence of Uniform Consistency and Convergence) For f(x) sufficiently smooth and u(x) defined by L[u] := ɛ u (x) + c u(x) f(x) = ; u() = u L, u() = u R, k th -order uniform consistency is equivalent to k th -order uniform convergence. That is to say, if and only if where L h [u h ] = u u h < Ch k (C independent of ɛ) L h [u] = L[u] + C l h l, l=k lim C l = O(h k l ). h ɛ=h There are a number of directions in which this work could be extended. For example, the derivation for a non-uniform mesh would be straightforward but tedious. It is anticipated that this would result in the loss of one order of formal accuracy. Also, the fitted operators presented herein could be combined with corresponding fitted uniform meshes []. Finally, treatment of multiple spatial dimensions, transients, and nonlinearities would each be of substantial interest. For the extension of these methods to variable coefficient reaction-diffusion problems with boundary layers, see [5]. The extension to the treatment of the highly oscillatory case, where the diffusion and reaction terms are of the same sign, will appear in [6]. Acknowledgements. As always, the assistance of Barbara A. McCartin in the production of this manuscript is most gratefully acknowledged.
11 Constant coefficient reaction-diffusion problem of boundary layer type 73 References [] G. Birkhoff and S. Gulati, Optimal Few-Point Discretizations of Linear Source Problems, SIAM J. Numerical Analysis, (974), no. 4, [] L. Collatz, The Numerical Treatment of Differential Equations, Springer- Verlag, New York, NY, [3] E. P. Doolan, J. J. H. Miller and W. H. A. Schilders, Uniform Numerical Methods for Problems with Initial and Boundary Layers, Boole Press, Dublin, 98. [4] P. A. Farrell, A. F. Hegarty, J. J. H. Miller, E. O Riordan and G. I. Shishkin, Robust Computational Techniques for Boundary Layers, Chapman & Hall/CRC, Boca Raton, Florida,. [5] B. J. McCartin, Higher-Order Accurate Schemes for a Variable Coefficient Singularly Perturbed Reaction-Diffusion Problem of Boundary Layer Type, Applied Mathematical Sciences, (6), no. 4, [6] B. J. McCartin, Higher-Order Accurate Schemes for a Constant Coefficient Singularly Perturbed Reaction-Diffusion Problem of Highly Oscillatory Type, Applied Mathematical Sciences, (6), to appear. [7] J. J. H. Miller, On the Convergence, Uniformly in ɛ, of Difference Schemes for a Two Point Boundary Singular Perturbation Problem, Chapter in Numerical Analysis of Singular Perturbation Problems, P. W. Hemker and J. J. H. Miller (Eds.), Academic Press, New York, NY, 979, [8] J. J. H. Miller, E. O Riordan and G. I. Shishkin, Fitted Numerical Methods for Singular Perturbation Problems, Revised Edition, World Scientific, Singapore,. [9] K. W. Morton, Numerical Solution of Convection-Diffusion Problems, Chapman & Hall, London, 996. [] H.-G. Roos, M. Stynes and L. Tobiska, Numerical Methods for Singular Perturbation Problems, Springer-Verlag, New York, NY, 996. [] G. I. Shishkin and L. P. Shishkina, Difference Methods for Singular Perturbation Problems, Chapman & Hall/CRC, Boca Raton, Florida, 8.
12 74 Brian J. McCartin Central Difference Scheme.5 error Figure 3: Central Difference Scheme (ɛ = 3, N = 3) ɛ N = 3 N = 64 N = 8 N = 56 N = ( 3) 9.594( 4).3478( 4) 5.88( 5).45495( 5) ( 3).8663( 3) ( 4).84( 4).9533( 5) ( ) 3.736( 3) ( 4).37636( 4) ( 5) ( ) ( 3).88546( 3) ( 4).945( 4) ( ).4647( ) 3.749( 3) 9.5( 4).3879( 4) ( ).4833( ) ( 3).88993( 3) ( 4) Table : Central Difference Scheme: Error and Convergence Rate
13 Constant coefficient reaction-diffusion problem of boundary layer type 75 Miller Scheme 4 x 4 error Figure 4: Miller Scheme (ɛ = 3, N = 3) ɛ N = 3 N = 64 N = 8 N = 56 N = ( 3) ( 4).9537( 4).747( 5) 6.854( 6) ( 3) 4.59( 4).333( 4).8363( 5) 7.99( 6) ( 3) ( 4).4886( 4).8788( 5) 7.3( 6) ( 3) 4.546( 4).55( 4).8943( 5) 7.448( 6) ( 3) ( 4).485( 4).8976( 5) 7.697( 6) ( 3) 4.598( 4).3566( 4).8943( 5) ( 6) Table : Miller Scheme: Error and Convergence Rate
14 76 Brian J. McCartin Mehrstellenverfahren Scheme 6 x 4 error Figure 5: Mehrstellenverfahren Scheme (ɛ = 3, N = 3) ɛ N = 3 N = 64 N = 8 N = 56 N = ( 5).4368( 6) ( 8) ( 9) () ( 5).8697( 6).79578( 7).337( 8) 7.35() ( 5) ( 6) 3.678( 7).655( 8).465( 9) ( 4).46( 5) ( 7) ( 8).84( 9) ( 4).44( 5).4459( 6) 9.769( 8) ( 9) ( 4) ( 5).8563( 6).87( 7).365( 8) Table 3: Mehrstellenverfahren Scheme: Error and Convergence Rate
15 Constant coefficient reaction-diffusion problem of boundary layer type 77 Fitted Mehrstellenverfahren Scheme 5 x 7 error Figure 6: Fitted Mehrstellenverfahren Scheme (ɛ = 3, N = 3) ɛ N = 3 N = 64 N = 8 N = 56 N = ( 7) 3.948( 8).47486( 9).54789() () ( 7) ( 8).69646( 9).68793().6() ( 7) 4.56( 8).8489( 9).78558().93() ( 7) 4.594( 8).94( 9).84999().588() ( 7) ( 8).9855( 9).88887().845() ( 7) ( 8).9787( 9).98().4() Table 4: Fitted Mehrstellenverfahren Scheme: Error and Convergence Rate
16 78 Brian J. McCartin McCartin Scheme 4 x 7 error Figure 7: McCartin Scheme (ɛ = 3, N = 3) ɛ N = 3 N = 64 N = 8 N = 56 N = ( 7) 5.939( 8) ( 9).38().4538() ( 6) ( 8) ( 9).535().5775() ( 6) ( 8) 4.79( 9).67934().679() ( 6) ( 8) ( 9).7775().7365() ( 6) 6.999( 8) 4.57( 9).83783().77683() ( 7) ( 8) ( 9).874().8() Table 5: McCartin Scheme: Error and Convergence Rate
17 Constant coefficient reaction-diffusion problem of boundary layer type 79 Blended McCartin Scheme 3 x 8 error Figure 8: Blended McCartin Scheme (ɛ = 3, N = 3) ɛ N = 3 N = 64 N = 8 N = 56 N = ( 9) 3.756() 5.833(3).599(4).78746(4) ( 9) 8.6().83().38778(4) 9.484(4) ( 8).734().734() (4) 6.74(4) ( 8) () () (4).465(4) ( 8) ().546().888(3) (5) ( 8).35956( 9).3() (3) (5) Table 6: Blended McCartin Scheme: Error and Convergence Rate Received: January 3, 6; Published: March 3, 6
An ε-uniform Initial Value Technique For Convection-Diffusion Singularly Perturbed Problems
An -Uniform Initial Value Technique For Convection-Diffusion Singularly Perturbed Problems RKBawa and Vinod Kumar Abstract In this paper, we have proposed an -uniform initial value technique for singularly
More informationUniformly convergent scheme for Convection Diffusion problem
c Copyright, Darbose International Journal of Applied Mathematics and Computation Volume 3(4),pp 290 299, 2011 http://ijamc.psit.in Uniformly convergent scheme for Convection Diffusion problem K. Sharath
More informationRobust Numerical Methods for Singularly Perturbed Differential Equations SPIN Springer s internal project number, if known
Hans-Görg Roos Martin Stynes Lutz Tobiska Robust Numerical Methods for Singularly Perturbed Differential Equations SPIN Springer s internal project number, if known Convection-Diffusion-Reaction and Flow
More informationA Quintic Spline method for fourth order singularly perturbed Boundary Value Problem of Reaction-Diffusion type
A Quintic Spline method for fourth order singularly perturbed Boundary Value Problem of Reaction-Diffusion type Jigisha U. Pandya and Harish D. Doctor Asst. Prof., Department of Mathematics, Sarvajanik
More informationA Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete Data
Applied Mathematical Sciences, Vol. 1, 16, no. 7, 331-343 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/1.1988/ams.16.5177 A Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete
More informationAn Adaptive Stencil Linear Deviation Method for Wave Equations
211 An Adaptive Stencil Linear Deviation Method for Wave Equations Kelly Hasler Faculty Sponsor: Robert H. Hoar, Department of Mathematics ABSTRACT Wave Equations are partial differential equations (PDEs)
More informationHyperbola for Curvilinear Interpolation
Applied Mathematical Sciences, Vol. 7, 2013, no. 30, 1477-1481 HIKARI Ltd, www.m-hikari.com Hyperbola for Curvilinear Interpolation G. L. Silver 868 Kristi Lane Los Alamos, NM 87544, USA gsilver@aol.com
More informationVARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung
POLYTECHNIC UNIVERSITY Department of Computer and Information Science VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung Abstract: Techniques for reducing the variance in Monte Carlo
More informationFOURTH ORDER COMPACT FORMULATION OF STEADY NAVIER-STOKES EQUATIONS ON NON-UNIFORM GRIDS
International Journal of Mechanical Engineering and Technology (IJMET Volume 9 Issue 10 October 2018 pp. 179 189 Article ID: IJMET_09_10_11 Available online at http://www.iaeme.com/ijmet/issues.asp?jtypeijmet&vtype9&itype10
More informationNonparametric Regression
Nonparametric Regression John Fox Department of Sociology McMaster University 1280 Main Street West Hamilton, Ontario Canada L8S 4M4 jfox@mcmaster.ca February 2004 Abstract Nonparametric regression analysis
More informationAsymptotic Error Analysis
Asymptotic Error Analysis Brian Wetton Mathematics Department, UBC www.math.ubc.ca/ wetton PIMS YRC, June 3, 2014 Outline Overview Some History Romberg Integration Cubic Splines - Periodic Case More History:
More informationPROOF OF THE COLLATZ CONJECTURE KURMET SULTAN. Almaty, Kazakhstan. ORCID ACKNOWLEDGMENTS
PROOF OF THE COLLATZ CONJECTURE KURMET SULTAN Almaty, Kazakhstan E-mail: kurmet.sultan@gmail.com ORCID 0000-0002-7852-8994 ACKNOWLEDGMENTS 2 ABSTRACT This article contains a proof of the Collatz conjecture.
More informationThe Total Variation Approach to Approximate Hyperbolic Wave Equations
2 CUMMING The Total Variation Approach to Approximate Hyperbolic Wave Equations Christine Cumming Faculty Sponsor: Dr. Robert H. Hoar, Department of Math ABSTRACT From ultrasound for medical purposes to
More informationFinite element solution of multi-scale transport problems using the least squares based bubble function enrichment
Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment A. Yazdani a, V. Nassehi b1 a Cranfield University, School of Applied Sciences, Cranfield,
More informationSolving a Two Dimensional Unsteady-State. Flow Problem by Meshless Method
Applied Mathematical Sciences, Vol. 7, 203, no. 49, 242-2428 HIKARI Ltd, www.m-hikari.com Solving a Two Dimensional Unsteady-State Flow Problem by Meshless Method A. Koomsubsiri * and D. Sukawat Department
More informationConditional Volatility Estimation by. Conditional Quantile Autoregression
International Journal of Mathematical Analysis Vol. 8, 2014, no. 41, 2033-2046 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.47210 Conditional Volatility Estimation by Conditional Quantile
More informationVertex Graceful Labeling of C j C k C l
Applied Mathematical Sciences, Vol. 8, 01, no. 8, 07-05 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ams.01.5331 Vertex Graceful Labeling of C j C k C l P. Selvaraju 1, P. Balaganesan,5, J. Renuka
More informationFinite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2
Finite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2 David W. Trott and Matthias K. Gobbert Department of Mathematics and Statistics, University of Maryland, Baltimore County,
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationAn Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)
An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationA New Approach to Meusnier s Theorem in Game Theory
Applied Mathematical Sciences, Vol. 11, 2017, no. 64, 3163-3170 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.712352 A New Approach to Meusnier s Theorem in Game Theory Senay Baydas Yuzuncu
More informationA Singular Example for the Averaged Mean Curvature Flow
To appear in Experimental Mathematics Preprint Vol. No. () pp. 3 7 February 9, A Singular Example for the Averaged Mean Curvature Flow Uwe F. Mayer Abstract An embedded curve is presented which under numerical
More informationCS321 Introduction To Numerical Methods
CS3 Introduction To Numerical Methods Fuhua (Frank) Cheng Department of Computer Science University of Kentucky Lexington KY 456-46 - - Table of Contents Errors and Number Representations 3 Error Types
More informationNumerical and theoretical analysis of shock waves interaction and reflection
Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for
More informationExtended Goursat s Hypergeometric Function
International Journal of Mathematical Analysis Vol. 9, 25, no. 3, 59-57 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.2988/ijma.25.532 Extended Goursat s Hypergeometric Function Daya K. Nagar and Raúl
More informationDomination Number of Jump Graph
International Mathematical Forum, Vol. 8, 013, no. 16, 753-758 HIKARI Ltd, www.m-hikari.com Domination Number of Jump Graph Y. B. Maralabhavi Department of Mathematics Bangalore University Bangalore-560001,
More informationComparison of different solvers for two-dimensional steady heat conduction equation ME 412 Project 2
Comparison of different solvers for two-dimensional steady heat conduction equation ME 412 Project 2 Jingwei Zhu March 19, 2014 Instructor: Surya Pratap Vanka 1 Project Description The purpose of this
More informationWhat is the Optimal Bin Size of a Histogram: An Informal Description
International Mathematical Forum, Vol 12, 2017, no 15, 731-736 HIKARI Ltd, wwwm-hikaricom https://doiorg/1012988/imf20177757 What is the Optimal Bin Size of a Histogram: An Informal Description Afshin
More informationA Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods
A Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods Scott A. Sarra, Derek Sturgill Marshall University, Department of Mathematics, One John Marshall Drive, Huntington
More informationGraph Sampling Approach for Reducing. Computational Complexity of. Large-Scale Social Network
Journal of Innovative Technology and Education, Vol. 3, 216, no. 1, 131-137 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/jite.216.6828 Graph Sampling Approach for Reducing Computational Complexity
More informationComputing Nearly Singular Solutions Using Pseudo-Spectral Methods
Computing Nearly Singular Solutions Using Pseudo-Spectral Methods Thomas Y. Hou Ruo Li January 9, 2007 Abstract In this paper, we investigate the performance of pseudo-spectral methods in computing nearly
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationA-posteriori Diffusion Analysis of Numerical Schemes in Wavenumber Domain
2th Annual CFD Symposium, August 9-1, 218, Bangalore A-posteriori Diffusion Analysis of Numerical Schemes in Wavenumber Domain S. M. Joshi & A. Chatterjee Department of Aerospace Engineering Indian Institute
More informationInclusion of Aleatory and Epistemic Uncertainty in Design Optimization
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Inclusion of Aleatory and Epistemic Uncertainty in Design Optimization Sirisha Rangavajhala
More informationThe Number of Fuzzy Subgroups of Cuboid Group
International Journal of Algebra, Vol. 9, 2015, no. 12, 521-526 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ija.2015.5958 The Number of Fuzzy Subgroups of Cuboid Group Raden Sulaiman Department
More informationSolve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method
Solve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method P.B. Choksi 1 and A.K. Pathak 2 1 Research Scholar, Rai University,Ahemdabad. Email:pinalchoksey@gmail.com 2 H.O.D.
More informationThe WENO Method in the Context of Earlier Methods To approximate, in a physically correct way, [3] the solution to a conservation law of the form u t
An implicit WENO scheme for steady-state computation of scalar hyperbolic equations Sigal Gottlieb Mathematics Department University of Massachusetts at Dartmouth 85 Old Westport Road North Dartmouth,
More informationISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 2, Issue 3, September 2012
Mitigation Curves for Determination of Relief Holes to Mitigate Concentration Factor in Thin Plates Loaded Axially for Different Discontinuities Shubhrata Nagpal, S.Sanyal, Nitin Jain Abstract In many
More informationA Moving Mesh Method for Time dependent Problems Based on Schwarz Waveform Relaxation
A Moving Mesh Method for Time dependent Problems Based on Schwarz Waveform Relaation Ronald D. Haynes, Weizhang Huang 2, and Robert D. Russell 3 Acadia University, Wolfville, N.S., Canada. ronald.haynes@acadiau.ca
More informationSolutions of Stochastic Coalitional Games
Applied Mathematical Sciences, Vol. 8, 2014, no. 169, 8443-8450 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.410881 Solutions of Stochastic Coalitional Games Xeniya Grigorieva St.Petersburg
More informationChapter 1 BACKGROUND
Chapter BACKGROUND. Introduction In many areas of mathematics and in applications of mathematics, it is often necessary to be able to infer information about some function using only a limited set of sample
More informationDiscrete Optimization. Lecture Notes 2
Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationChap.12 Kernel methods [Book, Chap.7]
Chap.12 Kernel methods [Book, Chap.7] Neural network methods became popular in the mid to late 1980s, but by the mid to late 1990s, kernel methods have also become popular in machine learning. The first
More information1 2 (3 + x 3) x 2 = 1 3 (3 + x 1 2x 3 ) 1. 3 ( 1 x 2) (3 + x(0) 3 ) = 1 2 (3 + 0) = 3. 2 (3 + x(0) 1 2x (0) ( ) = 1 ( 1 x(0) 2 ) = 1 3 ) = 1 3
6 Iterative Solvers Lab Objective: Many real-world problems of the form Ax = b have tens of thousands of parameters Solving such systems with Gaussian elimination or matrix factorizations could require
More informationA Moving Mesh Method for Time Dependent Problems based on Schwarz Waveform Relaxation
A Moving Mesh Method for Time Dependent Problems based on Schwarz Waveform Relaation Ronald D. Haynes, Weizhang Huang 2, and Robert D. Russell 3 Acadia University, Wolfville, N.S., Canada ronald.haynes@acadiau.ca
More informationSome Coupled Fixed Point Theorems on Quasi-Partial b-metric Spaces
International Journal of Mathematical Analysis Vol. 9, 2015, no. 6, 293-306 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.412388 Some Coupled Fixed Point Theorems on Quasi-Partial b-metric
More informationLocal nonobtuse tetrahedral refinements around an edge
Local nonobtuse tetrahedral refinements around an edge Sergey Korotov, Michal Křížek September, 00 8 9 0 8 9 0 M - asque enter for pplied Mathematics izkaia Technology Park, uilding 00 E 80, Derio, asque
More informationOn the Number of Tilings of a Square by Rectangles
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange University of Tennessee Honors Thesis Projects University of Tennessee Honors Program 5-2012 On the Number of Tilings
More informationDOWNLOAD PDF BIG IDEAS MATH VERTICAL SHRINK OF A PARABOLA
Chapter 1 : BioMath: Transformation of Graphs Use the results in part (a) to identify the vertex of the parabola. c. Find a vertical line on your graph paper so that when you fold the paper, the left portion
More informationNumerical Rectification of Curves
Applied Mathematical Sciences, Vol. 8, 2014, no. 17, 823-828 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.39500 Numerical Rectification of Curves B. P. Acharya, M. Acharya and S. B.
More informationc 2013 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 35, No. 5, pp. A2225 A2254 c 2013 Society for Industrial and Applied Mathematics ROBUST SOLUTION OF SINGULARLY PERTURBED PROBLEMS USING MULTIGRID METHODS SCOTT MACLACHLAN AND
More informationEnnumeration of the Number of Spanning Trees in the Lantern Maximal Planar Graph
Applied Mathematical Sciences, Vol. 8, 2014, no. 74, 3661-3666 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.44312 Ennumeration of the Number of Spanning Trees in the Lantern Maximal
More informationHarmonic Spline Series Representation of Scaling Functions
Harmonic Spline Series Representation of Scaling Functions Thierry Blu and Michael Unser Biomedical Imaging Group, STI/BIO-E, BM 4.34 Swiss Federal Institute of Technology, Lausanne CH-5 Lausanne-EPFL,
More informationBuilding the Graphics Memory of. the Stiffness Matrix of the Beam
Contemporary Engineering Sciences, Vol. 11, 2018, no. 92, 4593-4605 HIKARI td, www.m-hikari.com https://doi.org/10.12988/ces.2018.89502 Building the Graphics Memory of the Stiffness Matrix of the Beam
More informationApproximation of 3D-Parametric Functions by Bicubic B-spline Functions
International Journal of Mathematical Modelling & Computations Vol. 02, No. 03, 2012, 211-220 Approximation of 3D-Parametric Functions by Bicubic B-spline Functions M. Amirfakhrian a, a Department of Mathematics,
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationA Comparative Study of Data Mining Process Models (KDD, CRISP-DM and SEMMA)
International Journal of Innovation and Scientific Research ISSN 2351-8014 Vol. 12 No. 1 Nov. 2014, pp. 217-222 2014 Innovative Space of Scientific Research Journals http://www.ijisr.issr-journals.org/
More informationMATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES
MATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES Goals of research Dr. Uldis Raitums, Dr. Kārlis Birģelis To develop and investigate mathematical properties of algorithms
More informationNew Mexico Tech Hyd 510
Numerics Motivation Modeling process (JLW) To construct a model we assemble and synthesize data and other information to formulate a conceptual model of the situation. The model is conditioned on the science
More informationEdge Detection Free Postprocessing for Pseudospectral Approximations
Edge Detection Free Postprocessing for Pseudospectral Approximations Scott A. Sarra March 4, 29 Abstract Pseudospectral Methods based on global polynomial approximation yield exponential accuracy when
More informationHOUGH TRANSFORM CS 6350 C V
HOUGH TRANSFORM CS 6350 C V HOUGH TRANSFORM The problem: Given a set of points in 2-D, find if a sub-set of these points, fall on a LINE. Hough Transform One powerful global method for detecting edges
More informationIntroduction to Computational Mathematics
Introduction to Computational Mathematics Introduction Computational Mathematics: Concerned with the design, analysis, and implementation of algorithms for the numerical solution of problems that have
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationModeling with Uncertainty Interval Computations Using Fuzzy Sets
Modeling with Uncertainty Interval Computations Using Fuzzy Sets J. Honda, R. Tankelevich Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, CO, U.S.A. Abstract A new method
More informationHyperbolic structures and triangulations
CHAPTER Hyperbolic structures and triangulations In chapter 3, we learned that hyperbolic structures lead to developing maps and holonomy, and that the developing map is a covering map if and only if the
More informationA SEMI-ANALYTIC SPECTRAL METHOD FOR ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 43, pp. 1 11. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu A SEMI-ANALYTIC SPECTRAL METHOD FOR ELLIPTIC
More informationTruncation Errors. Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4.
Chapter 4: Roundoff and Truncation Errors Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4. 1 Outline Errors Accuracy and Precision
More informationReview Initial Value Problems Euler s Method Summary
THE EULER METHOD P.V. Johnson School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 INITIAL VALUE PROBLEMS The Problem Posing a Problem 3 EULER S METHOD Method Errors 4 SUMMARY OUTLINE 1 REVIEW 2 INITIAL
More informationFault Tolerant Domain Decomposition for Parabolic Problems
Fault Tolerant Domain Decomposition for Parabolic Problems Marc Garbey and Hatem Ltaief Department of Computer Science, University of Houston, Houston, TX 77204 USA garbey@cs.uh.edu, ltaief@cs.uh.edu 1
More informationDiffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012
Diffuse Optical Tomography, Inverse Problems, and Optimization Mary Katherine Huffman Undergraduate Research Fall 11 Spring 12 1. Introduction. This paper discusses research conducted in order to investigate
More informationEditorial Coupled Numerical Methods in Engineering Analysis
Mathematical Problems in Engineering Volume 2011, Article ID 436978, 4 pages doi:10.1155/2011/436978 Editorial Coupled Numerical Methods in Engineering Analysis Delfim Soares Jr., 1 Otto von Estorff, 2
More informationProximal Manifold Learning via Descriptive Neighbourhood Selection
Applied Mathematical Sciences, Vol. 8, 2014, no. 71, 3513-3517 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.42111 Proximal Manifold Learning via Descriptive Neighbourhood Selection
More informationConstructing Numerical Reference Modes from Sparse Time Series
Constructing Numerical Reference Modes from Sparse Series Peter S. Hovmand, Ph.D. Washington University George Warren Brown School of Social Work Campus Box 1196 One Brookings Drive St. Louis, MO 63130
More informationPARALLEL METHODS FOR SOLVING PARTIAL DIFFERENTIAL EQUATIONS. Ioana Chiorean
5 Kragujevac J. Math. 25 (2003) 5 18. PARALLEL METHODS FOR SOLVING PARTIAL DIFFERENTIAL EQUATIONS Ioana Chiorean Babeş-Bolyai University, Department of Mathematics, Cluj-Napoca, Romania (Received May 28,
More informationDisconnection Probability of Graph on Two Dimensional Manifold: Continuation
Applied Mathematical Sciences, Vol. 10, 2016, no. 40, 2003-2011 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.63123 Disconnection Probability of Graph on Two Dimensional Manifold: Continuation
More informationCollocation and optimization initialization
Boundary Elements and Other Mesh Reduction Methods XXXVII 55 Collocation and optimization initialization E. J. Kansa 1 & L. Ling 2 1 Convergent Solutions, USA 2 Hong Kong Baptist University, Hong Kong
More informationHigh Order Weighted Essentially Non-Oscillatory Schemes for Convection. Dominated Problems. Chi-Wang Shu 1
High Order Weighted Essentially Non-Oscillatory Schemes for Convection Dominated Problems Chi-Wang Shu Division of Applied Mathematics, Brown University, Providence, Rhode Island 09 ABSTRACT High order
More informationThe Method of Approximate Particular Solutions for Solving Elliptic Problems with Variable Coefficients
The Method of Approximate Particular Solutions for Solving Elliptic Problems with Variable Coefficients C.S. Chen, C.M. Fan, P.H. Wen Abstract A new version of the method of approximate particular solutions
More informationImplementation of the Continuous-Discontinuous Galerkin Finite Element Method
Implementation of the Continuous-Discontinuous Galerkin Finite Element Method Andrea Cangiani, John Chapman, Emmanuil Georgoulis and Max Jensen Abstract For the stationary advection-diffusion problem the
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationOn the order of accuracy and numerical performance of two classes of finite volume WENO schemes
On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationTable : IEEE Single Format ± a a 2 a 3 :::a 8 b b 2 b 3 :::b 23 If exponent bitstring a :::a 8 is Then numerical value represented is ( ) 2 = (
Floating Point Numbers in Java by Michael L. Overton Virtually all modern computers follow the IEEE 2 floating point standard in their representation of floating point numbers. The Java programming language
More informationShape Control of Cubic H-Bézier Curve by Moving Control Point
Journal of Information & Computational Science 4: 2 (2007) 871 878 Available at http://www.joics.com Shape Control of Cubic H-Bézier Curve by Moving Control Point Hongyan Zhao a,b, Guojin Wang a,b, a Department
More informationSome Algebraic (n, n)-secret Image Sharing Schemes
Applied Mathematical Sciences, Vol. 11, 2017, no. 56, 2807-2815 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.710309 Some Algebraic (n, n)-secret Image Sharing Schemes Selda Çalkavur Mathematics
More informationMODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
A. N. Johnson et al., Int. J. Comp. Meth. and Exp. Meas., Vol. 3, No. 3 (2015) 269 278 MODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
More informationGRADED MESHES FOR HIGHER ORDER FEM *
Journal of Computational Mathematics Vol.33, No.1, 2015, 1 16. http://www.global-sci.org/jcm doi:10.4208/jcm.1405-m4362 GRADED MESHES FOR HIGHER ORDER FEM * Hans-Görg Roos Institute of Numerical Mathematics,
More informationGraceful Labeling for Some Star Related Graphs
International Mathematical Forum, Vol. 9, 2014, no. 26, 1289-1293 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.4477 Graceful Labeling for Some Star Related Graphs V. J. Kaneria, M.
More informationCalibration of Nonlinear Viscoelastic Materials in Abaqus Using the Adaptive Quasi-Linear Viscoelastic Model
Calibration of Nonlinear Viscoelastic Materials in Abaqus Using the Adaptive Quasi-Linear Viscoelastic Model David B. Smith *, Uday Komaragiri **, and Romil Tanov ** ** * Ethicon Endo-Surgery, Inc., Cincinnati,
More informationarxiv: v1 [math.na] 26 Jun 2014
for spectrally accurate wave propagation Vladimir Druskin, Alexander V. Mamonov and Mikhail Zaslavsky, Schlumberger arxiv:406.6923v [math.na] 26 Jun 204 SUMMARY We develop a method for numerical time-domain
More informationResearch Article Data Visualization Using Rational Trigonometric Spline
Applied Mathematics Volume Article ID 97 pages http://dx.doi.org/.//97 Research Article Data Visualization Using Rational Trigonometric Spline Uzma Bashir and Jamaludin Md. Ali School of Mathematical Sciences
More informationSpline Solution Of Some Linear Boundary Value Problems
Applied Mathematics E-Notes, 8(2008, 171-178 c ISSN 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Spline Solution Of Some Linear Boundary Value Problems Abdellah Lamnii,
More informationVW 1LQH :HHNV 7KH VWXGHQW LV H[SHFWHG WR
PreAP Pre Calculus solve problems from physical situations using trigonometry, including the use of Law of Sines, Law of Cosines, and area formulas and incorporate radian measure where needed.[3e] What
More informationTravel-Time Sensitivity Kernels in Long-Range Propagation
Travel-Time Sensitivity Kernels in Long-Range Propagation Emmanuel Skarsoulis Foundation for Research and Technology Hellas Institute of Applied and Computational Mathematics P.O. Box 1385, GR-71110 Heraklion,
More informationModeling Ground Water Problems Using the Complex Polynomial Method
Modeling Ground Water Problems Using the Complex Polynomial Method A. W. Bohannon and T. V. Hromadka, AS-0020 Abstract Numerical methods for solving the governing partial differential equations involved
More informationLacunary Interpolation Using Quartic B-Spline
General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan
More information