A Meshless Simulations for 2D Nonlinear Reaction-diffusion Brusselator System
|
|
- Derick Waters
- 5 years ago
- Views:
Transcription
1 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 A Meshless Simulations for D Nonlinear Reaction-diffusion Brusselator Sstem Ahmad Shirzadi, Vladimir Sladek, Jan Sladek 3 Abstract: This paper is concerned with the development of a numerical approach based on the Meshless Local Petrov-Galerkin (MLPG method for the approimate solutions of the two dimensional nonlinear reaction-diffusion Brusselator sstems. The method uses finite differences for discretizing the time variable and the moving least squares (MLS approimation for field variables. The application of the weak formulation with the Heaviside tpe test functions supported on local subdomains (around the nodes used in MLS approimation to semi-discretized partial differential equations ields the finite-volume local weak formulation. A predictor-corrector scheme is used to handle the nonlinearit of the problem within each time step. Numerical test problems are given to verif the accurac of the proposed method. Under particular conditions, this sstem ehibits Turing instabilit which results in a pattern forming instabilit. This concept is studied and a test problem is given. Kewords: Meshless Local Petrov-Galerkin method; Moving least squares; Finite differences; Nonlinear problems; Brusselator equations; Turing instabilit; Pattern formation. Introduction In 95 Turing showed [Turing (95] that a simple mathematical model describing self-diffusing and reacting chemicals could give rise to stationar spatial concentration patterns of fied characteristic length from a random initial configuration and proposed that reaction-diffusion models might have relevance in describing morphogenesis, the growth of biological form. In the late 96s, IIa Department of Mathematics, Persian Gulf Universit, Bushehr, Iran. shirzadi@pgu.ac.ir, shirzadi.a@gmail.com Institute of Construction and Architecture, Slovak Academ of Sciences, 8453 Bratislava, Slovakia. vladimir.sladek@savba.sk 3 Institute of Construction and Architecture, Slovak Academ of Sciences, 8453 Bratislava, Slovakia. jan.sladek@savba.sk
2 6 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 Prigogine among others at the Universit of Brussels developed the Brusselator model [Prigogine and Lefever (968;Nicolis and Prigogine (997] which is one of the simplest reaction-diffusion models ehibiting Turing instabilit. From the mathematical point of view, the Brusselator is a sstem of two coupled nonlinear partial differential equations. The linear stabilit analsis can be used for predicting the parameter values that result in the Turing instabilit in a particular reaction-diffusion sstem. Based on the linear approimation one can also effectivel predict the characteristic length of the resulting pattern. Linear analsis cannot, however, predict the spatial characteristics of the resulting patterns, since the pattern selection is governed b comple nonlinear dnamics. Thus some nonlinear analsis is needed. Since analtical solutions are unavailable, finding efficient numerical methods for solving these kind of problems has been an active area of research. Meshless methods as a tool for numerical simulations of linear or nonlinear problems, especiall partial differential equations have become popular in recent ears. This paper aims to present a new meshless method for the numerical simulations of the Brusselator equations. The proposed method is based on the MLPG formulation [Atluri and Zhu (998a;Atluri and Zhu (998b;Atluri, Kim, and Cho (999;Atluri (4;Atluri and Shen (] which uses the local weak form equations and the MLS approimation to transform the model equations to the final sstem of algebraic equations. An application of the MLPG method for numericall solving coupled pair of nonlinear diffusion equations can be found in [Abbasband, Sladek, Shirzadi, and Sladek (]. For a review of applications of the MLPG method in engineering and sciences the readers are referred to [Sladek, Stanak, Han, Sladek, and Atluri (3]. Analsis of the convergence of the MLPG method is well studied in [Shirzadi and Ling (3]. The non-linear reaction-diffusion Brusselator sstem is given b: u t v t ( = A + u u v (B + u + D u + u ( = Bu u v v + D v + v, in the two-dimensional region Ω = [a,b], with initial condition, ( u(,, = f (,, ( v(,, = g(,, and Dirichlet boundar condition: u(,t = q (,t, Ω t, t >, v(,t = q (,t,
3 Meshless Simulations for Brusselator Sstem 6 or subject to Neumann s boundar conditions on the boundar Ω. Recall that A denotes the value of the source term and B is the bifurcation parameter setting the distance to the onset of instabilit. The chemical u is the activator concentration and v is the inhibitor. D u and D v are diffusion coefficients, f, g, q and q are known functions. Some eisting papers on the stud of this problem were concerned with theoretical concepts such as smmetr-breaking and bifurcations [Nicolis and Prigogine (997;Walgraef, Dewel, and Borckmans (98]. The theor b Turing [Turing (95] predicts that a spatiall uniform stationar state that is stable against perturbations in the absence of diffusion might become unstable against perturbations in the presence of diffusion resulting in a pattern forming instabilit. Some papers on numerical studies of this model emploed the computational approach and addressed the problem of pattern selection as a function of sstem parameters in two-dimensional and three-dimensional sstems [Verdasca, Wit, Dewel, and Borckmans (99;Borckmans, Wit, and Dewel (99;Leppänen (4]. For papers which deal with the numerical solution of the reactiondiffusion Brusselator sstem see Twizell et al. [Twizell, Gumel, and Cao (999] for a second order finite difference method, Adomian [Adomian (995] for decomposition method, Whe-Teong [Ang (3] for the dual-reciprocit boundar element method, S-u. Islam et al. [Siraj-ul, Arshed, and Sirajul (] for a numerical scheme based on RBF collocation method, Mittal and Jiwari [Mittal and Jiwari (] for polnomial based differential quadrature method, Kamranian et al. [Kamranian, Tatari, and Dehghan (] for the finite point method and so on. The remaining content of the paper is organized as follows. Section describes the basics of the MLS approimation used in this stud. The time discretization approach is presented in Section 3. Construction of local weak equations are illustrated in Section 4 and the obtained local weak equations are discretized in Section 5. The predictor-corrector scheme which is proposed to handle the nonlinear terms is also illustrated in Section 5. To simplif the evaluation of local integrals, certain techniques are described in Section 6. Two test problems are given in Section 7 to verif the accurac of the method. To address the problem of pattern selection as a function of sstem parameters, another test problem is given at the end of Section 7. Our conclusions are summarized in the last section. The MLS approimation Consider a sub-domain Ω, the neighborhood of a point which is located in the problem domain Ω. To approimate the distribution of function u in Ω over a number of randoml located nodes i,i =,,...n, the MLS approimant u h ( of u, Ω, is defined b u h ( = p T ( a( Ω, (3
4 6 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 where p T ( = [p (, p (,..., p m (] is a complete monomial basis of order m, and a( is a vector containing coefficients a j (, j =,,...m which are functions of the space coordinates. For eample, for a -D problem, p T ( = [,,] and p T ( = [,,,,, ], for linear basis (m = 3 and quadratic basis (m = 6, respectivel. The coefficient vector a( is determined b minimizing a weighted discrete L norm, defined as J( = n i= w i ([p T ( i a( û i ] (4 = [P.a( û] T.W.[P.a( û], where w i ( is the weight function associated with the node i, with w i ( > for all in the support of w i (, i denotes the value of at node i, n is the number of nodes in Ω for which the weight functions w i ( >, the matrices P and W are defined as p T ( P = p T (... p T ( n n m w (..., W = w n ( and û T = [û,û,...,û n ] where û i, i =,,...,n in (4 are the fictitious nodal values, and not the nodal values of the unknown trial function u h ( in general. In this work the Gaussian weight function is used: w i ( = ep[ ( d i c i ] ep[ ( r i ep[ ( r i c i ] c i ], d i r i,, d i r i, where d i = i, c i is a constant controlling the shape of the weight function w i and r i is the size of the support domain. The stationarit of J in (4 with respect to a( leads to the following linear relation between a( and û, A(a( = B(û, (5 where the matrices A( and B( are defined b A( = P T WP = B(P = n i= w i (p( i p T ( i, B( = P T W = [w (p(,w (p(,...,w n (p( n ].
5 Meshless Simulations for Brusselator Sstem 63 The MLS approimation is well defined if and onl if the rank of P equals m. A necessar condition for a well-defined MLS approimation is that at least m weight functions are non-zero (i.e. n > m for each sample point Ω and that the nodes in Ω will not be arranged in a special pattern such as on a straight line. Solving for a( from (5 and substituting it into (3 gives a relation which ma be written as the form of an interpolation function similar to that used in FEM, as u h ( = Φ T (.û = n i= and essentiall u i û i, where Φ T ( = p T (A (B(, or φ i ( = m p j ([A (B(] ji. j= φ i (û i ; u h ( i u i ; Ω, (6 φ i ( is called the shape function of the MLS approimation corresponding to nodal point i. The support of the nodal point i is usuall taken to be a circle of radius r i, centered at i. The fact that φ i ( =, for not in the support of nodal point i preserves the local character of the Moving Least Squares approimation. The partial derivatives of φ i ( are obtained as φ i,k = m j= [p j,k (A B ji + p j (A B,k + A,k B ji], (7 in which A,k = (A,k represents the derivative of the inverse of A with respect to k, kth coordinate of, which is given b A,k = A A,k A, where (,i denotes ( / i. 3 The finite difference approimation The finite-difference approimation of the time derivatives in the θ method is given as follows θ u k+ + ( θ u k = u k+ u k, θ. (8 t Considering Eq. ( at the time instants k t and (k + t, one obtains, respectivel θ u k+ = θd u u k+ + θa + θ(u k+ v k+ θ(b + u k+ ( θ u k = ( θd u u k + ( θa + ( θ(u k v k ( θ(b + u k
6 64 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 Hence and from (8, we have u k+ u k t ( = D u u k + θd u u k+ u k + A + θ(u k+ v k+ (9 + ( θ(u k v k (B + (θu k+ + ( θu k Using Crank-Nicholson scheme (θ = Eq. (9 becomes: u k+ u k t = D u ( u k + u k+ + A + ((u k+ v k+ + (u k v k ( (B + u k+ + u k or [ + t (B + D t u ] u k+ t (u k+ v k+ [ = t (B + + D t u ] u k + t (u k v k + ta. Similarl, we have [ t D v ] v k+ t Bu k+ + t (u k+ v k+ [ t = + D v ] v k + t Bu k t ( (u k v k Therefore, the time dependent PDEs are transformed to the semi-discrete PDEs of the elliptic tpe for the field variables u k+ and v k+, assuming the fields u k and v k being known from the computation in the previous time step. To treat the nonlinearit, a predictor-corrector scheme is proposed in this paper. 4 The local weak form In this paper, we construct the weak form over local sub-domains such as, located entirel inside Ω which is a circle of radius r and centered at node i. The local weak form of the equations ( and ( for i = ( i, i can be written as [[ + t t ] u k+ t (u k+ v k+ ] u d = (B + D u [[ t (B + + D t u ] u k + t ] (u k v k + ta u d, ( (
7 Meshless Simulations for Brusselator Sstem 65 [[ t D v ] v k+ t Bu k+ + t ] (u k+ v k+ u d = [[ t + D v ] v k + t Bu k t ] (u k v k u d where u and v are trial functions and u is a test function. If the Heaviside step function,, u ( =, /, (3 is chosen as the test function in each sub-domain, then u = and the local weak forms ( and (3 are transformed into the following equations ( + t t + D u t (B + u k+ d D u t (u k+ v k+ d = u k n ds + t ( t u k+ n ds (B + u k d ( (u k v k + A d, (4 + t t v k+ d D v t v k+ n ds ((u k+ v k+ Bu k+ d = ((u k v k Bu k d v k d + D v t v k n ds (5 5 Discretized equations and the predictor-corrector scheme We consider N nodal points, some of them located on the boundaries and some in the domain. Let u k ( = u(,k t be approimated according to (6 with denoting the nodal unknowns as û k,i. In each time step, our aim is to compute û k+,i, i =,,...,N assuming û k,i, i =,,...,N being known from the previous time step. Since there are nonlinear terms in the equations, we propose a predictor-corrector scheme in each time step to obtain a desired accurac. So we assume that û (l k+ is an approimation of û k+ after l correction step and û k+ is identified with the value û ( f k+ obtained in the final correction. We will begin from û( k+ = û k in nonlinear term. For internal nodes, from (4, (5 and using the MLS approimation (6, we
8 66 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 have the following linear equations: ( + t t + D u t (B + N φ j d û j=( (l+ k+, j D u t (ũ (l k+ ( ṽ (l k+ (d = ( t (B + N ( φ j j= n ds û (l+ k+, j N j=( φ j d û k, j N ( φ ( j j= n ds û k, j + t (ũ k( ṽ k ( + A d, N ( φ j d ˆv (l+ k+, j D v t j= + t N ( = j= t where ( (ũ (l k+ ( ṽ (l k+ ( d B t φ j d ˆv k, j + D v t ( (ũ k ( ṽ k ( Bũ k ( N ( φ j j= n ds ˆv (l+ k+, j N ( φ j j= d, N ( φ j (d û (l+ k+, j j= n ds ˆv k, j (6 (7 ũ k ( = N φ j (û k, j = j= N φ j (û ( f k, j, ṽ k( = j= N φ j ( ˆv k, j = j= N φ j ( ˆv ( f k, j, j= ũ (l N k+ = φ j (û (l k+, j, j= N ṽ(l k+ = φ j ( ˆv (l k+, j. j= The Dirichlet boundar conditions are imposed directl using the equations of the boundar conditions and the MLS approimation (6. To impose the Neumann s boundar conditions, we adopt a method proposed in Abbasband and Shirzadi (. 6 Computational techniques The MLS approimation with Gaussian weight function and quadratic basis (m = 6 are used in all test problems. We have to deal with two kinds of integrals in the discretized local integral equations (6 and (7. Basicall, both the integrations over the circular sub-domain and/or its boundar can be performed analticallsladek and Sladek (; Sladek, Sladek, and Zhang ( with saving the
9 Meshless Simulations for Brusselator Sstem 67 computational time substantiall because of eliminating the evaluation of the shape functions and/or its derivatives at the integration points. Then, we could write g(d = πr g( i + π 8 r4 g,kk ( i + O(r 6, g(ds = πr g( i + π r3 g,kk ( i + O(r 5, g n (ds = πr g,kk ( i + O(r 4. In order to avoid the inaccurac of the second and higher-order derivatives of the shape functions in meshless approimations, we prefer a numerical integration over. Since the number of the integration points in the domain integrals is much higher than in the contour integrals, the analtical integration in case of domain integrals ields more substantial savings in the computational time. Therefore we appl the analtical integrations over with omitting the terms O(r 4 involving the derivatives of the shape functions. Thus summarizing, the required integrals over are computed b using the regular Gauss-Legendre quadrature rule with seven integration points as π = φ j (ds φ j ( i + r cos(θ, i + r sin(θr dθ = πr φ j ( i + r cos(πθ + π, i + r sin(πθ + πdθ πr 7 p=w p φ j ( i + r cos(πθ p + π, i + r sin(πθ p + π where w p and θ p are the Gauss quadrature integration rule weights and points on [-, ]. The same technique can be used for evaluation of contour integrals of the normal derivatives of the shape functions n k (φ i,k ( with utilizing Eq. (7 for the gradients of the shape functions. Finall, the required integrals over sub-domains are approimated as follows: φ j (d πφ j ( i r.
10 68 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 and similarl, the domain integrals involving the nonlinear terms are approimated as: (ũ (l k+ ṽ (l k+d π(ũ(l k+ (i ṽ (l k+ (i r, ( (ũ k ṽ k + A d π( (ũ k( i ṽ k ( i + A r, ( (ũ (l k+ ( ṽ (l ( π (ũ (l k+ (i ṽ (l N ( k+ d ( B j= k+ (i B ((ũ k ṽ k Bũ k d π N j= φ j ( i û (l+ k+, j φ j (d û (l+ k+, j r ( (ũ k ( i ṽ k ( i Bũ k ( i r. Recall that the assumed approimations for the domain integrals are meaningful provided that the radius of the sub-domains is sufficientl small. 7 Test problems The domain and boundar integrals are evaluated as demonstrated in the previous section. Because of the computational techniques described in the previous section, the radius of each local sub-domains, r, should be small enough. A ver small r also causes much cancelation error and therefore, it is chosen as.5h r.h. The shape parameter of the Gaussian weight function is chosen as c i.6h, and the radius of the support of the weight function corresponding to node i, r i, is chosen as r i 7h where h is the minimum distance between nodes. For the predictor corrector scheme, 3 or 4 corrections are calculated in each time step. 7. Eample. For the first test problem consider the reaction diffusion Brusselator Sstem ( with D u = D v = 4, B = and A = in the region Ω = {(, :, } and the boundar condition etracted from the eact solution given b: u(,,t = ep( t, v(,,t = ep( t + +. In this eample, the relative error which will be reported is defined as: N i= e u R = (u i ū i N i= (u, i
11 Meshless Simulations for Brusselator Sstem 69 and the infinit-norm are defined as: e u = Ma{ u i ū i, i =,,...,N}, u i and ū i are the eact and approimate value of u at point i, respectivel, and N is the number of nodes. In order to show the convergence of the proposed method with respect to the spatial variables, Table presents the results obtained with different number of nodal points and fied time step t =. at time instant t =. B going through each column of Table, one can see increasing accurac with increasing the number of nodal points. Table : The results obtained with different number of nodal points, t =. at t =. N e u e v e u R e v R The results obtained at time instant t = 4 with using N = 44 nodal points and different size of time step t are presented in Table. B going through each column of Table, one can see increasing accurac with decreasing the size of the time step t. Table 3 presents the results obtained at different time instants with Table : Results obtained at t = 4 with using N = 44 nodal points and various time steps t t e u e v e u R e v R using N = 44 nodal points and the time step t =.. B going through each column of Table 3, we can see increasing accurac of u and decreasing accurac of v. B the finite differences applied to discretize the time variable, the behavior of v is natural. The reason of increasing the accurac of u b increasing time variable t is cancelation error; b increasing t, u and ver small values of u causes cancelation error.
12 7 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 Table 3: Results obtained with using time step t =., N = 44 nodal points at different time instants t e u e v e u R e v R Error profiles of u and v at t=4 obtained b using N = 44 number of nodal points and t =. are presented in Figures and. 5 Error of u Figure : Error profiles of u with N = 44 and t =. at t = Eample. For the second test problem, we consider the reaction diffusion Brusselator Sstem ( with D u = D v =., B = and A = in the region Ω = {(, :
13 Meshless Simulations for Brusselator Sstem 7 Error of v Figure : Error profiles of v withn = 44 and t =. at t = 4, } with the initial and Neumann s boundar conditions given b: u(,, = 3 3, v(,, = 3 3, u(,,t = u(,,t = u(,,t = u(,,t =, = = = = v(,,t = v(,,t = v(,,t = v(,,t =, = = = = for which the eact solution is unknown. For small values of the diffusion coefficients D u and D v, if B + A > then the stead state solution of the Brusselator sstem ( converges to equilibrium point (A, B/A. This is valid for Eample and is confirmed b the numerical solutions. Here, the relative error which will be reported is defined as: e u R = N i= (ūn i ū N i N, i= (ūn i
14 7 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 and the infinit-norm are defined as: e u = Ma{ ū N i ū N i, i =,,...,N }, where ū N i is the approimation of u at node i obtained with using N nodal points and { i, i =,,...,N } { i, i =,,...,N }. The numerical solution obtained at some points with using N = 44 nodal points and t =. at different time instants are presented in Table 4. This table confirms that (u, v (A, B/A when t. The results obtained at time instant t = 4 with using N = 44 and N = Table 4: The solution obtained at some points. (.,. (.4,.6 (.5,.5 (.8,.9 t u v u v u v u v nodal points and different time step t are presented in Table 5. From this table, it can be seen that b decreasing the size of the time step t, the accurac of the numerical solution increases. Table 5: Results obtained at t = 4, with different time step t, N = 44 and N = points. t e u e v e u R e v R Initial profiles of u and v are presented in Figures 3 and 4. Plots of the numerical solutions obtained for u and v at t = b using N = 44 nodal points and t =
15 Meshless Simulations for Brusselator Sstem 73. are presented in Figures 5 and 6. Figures 7 and 8 presents the numerical solution at point (.5,.5 versus time t. From these figures we can see that at ever point of the domain the solution (u,v tends to (A,B/A...5 u(,, Figure 3: Initial concentration profile of u...5 v(,, Figure 4: Initial concentration profile of v.
16 74 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 u(,, Figure 5: The solution u at t =, N = 44, t = v(,, Figure 6: The solution v solutions at t =, N = 44, t =..
17 Meshless Simulations for Brusselator Sstem u(.5,.5,t obtained b MLPG Figure 7: Plot of u(.5,.5, t versus time t with N = 44, t = v(.5,.5,t obtained b MLPG Figure 8: Plot of v(.5,.5, t versus time t with N = 44, t = Eample 3. Pattern formation in Brusselator sstem Consider the Brusselator equation ( with no flu boundar conditions. The threshold for the Turing instabilit is B T c = [ + A D u /D v ]. The wave number cor-
18 76 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 u, t= v, t= Figure 9: u and v at t= responding to the most unstable Turing mode is given b kc = A/ D u D v. The corresponding characteristic length of the pattern is defined b λ c = π k c. This test problem considers ( with A = 4.5 and D v D u = 8 resulting in B T c 6.7. The pa-
19 Meshless Simulations for Brusselator Sstem 77 u, t= v, t= Figure : u and v at t=5 rameter B is chosen as B = 6.75 resulting in spot patterns as shown in Leppänen (4. The domain of the problem is chosen as [,] and the diffusion coefficients are fied with D v =. D u =.5 resulting in λ c.79. The initial
20 78 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 u, t= v, t= Figure : u and v at t= configuration was random perturbations around the uniform stationar state. Figures 9- presents the process of the pattern formation for the considered model. In all contour graphs, coloration is determined b a constant threshold value, u s for
21 Meshless Simulations for Brusselator Sstem 79 u, t= v, t= Figure : u and v at t=4 u and v s for v, such that in the regions with white color, u < u s and v < v s while the regions which eperience the concentration u > u s and v > v s are colored with green. The results confirm that the profiles of the function v are alwas 8 out of phase to those of u.
22 8 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 8 Conclusions A numerical method for solving the two dimensional nonlinear reaction-diffusion Brusselator sstem was presented. The time variable was discretized b using onestep θ method. Then the resultant elliptic tpe PDEs were considered in weak forms on the local subdomains and the MLS approimation was used to approimate the field variables. The Heaviside step function was used as test function in each local subdomain. The nonlinear terms were treated b using the predictor corrector scheme. Test problems reveal that the method is of high accurac and stabilit. A model sstem ehibiting Turing instabilit was given as a test problem and the resulting patterns were presented. Acknowledgements The financial supports b APVV-3- (Slovak Research and Development Agenc is greatl acknowledged. References Abbasband, S.; Shirzadi, A. (: MLPG method for two-dimensional diffusion equation with Neumann s and non-classical boundar conditions. Appl. Numer. Math., vol. 6, pp Abbasband, S.; Sladek, V.; Shirzadi, A.; Sladek, J. (: Numerical simulations for coupled pair of diffusion equations b MLPG method. CMES Compt. Model. Eng. Sci., vol. 7, no., pp Adomian, G. (995: The diffusion-brusselator equation. Computers and Mathematics with Applications, vol. 9, no. 5, pp. 3. Ang, W.-T. (3: The two-dimensional reaction-diffusion brusselator sstem: a dual-reciprocit boundar element solution. Eng. Anal. Bound. Elem., vol. 7, no. 9, pp Atluri, S. (4: The meshless method (MLPG for domain and BIE discretizations. Tech Science Press. Atluri, S.; Kim, H.-G.; Cho, J. Y. (999: A critical assessment of the trul meshless local Petrov-Galerkin (MLPG and local boundar integral equation (LBIE methods. Comput. Mech., vol. 4, pp Atluri, S.; Shen, S. (: Method. Tech Science Press. The Meshless Local Petrov-Galerkin (MLPG Atluri, S.; Zhu, T. (998: A new meshless local Petrov-Galerkin (MLPG approach in computational mechanics. Comput. Mech., vol., no., pp. 7 7.
23 Meshless Simulations for Brusselator Sstem 8 Atluri, S.; Zhu, T. (998: A new meshless local Petrov-Galerkin (MLPG approach to nonlinear problems in computer modeling and simulation. Comput. Modeling Simulation in Engrg., vol. 3, pp Borckmans, P.; Wit, A. D.; Dewel, G. (99: Competition in ramped turing structures. Phsica A Statistical Mechanics and its Applications, vol. 88, no. -3, pp Kamranian, M.; Tatari, M.; Dehghan, M. (: The finite point method for reaction-diffusion sstems in developmental biolog. CMES Compt. Model. Eng. Sci., vol. 8, no., pp. 7. Leppänen, T. (4: Computational Studies of Pattern Formation in Turing Sstems. PhD thesis, Helsinki Universit of Technolog, Espoo, Finland, 4. Mittal, R.; Jiwari, R. (: Numerical solution of two-dimensional reactiondiffusion brusselator sstem. Appl. Math. Comp., vol. 7, no., pp Nicolis, G.; Prigogine, I. (997: Self-organization in nonequilibrium sstems: from dissipative structures to order through fluctuations. John Wile & Sons. Prigogine, I.; Lefever, R. (968: Smmetr breaking instabilities in dissipative sstems.ii. J. Chem. Phs., vol. 48, pp Shirzadi, A.; Ling, L. (3: Convergent overdetermined-rbf-mlpg for solving second order elliptic PDEs. Adv. Appl. Math. Mech., vol. 5, pp Siraj-ul, I.; Arshed, A.; Sirajul, H. (: A computational modeling of the behavior of the two- dimensional reaction-diffusion brusselator sstem. Appl. Math. Model., vol. 34, no., pp Sladek, J.; Stanak, P.; Han, Z. D.; Sladek, V.; Atluri, S. N. (3: Applications of the MLPG method in engineering & sciences: A review. CMES Compt. Model. Eng. Sci., vol. 9, no. 5, pp Sladek, V.; Sladek, J. (: Local integral equations implemented b MLSapproimation and analtical integrations. Eng. Anal. Bound. Elem., vol. 34, no., pp Sladek, V.; Sladek, J.; Zhang, C. (: On increasing computational efficienc of local integral equation method combined with meshless implementations. CMES Compt. Model. Eng. Sci., vol. 63, no. 3, pp Turing, A. (95: The chemical basis of morphogenesis. Phil. Trans. R. Soc. Lond. B, vol. 37, no. 64, pp Twizell, E. H.; Gumel, A. B.; Cao, Q. (999: A second-order scheme for the Brusselator reaction-diffusion sstem. J. Math. Chem., vol. 6, no. 4, pp (.
24 8 Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 Verdasca, J.; Wit, A. D.; Dewel, G.; Borckmans, P. (99: Reentrant heagonal turing structures. Phsics Letters A, vol. 68, no. 3, pp Walgraef, D.; Dewel, G.; Borckmans, P. (98: Fluctuations near nonequilibrium phase transitions to nonuniform states. Phs. Rev. A, vol., pp
A Meshless Method for Solving the 2D Brusselator Reaction-Diffusion System
Copright Tech Science Press CMES, ol., no., pp.-8, A Meshless Method for Soling the D Brusselator Reaction-Diffusion Sstem M. Mohammadi, R. Mokhtari, and R. Schaback Abstract: In this paper, the two-dimensional
More informationNumerical solution of Maxwell equations using local weak form meshless techniques
Journal of mathematics and computer science 13 2014), 168-185 Numerical solution of Maxwell equations using local weak form meshless techniques S. Sarabadan 1, M. Shahrezaee 1, J.A. Rad 2, K. Parand 2,*
More informationOn the Extension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models
On the Etension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models J. T. Oden, S. Prudhomme, and P. Bauman Institute for Computational Engineering and Sciences The Universit
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationNUMERICAL MLPG ANALYSIS OF PIEZOELECTRIC SENSOR IN STRUCTURES
DOI: 10.2478/sjce-2014-0009 NUMERICAL MLPG ANALYSIS OF PIEZOELECTRIC SENSOR IN STRUCTURES Peter STAŇÁK 1*, Ján SLÁDEK 1, Vladimír SLÁDEK 1, Slavomír KRAHULEC 1 Abstract The paper deals with a numerical
More information* τσ σκ. Supporting Text. A. Stability Analysis of System 2
Supporting Tet A. Stabilit Analsis of Sstem In this Appendi, we stud the stabilit of the equilibria of sstem. If we redefine the sstem as, T when -, then there are at most three equilibria: E,, E κ -,,
More informationAn improved ADI-DQM based on Bernstein polynomial for solving two-dimensional convection-diffusion equations
Mathematical Theor and Modeling ISSN -50 (Paper) ISSN 5-05(Online) Vol 3, No.1, 013 An improved ADI-DQM based on Bernstein polnomial for solving two-dimensional convection-diffusion equations A.S.J. Al-
More informationKrylov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations
DOI.7/s95-6-6-7 Krlov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations Dong Lu Yong-Tao Zhang Received: September 5 / Revised: 9 March 6 / Accepted: 8
More information520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationSolving Heat Conduction Problems by the Direct Meshless Local Petrov-Galerkin (DMLPG) method
Solving Heat Conduction Problems by the Direct Meshless Local Petrov-Galerkin (DMLPG) method Davoud Mirzaei, Robert Schaback Department of Mathematics, University of Isfahan, 8745-63 Isfahan, Iran. Institut
More informationPerformance of Multiquadric Collocation Method in Solving Lid-driven Cavity Flow Problem with Low Reynolds Number
Copyright c 2006 Tech Science Press CMES, vol.15, no.3, pp.137-146, 2006 Performance of Multiquadric Collocation Method in Solving Lid-driven Cavity Flow Problem with Low Reynolds umber S. Chantasiriwan
More informationOPTIMAL CONTROL FOR A PARABOLIC SYSTEM MODELLING CHEMOTAXIS
Trends in Mathematics Information Center for Mathematical Sciences Volume 6, Number 1, June, 23, Pages 45 49 OPTIMAL CONTROL FOR A PARABOLIC SYSTEM MODELLING CHEMOTAXIS SANG UK RYU Abstract. We stud the
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationUNCORRECTED. To recognise the rules of a number of common algebraic relations: y = x 1 y 2 = x
5A galler of graphs Objectives To recognise the rules of a number of common algebraic relations: = = = (rectangular hperbola) + = (circle). To be able to sketch the graphs of these relations. To be able
More informationSolving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells
Copyright 2010 Tech Science Press CMES, vol.57, no.2, pp.109-135, 2010 Solving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells E. J. Sellountos 1,
More informationAN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES *
Journal of Computational Mathematics, Vol.27, No.2-3, 29, 388 399. AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES * Yin Yang Hunan Ke Laborator for Computation and Simulation
More informationLESSON #42 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART 2 COMMON CORE ALGEBRA II
LESSON #4 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART COMMON CORE ALGEBRA II You will recall from unit 1 that in order to find the inverse of a function, ou must switch and and solve for. Also,
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
More informationSolving Variable-Coefficient Fourth-Order ODEs with Polynomial Nonlinearity by Symmetric Homotopy Method
Applied and Computational Mathematics 218; 7(2): 58-7 http://www.sciencepublishinggroup.com/j/acm doi: 1.11648/j.acm.21872.14 ISSN: 2328-565 (Print); ISSN: 2328-5613 (Online) Solving Variable-Coefficient
More informationDiscontinuous Galerkin method for a class of elliptic multi-scale problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 000; 00: 6 [Version: 00/09/8 v.0] Discontinuous Galerkin method for a class of elliptic multi-scale problems Ling Yuan
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More informationA high order Padé ADI method for unsteady convection-diffusion equations
Center for Turbulence Research Annual Research Briefs 2005 85 A high order Padé ADI method for unstead convection-diffusion equations B D. You 1. Motivation and objectives The unstead convection-diffusion
More informationSTABILIZED FEM SOLUTIONS OF MHD EQUATIONS AROUND A SOLID AND INSIDE A CONDUCTING MEDIUM
Available online: March 09, 2018 Commun. Fac. Sci. Univ. Ank. Ser. A1 Math. Stat. Volume 68, Number 1, Pages 197 208 (2019) DOI: 10.1501/Commua1_0000000901 ISSN 1303 5991 http://communications.science.ankara.edu.tr/index.php?series=a1
More informationSymmetry Arguments and the Role They Play in Using Gauss Law
Smmetr Arguments and the Role The la in Using Gauss Law K. M. Westerberg (9/2005) Smmetr plas a ver important role in science in general, and phsics in particular. Arguments based on smmetr can often simplif
More informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationRANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM 1
RANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM Jirka Roubal Vladimír Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Karlovo náměstí
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationResearch Article Chaotic Attractor Generation via a Simple Linear Time-Varying System
Discrete Dnamics in Nature and Societ Volume, Article ID 836, 8 pages doi:.//836 Research Article Chaotic Attractor Generation via a Simple Linear Time-Varing Sstem Baiu Ou and Desheng Liu Department of
More informationy R T However, the calculations are easier, when carried out using the polar set of co-ordinates ϕ,r. The relations between the co-ordinates are:
Curved beams. Introduction Curved beams also called arches were invented about ears ago. he purpose was to form such a structure that would transfer loads, mainl the dead weight, to the ground b the elements
More informationDynamics of multiple pendula without gravity
Chaotic Modeling and Simulation (CMSIM) 1: 57 67, 014 Dnamics of multiple pendula without gravit Wojciech Szumiński Institute of Phsics, Universit of Zielona Góra, Poland (E-mail: uz88szuminski@gmail.com)
More informationUnit 3 Notes Mathematical Methods
Unit 3 Notes Mathematical Methods Foundational Knowledge Created b Triumph Tutoring Copright info Copright Triumph Tutoring 07 Triumph Tutoring Pt Ltd ABN 60 607 0 507 First published in 07 All rights
More informationVibrational Power Flow Considerations Arising From Multi-Dimensional Isolators. Abstract
Vibrational Power Flow Considerations Arising From Multi-Dimensional Isolators Rajendra Singh and Seungbo Kim The Ohio State Universit Columbus, OH 4321-117, USA Abstract Much of the vibration isolation
More informationCOMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER PARABOLIC EQUATIONS. Lili Ju. Xinfeng Liu.
DISCRETE AND CONTINUOUS doi:13934/dcdsb214191667 DYNAMICAL SYSTEMS SERIES B Volume 19, Number 6, August 214 pp 1667 1687 COMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER
More informationNUMERICAL ANALYSIS OF AN ENERGY-CONSERVATION SCHEME FOR TWO-DIMENSIONAL HAMILTONIAN WAVE EQUATIONS WITH NEUMANN BOUNDARY CONDITIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 6, Number, Pages 39 339 c 9 Institute for Scientific Computing and Information NUMERICAL ANALYSIS OF AN ENERGY-CONSERVATION SCHEME FOR TWO-DIMENSIONAL
More informationMath 030 Review for Final Exam Revised Fall 2010 RH/ DM 1
Math 00 Review for Final Eam Revised Fall 010 RH/ DM 1 1. Solve the equations: (-1) (7) (-) (-1) () 1 1 1 1 f. 1 g. h. 1 11 i. 9. Solve the following equations for the given variable: 1 Solve for. D ab
More informationSimulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundary Element Method
10 th International LS-DYNA Users Conference Simulation Technolog (2) Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundar Element Method Yun Huang Livermore Software Technolog Corporation
More informationResearch Article Chaos and Control of Game Model Based on Heterogeneous Expectations in Electric Power Triopoly
Discrete Dnamics in Nature and Societ Volume 29, Article ID 469564, 8 pages doi:.55/29/469564 Research Article Chaos and Control of Game Model Based on Heterogeneous Epectations in Electric Power Triopol
More informationThe Hopf Bifurcation Theorem: Abstract. Introduction. Transversality condition; the eigenvalues cross the imginary axis with non-zero speed
Supercritical and Subcritical Hopf Bifurcations in Non Linear Maps Tarini Kumar Dutta, Department of Mathematics, Gauhati Universit Pramila Kumari Prajapati, Department of Mathematics, Gauhati Universit
More informationAdomian Decomposition Method for Solving the Kuramoto Sivashinsky Equation
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn:2319-765x. Volume 1, Issue 1Ver. I. (Jan. 214), PP 8-12 Adomian Decomposition Method for Solving the Kuramoto Sivashinsky Equation Saad A.
More informationarxiv: v1 [math.na] 19 Nov 2016
A new numerical algorithm for the solutions of hperbolic partial differential equations in ( + )-dimensional space Brajesh Kumar Singh arxiv:6697v [mathna] 9 Nov 6 Department of Applied Mathematics, School
More informationMeshless Local Petrov-Galerkin Mixed Collocation Method for Solving Cauchy Inverse Problems of Steady-State Heat Transfer
Copyright 2014 Tech Science Press CMES, vol.97, no.6, pp.509-533, 2014 Meshless Local Petrov-Galerkin Mixed Collocation Method for Solving Cauchy Inverse Problems of Steady-State Heat Transfer Tao Zhang
More informationVelocity Limit in DPD Simulations of Wall-Bounded Flows
Velocit Limit in DPD Simulations of Wall-Bounded Flows Dmitr A. Fedosov, Igor V. Pivkin and George Em Karniadakis Division of Applied Mathematics, Brown Universit, Providence, RI 2912 USA Abstract Dissipative
More informationImplicit Second Derivative Hybrid Linear Multistep Method with Nested Predictors for Ordinary Differential Equations
American Scientiic Research Journal or Engineering, Technolog, and Sciences (ASRJETS) ISSN (Print) -44, ISSN (Online) -44 Global Societ o Scientiic Research and Researchers http://asretsournal.org/ Implicit
More information3.6 Applications of Poisson s equation
William J. Parnell: MT3432. Section 3: Green s functions in 2 and 3D 69 3.6 Applications of Poisson s equation The beaut of Green s functions is that the are useful! Let us consider an example here where
More informationSAMPLE. A Gallery of Graphs. To recognise the rules of a number of common algebraic relationships: y = x 1,
Objectives C H A P T E R 5 A Galler of Graphs To recognise the rules of a number of common algebraic relationships: =, =, = / and + =. To be able to sketch the graphs and simple transformations of these
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationA MANIFOLD-BASED EXPONENTIALLY CONVERGENT ALGORITHM FOR SOLVING NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Journal of Marine Science and Technology Vol. No. pp. -9 () DOI:.69/JMST--- A MANIOLD-ASED EXPONENTIALLY CONVERGENT ALGORITHM OR SOLVING NON-LINEAR PARTIAL DIERENTIAL EQUATIONS Chein-Shan Liu Key words:
More information8.1 Exponents and Roots
Section 8. Eponents and Roots 75 8. Eponents and Roots Before defining the net famil of functions, the eponential functions, we will need to discuss eponent notation in detail. As we shall see, eponents
More informationPICONE S IDENTITY FOR A SYSTEM OF FIRST-ORDER NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS
Electronic Journal of Differential Equations, Vol. 2013 (2013), No. 143, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu PICONE S IDENTITY
More informationPeriodic Structures in FDTD
EE 5303 Electromagnetic Analsis Using Finite Difference Time Domain Lecture #19 Periodic Structures in FDTD Lecture 19 These notes ma contain coprighted material obtained under fair use rules. Distribution
More informationCoordinate geometry. + bx + c. Vertical asymptote. Sketch graphs of hyperbolas (including asymptotic behaviour) from the general
A Sketch graphs of = a m b n c where m = or and n = or B Reciprocal graphs C Graphs of circles and ellipses D Graphs of hperbolas E Partial fractions F Sketch graphs using partial fractions Coordinate
More informationA Numerical Algorithm for Ambrosetti-Prodi Type Operators
A Numerical Algorithm for Ambrosetti-Prodi Tpe Operators José Teieira Cal Neto and Carlos Tomei March 9, Abstract We consider the numerical solution of the equation u f(u) = g, for the unknown u satisfing
More informationTowards Higher-Order Schemes for Compressible Flow
WDS'6 Proceedings of Contributed Papers, Part I, 5 4, 6. ISBN 8-867-84- MATFYZPRESS Towards Higher-Order Schemes for Compressible Flow K. Findejs Charles Universit, Facult of Mathematics and Phsics, Prague,
More informationDATA STRUCTURES FOR THE NONLINEAR SUBGRID STABILIZATION METHOD
DATA STRUCTURES FOR TE NONLINEAR SUBGRID STABILIZATION METOD Isaac P. Santos Lucia Catabriga Departamento de Informática Universidade Federal do Espírito Santo - UFES Av. Fernando Ferrari, 54 - Vitória,
More informationDISCONTINUOUS GALERKIN DISCRETIZATION WITH EMBEDDED BOUNDARY CONDITIONS
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vol.3(23), No., pp.35 58 c 23 Editorial Board of the Journal Computational Methods in Applied Mathematics Joint Co. Ltd. DISCONTINUOUS GALERKIN DISCRETIZATION
More informationA FLOW-CONDITION-BASED INTERPOLATION MIXED FINITE ELEMENT PROCEDURE FOR HIGHER REYNOLDS NUMBER FLUID FLOWS
Mathematical Models and Methods in Applied Sciences Vol. 12, No. 4 (2002) 525 539 c World Scientific Publishing Compan A FLOW-CONDITION-BASED INTERPOLATION MIXED FINITE ELEMENT PROCEDURE FOR HIGHER REYNOLDS
More informationTwo-sided Mullins-Sekerka flow does not preserve convexity
Differential Equations and Computational Simulations III J. Graef, R. Shivaji, B. Soni, & J. Zhu (Editors) Electronic Journal of Differential Equations, Conference 01, 1997, pp. 171 179. ISSN: 1072-6691.
More informationNonlinear Systems Examples Sheet: Solutions
Nonlinear Sstems Eamples Sheet: Solutions Mark Cannon, Michaelmas Term 7 Equilibrium points. (a). Solving ẋ =sin 4 3 =for gives =as an equilibrium point. This is the onl equilibrium because there is onl
More informationMATHEMATICAL MODELING AND THE STABILITY STUDY OF SOME CHEMICAL PHENOMENA
HENRI COND IR FORCE CDEM ROMNI INTERNTIONL CONFERENCE of SCIENTIFIC PPER FSES rasov -6 Ma GENERL M.R. STEFNIK RMED FORCES CDEM SLOVK REPULIC MTHEMTICL MODELING ND THE STILIT STUD OF SOME CHEMICL PHENOMEN
More information7. Choice of approximating for the FE-method scalar problems
7. Choice of approimating for the FE-method scalar problems Finite Element Method Differential Equation Weak Formulation Approimating Functions Weighted Residuals FEM - Formulation Weak form of heat flow
More informationRELATIONS AND FUNCTIONS through
RELATIONS AND FUNCTIONS 11.1.2 through 11.1. Relations and Functions establish a correspondence between the input values (usuall ) and the output values (usuall ) according to the particular relation or
More informationChapter 4 Analytic Trigonometry
Analtic Trigonometr Chapter Analtic Trigonometr Inverse Trigonometric Functions The trigonometric functions act as an operator on the variable (angle, resulting in an output value Suppose this process
More informationEntropy ISSN
Entrop 5, 7[4], 3 37 3 Entrop ISSN 99-43 www.mdpi.org/entrop/ Full Paper Second law analsis of laminar flow in a channel filled with saturated porous media: a numerical solution Kamel Hooman, Arash Ejlali
More informationStability Analysis of Laminated Composite Thin-Walled Beam Structures
Paper 224 Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures echnolog, B.H.V. opping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Stabilit nalsis
More information39. (a) Use trigonometric substitution to verify that. 40. The parabola y 2x divides the disk into two
35. Prove the formula A r for the area of a sector of a circle with radius r and central angle. [Hint: Assume 0 and place the center of the circle at the origin so it has the equation. Then is the sum
More informationGreen s Theorem Jeremy Orloff
Green s Theorem Jerem Orloff Line integrals and Green s theorem. Vector Fields Vector notation. In 8.4 we will mostl use the notation (v) = (a, b) for vectors. The other common notation (v) = ai + bj runs
More informationEFFECT OF DAMPING AND THERMAL GRADIENT ON VIBRATIONS OF ORTHOTROPIC RECTANGULAR PLATE OF VARIABLE THICKNESS
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 1, Issue 1 (June 17), pp. 1-16 Applications and Applied Mathematics: An International Journal (AAM) EFFECT OF DAMPING AND THERMAL
More informationA NEW CHAOTIC ATTRACTOR FROM 2D DISCRETE MAPPING VIA BORDER-COLLISION PERIOD-DOUBLING SCENARIO
A NEW CHAOTIC ATTRACTOR FROM D DISCRETE MAPPING VIA BORDER-COLLISION PERIOD-DOUBLING SCENARIO ZERAOULIA ELHADJ Received 1 April 5 The following map is studied: (, ) (1 + a( )+,b). It is proved numericall
More informationA new operational matrix for solving twodimensional nonlinear integral equations of fractional order
Jabari Sabeg et al. Cogent Mathematics (07) 4: 34707 https://doi.org/0.080/33835.07.34707 APPLIED & INTERDISCIPLINARY MATHEMATICS RESEARCH ARTICLE A new operational matri for solving twodimensional nonlinear
More informationStability Of Structures: Continuous Models
5 Stabilit Of Structures: Continuous Models SEN 311 ecture 5 Slide 1 Objective SEN 311 - Structures This ecture covers continuous models for structural stabilit. Focus is on aiall loaded columns with various
More information4.7. Newton s Method. Procedure for Newton s Method HISTORICAL BIOGRAPHY
4. Newton s Method 99 4. Newton s Method HISTORICAL BIOGRAPHY Niels Henrik Abel (18 189) One of the basic problems of mathematics is solving equations. Using the quadratic root formula, we know how to
More informationElliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations
5 4. Properties of Laplace s and Poisson s Equations Chapter 4 Elliptic Equations Contents. Neumann conditions the normal derivative, / = n u is prescribed on the boundar second BP. In this case we have
More informationThe Force Table Introduction: Theory:
1 The Force Table Introduction: "The Force Table" is a simple tool for demonstrating Newton s First Law and the vector nature of forces. This tool is based on the principle of equilibrium. An object is
More informationEstimation Of Linearised Fluid Film Coefficients In A Rotor Bearing System Subjected To Random Excitation
Estimation Of Linearised Fluid Film Coefficients In A Rotor Bearing Sstem Subjected To Random Ecitation Arshad. Khan and Ahmad A. Khan Department of Mechanical Engineering Z.. College of Engineering &
More informationA Method for Geometry Optimization in a Simple Model of Two-Dimensional Heat Transfer
A Method for Geometr Optimization in a Simple Model of Two-Dimensional Heat Transfer X. Peng, K. Niakhai B. Protas Jul, 3 arxiv:37.48v [math.na] 4 Jul 3 Abstract This investigation is motivated b the problem
More informationStrain Transformation and Rosette Gage Theory
Strain Transformation and Rosette Gage Theor It is often desired to measure the full state of strain on the surface of a part, that is to measure not onl the two etensional strains, and, but also the shear
More informationPolynomial and Rational Functions
Polnomial and Rational Functions Figure -mm film, once the standard for capturing photographic images, has been made largel obsolete b digital photograph. (credit film : modification of work b Horia Varlan;
More informationAn Error Minimized Pseudospectral Penalty Direct Poisson Solver
An Error Minimized Pseudospectral Penalt Direct Poisson Solver b Horng Tz-Leng Dept. of Applied Mathematics Feng Chia Universit Taichung Taiwan tlhorng123@gmail.com http://newton.math.fcu.edu.tw/~tlhorng
More informationExponential time differencing methods for nonlinear PDEs
XXI Congreso de Ecuaciones Diferenciales Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, 1- septiembre 9 (pp. 1 8) Exponential time differencing methods for nonlinear PDEs F. de la Hoz 1,
More informationStability Analysis of a Geometrically Imperfect Structure using a Random Field Model
Stabilit Analsis of a Geometricall Imperfect Structure using a Random Field Model JAN VALEŠ, ZDENĚK KALA Department of Structural Mechanics Brno Universit of Technolog, Facult of Civil Engineering Veveří
More informationA Direct Derivation of the Griffith-Irwin Relationship using a Crack tip Unloading Stress Wave Model.
A Direct Derivation of the Griffith-Irwin Relationship using a Crack tip Unloading Stress Wave Model. C.E. Neal-Sturgess. Emeritus Professor of Mechanical Engineering, The Universit of Birmingham, UK.
More informationImproved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations
American Scientific Researc Journal for Engineering, Tecnolog, and Sciences (ASRJETS) ISSN (Print) 33-44, ISSN (Online) 33-44 Global Societ of Scientific Researc and Researcers ttp://asrjetsjournal.org/
More informationEstimators in simple random sampling: Searls approach
Songklanakarin J Sci Technol 35 (6 749-760 Nov - Dec 013 http://wwwsjstpsuacth Original Article Estimators in simple random sampling: Searls approach Jirawan Jitthavech* and Vichit Lorchirachoonkul School
More informationCONTINUOUS SPATIAL DATA ANALYSIS
CONTINUOUS SPATIAL DATA ANALSIS 1. Overview of Spatial Stochastic Processes The ke difference between continuous spatial data and point patterns is that there is now assumed to be a meaningful value, s
More informationUNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD
Proceedings of the 3rd ASME/JSME Joint Fluids Engineering Conference Jul 8-23, 999, San Francisco, California FEDSM99-8 UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX
More informationMeshless Local Petrov-Galerkin (MLPG) Methods in
Meshless Local Petrov-Galerkin MLPG) Methods in Quantum Mechanics *Williams L Nicomedes Renato C Mesquita *Ferno J S Moreira * Federal University of Minas Gerais Dept of Electronic Engineering Belo Horizonte
More informationSubband Coding and Wavelets. National Chiao Tung University Chun-Jen Tsai 12/04/2014
Subband Coding and Wavelets National Chiao Tung Universit Chun-Jen Tsai /4/4 Concept of Subband Coding In transform coding, we use N (or N N) samples as the data transform unit Transform coefficients are
More informationControl Schemes to Reduce Risk of Extinction in the Lotka-Volterra Predator-Prey Model
Journal of Applied Mathematics and Phsics, 2014, 2, 644-652 Published Online June 2014 in SciRes. http://www.scirp.org/journal/jamp http://d.doi.org/10.4236/jamp.2014.27071 Control Schemes to Reduce Risk
More informationChaos Control and Synchronization of a Fractional-order Autonomous System
Chaos Control and Snchronization of a Fractional-order Autonomous Sstem WANG HONGWU Tianjin Universit, School of Management Weijin Street 9, 37 Tianjin Tianjin Universit of Science and Technolog College
More informationOn the Missing Modes When Using the Exact Frequency Relationship between Kirchhoff and Mindlin Plates
On the Missing Modes When Using the Eact Frequenc Relationship between Kirchhoff and Mindlin Plates C.W. Lim 1,*, Z.R. Li 1, Y. Xiang, G.W. Wei 3 and C.M. Wang 4 1 Department of Building and Construction,
More informationNew approach to study the van der Pol equation for large damping
Electronic Journal of Qualitative Theor of Differential Equations 2018, No. 8, 1 10; https://doi.org/10.1422/ejqtde.2018.1.8 www.math.u-szeged.hu/ejqtde/ New approach to stud the van der Pol equation for
More informationInterspecific Segregation and Phase Transition in a Lattice Ecosystem with Intraspecific Competition
Interspecific Segregation and Phase Transition in a Lattice Ecosstem with Intraspecific Competition K. Tainaka a, M. Kushida a, Y. Ito a and J. Yoshimura a,b,c a Department of Sstems Engineering, Shizuoka
More informationCopyrighted by Gabriel Tang B.Ed., B.Sc. Page 1.
Chapter : Linear and Quadratic Functions Chapter : Linear and Quadratic Functions -: Points and Lines Sstem of Linear Equations: - two or more linear equations on the same coordinate grid. Solution of
More informationAnti-diffusive Methods for Hyperbolic Heat Transfer
Anti-diffusive Methods for Hperbolic Heat Transfer Wensheng Shen and Leigh Little Department of Computational Science, SUNY Brockport Brockport, NY 44, USA Liangjian Hu Department of Applied Mathematics,
More informationSHADOW AND BACKGROUND ORIENTED SCHLIEREN INVESTIGATION OF SHOCK WAVES IN GAS-DISCHARGE MEDIUM
SHADOW AND BACKGROUND ORIENTED SCHLIEREN INVESTIGATION OF SHOCK WAVES IN GAS-DISCHARGE MEDIUM J. JIN 1, I.V. MURSENKOVA 1,c, N. N. SYSOEV 1 and I.A. ZNAMENSKAYA 1 1 Department of Phsics, Lomonosov Moscow
More informationAn Optimization Method for Numerically Solving Three-point BVPs of Linear Second-order ODEs with Variable Coefficients
An Optimization Method for Numericall Solving Three-point BVPs of Linear Second-order ODEs with Variable Coefficients Liaocheng Universit School of Mathematics Sciences 252059, Liaocheng P.R. CHINA hougbb@26.com;
More informationSimulation and Experimental Validation of Chaotic Behavior of Airflow in a Ventilated Room
Simulation and Eperimental Validation of Chaotic Behavior of Airflow in a Ventilated Room Jos van Schijndel, Assistant Professor Eindhoven Universit of Technolog, Netherlands KEYWORDS: Airflow, chaos,
More informationMeshless equilibrium on line method (MELM) for linear elasticity
Structural Engineering and Mechanics, Vol. 35, No. 4 (2010) 000-000 1 Meshless equilibrium on line method (MELM) for linear elasticity A. Sadeghirad*, S. Mohammadi and I. Mahmoudzadeh Kani School of Civil
More informationNEW IMPLEMENTATION OF MLBIE METHOD FOR HEAT CONDUCTION ANALYSIS IN FUNCTIONALLY GRADED MATERIALS
NEW IMPLEMENTATION OF MLBIE METHOD FOR HEAT CONDUCTION ANALYSIS IN FUNCTIONALLY GRADED MATERIALS DAVOUD MIRZAEI,,, MEHDI DEHGHAN Abstract. The meshless local boundary integral equation (MLBIE) method with
More information