CHAPTER-III GENERAL GROUP THEORETIC TRANSFORMATIONS FROM BOUNDARY VALUE TO INITIAL VALUE PROBLEMS
|
|
- Adam Skinner
- 5 years ago
- Views:
Transcription
1 CHAPTER-III GENERAL GROUP THEORETIC TRANSFORMATIONS FROM BOUNDARY VALUE TO INITIAL VALUE PROBLEMS 3.1 Introduction: The present chapter treats one of the most important applications of the concept of continuous transformation groups: the numerical solution of boundary value problems as related to the class of transformation from boundary value to initial value problems. A completely new and general method will be developed based on S. Lie s continuous and compact transformation group. The concept that already discussed in detailed with some counter examples in Chapter 2 of searching for all possible groups of transformations in a similarity analysis of partial differential equations is applied here too to develop method of transformation of Boundary Value Problem to Initial Value Problem. The advancement in the analytic theories of fluid mechanics have contributed greatly to our knowledge of non-linear differential equations, both ordinary and partial. One of the classes of differential equations resulted from fluid mechanics problem is the nonlinear two point ordinary differential equations. As in other types of nonlinear equations, analytical solutions to these equations generally do not exist except for some very special cases and numerical integration methods, such as the Runge-Kutta Method, usually have to be used. The main difficulty in the numerical integration of this type of equation is due to the fact that the boundary conditions are given at two points. Thus the missing boundary condition at the initial point has to be assumed and the equation is then integrated as an initial value problem. If the assumed initial condition is correct, the solution will satisfy the boundary condition at the second point. Otherwise, another value has to be chosen and the process is repeated. This trial and error method usually takes a long computing time and the solution frequency is very sensitive to even a very small change of the assumed initial condition. A major simplification can be achieved if a transformation can be reduced to an initial value problem with all boundary conditions specified at a single point Runge-Kutta method. 70
2 A well-known example of this class of transformation is found in the solution of Blasius boundary layer equations. After a similarity transformation is made, the boundary layer equations are reduced to nonlinear third order ordinary differential equation with two boundary conditions given at the surface of the plate and third at infinity. A method of transformation was introduced by Toepfer (1912) and the problem was transformed to an initial value problem. The successful application of this type of transformations in Blasius equation makes one wonder whether it can be extended to other equations. A start in this direction was made by Klamkin (1962) in which Topfer s method was extended to many similar type of equations. The method was further extended [Na (1968)] by using the theory of transformation groups to problems whose boundary conditions are given at finite intervals. From this point of view, Toepfer (1912) and Klamkin s (1962) transformations belong to the so-called linear group of transformation. By introducing a spiral group of transformation, the method was seen to be applicable to many new types of equations. One example is found in the solution of conduction heat transfer with nonlinear heat generation, Na and Tang (1969). Although the method developed so far has been successful, the class of equations, which can be treated in this manner, is still limited. Many attempts have been made recently to extend the method and present work represents one successful result. The extended method to be given here can be applied to problems in which certain physical parameters appear in the differential equation and/or in the boundary conditions. The well-known example of fluid mechanics namely, Blasius equation is used to illustrate the entire development of present technique. As pointed out earlier, the initial value method was first introduced by Toepfer (1912), in his attempt to solve Balsius s equation in the boundary layer theory by a series expansion method. No progress was made for half a century, until 1962 when Klamkin (1962), following the same reasoning, extended the method to a wider class of problems. Major extensions were made possible only when the concept was reexamined and interpreted by Na ( ) in terms of transformation groups. 71
3 Even though the method is sometimes referred to as the method of transformation groups, perhaps the only group concept used in the method. Attempts have been made to derive groups of transformations, which can reduce boundary value problems to initial value problems by using the invariant properties of a differential equation under a socalled infinitesimal transformation group. The theory is, however, incomplete at this time and considerable work is needed before the method can be established. Until that happens, it is much simpler to view the method as a straightforward transformation. For a complete review of the various aspects of the group-theoretic method, one has to refer Ames (1972), Na and Hansen (1971), Bluman and Cole (1974), Ovsjannikov (1967), and Na (1974). Basically, method starts by defining a group of transformations. The particular transformation within this group of transformations which can reduce the boundary value problem to an initial value problem is identified by the requirement of invariance of the given differential equation under the applied group of transformation at the missing boundary conditions. The method is simple and its application is straight forward Blasius Solution in the Boundary Layer Theory: To introduce the background, which led to the method discussed in this chapter, let us consider the flow of an incompressible fluid over a semi-infinite flat plate, as shown in figure 3.1. Such a flow is usually called the boundary layer flow since the viscous effects are limited to a thin layer near the surface, across which the velocity changes from zero at the wall to the main-stream velocity U at the edge of this layer. The analysis of such a flow requires the application of the Navier-Stokes equations, which are derived from the laws of conservation of mass and of momentum. By introducing the boundary layer assumptions (i.e. the existence of a very thin layer and that u>> v), the Navier-Stokes equation become, u v x + y 0 (3.1) u u x + v u y ν 2 u y 2 (3.2) 72
4 Figure 3.1 Boundary layer over a plate 73
5 subject to the boundary conditions y 0: u v 0 y : u U Where u and v are the velocity components in the x and y direction, respectively, and ν is the viscosity of the fluid. The boundary conditions at the wall (y 0) are based o the physical condition that neither slip nor mass transfer occurs on the surface. The boundary condition at the edge of the boundary layer simply means that, mathematically, the velocity component u approaches the mainstream velocity U asymptotically. Equations (3.1) and (3.2) were solved by Blasius in As the next step, a stream function ψ is introduced such that u ψ y, v - ψ (3.3) x Besides the physical reason of introducing this function (i.e. constant ψ lines are stream lines), the mathematical significance is that the equation of continuity, Equation (3.1), is satisfied identically and the momentum equation becomes, ψ y 2 ψ - ψ 2 ψ x y x y 2 ν 3 ψ y 3 (3.4) subject to boundary conditions, y 0: ψ x ψ y 0 y : ψ y U In other words, the solution of two equations with two variables (u, v) is reduced to the solution of one equation with one variable (ψ). Blasius introduced the following transformations: η y U νx (3.5) 74
6 and f(η) ψ νx.u (3.6) This transforms equation (3.4) to the following equation: d 3 f 3 d 2 f + 1/2 f 2 0 (3.7) subject to the boundary conditions f(0) df(0) 0, df( ) 1 The class of transformations in which a nonlinear partial differential equation is transformed to a nonlinear ordinary differential equation is known as the class of similarity transformations. The foundation of the method is contained in the general theories of continuous transformation groups that were introduced and related extensively by Lie in the latter part of the last century (1884). Equation (3.7) is a boundary value problem. In order to follow the reasoning which led Toepfer (1912) to introduce a transformation of variables and transform the problem to an initial value problem, solution of equation (3.7) by series expansion should be given in detail. Let us assume the solution of equation (3.7) to be f(η) C n η n (3.8) n0 Before substituting the solution in to the differential equation, some of the coefficients can be evaluated from the boundary conditions. From the two boundary conditions on the surface, η0: f(0) 0, df(0) 0 75
7 We get C 0 C 1 0 (3.9) Next, let us introduce the notation d 2 f(0) 2 λ (3.10) Where, of course the constant λ is unknown quantity. In terms of λ, the coefficient C 2 is given by C 2 λ 2! (3.11) Equation (3.8) can now be substituted in to equation (3.7), and we get n(n-1)(n-2) C n η n-3 +1/2 C n η n n (n-1) C n η n-2 0 n3 n0 n2 or, in expanded form, [(3)(2)(1) C 3 +1/2 C 0 (2)(1) C 2 ] +[(4)(3)(2) C 4 +1/2 C 0 (3)(2) C 3 +1/2 (2)(1) C 2 C 1 ] η +[(5)(4)(3) C 5 +1/2 C 0 (4)(3) C 4 +1/2 (2)(1) C 2 C 2 + 1/2 C 1(3)(2) C 3 ] η 2 +[(6)(5)(4) C 6 +1/2 C 0 (5)(4) C 5 +1/2 (2)(1) C 2 C 3 + 1/2 C 1(4)(3) C 4 + 1/2 C 2 (3)(2) C 3 ] η (3.13) By setting the coefficients of equation (3.13) equal to zero and employing equations (3.9) and (3.11), we get C 3 C 4 0, C 7 C 6 0, C 5 λ2 5!(2) 11λ 3 C 8 8!(2) 2 (3.14) 76
8 The solution, equation (3.8), can be therefore be written as, f(η) λη2 λ 2 η 5 11λ 3 η λ 4 η 11 2! + 5!(2) 8!(2) 2 11!(2) 3 (3.15) Which still involves the unknown constant λ. In principle, the boundary condition at the second point should be applied to determine the constantλ. However, this can not be done directly here due to the fact that boundary condition at the second point is given at η. Blasius (1950) used a matching method, in which the solution given by equation (3.15) applied to the region close to the wall is matched to an asymptotic expansion for largeη. The process is tedious and is of no interest here, By arranging the variables, equation (3.15) can be written as, (λ -1/3 f) (λ 1/3 η) 2 (λ 1/3 η) 5 11(λ 1/3 η) 8 375(λ 1/3 η) 11 2! 5!(2) + 8!(2) 2 11!(2) (3.16) This suggests a transformation of the form F λ -1/3 f, ξλ 1/3 η (3.17) Equation (3.16) then becomes, F ξ 5 11ξ ξ 11 2! +.. 5!(2) 8!(2) 2 11!(2) 3 ξ 2 (3.18) The significance of the transformation of variables will become obvious of the differential equation (3.7) and its boundary conditions are transformed. From equation (3.7), we get d 3 F dξ 3 + 1/2 F d 2 F dξ 2 0 (3.19) The boundary conditions at the wall are transformed to ξ 0: F(0) 0, df(0) dξ 0 (3.20) 77
9 Using the transformation, equation (2.20) gives d 2 F(0) dξ 2 1 (3.21) The transformed equation (3.19), subject to the boundary conditions (3.20) and (3.21), is an initial value problem. The solution of equation (3.19) can be obtained by a forward integration scheme. The relation between (η, f) and (ξ, F) is defined in equation (3.17). If the constant λ can be found, the solution of (η, f) can be calculated from the solution of (ξ,f). To this end, the boundary condition at infinity df( ) 1 is transformed according to equation (3.17) and we get, λ -2/3 df( ) dξ 1 or λ df( ) dξ - 3/2 (3.22) Thus, the solution of equation (3.7) consists of three steps. 1. Solve equation (3.19), subject to the boundary conditions (3.20) and (3.21), as an initial value problem by forward integration. In particular, we get the value of df( ) /dξ. 2. Compute λ by equation (3.22) 3. Compute f(η) by using equation (3.17). The above demonstrates an important idea of solving boundary value problems, that of solving a complementary initial value problem. In this way the solution becomes non iterative. 78
10 While the value of such a transformation was obvious, no further work was done to extend the idea to other types of problems until 1962, when Klamkin (1962) published an extension of this method by following the same reasoning as given in Toepfer s method. For example, it was found that for a general second order ordinary differential equation of the form N A i i 1 d 2 m f i 2 df n i f r i η s i 0 (3.23) Subject to the boundary conditions f(0) 0, df( ) k a transformation F λ α -1 f(η), ξ λ α η (3.24) enables the solution of equation(3.23) to be obtained by first solving the initial value problem, N A i i 1 d 2 F dξ 2 m i df dξ n i ξ s i F r i 0. (3.25) subject to the boundary conditions, df(0) F(0) 0, 1 dξ and then calculating the constant λ from λ k / F ' ( ) (3.26) Finally, the solution of equation (3.23) can be calculated by using equation (3.24). 79
11 The constant α in the transformation (3.24) is determined by requiring that, upon transformation, equation (3.25) should be independent of A. Thus, from the above theory, we can conclude that although Klamkin s method represents an important extension of the technique, the method is still quite limited for the following three reasons. First, only the types of transformations similar to those given above are known to be applicable. Second, the boundary conditions at the second point have to be given at infinity. Third, the boundary conditions at the initial point must be homogeneous. To remove some of the limitations, the present method has to be interpreted in terms of the concepts of transformation groups. All these limitations are still severe. T.Y. Na ( ) has tried to remove above limitations by introducing group transformations. He found partly successful. This is because he has explained his views through same particular examples and he couldn t derive same general mathematical theory. Actually, the method of Klamkin (1962) can be applicable to ordinary differential equations or system of ordinary differential equations, which were invariant under certain groups of homogeneous linear transformations. Additionally, the boundary conditions were specified at origin and at infinity were found homogeneous at origin. Subsequently in two papers by Na ( ), it was shown that the method was applicable to finite intervals and also to equations, which were invariant under other groups of transformations. It was also stated by Na (1967) that the boundary conditions at the initial point must be homogeneous at the initial point for the method to be applicable. Although all the boundary conditions considered by Klamkin (1962), Na ( ) were homogeneous at the initial point, it will be shown subsequently that this condition is not always necessary. 80
12 In all the boundary value problems and initial value problems to be considered, we are tacitly assuming the existence and uniqueness of the latter, which in turn will have implications for the former. In the next section, the general theory of reduction of Boundary Value Problem to Initial Value Problem is given Second-Order Boundary Value Problems: First we consider general second order differential equation: A mnrs y'' m y' n y r x s 0 (3.27) m,n,r,s subject to the boundary conditions, y ' (0) ay(0) + b, y e ( ) k (3.28) m, n, s,r are arbitrary indices, A mnrs are arbitrary constants and (e) is an arbitrary integer. If (3.27) is multi valued for y'', we assume a particular branch is specified. We now assume that y can be expressed in the form y λ F(μ x) (3.29) Where F(x) also satisfies (3.27) for arbitrary constants λ, μ, but subject to the initial conditions F(0) F (0) 1 In order that both y and F(x) satisfy (3.27), the D.E must be invariant under the twoparameter group of homogeneous linear transformations. x 1 μ x, y 1 y λ The condition that this implies on the indices m,n,r,s, got by substituting (3.29) in to (3.27) i.e. A mnrs F n (μ x) m F (μ x) n F (μ x) r (μ x) s λ c μ d 0 m,n,r,s where, c m + n + r d 2m + n s (3.30) 81
13 Consequently, c & d must be constant for all sets of indices m, n, r, s and then (3.27) reduces to A mn y m y n y r x s 0 (3.31) m,n Where, r and s are given by (3.30) It now follows that y(0) λ, y (0) aλ + b λμ (3.32) K λ μ e F (e) ( ) We now solve the initial value problem for F and determine F (e) ( ). Then λ and μ are determined from the simultaneous equations in (3.32) which then give y(0) and y (0), Thus we have converted the original Boundary value problem in to similar Initial value problems or else as noted earlier, y can be determined from y λf (μ x). Having tacitly assumed the existence and uniqueness of F(x) in (0, ), the existence of y depends on the existence and uniqueness of μ and λ. Eliminating λ in (3.32), we have μ e k' (μ -a) Thus, depending on the relative values of e, k and a, there can be zero, one, two or three solutions for μ and the same corresponding for y. The case for a finite interval instead of an infinite internal for the preceding problem can be treated in a similar fashion. If the boundary condition was given by y (e) (L) k We would then have y ' (0) aλ + b λ μ, K λ μ e F (e) (μl) again there are two equations for λ and μ but not quite so simple as before. Although these equations can be easily solved numerically, it may be easier to start out with the original system, pick a trial value for λ and then refine it by interpolation techniques. For higher order equations, the method will more useful. 82
14 If instead of boundary conditions (3.28), we had y(0) a, y (e) ( ) k, We would assume that y could be expressed in the form y F (μ x), where, F(x) also satisfy (3.27) for arbitraryμ. This entails that in (3.27) 2m + n s Constant. Since, y(0) a, F(0) a If we now let F (0) 1, then y (0) k, which is determined from K μ e F (e) ( ) Provided e 0. For the case e0, we have an anomaly that could be due to impossible boundary conditions. An example, consider the differential equation (xd-1) (xd +1) y 0 Whose solution is yax + (B/x) For a finite interval, with boundary conditions y (e) (L) k instead of y (e) ( ) k, we can proceed as before. We could also replace boundary condition at the terminal point by either of the more general ones or A i y ( i ) (L) K i lim x A i x (e i ) ( i ) y (x) K i 3.4 Transformations of Boundary Conditions: There are Boundary Value Problems where the previous method will not apply directly unless the Boundary Conditions are first transformed into a suitable form. As an example, consider the Boundary Value Problems: x n / (n+1) y y (2n+1) / (n+1) y (0) 1 (3.33) y ( ) 0 83
15 For n1, we get the Thomas-Fermi equation (1960) which had arisen in the determination of the effective nuclear change in heavy atoms. It follows that if F(x) is a solution of the above Differential Equation, so also is y λ (n+2)/n F(λ x) even though the equation has a one-parameter group of transformations, the method does not work. It would work, if the boundary conditions were interchanged to y(0) 0, y ( ) 1 This can be done by letting x 1/x. Since y '' x 4 y '' + 2xy ' It might appear as if we will lose the one-parameter group of transformations. Fortunately, if for any differential equation where F(x) and λ α F(λx) are both solutions, one makes the transformation x x -r, then if G(x) is a solution so also is λ α G(x/λ 1 /r ). Similar results apply if we transform y y s. Carrying out the transformation x 1/x on (3.33), We get x -n/ (n+1) { x 4 y '' + 2x 3 y ' } y (2n+1)/ (n+1) (3.34) y(0) 0, y( ) 1 Now, letting y λ (n+2)/ n F(x/λ) with F(0) 0, F ' (0)1, We get λ {F( )} -n/(n+2) It should be noted that the differential equation is singular at the origin. Consequently, in order to start the numerical solution for the Initial Value Problem, first one has to find the asymptotic solution in the neighbourhood of the origin. Now in the next section using simple group transformation, we show thoroughly how to transform above boundary value problem in to initial value problem. 84
16 3.5 Initial Value method by group of Transformation: The Blasius Equation in Boundary Layer Flow: The Blasius Equation, equation (3.7), which Toepfer (1912) treated, will now be reconsidered from the point of view of transformation groups. The equation is d 3 f 3 + 1/2 f d 2 f 2 0 (3.35) Subject to the boundary conditions, f(0) df(0) 0, df( ) 1 Let the given group of transformations be the linear group defined by η A α 1 η, f A α 2 f (3.36) Where, A is the parameter of transformation and α 1 and α 2 are constants to be determined. Under this group of transformations, equation (3.1) becomes d 3 f 3 + 1/2 f d 2 f 2 + 1/2 A(2α 2-2 α A α 2-3α 1 d 3 f 1) 3 d 2 f f 2 0 (3.37) We first require that the transformed differential equation be independent of the parameter of transformation A. This requires that the powers of A in the two terms be equal, i.e. α 2-3α 1 2α 2-2α 1 (3.38) Equation (3.3) gives d 3 f + f d 2 f d 3 f ½ /2 f d 2 f Aα2-3α From equation (3.39), we get d 3 f + 1/2 0 3 f d 2 f 2 (3.39) (3.40) 85
17 Fig. 3.2 Solution of equation
18 In the second step, the missing boundary condition is set equal to the parameter of transformation, i.e., d 2 f(0) (3.41) 2 A This step is somewhat arbitrary, but it does not place any limitation on the method since A is an unknown constant to begin with. Under the linear group of transformations, this condition becomes A α 2-2α 1 A d 2 f(0) 2 (3.42) and we further require that the transformed boundary condition be independent of A. This is possible if α 2-2α 1 1 (3.43) Equation (3.42) then gives, d 2 f(0) 1 2 Equations (3.38) and (3.43) give (3.44) α 1 - α 2-1/3 (3.45) Finally, the parameter A can be evaluated from the boundary condition at infinity, which gives df( ) A α 2-α 1 1 (3.46) Or, using equation (3.45), -3/2 A df( ) 1 (3.47) The other boundary conditions are transformed simply to 87
19 f (0) df(0) 0 (3.48) Equation (3.40), together with the boundary conditions (3.44) and (3.48), constitutes an initial value problem, the numerical solution of which is shown in figure 3.2. df It is seen that as η approaches infinity approaches the value of , which is taken as df( ). Substituting in to equation (3.47), we get -3/2 A df( ) Which is the missing boundary condition at η 0, as defined by equation (3.41). Since α 1, α 2 and A are all known, the solution of the original equation, equation(3.35) can be completed by using equation(3.36). It will be of interest to give a physical interpretation of fig-3.3, which shows df/ as a function of η. From equation (3.3) and (3.6), we get u ψ y y νxu f νxu df η y ν xu U df νx U df or u U df (3.49) 88
20 Fig. 3.3 Solution of equation
21 Thus, the curve given in figure 3.3 is a plot of the velocity distribution as a function of η. Referring to figure 3.1, if the fluid and the mainstream velocity U are given and the velocity distribution at a distance x x 0 is needed, equation (3.5) then gives the relation between η and y, namely η y U νx 0 (3.50) Since, U, v and x 0 are known; equation (3.50) gives the corresponding value of η for any y. The curve in fig.3.3 can then be used to read u/ U at this value of y. The method can now be summarized as consisting of the following steps. 1. A group of transformations is defined. From the point of view of transformation groups, it is immediately recognized that the class of transformation treated by Toepfer (1912) and Klamkin (1962) is the familiar linear group of transformations. The immediate conclusion is that there must be other groups. Questions as to whether there is a deductive way to derive the particular group of transformations may also arise (1968). At the present time, these theories are still incomplete and we have to proceed with the method based on a given group of transformations. 2. The differential equation is required to be independent of the parameter of transformation, which should lead to only one relation between α 1 and α 2. This step offers the first test as to whether or not the method is applicable. For example, the Falkner-Skan equation, d 3 f 3 + 1/2 f d 2 f + 2 β 1 - df 2 0 is independent of the parameter of transformation in the linear group defined in equation (3.36) only for α 1 α 2 0. The method thus fails to apply. 90
22 3. The unknown boundary condition is set equal to the parameter of transformation. Another relation connecting α 1 and α 2 can be written from which α 1 and α 2 can be solved, leaving A as the only unknown in the transformation. 4. The parameter of transformation is determined from the boundary condition at the second point such as in equation (3.47). One interesting case will be when the boundary condition at infinity is homogeneous, i.e. when the right hand side of equation (3.46) is zero. One example is found in the flow of a laminar Non- Newtonian jet (1968). Another case is when the boundary condition at the second point is given at a finite distance instead of at infinity. 5. As a last step, the boundary conditions at the initial point are transformed. Here, it is required that they be homogeneous. Otherwise, the method can not be applied. For example, if f(0) k, then upon transformation. A α2 f (0) k Which is independent of A only if α 2 is equal to zero. At this point, a work by Klamkin (1970) should be cited. In his work, Klamkin tried to remove the requirement that the boundary condition at the initial point be homogeneous. To illustrate his reasoning, let us consider one example given in Klamkin (1970). d 2 f 2 + 1/η df α + β 1 f 2 We will, however, write the boundary conditions as df(0) af(0) + b, f( ) 1 After a linear transformation, equation (3.36), the differential becomes A α 2-2α 1 d 2 f 2 + 1/ η df α + β A -2α 2 1 f 2 We will consider three cases. 91
23 Case:1 If neither α nor β is zero, invariance of the differential equation gives two equations, namely α 2-2α 1 0, 2α 2 0 Which means the problem can not be transformed to an initial value problem since now both α 1 and α 2 are zero. Case:2 If α equals to zero (i.e., the differential equation does not contain constant term), then one equation results from the invariance of the differential equation, namely α 2-2α 1-2α 2 and the method presented above can be applied. Since we now have one equation for the solution of α 1 and α 2, only one equation relating α 1 and α 2 should come from the boundary conditions. This is the reason why the boundary condition at the initial point should be homogeneous, since the boundary condition at the second point always gives one such relation. Case:3 If both α and β are zero (i.e. for differential equations whose transformed form, like the two terms on the left-side of the transformed equation above, has exactly the same powers of A, which can then be factored out for all values of α 1 and α 2 without the necessity of imposing any relation between α 1 and α 2 ), then no relation between α 1 and α 2 will emerge from the differential equation. This leaves the two boundary conditions to provide the two equations necessary for the solution of α 1 and α 2. We now set f (0) A The first boundary condition gives By outing A α 2-α 1 df(0) aa + b f(0) 1 and df(0) 1 We get two equations, α 2 1, A α2-α1 Aa + b 92
24 Which together with the transformed boundary condition at the second point, A α 2 1 f( ) are the equations required for the solution of α 1, α 2 and A. The above cases show clearly that by placing limitations on the form of differential equation, some of the restrictions imposed on the method can be removed. 3.6 Development of General Procedure: Consider general second-order differential equation N i 1 A i d 2 y 2 m i dy n i y r i x s i 0 (3.51) Subject to the two cases of boundary conditions Case I d d y( ) y(0) 0, k d Case II y(0) 0, d d y(l) d 0 Consider the linear group of transformations x B β 1 x ; y B β 2 y (3.52) Under this group of transformation, Eq. (3.51) becomes N B m i (β 2-2β 1 ) + n i (β 2-2β 1 ) + r i β 2 + s i β 1 i 1 N d x m y i dy A i 2 n i x s i y r i 0 (3.53) i 1 93
25 Equation (3.53) will be independent of the parameter of transformation, B, if the powers of B in each term are equal, i.e. m 1 (β 2-2β 1 ) + n 1 (β 2 - β 1 ) + r 1 β 2 + s 1 β 1 m i (β 2-2β 1 ) + n i (β 2 - β 1 ) + r i β 2 + s i β 1 (3.54) Where i 2,..., N. In general, (3.54) gives (N-1) equations with only two unknowns, β 1 and β 2. The method is applicable only if the (N-1) equations actually reduce to one independent equation. To illustrate the problem that may arise, the Falker-Skan equation may be cited. d 3 f 3 + 1/2 f d 2 f + 2 β 1 - df 2 0 (3.55) Under the linear group of transformation defined by Eq. (3.52), Eq. (3.55) becomes d 3 f d 2 f B β β 0 2-3β B 2β 2-2β 1 1/2 f B 2β 2-2β 1 df 2 (3.56) Which is independent of the parameter of transformation if β 2-3β 1 2 β 2-2β 1 0 2β 2-2β 1 (3.57) Two independent equations are obtained from (3.57). As a result, we get β 1 β 2 0 which means the method is inapplicable. Assuming for now that such a situation does not exist, Eq. (3.53) becomes N i 1 A i d 2 y 2 m i dy n i x s i y r i 0 (3.58) With one relation between β 1 and β 2 obtained from Eq. (3.54), the other relation required for the solution of β 1 and β 2 can be obtained by putting the slope at x 0 equal to the parameter of transformation, B, i.e., dy(0) B After transformation, we have: (3.59) β2- β1 B dy (0) B 94
26 Which is independent of B if β 2 - β 1 1 (3.60) The transformed boundary conditions are therefore, y (0) 0, dy(0) 1 (3.61) Equations (3.54) and (3.60) give solutions to β 1 and β 2. Finally, the value of B can be found by applying the boundary condition at the second point. Thus: Case I: or, B β 2 - dβ 1 d d y ( ) d k B k d d y ( ) d 1/ (β 2 - dβ 1 ) (3.62) Case II: d d y d 0 at B β1 x L Or, L 1/ β 1 B x where d d y d 0 (3.63) Thus, Eq. (3.58) is solved with the boundary conditions given in Eq. (3.61) and the value of x where d d y 0 can be found from the solution of the transformed equation. This d result is then substituted in to Eq. (3.63) and the value of B computed. It should be noted that there are cases where additional problems may arise. 95
27 To illustrate the practical application of the above theory, we now consider the following boundary value problem: d 2 T 2 - λ 2 0 T With the boundary conditions T (0) 0, The solution is dt(0) 1 T 1/λ sin h λ x This is never zero. This places another limitation on the method. The boundary conditions at x L in case II need not be homogeneous. For example, one may have at x L, d d y d K Thus, B β 2 - dβ 1 d d y d K at B β 1 x L (3.64) Since, K, L, β 1 and β 2 are known constants, the value of B can be found by searching for d d y values of x and in the solution of Eq. (3.58) which give the same value of B in d both equations of Eq (3.64). One way of doing this is by eliminating B in Eq. (3.64) which leads to d d y K d x d L x d- β 2 /β 2 (3.65) Next, (d d y/d x d ) Vs. x can be plotted as a curve. Another curve from the solution to Eq. (3.58) can be plotted with the same co-ordinates. The intersection of these two curves will give the required value of x and d d y/d x d which in turn can be used to compute B from Eq. (3.64). 96
28 One final remark about the method is necessary. Suppose the boundary condition at the initial point is dy(0)/ 0. In this case, we merely have to put Thus, y(0) B B β2 y (0) B, And if the result is to be independent of B, β 2 must be equal to 1. Under no circumstances, however, should the boundary condition at the initial point be non homogeneous. If it is, one more equation relating β 1 and β 2 will result. The method then cannot be applied Spiral Group of Transformations: We now consider a class of nonlinear ordinary differential equations in which a spiral group of transformation rather than a linear group is needed for the method to apply. We consider here the class of equations N i 1 C i d 2 y 2 m i dy n i x q i e p iy 0 (3.66) With the boundary conditions Case I: x 0 : Case II: x 0 : dy dy 0 ; x 1 : y 0 0 ; x : y k 1 Where C i, m i, n i, p i and q i are constants and N is the number of terms in Eq. (3.66) Let us define the one-parameter spiral group of transformations x e β 1A x ; y y + α2 A (3.67) 97
29 where A is the parameter of transformation and α 1 and α 2 are constants to be determined. Under this group of transformation, Eq. (3.66) becomes N i 1 C i e (-2m iα 1 n i α 1 + p i α 2 + qiα 1 )A x d 2 y 2 m i dy n i e piy i x qi 0 (3.68) The equation is seen to be independent of the parameter of transformation, A, if the powers of e in each term are equal, i.e., (-2m i n i + q i ) α 1 + p i α 2 (-2m 1 n 1 + q 1 ) α 1 + p 1 α 2 (3.69) where i 2,..., N. The transformed equation becomes N i 1 C i d2 y 2 m i dy n i e piy x qi 0 (3.70) Equation (3.69) represents (N-1) equations. In general, the method can be applied only if one independent equation results from these (N-1) equations. For example, if Eq. (3.66) takes the form d 2 y x dy + e y 0 (3.71) Equation (3.69) then gives one independent equation for α 1 and α 2 as -2α 1 α 2 (3.72) Physically, Eq. (3.71) may be interpreted as the equation for heat conduction in spheres with exponential heat generation. To determine the second relation for the solution of α 1 and α 2, we put y(0) A (3.73) Upon transformation, this condition becomes: y (0) + α 2 A A 98
30 Which is seen to be independent of A if α 2 1 (3.74) The transformed boundary conditions become dy (0) y (0) 0 and 0 (3.75) For the example given in Eq. (3.71), α 1 and α 2 can be found from Eqs. (3.72) and (3.74): α 1-1/2 and α 2 1 Finally, to get the parameter of transformation, A, the boundary condition at the second point is used. The two cases are considered separately. Case I. The boundary condition at x 1 becomes y + A 0 at e α 1 A x 1 (3.76) Eliminating A, we get e - α1 y x 1 (3.77) Case II. The boundary condition at x becomes Or x : y + α 2 A k 1 A 1/α 2 [ k 1 - y ( ) ] (3.78) Therefore, the method proceeds as follows: First, Eq. (3.70) is solved with the boundary conditions (3.75). In case I, the solution curve to Eq. (3.70) can be plotted on y vs. x co-ordinates. Equation (3.77) is plotted on the same co-ordinates. The intersection of these two curves gives the values of y and x which give the same value of A from Eq. (3.76). The value of A is then determined. In case II, the value of A can be computed from Eq. (3.78). 99
31 The solution to the original equation, Eq. (3.66) can be obtained from Eq. (3.67) since now α 1, α 2 and A are known constants Two Parameter Method: Transformation of Two Boundary Conditions: The method developed in the preceding section can be extended to higher order differential equations as long as only one boundary condition is required to be transformed. In this section, the method will be extended to higher-order differential equation in which more than one boundary condition need to be transformed. For such cases, multi parameter groups of transformations are required. Consider now the third order differential equation N i 1 Subject to the boundary conditions Case I. Case II. A i d 3 y 3 m i d 2 y 2 y s i x t i We now define a two-parameter group of transformation n i dy d d 1 y( ) k 1, d d d 2 y( ) y(0) 0, k 1 2 d 2 d d 1 y(l) 0, d d d 2 y(l) y(0) 0, k 1 2 d 2 r i 0 (3.79) x B β 1 C γ1 x ; y B β2 C γ2 y (3.80) Under this group of transformation, Eq. (3.79) becomes N i 1 A i B m i (β-3β 1) + n i (β 2-2β 1 ) + r (β i 2 -β 1 ) + s β i 2 + t β i 1 x C m i (γ 2-3γ 1 ) + n i (γ 2-2γ 1 ) + r i (γ 2-γ 1 ) + s i γ 2 + t i γ 1 X d 3 y 3 m n i d 2 r y i dy i y s i t x i 0 2 (3.81) 100
32 The method can be applied if, for all i s, m i (β 2-3β 1 ) + n i (β 2-2β 1 ) + r i (β 2 -β 1 ) + s i β 2 + t i β 1 C 1 m i (γ 2-3γ 1 ) + n i (γ 2-2γ 1 ) + r i (γ 2 -γ 1 ) + s i γ 2 + t i γ 1 C 2 (3.82a) (3.82b) Where C 1 and C 2 are two arbitrary constants. Equation (3.81) then becomes N i 1 A i d 3 y 3 For example, the equation m n i d 2 r y i dy i y s i t x i 0 2 (3.83) d 3 y x 3 dy d 2 y 2 + y x (3.84) belongs to this class. The boundary condition at the initial point y(0) 0, can be transformed to y (0) 0 (3.85) To get the other boundary conditions at the initial point, let us put dy(0) B, and d 2 y(0) C 2 (3.86) Upon transformation, (3.86) becomes B β 2- β 1 C γ 2- γ 1 dy (0) B (3.87a) and B β 2-2β 1 C γ 2-2γ 1 d 2 y (0) 2 C (3.87b) Which are seen to be independent of B and C if, from (3.87a), 101
33 β 2 -β 1 1, γ 2 - γ 1 0 and, from (3.87b), β 2-2β 1 0, γ 2-2γ 1 1 (3.87a) and (3.87b) are then transformed to dy (0) 1, d 2 y (0) 1 2 (3.88) and the values of β 1, β 2, γ 1 and γ 2 β 1 1, β 2 2, γ 1-1, γ 2-1 To get the parameters of transformation, B and C, the boundary conditions at the second point are used. Case I B 2-d 1 C d 1-1 k 1 d d 1 y ( ) d 1 (3.89a) B 2-d 2 C d 2-1 k 2 d d 2 y ( ) d 2 (3.89b) Case II BC -1 L x where d d 1 y d 1 0 (3.90a) B 2-d 2 C -1+d 2 k 2 d d 2 y (3.90b) d 2 102
34 Therefore, B and C can be solved from Eqs. (3.89) or (3.90). The method can easily be extended to equation of the type: N i 1 A i d 3 y 3 m i d 2 y 2 n i e s iy x t i The only difference here is that one assumes a transformation group defined by dy r i 0 (3.91) x e β 1B+γ1C x ; y y + β 2 B + γ 2 C (3.92) Other steps remain the same Simutaneous Differential Equations: Application of the method to simutaneous differential equations again involves multi-parameter groups. Consider now the following system of two simutaneous equations: N i 1 A i d 2 y 2 m i dy n i y p i d 2 z 2 r s i dz i 0 z t i x q i (3.93a) M i 1 Bj d 2 y 2 m j dy n j y p j d 2 z 2 r s j dz j z t j x q j 0 (3.93b) Subject to the boundary conditions Case I. y(0) 0, z(0) 0, d d y( ) d k 1, d d z( ) d k 2 Case II. y(0) 0, z(0) 0, d d y(l) d k 1, d d z(l) d k 2 103
35 Let us now define a two-parameter transformation group x λ β 1 x, y λ β2 y, z λ β 3 μ δ z (3.94) Again, the differential equation, (3.93), are independent of the parameters of transformation, λ and μ, if the powers of λ and μin each term are respectively the same. This leads to the following system of equations for the solution of β 1, β 2, β 3 and δ: (β 2-2β 1 ) m i + (β 2 - β 1 ) n i + β 2 p i + r i (β 3-2β 1 ) + s i (β 3 -β 1 ) + t i β 3 + q i β 1 (β 2-2β 1 ) m 1 + (β 2 - β 1 ) n 1 + β 2 p 1 + r 1 (β 3-2β 1 ) + s 1 (β 3 -β 1 ) + t 1 β 3 + q 1 β 1 (3.95) r i + s i + t i r 1 + s 1 + t 1 (3.96) (β 2-2β 1 ) m j + (β 2 - β 1 ) n j + β 2 p j + r j (β 3-2β 1 ) + s j (β 3 -β 1 ) + t j β 3 + q j β 1 (β 2-2β 1 ) m 1 + (β 2 - β 1 ) n 1 + β 2 p 1 + r 1 (β 3-2β 1 ) + s 1 (β 3 -β 1 ) + t 1 β 3 + q 1 β 1 (3.97) r j + s j + t j r 1 + s 1 + t 1 (3.98) where i 2,..., N and j 2,...,M. Substitution of Eqs. (3.96) and (3.98) into Eq. (3.95) and (3.97), respectively, gives (m i + n i + p i ) β 2 - (2m i + n i + 2r i + s i - q i ) β 1 (m 1 + n 1 + p 1 ) β 2 - (2m 1 + n 1 + 2r 1 + s 1 - q 1 ) β 1 (3.99) ( m j + n j + p j ) β 2 - ( 2 m j + n j + 2 r j + s j - q j ) β 1 ( m 1 + n 1 + p 1 ) β 2 - (2 m 1 + n r 1 + s 1 - q 1 ) β 1 (3.100) 104
36 The method is applicable if Eqs. (3.99) and (3.100) each represent only one independent equation and that both give the same ratio of β 2 /β 1. If these conditions are satisfied, the ratio of β 2 /β 1 is known. Next, the required boundary conditions are defined to be equal to λ and μ respectively, i.e., y ' (0) λ and z ' (0) μ upon transformation, λ β2- β1 y ' (0) λ and μ δ λ β3- β1 z ' (0) μ Which then give, β 2 -β 1 1, δ 1 and β 3 -β 1 0 (3.101) The ratio of β 2 /β 1 obtained from Eqs. (3.99) and (3.100), together with Eq. (3.101), gives solutions of β 1, β 2, β 3 and δ. To get the parameters of transformation, the same method discussed in previous paragraphs can be applied. We will not repeat here. The method can be easily generalized to include cases with exponentials of y or z or both in Eqs. (3.93). 3.7 More General Types of Equations: The method developed above can be extended to more general types of equation. Two cases are considered here. Consider now the general second order differential equations: G M j 1 b j d 2 y 2 t j dy s j y r j x u j 0 (3.102) Where G represents an arbitrary function of the argument indicated. Case I. The method can be applied if, under the linear transformation group x A α 1 x, y A α 2 y, 105
37 only one relation between α 1 and α 2 is obtained from the condition that Eq. (3.102) is invariant under this group of transformation. As an example, the equation y 3 d 2 y 2 + sin y dy will give one relation between α 1 and α 2, namely 2α 2 - α 1 0 Case II. If y is absent in all terms in Eq. (3.102), the spiral group of transformation can always be applied for any arbitrary function in Eq. (3.102). As an example, consider d 2 y 2 + f 1 dy (3.103) where f 1 is any arbitrary function of dy/. Under the spiral group of transformation. x e α 1 a x, y y + α 2 a Eq. (3.103) becomes e -2α1 a d 2 y 2 + f 1 e -α1 a dy which is independent of a if α 1 0 for any arbitrary function f 1. The remaining steps remain the same. 3.8 Conclusion : The application of a linear or spiral group of transformation to the class of transformation from a boundary value to an initial value problem is treated. The method consists of three basic steps. First, a transformation group is defined and the given differential equation is required to be invariant, i.e. independent of the parameter of transformation, under this group of transformation. In step 2, the required boundary 106
38 condition is set to be equal to the parameter of transformation. Finally, the parameter of transformation is found by using the boundary condition at the second point. Knowing this general concept, the method treated can be applied to higher-order equations or other types of equations. It is simple to apply and only algebraic solutions are required to get the transformation. The main disadvantage of this method lies, however, on the arbitrariness in the selection of a proper group for a given differential equation. 107
Chapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationChapter 2 Non-linear Equations
Chapter 2 Non-linear Equations 2.1 First Examples of Non-linear Ordinary Differential Equations Before describing some general properties of non-linear equations, we find it useful to present some well-known
More informationSymmetry Solutions of a Third-Order Ordinary Differential Equation which Arises from Prandtl Boundary Layer Equations
Journal of Nonlinear Mathematical Physics Volume 15, Supplement 1 (2008), 179 191 Article Symmetry Solutions of a Third-Order Ordinary Differential Equation which Arises from Prandtl Boundary Layer Equations
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationBoundary value problem with integral condition for a Blasius type equation
ISSN 1392-5113 Nonlinear Analysis: Modelling and Control, 216, Vol. 21, No. 1, 114 12 http://dx.doi.org/1.15388/na.216.1.8 Boundary value problem with integral condition for a Blasius type equation Sergey
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationFURTHER SOLUTIONS OF THE FALKNER-SKAN EQUATION
FURTHER SOLUTIONS OF THE FALKNER-SKAN EQUATION LAZHAR BOUGOFFA a, RUBAYYI T. ALQAHTANI b Department of Mathematics and Statistics, College of Science, Al-Imam Mohammad Ibn Saud Islamic University (IMSIU),
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
More informationLECTURE 10: REVIEW OF POWER SERIES. 1. Motivation
LECTURE 10: REVIEW OF POWER SERIES By definition, a power series centered at x 0 is a series of the form where a 0, a 1,... and x 0 are constants. For convenience, we shall mostly be concerned with the
More informationCHAPTER-7 SIMILARITY SOLUTIONS OF THE THREE DIMENSIONAL BOUNDARY LAYER EQUATIONS OF A CLASS OF GENERAL NON-NEWTONIAN FLUIDS
CHAPTER-7 SIMILARITY SOLUTIONS OF THE THREE DIMENSIONAL BOUNDARY LAYER EQUATIONS OF A CLASS OF GENERAL NON-NEWTONIAN FLUIDS 7.0. INTRODUCTION: The subject of boundary layer flows of non-newtonian fluids
More informationSolution of Partial Differential Equations
Solution of Partial Differential Equations Introduction and Problem Statement Combination of Variables R. Shankar Subramanian We encounter partial differential equations routinely in transport phenomena.
More information4.2 Concepts of the Boundary Layer Theory
Advanced Heat by Amir Faghri, Yuwen Zhang, and John R. Howell 4.2 Concepts of the Boundary Layer Theory It is difficult to solve the complete viscous flow fluid around a body unless the geometry is very
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationMathCAD solutions for problem of laminar boundary-layer flow
MathCAD solutions for problem of laminar boundary-layer flow DANIELA CÂRSTEA Industrial High-School Group of Railways Transport, Craiova, ROMANIA danacrst@yahoo.com Abstract: -The problem of laminar boundary-layer
More informationMarching on the BL equations
Marching on the BL equations Harvey S. H. Lam February 10, 2004 Abstract The assumption is that you are competent in Matlab or Mathematica. White s 4-7 starting on page 275 shows us that all generic boundary
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More informationChapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More information9. Boundary layers. Flow around an arbitrarily-shaped bluff body. Inner flow (strong viscous effects produce vorticity) BL separates
9. Boundary layers Flow around an arbitrarily-shaped bluff body Inner flow (strong viscous effects produce vorticity) BL separates Wake region (vorticity, small viscosity) Boundary layer (BL) Outer flow
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationDouble reduction of partial differential equations with applications to laminar jets and wakes
Double reduction of partial differential equations with applications to laminar jets and wakes LADY NOMVULA KOKELA A Dissertation submitted to the Faculty of Science, University of the Witwatersrand, Johannesburg,
More informationCHAPTER 4 ANALYTICAL SOLUTIONS OF COUPLE STRESS FLUID FLOWS THROUGH POROUS MEDIUM BETWEEN PARALLEL PLATES WITH SLIP BOUNDARY CONDITIONS
CHAPTER 4 ANALYTICAL SOLUTIONS OF COUPLE STRESS FLUID FLOWS THROUGH POROUS MEDIUM BETWEEN PARALLEL PLATES WITH SLIP BOUNDARY CONDITIONS Introduction: The objective of this chapter is to establish analytical
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationTHE LIE-GROUP SHOOTING METHOD FOR BOUNDARY LAYER EQUATIONS IN FLUID MECHANICS *
Conference of Global Chinese Scholars on Hydrodynamics HE LIE-GROUP SHOOING MEHOD FOR BOUNDARY LAYER EQUAIONS IN FLUID MECHANICS * CHANG Chih-Wen CHANG Jiang-Ren LIU Chein-Shan Department of Systems Engineering
More informationThe Gauss-Jordan Elimination Algorithm
The Gauss-Jordan Elimination Algorithm Solving Systems of Real Linear Equations A. Havens Department of Mathematics University of Massachusetts, Amherst January 24, 2018 Outline 1 Definitions Echelon Forms
More informationUNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence
UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: Day of examination: Friday 15. June 212 Examination hours: 9. 13. This problem set consists of 5 pages. Appendices: Permitted
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationLecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)
Lecture 0: A (Brief) Introduction to Group heory (See Chapter 3.3 in Boas, 3rd Edition) Having gained some new experience with matrices, which provide us with representations of groups, and because symmetries
More informationUNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of boundary layer Thickness and classification Displacement and momentum thickness Development of laminar and turbulent flows in circular pipes
More informationComplex Analysis MATH 6300 Fall 2013 Homework 4
Complex Analysis MATH 6300 Fall 2013 Homework 4 Due Wednesday, December 11 at 5 PM Note that to get full credit on any problem in this class, you must solve the problems in an efficient and elegant manner,
More informationNotes for Expansions/Series and Differential Equations
Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationA SEMI-ANALYTICAL ANALYSIS OF A FREE CONVECTION BOUNDARY-LAYER FLOW OVER A VERTICAL PLATE
A SEMI-ANALYTICAL ANALYSIS OF A FREE CONVECTION BOUNDARY-LAYER FLOW OVER A VERTICAL PLATE Haldun Alpaslan PEKER and Galip OTURANÇ Department of Mathematics, Faculty of Science, Selcu University, 475, Konya,
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
More informationENGI 3424 First Order ODEs Page 1-01
ENGI 344 First Order ODEs Page 1-01 1. Ordinary Differential Equations Equations involving only one independent variable and one or more dependent variables, together with their derivatives with respect
More informationTHE DIFFERENTIAL TRANSFORMATION METHOD AND PADE APPROXIMANT FOR A FORM OF BLASIUS EQUATION. Haldun Alpaslan Peker, Onur Karaoğlu and Galip Oturanç
Mathematical and Computational Applications, Vol. 16, No., pp. 507-513, 011. Association for Scientific Research THE DIFFERENTIAL TRANSFORMATION METHOD AND PADE APPROXIMANT FOR A FORM OF BLASIUS EQUATION
More informationTransport by convection. Coupling convection-diffusion
Transport by convection. Coupling convection-diffusion 24 mars 2017 1 When can we neglect diffusion? When the Peclet number is not very small we cannot ignore the convection term in the transport equation.
More informationFALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
Proceedings of the International Conference on Mechanical Engineering 2 (ICME2) 8-2 December 2, Dhaka, Bangladesh ICME-TH-6 FALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
More informationThe linear equations can be classificed into the following cases, from easier to more difficult: 1. Linear: u y. u x
Week 02 : Method of C haracteristics From now on we will stu one by one classical techniques of obtaining solution formulas for PDEs. The first one is the method of characteristics, which is particularly
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationSome Notes on Linear Algebra
Some Notes on Linear Algebra prepared for a first course in differential equations Thomas L Scofield Department of Mathematics and Statistics Calvin College 1998 1 The purpose of these notes is to present
More informationScaling and Spiral Equivalence from a Numerical Viewpoint
Scaling and Spiral Equivalence from a Numerical Viewpoint Riccardo Fazio and Salvatore Iacono Abstract The non-iterative numerical solution of nonlinear boundary value problems is a subject of great interest.
More informationFACTORIZATION AND THE PRIMES
I FACTORIZATION AND THE PRIMES 1. The laws of arithmetic The object of the higher arithmetic is to discover and to establish general propositions concerning the natural numbers 1, 2, 3,... of ordinary
More informationρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More informationChapter 9 Flow over Immersed Bodies
57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 009 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,
More informationBoiling Heat Transfer and Two-Phase Flow Fall 2012 Rayleigh Bubble Dynamics. Derivation of Rayleigh and Rayleigh-Plesset Equations:
Boiling Heat Transfer and Two-Phase Flow Fall 2012 Rayleigh Bubble Dynamics Derivation of Rayleigh and Rayleigh-Plesset Equations: Let us derive the Rayleigh-Plesset Equation as the Rayleigh equation can
More informationA SYSTEM OF AXIOMATIC SET THEORY PART VI 62
THE JOOBNAL OF SYMBOLIC LOGIC Volume 13, Number 2, June 1948 A SYSTEM OF AXIOMATIC SET THEORY PART VI 62 PAUL BEKNAYS 16. The r61e of the restrictive axiom. Comparability of classes. Till now we tried
More informationA SIMILARITY SOLUTION OF FIN EQUATION WITH VARIABLE THERMAL CONDUCTIVITY AND HEAT TRANSFER COEFFICIENT
Mathematical and Computational Applications, Vol. 11, No. 1, pp. 25-30, 2006. Association for Scientific Research A SIMILARITY SOLUTION OF FIN EQUATION WITH VARIABLE THERMAL CONDUCTIVITY AND HEAT TRANSFER
More informationB. Differential Equations A differential equation is an equation of the form
B Differential Equations A differential equation is an equation of the form ( n) F[ t; x(, xʹ (, x ʹ ʹ (, x ( ; α] = 0 dx d x ( n) d x where x ʹ ( =, x ʹ ʹ ( =,, x ( = n A differential equation describes
More informationON THE EFFECTIVENESS OF HEAT GENERATION/ABSORPTION ON HEAT TRANSFER IN A STAGNATION POINT FLOW OF A MICROPOLAR FLUID OVER A STRETCHING SURFACE
5 Kragujevac J. Sci. 3 (29) 5-9. UDC 532.5:536.24 ON THE EFFECTIVENESS OF HEAT GENERATION/ABSORPTION ON HEAT TRANSFER IN A STAGNATION POINT FLOW OF A MICROPOLAR FLUID OVER A STRETCHING SURFACE Hazem A.
More informationReview of Vector Analysis in Cartesian Coordinates
Review of Vector Analysis in Cartesian Coordinates 1 Scalar: A quantity that has magnitude, but no direction. Examples are mass, temperature, pressure, time, distance, and real numbers. Scalars are usually
More informationSOLUTIONS TO ADDITIONAL EXERCISES FOR II.1 AND II.2
SOLUTIONS TO ADDITIONAL EXERCISES FOR II.1 AND II.2 Here are the solutions to the additional exercises in betsepexercises.pdf. B1. Let y and z be distinct points of L; we claim that x, y and z are not
More informationWe are going to discuss what it means for a sequence to converge in three stages: First, we define what it means for a sequence to converge to zero
Chapter Limits of Sequences Calculus Student: lim s n = 0 means the s n are getting closer and closer to zero but never gets there. Instructor: ARGHHHHH! Exercise. Think of a better response for the instructor.
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 7
Numerical Fluid Mechanics Fall 2011 Lecture 7 REVIEW of Lecture 6 Material covered in class: Differential forms of conservation laws Material Derivative (substantial/total derivative) Conservation of Mass
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationApplications of parabolized stability equation for predicting transition position in boundary layers
Appl. Math. Mech. -Engl. Ed., 33(6), 679 686 (2012) DOI 10.1007/s10483-012-1579-7 c Shanghai University and Springer-Verlag Berlin Heidelberg 2012 Applied Mathematics and Mechanics (English Edition) Applications
More informationReview for Chapter 2 Test
Review for Chapter 2 Test This test will cover Chapter (sections 2.1-2.7) Know how to do the following: Use a graph of a function to find the limit (as well as left and right hand limits) Use a calculator
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationLECTURE 5, FRIDAY
LECTURE 5, FRIDAY 20.02.04 FRANZ LEMMERMEYER Before we start with the arithmetic of elliptic curves, let us talk a little bit about multiplicities, tangents, and singular points. 1. Tangents How do we
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationSYMMETRY REDUCTION AND NUMERICAL SOLUTION OF A THIRD-ORDER ODE FROM THIN FILM FLOW
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 709-719, 2010. c Association for Scientific Research SYMMETRY REDUCTION AND NUMERICAL SOLUTION OF A THIRD-ORDER ODE FROM THIN FILM FLOW E.
More informationUNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH *
4.4 UNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH 19 Discussion Problems 59. Two roots of a cubic auxiliary equation with real coeffi cients are m 1 1 and m i. What is the corresponding homogeneous
More informationA Solution Method for the Reynolds-Averaged Navier-Stokes Equation
A Solution Method for the Reynolds-Averaged Navier-Stokes Equation T.-W. Lee Mechanical and Aerospace Engineering, Arizona State University, Tempe, AZ, 8587 Abstract- Navier-Stokes equations are difficult
More information18.303: Introduction to Green s functions and operator inverses
8.33: Introduction to Green s functions and operator inverses S. G. Johnson October 9, 2 Abstract In analogy with the inverse A of a matri A, we try to construct an analogous inverse  of differential
More informationPhysics 221B Spring 2018 Notes 30 The Thomas-Fermi Model
Copyright c 217 by Robert G. Littlejohn Physics 221B Spring 218 Notes 3 The Thomas-Fermi Model 1. Introduction The Thomas-Fermi model is a relatively crude model of multi-electron atoms that is useful
More informationChapter 2. Mathematical Reasoning. 2.1 Mathematical Models
Contents Mathematical Reasoning 3.1 Mathematical Models........................... 3. Mathematical Proof............................ 4..1 Structure of Proofs........................ 4.. Direct Method..........................
More informationLECTURE 3. RATIONAL NUMBERS: AN EXAMPLE OF MATHEMATICAL CONSTRUCT
ANALYSIS FOR HIGH SCHOOL TEACHERS LECTURE 3. RATIONAL NUMBERS: AN EXAMPLE OF MATHEMATICAL CONSTRUCT ROTHSCHILD CAESARIA COURSE, 2011/2 1. Rational numbers: how to define them? Rational numbers were discovered
More informationResearch Article Lie Group Analysis of Unsteady Flow and Heat Transfer over a Porous Surface for a Viscous Fluid
Journal of Applied Mathematics Volume 202, Article ID 675287, 7 pages doi:0.55/202/675287 Research Article Lie Group Analysis of Unsteady Flow and Heat Transfer over a Porous Surface for a Viscous Fluid
More informationA Note on the Differential Equations with Distributional Coefficients
MATEMATIKA, 24, Jilid 2, Bil. 2, hlm. 115 124 c Jabatan Matematik, UTM. A Note on the Differential Equations with Distributional Coefficients Adem Kilicman Department of Mathematics, Institute for Mathematical
More informationA NOTE ON BLASIUS TYPE BOUNDARY VALUE PROBLEMS. Grzegorz Andrzejczak, Magdalena Nockowska-Rosiak, and Bogdan Przeradzki
Opuscula Math. 33, no. 1 213, 5 17 http://dx.doi.org/1.7494/opmath.213.33.1.5 Opuscula Mathematica A NOTE ON BLASIUS TYPE BOUNDARY VALUE PROBLEMS Grzegorz Andrzejczak, Magdalena Nockowska-Rosiak, and Bogdan
More informationLocal properties of plane algebraic curves
Chapter 7 Local properties of plane algebraic curves Throughout this chapter let K be an algebraically closed field of characteristic zero, and as usual let A (K) be embedded into P (K) by identifying
More informationLagrangian analysis of the laminar flat plate boundary layer
Lagrangian analysis of the laminar flat plate boundary layer Mohammad Gabr* EgyptAir Maintenance & Engineering, Cairo, Egypt 1 The flow properties at the leading edge of a flat plate represent a singularity
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationContents of today s lecture
Contents of today s lecture Blasius solution for laminar flat-plate boundary layer where external velocity is constant ( ) About solution methods for laminar boundary layers Thwaites method as an example
More informationPROBLEMS In each of Problems 1 through 12:
6.5 Impulse Functions 33 which is the formal solution of the given problem. It is also possible to write y in the form 0, t < 5, y = 5 e (t 5/ sin 5 (t 5, t 5. ( The graph of Eq. ( is shown in Figure 6.5.3.
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationSIMILARITY SOLUTION FOR MHD FLOW THROUGH VERTICAL POROUS PLATE WITH SUCTION
Journal of Computational and Applied Mechanics, Vol. 6., No. 1., (2005), pp. 15 25 SIMILARITY SOLUTION FOR MHD FLOW THROUGH VERTICAL POROUS PLATE WITH SUCTION Mohammad Ferdows, Masahiro Ota Department
More informationTopic 5: The Difference Equation
Topic 5: The Difference Equation Yulei Luo Economics, HKU October 30, 2017 Luo, Y. (Economics, HKU) ME October 30, 2017 1 / 42 Discrete-time, Differences, and Difference Equations When time is taken to
More informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 1982 NOTES ON MATRIX METHODS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 198 NOTES ON MATRIX METHODS 1. Matrix Algebra Margenau and Murphy, The Mathematics of Physics and Chemistry, Chapter 10, give almost
More informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
More informationGetting Started with Communications Engineering
1 Linear algebra is the algebra of linear equations: the term linear being used in the same sense as in linear functions, such as: which is the equation of a straight line. y ax c (0.1) Of course, if we
More informationBoundary Layer Theory. v = 0, ( v)v = p + 1 Re 2 v. Consider a cylinder of length L kept in a inviscid irrotational flow.
1 Boundary Layer Theory 1 Introduction Consider the steady flow of a viscous fluid. The governing equations based on length scale L and velocity scale U is given by v =, ( v)v = p + 1 Re 2 v For small
More informationA curvature-unified equation for a non-newtonian power-law fluid flow
Int. J. Adv. Appl. Math. and Mech. 2(3) (215) 72-77 (ISSN: 2347-2529) Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics A curvature-unified equation
More informationApplication of the perturbation iteration method to boundary layer type problems
DOI 10.1186/s40064-016-1859-4 RESEARCH Open Access Application of the perturbation iteration method to boundary layer type problems Mehmet Pakdemirli * *Correspondence: mpak@cbu.edu.tr Applied Mathematics
More informationSymmetry Reduction of Chazy Equation
Applied Mathematical Sciences, Vol 8, 2014, no 70, 3449-3459 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ams201443208 Symmetry Reduction of Chazy Equation Figen AÇIL KİRAZ Department of Mathematics,
More informationu = λ u 2 in Ω (1) u = 1 on Ω (2)
METHODS AND APPLICATIONS OF ANALYSIS. c 8 International Press Vol. 15, No. 3, pp. 37 34, September 8 4 SYMMETRY ANALYSIS OF A CANONICAL MEMS MODEL J. REGAN BECKHAM AND JOHN A. PELESKO Abstract. A canonical
More informationHeat and Mass Transfer over an unsteady Stretching Surface embedded in a porous medium in the presence of variable chemical reaction
Vol:5, No:2, 20 Heat and Mass Transfer over an unsteady Stretching Surface embedded in a porous medium in the presence of variable chemical reaction TGEmam International Science Index, Mathematical and
More informationMath 4381 / 6378 Symmetry Analysis
Math 438 / 6378 Smmetr Analsis Elementar ODE Review First Order Equations Ordinar differential equations of the form = F(x, ( are called first order ordinar differential equations. There are a variet of
More information