AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES *
|
|
- Sabina Willis
- 6 years ago
- Views:
Transcription
1 Journal of Computational Mathematics, Vol.27, No.2-3, 29, AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES * Yin Yang Hunan Ke Laborator for Computation and Simulation in Science and Engineering, School of Mathematics and Computational Science, Xiangtan Universit, Xiangtan 41115, China angintu@tu.edu.cn Yanping Chen School of Mathematical Sciences, South China Normal Universit, Guangzhou, China anpingchen@scnu.edu.cn Yunqing Huang Hunan Ke Laborator for Computation and Simulation in Science and Engineering, School of Mathematics and Computational Science, Xiangtan Universit, Xiangtan 41115, China huangq@tu.edu.cn Abstract This paper presents an efficient moving mesh method to solve a nonlinear singular problem with an optimal control constrained condition. The phsical problem is governed b a new model of turbulent flow in circular tubes proposed b Luo et al. using Prandtl s miing-length theor. Our algorithm is formed b an outer iterative algorithm for handling the optimal control condition and an inner adaptive mesh redistribution algorithm for solving the singular governing equations. We discretize the nonlinear problem b using a upwinding approach, and the resulting nonlinear equations are solved b using the Newton- Raphson method. The mesh is generated and the grid points are moved b using the arc-length equidistribution principle. The numerical results demonstrate that proposed algorithm is effective in capturing the boundar laers associated with the turbulent model. Mathematics subject classification: 65L1, 65L12. Ke words: Edd viscosit, Turbulent pipe flow, Boundar laer, Optimal control, Moving mesh. 1. Introduction The modeling of turbulent flows still plas an important role in computational fluid dnamics because direct simulation of flows are restricted to ver simple geometries and low Renolds number [6,9,18,19]. Development of turbulence model is therefore still an important task and even some semi-empirical means such as the edd viscosit or Prandtl s miing length are ver helpful to deal with man problems in engineering practice due to their simplicit. Luo et al. [15] established a new model of turbulent flow in circular tubes which is an application and improvement of Prandtl s miing-length theor. The model epresses the single phase flow in circular tubes, which is an optimal parameter control problem governed b a nonlinear singular equation. The model ields man comple mathematical characters such as strong boundar laer. The computational results resulting from the new model are found in good agreement with the eperimental results on fluid velocit distribution, edd viscosit distribution and friction factor. On the mathematical side, the governing equations associate with this model * Received Januar 5, 28 / Revised version received August 5, 28 / Accepted September 5, 28 /
2 Moving Mesh Method for a Model of Turbulent Flow 389 are quite complicated, and an effective numerical scheme for finding numerical approimations seems useful. The difficulties of this problem include the eistence of the boundar laers and an optional control condition enforced on the governing equations. To resolve the laers, numerical simulations require etremel fine meshes on the small localized portions of the phsical domain. It will become ver epensive if a uniform fine mesh is emploed. Recent stud has demonstrated that moving mesh methods are powerful in resolving large solution variations b increasing the solution accurac and decreasing the cost of computations, see, e.g., Adaptive moving mesh methods have important applications in solving partial differential equations (PDEs). Up to now, there have been important progresses [1, 3, 2]. Harten and Hman [11] began the earliest stud of the adaptive methods to improve resolution of discontinuous solutions of hperbolic equations. After their work, man other moving mesh methods on this direction have been proposed based on combining the variational grid methods with high resolution shock capturing methods, including the so-called moving mesh PDE (MMPDE) approach of moving mesh methods of W. Huang [1], moving finite element methods of Miller [17], and moving finite volume methods [22]. Recentl, there have been works on moving mesh methods based on Harmonic maps [8,16,22]. Theoretical results on adaptive mesh arising from equidistribution of a monitor function can be found in [3,4,12 14,21]. In particular, Kopteva [13] derived certain maimum norm a posterior error estimates for one-dimensional singularl perturbed convection-diffusion problems, see also a recent paper [14] for a similar posterior error estimate. The aim of this paper is to present an efficient and fast numerical method for the turbulent model. The proposed numerical algorithm includes two parts: (i) the outer iterative algorithm is used to solve the optimal control condition and (ii) the inner adaptive mesh redistribution algorithm is used to solve the singular problem. We discretize this nonlinear problem b using upwinding scheme. The discretized nonlinear equations is solved b Newton-Raphson method. The arc-length equidistribution principle is used in the part (ii) above. The numerical eamples will be provided to demonstrate the effectiveness of the proposed algorithm. This paper is organized as follows. In Section 2, we briefl review the model of turbulent flow in circular tubes b emploing Prandtl s miing-length theor. In Section 3, we will present the discrete schemes and algorithms. Numerical eperiments will be carried out in Section 4. Some concluding remarks will be presented in the final section. 2. A New Model of Turbulent Flow in Circular Tubes In this section, we briefl review the background of the model of turbulent flow in circular tubes which was proposed b Luo et al. [15]. Moreover, using dimensionless analsis we will derive a complete mathematical description for this model. Note that the shearing stress of Newtonian fluid for turbulent flow can be described b edd viscosit with dimensionless analsis [9]. We then have following epression: dũ dφ = ˆRφ 1 + µ t /µ L, (2.1) where ũ is dimensionless time-smoothed velocit, µ t is edd viscosit, µ L is kinematic viscosit, ˆR is dimensionless radius of a circular pipe and ˆR = ρûr/µ with R is the tube radius, ρ is the liquid densit, û is friction velocit, µ is the molecule viscosit, φ is dimensionless radial position in a circular pipe, φ = is the center of the tube and φ = 1 corresponds to the wall of
3 39 Y. YANG, Y. CHEN AND Y. HUANG the tube. As the velocit is disappeared on the wall, we have the boundar condition ũ(φ) φ=1 =. (2.2) Based on Prandtl s miing-length theor, the edd viscosit can be written as µ t = µ ˆR L L dũ L R R dφ, (2.3) and the constrained condition is 2 1 ũφdφ = ṷ u = Re ˆR, (2.4) where L is the miing length, u is average velocit in cross section of tubes, Re is the Renolds number defined b Re = ρur/µ. Let λ = L/R be the dimensionless miing length. Then (2.3) can be rewritten as µ t = µ ˆRλ L λdũ dφ. (2.5) Coupling (2.1) and (2.5), we have the following equation which is the application of the classical Prandtl s miing-length, dũ dφ = ˆRφ 1 + ˆRλ λ dũ. (2.6) Set ũ (φ) = dũ/dφ. Then (2.6) can be written as ũ (φ) = Using the smmetr of turbulent flow, we have dũ dφ =. φ= dφ ˆRφ 1 + ˆRλ λũ (φ). (2.7) It can be seen from Eq. (2.5) that the edd viscosit µ t at the center of a circular tube will be zero. However, the edd viscosit at the center is non-zero according to either the conception or the eperiment. This inconsistenc ma be caused mainl b the fact that the Prandtl theor onl took the first-order derivative of velocit to approimate the edd velocit (2.5). Thus, Luo et al. [15] introduced a simple wa to modif the Prandtl s miing length theor b utilizing the second-order derivative to approimate the edd velocit. B doing this, Eq. (2.5) can be changed to µ t = µ ˆRλ L λũ (φ) λ2 ũ (φ). (2.8) Combining (2.1) and (2.8) ields ũ (φ) = ˆRφ 1 + ˆRλ λũ (φ) λ2 ũ (φ). (2.9) The dimensionless miing length λ in (2.9) is found to be not a universal constant but a function of position, λ = L/R = kf(φ), f(φ), f() = 1, f(1) =, (2.1)
4 Moving Mesh Method for a Model of Turbulent Flow 391 where k is a dimensionless positive constant. For single phase flow in circular tubes, so ũ and ũ are nonpositive, the absolute sign in (2.9) can be dropped. Then (2.9) can be rewritten as [ ũ (φ) 1 ˆR(kf(φ)) 2 ũ (φ) ˆR(kf(φ)) ] 3 ũ (φ) = ˆRφ. (2.11) Since ( 1 ˆR(kf(φ)) 2 ũ (φ) ˆR(kf(φ)) 3 ũ (φ)) φ= >, and ( ũ () 1 ˆR(kf()) 2 ũ () ˆR(kf()) ) 3 ũ () =, we can derive that (2.11) satisfies the smmetr of the turbulent flow: We then obtain from (2.11) that ũ () =. (2.12) 1 ˆRk 2 f(φ) 2 ũ (φ) ˆRk 3 f(φ) 3 ũ (φ) = ˆRφ ũ, φ (, 1]. (2.13) (φ) Set a(φ) = ˆRk 2 f(φ) 2 and b(φ) = ˆRk 3 f(φ) 3. Letting φ in (2.13) gives Then we obtain that lim (1 φ a(φ)ũ (φ) b(φ)ũ (φ)) = lim ˆRφ φ ũ (φ) = (1 b()ũ ()) = ˆR ũ () ũ () = b() ˆR 2b() (using L Hospital rule). = ˆR 2 k 3 ˆRk 3. (2.14) For the turbulent flow, Re (4, 1). In (2.1), k and f(φ) have man possible choices. Because we are aimed at showing how to solve the non-linear singular model (2.11) with an optimal control condition (2.4), we use the epression proposed b Prandtl. When the Renolds number is sufficientl high, the miing length in the wall turbulence can be written as f(φ) = 1 φ, k =.2. (2.15) Then we change the problem to the following mathematical model. For an given α (4, 1), find constant c such that satisf ( () 1 k 2 c(1 2 ) 2 () 1 ) 2 k3 c(1 2 ) 3 () = c, [, 1] (2.16) and the constrained condition together with three boundar conditions 2 1 d = α c, (2.17) (1) =, () =, () = c 2 k 3 ck 3. (2.18)
5 392 Y. YANG, Y. CHEN AND Y. HUANG It is eas to show that c >. Set = 1 in (2.16), it can be verified that (1) = c. (2.19) Since () and () are nonpositive, k and c are positive, we have (), ( 1 k 2 c(1 2 ) 2 () 1 ) 2 k3 c(1 2 ) 3 (), c. (2.2a) (2.2b) (2.2c) It is then epected that we have the boundar condition () = and () has no other zero point. Moreover, = 1 is the singular point of the equation, which ma ield a boundar laer on the right end. 3. A Moving Mesh Algorithm In this section, we present an efficient and fast approach for the model (2.16)-(2.18), which includes an iterative method and an adaptive method Analsis of the model We consider the model (2.16)-(2.17). The difficult of this problem involves two parts. Firstl, c is unknown as well as ; this problem is an optimal control problem and (2.17) is a control equation. For a given α, we want to obtain the solution and c which not onl satisf the differential equation (2.16) but also the control (2.17). Secondl, (2.16) is not onl a nonlinear differential equation, but also a singular problem. In order to solve the optimal control problem, we use an iterative algorithm to compute and c. For an initial c, solve (2.16) to obtain. Then compute c b solving (2.17) with the obtained. Compare the computed c and the initial c, if the discrepanc satisfies an preassigned tolerance, stop computing, else, the obtained c becomes the initial value and repeat the process. This is an outer iteration for solving the model. A detailed algorithm will be described in the net subsection. Then we consider Eq. (2.16) with boundar conditions in (2.18). As a singular problem, there is a region in which the solution of the differential equation is steep. For (1 2 ) 1 when 1, the solution has a boundar laer near the boundar = 1. It is well known that central or upwinding differential scheme on an uniform mesh will not give a satisfactor numerical solution in this case. To obtain a reliable numerical solution for (2.16) and (2.18), it is better to use a mesh that concentrates nodes in the boundar laer. Ideall, the mesh should be generated b adapting it to the features of the computed solution and this is usuall done b equidistributing a monitor function over the domain of the problem. This adaptive methods can handle not onl boundar laer but also interior laer problems. In fact, this is a moving mesh approach that can capture the laer effectivel Numerical method and algorithm Given a partition of [, 1]: Ω N = { j = < 1 < < N = 1}.
6 Moving Mesh Method for a Model of Turbulent Flow 393 On Ω N, we discretize (2.16) and (2.18) b using a upwinding scheme as follows: D + i (1 a i D + i b i DD i ) = c i for 1 i N 1 (3.1) with N =, where the operators used are given b D v i = v i v i 1 h i, D + v i = v i+1 v i h i+1, Dv i = v i+1 v i h i, h i = i i 1, hi = (h i+1 + h i )/2, a i = k 2 c(1 i 2 ) 2, b i = 1 2 k3 c(1 i 2 ) 3. At =, since () () + () = (), we set () 1 h 1 1 h 1 = 2 1, h 1 where 1 = 1 is used. Consequentl, the discretization of (2.16) at = becomes h ( h 2 1 b 2 1 ) 1 h 2 = c. (3.2) 1 As Eq. (3.1) is a nonlinear equation, we use Newton-Raphson iteration method to solve it. The integration in (2.17) is approimated b the following quadrature formula: 2 1 d N h i ( i i + i 1 i 1 ). (3.3) Denote i () C[, 1] the piecewise linear interpolant through the knots ( i, i ). We choose a monitor function: w() = 1 + i () 2 to equidistribute mesh which is the discrete analogue of the standard arc-length monitor function w() = 1 + () 2. Note that ( i ()) = D i, I i = ( i 1, i ), 1 i N 1. In an other words, we construct a mesh b solving the following equations: ( i+1 i ) 2 + ( i+1 i ) 2 = ( i i 1 ) 2 + ( i i 1 ) 2, 1 i N 1, =, N = 1. (3.4) To equidistribute the mesh and to capture the boundar laer, we move mesh iterativel to generate a fine mesh.
7 394 Y. YANG, Y. CHEN AND Y. HUANG Main Algorithm 1. Initialize c and give mesh Ω N ; 2. Use Newton-Raphson method and call adaptive algorithm to compute the discrete solution { i } b (3.1) on Ω N ; 3. Solve c = α 2 1 d b (3.3), if c c η, goto step 5, otherwise continue; 4. Set c = c c, return to step 2; 5. Set c = c, then stop. Adaptive Algorithm 1. Initialize mesh: Ω () N 2. Solve (3.1) for { (k) i }, let l (k) i = = {() j = () < () 1 < < () N = 1}, k = ; ( (k) ( (k) i ) (k) ( (k) i 1 ))2 + (h (k) i ) 2, L (k) = 3. Test mesh: let τ be a user chosen constant with τ > 1. If mal (k) /L (k) τ/n, then goto step 5. N i=1 l (k) i. 4. Generate new mesh: choose point = (k+1) < (k+1) 1 < < (k+1) N that for each i, the distance from ( (k+1) i 1 L (k) /N. Let k = k + 1, return to step Set i = (k) i, i = (k) i, then stop., (k) ( (k+1) i 1 ) to ( (k+1) i, (k) ( (k+1) i = 1, such )) equals 4. Numerical Results This section considers numerical solutions of the model (2.16)-(2.17) to demonstrate the performance of the iterative method and the moving mesh method proposed in last section. The parameters α in (2.17) are taken as 2, 4, 6, 8, 1 respectivel. For a given α, there is a corresponding c satisfing (2.16)-(2.17). Table 4.1 and Fig. 4.1 show the relation between α and c. From Table 4.1, we conclude that c is monotonicall increasing with respect to α for a fied N and c is convergent with increasing N. Fig. 4.1 shows that the relation of α and c is almost linear. So using least-square methods we can write α and c in the form of c = kα + b with k =.79 and α = Table 4.1: The value of c for different α and different N. N α = 2 α = 4 α = 6 α = 8 α =
8 Moving Mesh Method for a Model of Turbulent Flow c Fig The relation of dimensionless radius c and the Renold number α. In this model () is the velocit distribution. Eq. (2.16) is a nonlinear convection-dominated stationar convection-diffusion problem. The solution () to (2.16) and (2.18) has a boundar laer near = 1. Our numerical solution { i } showing b Fig. 4.2 confirms that { i } is ver steep near = =2, c= =6, c= =1, c= =4, c= =8, c= =5, c= Fig Liquid velocit distribution for different Renold number α.
9 396 Y. YANG, Y. CHEN AND Y. HUANG With the boundar condition () = and since () ma have no other zero point, we change Eq. (2.16) to (1 k 2 c(1 2 ) 2 12 ) k3 c(1 2 ) 3 = c, (, 1]. Set z() = (1 k 2 c(1 2 ) 2 12 ) k3 c(1 2 ) 3 At =, b (2.18), we have z() = c () = k3 c 2. 2 = c, (, 1]. (4.1) Note that z() is the edd viscosit distribution. Let z i be the discrete form of z(). Then we have z = k 3 c 2, 2 z i = c ih i+1, for i = 1,, N 1, i+1 i z N = z(1) = c () = c =1 c = z Fig The edd viscosit distribution z for Renold number α = 6. Denote z i () the piecewise linear interpolation though the knots ( i, z i ). Fig. 4.3 presents the value of z() in the case of α = 6 with N = 124. Fig. 4.4 shows the whole histor of the mesh moving from its initial uniform mesh and the convergence histor of c from its initial value c = 1.. It is clearl demonstrated that the mesh concentrates nodes near the boundar laer. From Fig. 4.4, it is observed that for c = 1. the mesh move from the uniform mesh to the 1st level mesh; for c = 156 the mesh moves from 2nd to 3rd mesh while the 3rd mesh is the final adaptive mesh corresponding to c = 156. The mesh of 6th level is good enough for other
10 Moving Mesh Method for a Model of Turbulent Flow 397 c=66.6, 13 th N=64, k th moving mesh c=66.5, 12 th c=66.23, 11 th c=65.55, 1 th c=63.82, 9 th c=599.45, 8 th c=588.48, 7 th c=561.44, 6 th c=497.8, 5 th c=364.34, 4 th c=156.3, 3 th c=156.3, 2 th c=1., 1 th c=1., th initial mesh: uniform mesh 1 Fig Moving mesh histor. choices of c, so onl outer iterations are needed to compute the convergent c while the mesh does not need to move an more. In our algorithm, for a fied α, c can be convergent from an arbitrar initial value. Taking α = 6 for eample, the following datum shows how c = 1 and c = 1 convergence respectivel with N = 124: c = c = Our outer iterative is mainl to compute c and. Take α = 6 and N = 124 with initial c = 1. as an eample. Fig. 4.5 gives the convergent velocit of c b plotting log( c c ), where c is computed b the present iterative value, c the last iterative value. We can find that its convergence is almost linear. We close this section b making several remarks. First, in the main algorithm, we have another choice of step 4, i.e., set c = c, return to step 2. Secondl, in the Newton-Raphson iteration, we choose the initial guess for () i = and () i = 1; both give the same convergent value. Thirdl, in step 3 of the adaptive algorithm, τ can be chosen that to satisf τ 1. In this paper, we set τ = Conclusions and Future Works In this paper, following a new proposal of Luo et al. [15] the mathematical governing equations for a model for turbulent flow in circular tubes are derived. The mathematical setting involves a nonlinear singular problem with an optimal control condition. An efficient numerical strateg based on the moving mesh method is proposed. The numerical results show that
11 398 Y. YANG, Y. CHEN AND Y. HUANG log( c c ) i th outer iterative Fig Variation of log( c c ) vs. i-th outer iteration. our proposed numerical scheme is efficient in solving the nonlinear singular problem with constraints. In our future work, we will consider the theoretical aspects of the proposed method including issues on stabilit and accurac. Acknowledgments. We thank Professors Hean Luo, Pingle Liu and Xiao Luo for man helpful discussions during the preparation of this work. We would also like to thank the referees for ver constructive comments. This work is supported b NSFC Outstanding Youth Scientist Grant and NSFC grant , the National Basic Research Program of China under the grant 25CB3217(1,3), and Scientific Research Fund of Hunan Provincial Education Department. Special thanks go to Guangdong Provincial Zhujiang Scholar Award Project awarded to the second author. References [1] J.U. Brackbill and J.S. Saltzman, Adaptive gid with directional control, J. Comput. Phs, 18 (1993), [2] M. Beckett and J.A. Mackenzie, Convergence analsis of finite difference approimations on equidistributed grids to a boundar value problem, Appl. Numer. Math., 35 (2), [3] H.D. Ceniceros and T.Y. Hou, An efficient dnamicall adaptive mesh for potentiall singular solutions, J. Comput. Phs., 172 (21), [4] Y. Chen, Uniform pointwise convergence for a songularl perturbed problem using arc-length equidistribution, J. Comput. Appl. Math., 159 (23), [5] Y. Chen, Uniform convergence analsis of finite difference approimations for songular perturbations on an adaptive grid, Adv. Comput. Math., 24 (26), [6] R. Datta, Edd viscosit and velocit distribution in turbulent pipe flow revisited, AICHE J., 39 (1993), [7] S.F. Davis and J.E. Flahert, An adaptive finite element method for initial boundar value problems for partial differential equations, SIAM J. Sci. Comput., 3 (1982), [8] Y. Di, R. Li and T. Tang, A general moving mesh framework in 3D and its application for simulating the miture of multi-phase flows. Commun. Comput. Phs., 3 (28), [9] U. Frisch, Turbulence, Cambridge Universit Press, (1995).
12 Moving Mesh Method for a Model of Turbulent Flow 399 [1] W.Z. Huang, Y. Ren and R. Russell, Moving mesh partial differential equations (MMPDEs) based on the equidistribution principle, SIAM J. Numer Anal., 31 (1994), [11] A. Harten and J.M. Hman, Self-adjusting grid methods for one-dimensional hperbolic conservation laws, J. Comput. Phs., 5 (1983), [12] N. Kopteva, Maimum norm a posteriori error estimates for a one-dimensional convection-diffusion problem, SIAM J. Numer. Anal., 39 (21), pp [13] N. Kopteva and M. Stnes, A robust adaptive method for a quasi-linear one-dimensional convection-diffusion problem, SIAM J. Numer. Anal, 39 (21), [14] T. Linβ, Maimum-norm error analsis of a non-monotone FEM for a singularl perturbed reaction-diffusion problem, BIT, 47 (27), [15] X. Luo, P. Liu and H. Luo, Modeling of turbulent flow in circular tubes the application and improvement of Prandtl s miing-length theor, J. Cent. South Univ. T. (English Edition), submitted. [16] R. Li, T. Tang and P.W. Zhang, Moving mesh methods in multiple dimensions based on Harmonic maps, J. Comput. Phs., 17 (21), [17] K. Miller and R.N. Miller, Moving finite element I, SIAM J. Numer. Anal., 18 (1981), [18] Y.S.M. Morsi, P.G. Holland, B.R. Claton, Prediction of turbulent swirling flows in aismmetric annuli, Appl. Math. Model., 19 (1996), [19] T. Mizushina, F. Ogino, Edd viscosit and universal velocit profile in turbulent flow in a straight pipe, Journal of Chemical Engineering of Japan, 35 (197), [2] Z.-H. Qiao, Numerical investigations of the dnamical behaviors and instabilities for the Gierer- Meinhardt sstem. Commun. Comput. Phs., 3 (28), [21] Y. Qiu, D.M. Sloan and T. Tang, Numerical solution of a singularl perturbed two-point boundar value problem using equidistribution:analsis of convergence, J. Comput. Appl. Math., 116 (2), [22] H.Z. Tang, T. Tang and P.W. Zhang, An adaptive mesh redistribution methods for non-linear Hamilton-Jacobi equations in two- and three-dimensions, J. Comput. Phs., 188 (23),
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More informationCONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 742-761, 21. c Association for Scientific Research CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL R. Naz 1,
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 informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol., No. 6, pp. 978 994 c 999 Societ for Industrial and Applied Mathematics A STUDY OF MONITOR FUNCTIONS FOR TWO-DIMENSIONAL ADAPTIVE MESH GENERATION WEIMING CAO, WEIZHANG HUANG,
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 information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
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 information7. TURBULENCE SPRING 2019
7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport
More informationJoule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder
Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 7-3, 04. Available online at www.jafmonline.net, ISSN 735-357, EISSN 735-3645. Joule Heating Effects on MHD Natural Convection Flows in Presence of
More informationSolitary Wave Solutions of KP equation, Cylindrical KP Equation and Spherical KP Equation
Commun. Theor. Phs. 67 (017) 07 11 Vol. 67 No. Februar 1 017 Solitar Wave Solutions of KP equation Clindrical KP Equation and Spherical KP Equation Xiang-Zheng Li ( 李向正 ) 1 Jin-Liang Zhang ( 张金良 ) 1 and
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 informationThree-dimensional numerical simulation of a vortex ring impacting a bump
HEOREICAL & APPLIED MECHANICS LEERS 4, 032004 (2014) hree-dimensional numerical simulation of a vorte ring impacting a bump Heng Ren, 1, a) 2, Changue Xu b) 1) China Electronics echnolog Group Corporation
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 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 informationUNIT 4 HEAT TRANSFER BY CONVECTION
UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional
More informationFrom bell-shaped solitary wave to W/M-shaped solitary wave solutions in an integrable nonlinear wave equation
PRAMANA c Indian Academ of Sciences Vol. 74, No. journal of Januar 00 phsics pp. 9 6 From bell-shaped solitar wave to W/M-shaped solitar wave solutions in an integrable nonlinear wave equation AIYONG CHEN,,,
More informationConvex ENO Schemes for Hamilton-Jacobi Equations
Convex ENO Schemes for Hamilton-Jacobi Equations Chi-Tien Lin Dedicated to our friend, Xu-Dong Liu, notre Xu-Dong. Abstract. In one dimension, viscosit solutions of Hamilton-Jacobi (HJ equations can be
More informationConsideration of Shock Waves in Airbag Deployment Simulations
Consideration of Shock Waves in Airbag Deploment Simulations Doris Rieger BMW Group ABSTRACT When the inflation process of a simple flat airbag was simulated with the MADYMO gas flow module, the resulting
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 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 informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationChapter 4 Transport of Pollutants
4- Introduction Phs. 645: Environmental Phsics Phsics Department Yarmouk Universit hapter 4 Transport of Pollutants - e cannot avoid the production of pollutants. hat can we do? - Transform pollutants
More informationDIGITAL CORRELATION OF FIRST ORDER SPACE TIME IN A FLUCTUATING MEDIUM
DIGITAL CORRELATION OF FIRST ORDER SPACE TIME IN A FLUCTUATING MEDIUM Budi Santoso Center For Partnership in Nuclear Technolog, National Nuclear Energ Agenc (BATAN) Puspiptek, Serpong ABSTRACT DIGITAL
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 informationFluid Mechanics II. Newton s second law applied to a control volume
Fluid Mechanics II Stead flow momentum equation Newton s second law applied to a control volume Fluids, either in a static or dnamic motion state, impose forces on immersed bodies and confining boundaries.
More informationLarge Eddy Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes
Int. J Heat Mass Transfer, Vol. 76 pp. 432-446, 214 Large Edd Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes Maram Mirzaei, Ahmad Sohankar, Lars Davidson,
More informationLECTURE NOTES - VIII. Prof. Dr. Atıl BULU
LECTURE NOTES - VIII «LUID MECHNICS» Istanbul Technical Universit College of Civil Engineering Civil Engineering Department Hdraulics Division CHPTER 8 DIMENSIONL NLYSIS 8. INTRODUCTION Dimensional analsis
More informationNumerical Analysis of Heat Transfer in the Unsteady Flow of a non- Newtonian Fluid over a Rotating cylinder
Bulletin of Environment, Pharmacolog and Life Sciences Bull. Env.Pharmacol. Life Sci., Vol 4 [Spl issue 1] 215: 318-323 214 Academ for Environment and Life Sciences, India Online ISSN 2277-188 Journal
More information2015 AASRI International Conference on Industrial Electronics and Applications (IEA 2015) Study on Optimization for Grey Forecasting Model.
05 AASRI International Conference on Industrial Electronics and Applications (IEA 05 Stud on Optimiation for Gre Forecasting Model Ying Li Min-an Tang Tao Liu Min-an Tang School of Automation and Electrical
More informationSubmarine sand ripples formation in a viscous fluid: 2D and 3D linear stability analysis
Marine Sandwave and River Dune Dnamics 1 & April 4 - Enschede, the Netherlands Submarine sand ripples formation in a viscous fluid: D and 3D linear stabilit analsis V. Langlois (1) and A. Valance (1) Groupe
More informationThe Gas-assisted Expelled Fluid Flow in the Front of a Long Bubble in a Channel
C. H. Hsu, P. C. Chen,. Y. ung, G. C. uo The Gas-assisted Epelled Fluid Flow in the Front of a Long Bubble in a Channel C.H. Hsu 1, P.C. Chen 3,.Y. ung, and G.C. uo 1 1 Department of Mechanical Engineering
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics http://jipam.vu.edu.au/ Volume 2, Issue 2, Article 15, 21 ON SOME FUNDAMENTAL INTEGRAL INEQUALITIES AND THEIR DISCRETE ANALOGUES B.G. PACHPATTE DEPARTMENT
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 informationLift Enhancement on Unconventional Airfoils
Lift Enhancement on nconventional Airfoils W.W.H. Yeung School of Mechanical and Aerospace Engineering Nanang Technological niversit, North Spine (N3), Level 2 50 Nanang Avenue, Singapore 639798 mwheung@ntu.edu.sg
More informationMechanics Departmental Exam Last updated November 2013
Mechanics Departmental Eam Last updated November 213 1. Two satellites are moving about each other in circular orbits under the influence of their mutual gravitational attractions. The satellites have
More informationResearch Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free Surface Flow
Journal of Applied Mathematics Volume 2, Article ID 793653, 6 pages doi:.55/2/793653 Research Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free
More informationIntroduction to Differential Equations. National Chiao Tung University Chun-Jen Tsai 9/14/2011
Introduction to Differential Equations National Chiao Tung Universit Chun-Jen Tsai 9/14/011 Differential Equations Definition: An equation containing the derivatives of one or more dependent variables,
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 information16.5. Maclaurin and Taylor Series. Introduction. Prerequisites. Learning Outcomes
Maclaurin and Talor Series 6.5 Introduction In this Section we eamine how functions ma be epressed in terms of power series. This is an etremel useful wa of epressing a function since (as we shall see)
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 informationSound Propagation in Ducts
Sound Propagation in Ducts Hongbin Ju Department of Mathematics Florida State Universit, Tallahassee, FL.3306 www.aeroacoustics.info Please send comments to: hju@math.fsu.edu In this section we will discuss
More informationApplications of Proper Orthogonal Decomposition for Inviscid Transonic Aerodynamics
Applications of Proper Orthogonal Decomposition for Inviscid Transonic Aerodnamics Bui-Thanh Tan, Karen Willco and Murali Damodaran Abstract Two etensions to the proper orthogonal decomposition (POD) technique
More informationIntegrals along a curve in space. (Sect. 16.1)
Integrals along a curve in space. (Sect. 6.) Line integrals in space. The addition of line integrals. ass and center of mass of wires. Line integrals in space Definition The line integral of a function
More informationAnalysis of condenser shell side pressure drop based on the mechanical energy loss
Article Calorifics December 01 Vol.57 No.36: 4718475 doi: 10.1007/s11434-01-5558-4 Analsis of condenser shell side pressure drop based on the mechanical energ loss ZENG Hui, MENG JiAn & LI ZhiXin * Ke
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at.sciencedirect.com ScienceDirect Procedia Engineering 9 (14 383 388 1th International Conference on Mechanical Engineering, ICME 13 Effects of volumetric heat source and temperature
More informationStudy on Nonlinear Vibration and Crack Fault of Rotor-bearing-seal Coupling System
Sensors & Transducers 04 b IFSA Publishing S. L. http://www.sensorsportal.com Stud on Nonlinear Vibration and Crack Fault of Rotor-bearing-seal Coupling Sstem Yuegang LUO Songhe ZHANG Bin WU Wanlei WANG
More informationHigh-order finite difference methods with subcell resolution for 2D detonation waves
Center for Turbulence Research Annual Research Briefs 7 High-order finite difference methods with subcell resolution for D detonation waves B W. Wang, C.-W. Shu, H. C. Yee AND B. Sjögreen. Motivation and
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
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 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 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 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 informationA Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
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 informationBASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS
BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS J. Roubal, V. Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Abstract The distributed
More informationClosed form expressions for the gravitational inner multipole moments of homogeneous elementary solids
Closed form epressions for the gravitational inner multipole moments of homogeneous elementar solids Julian Stirling 1,2, and Stephan Schlamminger 1 1 National Institute of Standards and Technolog, 1 Bureau
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationIncompressible Viscous Flows
Incompressible Viscous Flows Q. Choose the crect answer (i) The maximum velocit of a one-dimensional incompressible full developed viscous flow between two fixed parallel plates is 6m/s. The mean velocit
More informationNUMERICAL COMPUTATION OF PARTICLE TRANSPORT IN TURBULENT FLOW FIELD IN ELECTROSTATIC PRECIPITATOR
NUMERICAL COMPUTATION OF PARTICLE TRANSPORT IN TURBULENT FLOW FIELD IN ELECTROSTATIC PRECIPITATOR Jenő SUDA, István BERTA, Gergel KRISTÓF +, Tamás LAJOS # Ph.D. Student, + Assistant, # Professor, Department
More informationL O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S
L O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S Yoël FORTERRE, Olivier POULIQUEN Institut Universitaire des Sstèmes Thermiques et Industriels (UMR 6595 CNRS) 5, rue Enrico Fermi 13453
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 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 informationThe Numerical Analysis and Experiment of Shock Processing for Bouef
Int. Jnl. of Multiphsics Volume 4 Number 4 2010 329 The Numerical Analsis and Eperiment of Shock Processing for Bouef Yusuke YAMASHITA*, Asuka ODA*, Takeshi HUJII*, Hironori MAEHARA*, M Moatamedi**, Shigeru
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 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 informationMath 21a Homework 07 Solutions Spring, 2014
Math a Homework 7 Solutions Spring, 4. valuate the iterated integral. a) Stewart.7 # 6 ) e d d d We perform the iterated integral: e d d d e d d e d [ e [ ] 4 e + 4e. Note that we ve twice done an integral
More informationExperimental Calibration and Head Loss Prediction of Tuned Liquid Column Damper
Tamkang Journal of Science and Engineering, Vol. 8, No 4, pp. 319-35 (005) 319 Experimental Calibration and Head Loss Prediction of Tuned Liquid Column Damper Jong-Cheng Wu Department of Civil Engineering,
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 informationConservation of Linear Momentum for a Differential Control Volume
Conservation of Linear Momentum for a Differential Control Volume When we applied the rate-form of the conservation of mass equation to a differential control volume (open sstem in Cartesian coordinates,
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 informationSurvey of Wave Types and Characteristics
Seminar: Vibrations and Structure-Borne Sound in Civil Engineering Theor and Applications Surve of Wave Tpes and Characteristics Xiuu Gao April 1 st, 2006 Abstract Mechanical waves are waves which propagate
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 informationInternational Journal of Scientific Research and Reviews
Research article Available online www.ijsrr.org ISSN: 2279 543 International Journal of Scientific Research and Reviews Soret effect on Magneto Hdro Dnamic convective immiscible Fluid flow in a Horizontal
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
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 informationAn ellipsoidal particle in tube Poiseuille flow
J. Fluid Mech. (217), vol. 822, pp. 664 688. c Cambridge Universit Press 217 doi:1.117/jfm.217.298 An ellipsoidal particle in tube Poiseuille flow 664 Haibo Huang 1, and Xi-Yun Lu 1 1 Department of Modern
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 informationVibration Analysis of Isotropic and Orthotropic Plates with Mixed Boundary Conditions
Tamkang Journal of Science and Engineering, Vol. 6, No. 4, pp. 7-6 (003) 7 Vibration Analsis of Isotropic and Orthotropic Plates with Mied Boundar Conditions Ming-Hung Hsu epartment of Electronic Engineering
More informationA Meshless Simulations for 2D Nonlinear Reaction-diffusion Brusselator System
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
More informationBubble curtain system for prevention of seawater intrusion in coastal aquifers
Available online at www.pelagiaresearchlibrar.com Advances in Applied Science Research, 013, 4(4):81-90 ISSN: 0976-8610 CODEN (USA): AASRFC Bubble curtain sstem for prevention of seawater intrusion in
More informationCHAPTER 2: Partial Derivatives. 2.2 Increments and Differential
CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial
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 informationSpeeding chemical reactions by focusing
Speeding chemical reactions b focusing A. M. Lacasta, L. Ramírez Piscina, J. M. Sancho, and K. Lindenberg Citation: J. Chem. Phs. 138, 144502 2013); doi: 10.1063/1.4799875 View online: http://d.doi.org/10.1063/1.4799875
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 informationThe Plane Stress Problem
. 4 The Plane Stress Problem 4 Chapter 4: THE PLANE STRESS PROBLEM 4 TABLE OF CONTENTS Page 4.. INTRODUCTION 4 3 4... Plate in Plane Stress............... 4 3 4... Mathematical Model.............. 4 4
More informationCombined tidal and wind driven flows and residual currents
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Combined tidal and wind driven flows and residual currents 1 Lars Erik Holmedal, Hong Wang Abstract : The effect of a residual current on the combined
More informationDiffusion of Radioactive Materials in the Atmosphere
American Journal of Environmental Sciences (): 3-7, 009 ISSN 3-34X 009 Science Publications Diffusion of Radioactive Materials in the Atmosphere Abdul-Wali Ajlouni, Manal Abdelsalam, Hussam Al-Rabai ah
More informationInternational Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India
International Journal of Applied Mathematics and Phsics, 3(), Jul-December 0, pp. 55-67 Global Research Publications, India Effects of Chemical Reaction with Heat and Mass Transfer on Peristaltic Flow
More informationME615 Project Report Aeroacoustic Simulations using Lattice Boltzmann Method
ME615 Project Report Aeroacoustic Simulations using Lattice Boltzmann Method Kameswararao Anupindi School of Mechanical Engineering, Purdue Universit, West Lafaette, IN, 4796, USA Lattice Boltzmann method
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationA numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 6 (2016, 152 161 Research Article A numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows Lei Hou a, Yunqing Feng a,, Lin
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 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 informationSimultaneous Orthogonal Rotations Angle
ELEKTROTEHNIŠKI VESTNIK 8(1-2): -11, 2011 ENGLISH EDITION Simultaneous Orthogonal Rotations Angle Sašo Tomažič 1, Sara Stančin 2 Facult of Electrical Engineering, Universit of Ljubljana 1 E-mail: saso.tomaic@fe.uni-lj.si
More informationResearch Article Effect of Baffles on the Sloshing in Road Tankers Carrying LPG: A Comparative Numerical Study
Mathematical Problems in Engineering Volume 2015, Article ID 359470, 9 pages http://d.doi.org/10.1155/2015/359470 Research Article Effect of Baffles on the Sloshing in Road Tankers Carring LPG: A Comparative
More informationProper Orthogonal Decomposition Extensions For Parametric Applications in Transonic Aerodynamics
Proper Orthogonal Decomposition Etensions For Parametric Applications in Transonic Aerodnamics T. Bui-Thanh, M. Damodaran Singapore-Massachusetts Institute of Technolog Alliance (SMA) School of Mechanical
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 informationVisualizing Non-Equilibrium Flow Simulations using 3-D Velocity Distribution Functions
Visualizing Non-Equilibrium Flow Simulations using 3-D Velocit Distribution Functions A. Venkattraman and A.A. Alexeenko School of Aeronautics & Astronautics, Purdue Universit, West Lafaette IN 797 Abstract.
More informationAircraft Structures Structural & Loading Discontinuities
Universit of Liège Aerospace & Mechanical Engineering Aircraft Structures Structural & Loading Discontinuities Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/
More informationHigh resolution of 2D natural convection in a cavity by the DQ method
Journal of Computational and Applied Mathematics 7) 9 6 www.elsevier.com/locate/cam High resolution of D natural convection in a cavit b the DQ method D.C. Lo a,, D.L. Young b, C.C. Tsai c a Department
More information