Numerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
|
|
- Allen Hall
- 5 years ago
- Views:
Transcription
1 Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., nd Ave. N. E., Seattle, WA ; 6.5. Orthogonal Collocation on Finite Elements The method of orthogonal collocation on finite elements provides a collocation method that is high order. It is applied here to the convective diffusion equation to illustrate its success. When the Peclet number is large, the orthognal collocation method will not provide good results since the number N must be so large to approximate the solution. In that case we are better off using the method of orthogonal collocation on finite elements, as developed by Carey and Finlayson [975]. The domain is divided into elements and the method of orthogonal collocation is applied within each element. The numbering scheme is shown in Figure 6E.. There are three types of nodes: those at which a boundary condition is applied (2 of these), those at which the residual is set to zero (NCOL x NE of these), and those at which the first derivative is made continuous across element boundaries ( NE of these). The total number of points is then NT = (NCOL+) NE +. Within an element we transform the domain using ξ = x x k x k, x k = x k+ x k m2 so that the variable ξ goes from to within an element (rather than - to as we have used for the Galerkin methods). The differential equation can be transformed to give c t + Pe c x k ξ = x 2 k 2 c ξ 2 m2 and the boundary conditions are c(,t) = c x NE ξ (,t) = c(x,) = f(x) m22 Now we are going to use four different numbering systems depending on what we are solving for. The local numbering systems within an element is shown in Figure 6E. and the global number system is shown in Figure 6E.2. The index i represents a global number and the index I represents a local number; in the later case the element must also be identified. Thus we have the correspondence i = (k ) (NCOL + ) + I m23 x i = x (k) + ξ I x k
2 I = = NP = NCOL+2 ξ = Figure 6E.. Local numbering On an element we set the residuals to zero. c t (ξ I, t) + Pe c x k ξ (ξ I, t) = x 2 k 2 c ξ 2(ξ I, t), I =,..., NCOL m24 We also apply the boundary conditions x NE c i= = or c I= () = NCOL+2 A NCOL+2,I c I (NE) = I= m25 c (k) I () = f(x i ) Next we make the first derivatives continuous between the elements by taking Evaluated in the collocation method this is dc (x(k ) NCOL+2 ) = dc (x(k) ) m26 x k NCOL+2 A NCOL+2,I c I (k ) = I= x k NCOL+2 A,I c I (k) I= m27 The method of solution is to use an explicit method to step forward in time. c J (k),n+ c J (k),n t + Pe x k N+2 I= A JI c (k),n I = = N+2 xk 2 I= B JI c (k),n I, J = 2,...,NCOL+,k=,NE m28 If we know the solution at time level n, c i n we can find c I (k),n. Thus we can obtain the value of c i n+ for
3 x k x (k) x (k+) Residual Condition Boundary Conditions Element boundary, continuity x i Global numbering = NT = 2 NE Numbering for solving the ordinary differential equations = NE+ Figure 6E.2. Nodal numbering system each of the points interior to the element, marked by in Figure 6E.2. Implicit methods are also possible, and the linear algebra is discussed by Finlayson [98, p. 6]. Suppose we have made the calculation Eq.(m28) for all of the collocation points interior to each element. We still need the values of c i n+ for each of the points between elements and the two end points. Rearrange the equations Eq.(m25) and Eq.(m27) to have the unknowns on the left-hand side.
4 x k A NP, c (k ),n+ I= + = NP x k I=2 x k A NP,NP c (k ),n+ NP A Np,I c (k ),n+ I + NP x k I=2 x k A, c (k),n+ I= A J c (k),n+ I x k A,NP c (k),n+ NP = m29 c (),n+ I= = x NE A NP, c (NE),n+ + x NE A NP,NP c (NE),n+ NP = NP x NE I=2 A NP,I c (NE),n+ m3 We make the function continuous by taking c (k ),n+ NP = c (k),n+ m3 This system of equations can be represented as a tri-diagonal form (see Eq. 2.5). The values of x i in Eq. 2.5 are the values of c i n+ at the two end points and the points between the elements, using the third numbering system shown in Figure 6E.2. This tri-diagonal system can be solved and the results rearranged to provide the solution for each and every collocation point shown in Figure 6E.2, first numbering system, at the n+-th time level. Results are given for the method, but using the implicit method (which was readily available). For the convective diffusion equation the collocation results for 25 elements and 2 interior points are shown in Figure 6.33a. There are small oscillations. Jensen and Finlayson [98] have shown that oscillations do not occur in the steady state problem provided Pe x m35 This is the number that is comparable with Pe x/2 in the Galerkin method. For this calcualtion the quantity in Eq. (m35) takes the value.5. The Galerkin method with Pe x/2 = gives the results shown in Figure 5.8a; the Taylor-Galerkin method with Pe x/2 = gives the results shown in Figure 5.8e. There are small oscilllations in those solutions, too. When more elements are used, 5 elements and 2 interior collocation points, the method of orthogonal collocation on finite elements gives the results shown in Figure 6.33b. In this case there are no oscillations. All of these methods provide good solutions. The collocation method can also be applied with higher degree polynomials. A solution with 25 elements and 5 interior collocation points (a sixth degree polynomial within each element) gives a solution with oscillations that are very small. Here the key value is Pe x 2 5 m36 and it takes the value 4. The method of orthogonal collocation on finite elements can also be applied to non-linear prob-
5 lems as well. If Burger s equation is solved in the form u t + u u x = ν 2 u x 2 m37 we just replace Eq.(m24) with u t (ξ I, t) + u(ξ I, t) x k u ξ (ξ I, t) = ν xk 2 2 u ξ 2(ξ I, t), I =,..., NCOL m38 An explicit method can be used with this set, as well. Solutions for such a method are not shown here.
6 7.4. Special Petrov-Galerkin Weighting Functions Because S is an even function and satisfies Eq. (7.44) we have ds ds ds = h Since S() = this also says S() =. The integral s ds ds ds h2 can be evaluated by regarding s as the derivative of s 2 /2 and then integrating by parts. s ds ds ds = d ds (s S ) ds S ds = S ds h3 Since S is an even function we have and hence S ds = 2 s ds ds ds = 2 S ds = 2 h4 h5 The following relation can be obtained by integrating by parts the left-hand side. s 2 ds ds ds = 2 s S ds h6 With similar manipulations, (and using the fact that T is an odd function), we obtain dt ds ds =, s dt ds ds = T ds, s 2 dt ds ds = 2 s T ds h7 In the Petrov-Galerkin method the weighted residual is set to zero, Eq. (7.4). Integrate this by parts to obtain ν dw j du dw j V(u) = h8
7 The first term is ν dw j du = ν du + ν du h9 Consider each part in turn. du = NT i= u i dn i h2 But the derivative of the trial function N j is constant within an element so that we get dn j = h ds ds dn j = h ds ds + dn j+ = h ds ds h ds ds h2 Combining all results gives du = u j h ds ds + u j h ds ds ds ds + u j + h ds ds h22 The terms are evaluated using Eq. (7.42) and (7.44). du = h ( u j + 2 u j u j+ ) = u j u j u j+ u j = m h h j m j+ h23 Similarily for the second part of Eq. (h9). du = u j h ds ds + u j h ds ds ds ds + u j + h ds ds h24 By Eq. (7.42) we have
8 ds ds = ds ds h25 and by Eq. (7.44) both terms are zero. Thus there is not contribution due to the T j function. The first term of the weighted residual is then ν dw j du = ν (m j m j+ ) h26 The second term of the weighted residual is dw j V(u) = V(u) + V(u) = ds V(u) ds + ds V(u) ds + ds V(u) ds + ds V(u) ds h27 In the left integral (s = - to ) we use u = u j + ( u j u j ) ( s + ) = u j + m j h s h28 In the right integral (s = to ) we use u = u j + ( u j+ u j ) s = u j + m j+ h s h29 Examine each of the terms in Eq. (h27). V(u) = ds V(u j + m j h s) ds + ds V(u j + m j+ h s) ds h3 Since /ds is an odd function ds V(u j + m j h s) ds = ξ= ξ= dξ V( u j m j h ξ ) dξ = ds V(u j m j h s) ds h3 and we get
9 V(u) = ds V(u j m j h s) ds + ds V(u j + m j+ h s) ds h32 The same steps applied to the term involving /ds, which is an even function, gives V(u) = ds V(u j m j h s) ds + ds V(u j + m j+ h s) ds h33 Thus the weighted residual becomes ν (m j+ m j ) ds V(u j m j h s) ds + ds V(u j + m j+ h s) ds + + ds V(u j m j h s) ds + dt h34 j ds V(u j + m j+ h s) ds = Now for this application we have V(u j m j h s) = 2 u2 j m j u j h s + 2 m2 j h 2 s 2 V(u j + m j+ h s) = 2 u2 j + m j+ u j h s + 2 m2 j+ h2 s 2 h35 Substitution into Eq. (h34) gives ν + u j h α 2ν ( m j+ m j ) u j 2 ( u j+ u j ) = = h 2 ( mj+ 2 m2 j ) s S(s) ds + h 2 ( mj+ 2 + m2 j ) s T (s) ds h36 In this equation we have let α = 2 T (s) ds h37 The parameter α thus determines the amount of upstream weighting; if both α and u j have the same sign then the inclusion of the function T j increases the coefficient of the diffusion term. The presence of the function T j also decreases the truncation error. This development shows how to get an upstream method using the weighting function (7.39).
10 .3. Combustion in Flames The equations governing the propagation of flames in one dimension are given by Ramos [987]: ρ t + x ( ρ u ) = ρ Y i t + u Y i = R x i x ( ρ Y i V i ) b ρc T p t + u T x = p t N i = h i R i q x N i= ρ Y i V i C pi T x b2 Y N = N i= p = ρ R T N i= h i = h i + T Y i Y i Mˆ i T C pi dt,, C p = N i= q = λ T x Y i C pi b3 These are the continuity equation, conservation of mass, the energy equation, an equation saying the mass fractions add to one, the ideal gas law, the equation for heat capacity of a mixture, the equation for enthalpy of a species in an ideal gas, and q is the heat flux. The variables are ρ the density, u the velocity, t the time, x the distance, Y i the mole fraction of the i-th species, w i dot the reaction rate for the i-th species, V i the diffusion velocity of the i-th species, p the pressure, T the absolute temperature, R squiggle the gas constant, C p the heat capacity of the mixture, h i the enthalpy of the i-th species, and λ the thermal conductivity. In these equations we have assumed that the flame speed is slow enough that the Mach number is low; in that case the pressure is constant in space. We have also neglected the Soret and Dufour effects (or thermal diffusion). One form of the equations is for an infinite premixed combustible fluid mixture. The x-position goes from - < x < and the flame propagates from one boundary (x = - ) to the other. If the mixture is stationary, the fluid velocity ahead of the flame is zero. We transform the equations into a Lagrangian coordinate system in order to eliminate the convection term. This is possible because of the infinite extent of the physical domain. We go from where ( t, x ) ( τ, ψ ) b4
11 τ = t x = ρ t = x ρ t = x (ρ u ) x = ρ u b5 The derivatives are then calculated using ρ t = ρ τ (ρu) x τ t + ρ = (ρu) τ t τ x + (ρu) x b6 Since τ t =, t = ρ u τ x =, x = ρ b7 we obtain ρ t = ρ τ (ρu) x ρu ρ = ρ (ρu) b8 The continuity equation is then ρ τ ρu ρ + ρ (ρu) = b9 which becomes ρ τ + u ρ2 = b Likewise the derivatives of temperature and mass fraction are given by T x = T τ τ x + T x = ρ T b Thus the convective derivatives are
12 T t + u T x = T T ρu τ + u ρ T = T τ b2 With these manipulations the equations become Y i τ = ρ R i ( ρ Y i V i ) T τ = ρc p N i = h i R i ρ q N i= ρ 2 Y i V i C T pi b3 The boundary conditions are T T (,τ) = (,τ) = Y i (,τ) = Y i (,τ) = T(,τ) = T un, Y i (,τ) = Yi un b4 and the initial conditions are T(ψ,) = T un (ψ) Y i (ψ,) = Yi un (ψ) b5 where the superscript un represents the value in the unburned region. The Fickian model is obtained by making further assumptions. Assume that the binary diffusion coefficients of all pairs of chemicals are the same, the specific heats of all species are the same and that ρ 2 D = β = constant The diffusion velocity is given by Fick s law b6 V i = D x ( ln Y i ) = D ρ ( ln Y i ) = Dρ Y i ρ Y i V i = D ρ 2 Y i = β Y i Y i b8 The Lewis number is taken as one.
13 λ Le = = ρdc p b45 The equations are then Y i τ = β 2 Y i 2 + R i ρ T τ = β 2 T 2 N i= h i ρ R i ρc p b7 Another simplification of Eq.(b3) applies when a burner-stabilized flame is achieved. The premixed gas enters the burner at the flame holder at x = and flows to the flame zone and on to the burned zone (see Figure E.3). In this case we make the same transformation in Eq. (b.5) except that we replace one equation with: t = x ρ t = x (ρ u ) x = ρ u + ṁ(t) b9 where ṁ (t) = ρ(,t) u(,t) b2 is the mass flux at the flame holder. Eq. (b3) are then T b Y b flame T un Y un flow Figure E.3. Burner-stabilized Flame Geometry
14 Y i τ + ṁ Y i = ρ R i T τ + ṁ T = ρc p N ( ρ Y i V i ) i = h i R i + ρ q N i= ρ 2 Y i V i C T pi b2 We consider only steady-state solutions of these equations ṁ dy i dψ = R i ṁ dt dψ = C p N d dψ ( ρ Y i V i ) i = h i R i dq dψ N i= ρ Y i V i C dt pi dψ b22 together with boundary conditions T() = T un, Y i () = Y un i dt () = dt ( ) = dy i () = dy i ( ) = b23 We must solve these equations for the distribution of mole fractions, Y i, and temperature, T; we must also solve for the mass flux, m dot.
15 .4. Transient Catalytic Muffler The unsteady problem uses the full equations Eq. (.5)-(.5). A finite difference formulation of this equation would give α 5 α dts i 3 dt r hs L 2 s T i + Ti s + Ti s z Nu ( Ts i T i ) + β 5 < r i > = c which is to be solved subject to the steady equations Eq. (.35)-(.36) and Eq. (.38). Young and Finlayson [976] did this using orthogonal collocation on finite elements for the solid temperature and solving the integro-differential equations, Eq. (.5)-(.5). Here we consider several other options. If we evaluate the integral in Eq. Eq. (.48)-(.49) from one grid point to the next we get T i+ = e 2 α Nu z T i + 2 α Nu z i+ If the integral is evaluated using the trapezoid rule, the result is z i e 2 α Nu ( z i+ z') T s (z' ) dz', z i+ = z i + z c2 T i+ = e 2 α Nu z T i + 2 α Nu z T s 2 α Nu z 2 i + T s 2 i+ c3 This equation is of the form T i+ = γ T i + γ 2 T s i+ + γ 3 T s i c4 If the Euler method is applied to Eq. (.35)-(.36) we get T i+ T i z + 2 α Nu T i = 2 α Nu T s i c5 or if a second-order Runge-Kutta method is applied we get instead T i+ T i z + 2 α Nu 2 ( T i + T i+ ) = 2 α Nu 2 ( T s i + T s i+ ) c6 Both Eq. (c5-6) can be rearranged to give or T i+ = T i ( 2 α Nu z ) + 2 α Nu z T s i α Nu z T i+ = T + α Nu z i + z 2 2 α Nu + α Nu z ( Ts i + T s i+ ) c7
16 Notice that these are appoximations in various degrees to Eq. (c3). They can also be put into the form of Eq. (c4). Thus the temperature equations are Eq. (c) and Eq. (c4). The equations are similar for the mole fraction Y and Y s. The equation (.47) or Eq. (.49) can be written in finite difference form as where the Euler method uses Y i+ = δ Y i + δ 2 Y s i+ + δ 3 Y s i c8 δ = 2 α 2 Sh z, δ 2 =, δ 3 = 2 α 2 Sh z c9 and the Runge-Kutta method uses δ = α 2 Sh z + α 2 Sh z, δ 2 = α 2 Sh z + α 2 Sh z, δ 3 = α 2 Sh z + α 2 Sh z c2 The Eq. (.38) is then α 6 β 4 < r > i + 2 Sh (Ys i Y i ) = c2 Consider the case when the variables T i, Y i, T i s, Y i s are known at some time t. Then Eq. (c) can be used to step forward to determine T i s at time t + t for all i. Eq. (c4) can then be used to step forward to determine T i at time t + t for all i. Next look at the concentration equation for the first node. α 6 β 4 < r > + 2 Sh (Ys Y ) = c22 The value of Y would be known from the boundary condtions at every time t, including time t + t. Thus we must evaluate the reaction rate. Use y = Y in Eq. (.39), then y = Y in Eq. (.42), and then calculate < r > from Eq. (.43). Then Eq. (c22) gives Y s. At this point in the calculations we have the value of T i and T i s at time t + t for all i, and we know Y, Y s at time t + t. We then solve Eq. (c8) and (c2) for Y 2, Y 2 s at time t + t. These are written in terms of the unknowns Y i+ n+, Y i+ s,n+ and the known quantities Y i n+, Y i s,n+. Yi+ n+ δ 2 Yi+ s,n+ = δ Yi n+ + + δ 3 Yi s,n+ c23 α 6 β 4 < r > n+ i+ + Sh (Ys,n+ 2 i+ Yi+ n+ ) = c24 This set of two equations must be solved simultaneously for i=. Then we can increase i and repeat until we reach an i representing the end of the reactor. Then we know T i, Y i, T i s, Y i s at time t + t for all i. We then repeat the process. Of course, all the techniques discussed in the book can be applied. In particular the nodes can be moved as the front moves and concentrated near the reaction zone. Travelling waves could also be evaluated.
17 2.3. Compressibility Method The first method originated with Chorin [967] and is called the artificial compressibility method because the incompressible fluid is allowed to have a small compressibility. The chief problem with the Navier-Stokes equations and the continuity equation is that the Navier-Stokes equations are equations for velocities, and they include the pressure gradient. The continuity equation is a constraint on the velocity but is in essence the equation for pressure. The continuity equation for an incompressible fluid does not have a time derivative in it, so in the artificial compressibility method a time derivative is added. The continuity equation is expanded to be p t + c2 u = u The constant parameter c 2 is chosen to aid in convergence and the transient itself does not have any physical significance. The method is termed artificial compressibility because Eq. (u) is the continuity equation for a compressible fluid with an equation of state p = c 2 ρ u2 The method itself is listed in Table 2.2. Now we have three (or four) parabolic equations in time for the pressure and the two (or three) velocity components. Standard algorithms for parabolic equations can be applied. When the algorithm is applied each variable can be identified with the same node, as shown in Figure 2.E(a), or different grids can be used for different variables, as in Figure 2.E(b). u v p u u p u v p u p v p (a) regular mesh v v (b) staggered mesh Figure 2.E. Finite difference meshes The first grid, Figure 2.E(a) is not a preferred one because the boundary conditions for pressure are inconvenient to derive. We discuss here the equations involving boundary terms for the staggered grid, Figure 2.E(b), which is known as the MAC mesh. Consider the problem shown in Figure 2.E2, an entry flow problem with a solid boundary, a line of symmetry, an entry condition and an exit condition. For side A the u nodes occur on the boundaries, so that we set the u velocity there. The same thing is true for the v velocity on B and C. However, the v velocity nodes are not on the boundary A, nor are the u velocity nodes on the boundary B. On A we force the boundary condition by taking the v-
18 value x/2 outside the mesh to equal minus the v-value x/2 inside the mesh. Then the average value is zero, and the average value represents the v-value on A. u = v = u = 2(-y ) v = 2 A B C D u/ x = v/ x = v =, u/ y = Figure 2.E2. Entry-length problem v,j+/2 = v,j+/2 In an analogous fashion we use u i+/2,m+ = u i+/2,m for the u-boundary conditions on surface B. To force the first derivative of u to zero on boundary C we use u i+/2,m+ = u i+/2,m This represents the value of the first derivative centered on y=. On the exit boundary D the v derivative is centered on the exit. v N+,j+/2 = v N,j+/2 The u-velocity is expressed at nodes that are on the boundary. To set the first derivative to zero we take the value x from the exit to equal to the value + x from the exit. u N+3/2,j = u N /2,j Following these rules, and using an explicit method for integrating in time, gives the following algorithm for this case. u n+ i+/2,j un i+/2,j = pn i+,j pn i,j u n u n i+3/2,j un i /2,j t x i+/2,j + 2 x 4 (vn i+,j+/2 + vn i,j+/2 + vn i,j /2 + vn i+,j /2 )un i+/2,j+ un i+/2,j 2 y + + Re u n i + 3/2,j 2un i+/2,j + un i /2,j x 2 i =,..., N; j=,..., M u n i+/2,j+ + i+/2,j + un i+/2,j y 2 u5
19 v n+ i,j+/2 vn i,j+/2 = pn i,j+ pn i,j t y v n i,j+/2 v n i,j+3/2 vn i,j /2 2 y 4 (un i+/2,j + un i+/2,j+ + un i /2,j + un i /2,j+ ) vn i+,j+/2 vn i,j+/2 2 x + + v n i +,j+/2 2vn i,j+/2 + vn i,j+/2 Re x 2 + v n i,j+3/2 vn i,j+/2 + vn i,j /2 y 2 i =,..., N; j=,..., M- p n+ i,j p n i,j t = c 2 u n i + + /2,j u n+ x i /,j + vn+ i,j+/2 v n+ i,j /2 y i =,..., N; j=,..., M surface A: u / 2, j = 2 ( y 2 j ), v,j+/2 = v,j+/2 surface B: u i+/2,m+ = u i+/2,m, v i,m+/2 = surface C: u i+/2, = u i+/2,, v i,/2 = surface D: u N+3/2,j = u N /2,j, v N+,j+/2 = v N,j+/2
20 3.4. u-v Formulation Another way to solve this problem is to tackle the second-order problem, Eq. (3.9), but with the new variable v = u t f3 Then Eq. (3.9) is v t = v E + E 2 u y 2 f4 coupled with u t = v f2 We now have a system of two equations to solve for u(y,t) and v(y,t). Note that the two equations look like a diffusion Eq.(f4) with a source term and an equation with only a source term, Eq. (f2). The convective nature of the problem is not apparent with this formulation. Only a centered finite difference method is suitable here, since upstream derivatives cannot be used on the non-existent convection terms. dv i dt = v i E + u i+ 2 u i + u i E y 2 du i dt = v i i = 2,..., N f5 This set is applied at the points i = 2,..., N-. The two end points must be handled with boundary conditions. u(,t) = or u (t) = u(,t) = or u N (t) = f6 Thus we know the value of u and u N that are needed in the first of Eq. (f5) when i = 2 or i = N-. We do not know the corresponding values of v and v N, but fortunately we do not need them.. Thus we have enough equations for the solution. The initial conditions are u(y,) = or u i () =, i =,...N u t (y,) = or v i () =, i=2,...,n f7
21 We then integrate the 2(N-2) equations, Eq.(f5), together with the boundary conditions, Eq.(f6), beginning with the initial conditions, Eq.(f7). The stability of this set of equations can be derived using the Fourier analysis, which gives t y2 4 f2 or a time step four times smaller than for the other set of equations.
Numerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
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 informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
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 informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationEULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS
EULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS ML Combrinck, LN Dala Flamengro, a div of Armscor SOC Ltd & University of Pretoria, Council of Scientific and Industrial Research & University
More informationTheoretical Developments in Group Combustion of Droplets and Sprays
Theoretical Developments in Group Combustion of Droplets and Sprays William A. Sirignano University of California, Irvine Collaborations: Guang Wu, current student; Randall Imaoka, former student, US Navy;
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationNicholas Cox, Pawel Drapala, and Bruce F. Finlayson Department of Chemical Engineering, University of Washington, Seattle, WA, USA.
Transport Limitations in Thermal Diffusion Nicholas Cox, Pawel Drapala, and Bruce F. Finlayson Department of Chemical Engineering, University of Washington, Seattle, WA, USA Abstract Numerical simulations
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More information5. Coupling of Chemical Kinetics & Thermodynamics
5. Coupling of Chemical Kinetics & Thermodynamics Objectives of this section: Thermodynamics: Initial and final states are considered: - Adiabatic flame temperature - Equilibrium composition of products
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationChapter 2 Finite-Difference Discretization of the Advection-Diffusion Equation
Chapter Finite-Difference Discretization of the Advection-Diffusion Equation. Introduction Finite-difference methods are numerical methods that find solutions to differential equations using approximate
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More information10.34 Numerical Methods Applied to Chemical Engineering. Quiz 2
10.34 Numerical Methods Applied to Chemical Engineering Quiz 2 This quiz consists of three problems worth 35, 35, and 30 points respectively. There are 4 pages in this quiz (including this cover page).
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationAME 513. " Lecture 7 Conservation equations
AME 51 Principles of Combustion Lecture 7 Conservation equations Outline Conservation equations Mass Energy Chemical species Momentum 1 Conservation of mass Cubic control volume with sides dx, dy, dz u,
More informationFundamentals of Transport Processes Prof. Kumaran Indian Institute of Science, Bangalore Chemical Engineering
Fundamentals of Transport Processes Prof. Kumaran Indian Institute of Science, Bangalore Chemical Engineering Module No # 05 Lecture No # 25 Mass and Energy Conservation Cartesian Co-ordinates Welcome
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationChapter 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 informationOn Computing Reactive Flows Using Slow Manifolds
On Computing Reactive Flows Using Slow Manifolds by S. Paolucci 1, J. M. Powers 2, and S. Singh 3 Department of Aerospace and Mechanical Engineering University of Notre Dame Notre Dame, Indiana 46556-5637,
More informationarxiv: v1 [physics.chem-ph] 6 Oct 2011
Calculation of the Minimum Ignition Energy based on the ignition delay time arxiv:1110.1163v1 [physics.chem-ph] 6 Oct 2011 Jens Tarjei Jensen a, Nils Erland L. Haugen b, Natalia Babkovskaia c a Department
More informationTAU Solver Improvement [Implicit methods]
TAU Solver Improvement [Implicit methods] Richard Dwight Megadesign 23-24 May 2007 Folie 1 > Vortrag > Autor Outline Motivation (convergence acceleration to steady state, fast unsteady) Implicit methods
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More information2. Conservation of Mass
2 Conservation of Mass The equation of mass conservation expresses a budget for the addition and removal of mass from a defined region of fluid Consider a fixed, non-deforming volume of fluid, V, called
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationIntroduction to Physical Acoustics
Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with
More information2 Law of conservation of energy
1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion
More informationSIMPLE Algorithm for Two-Dimensional Channel Flow. Fluid Flow and Heat Transfer
SIMPLE Algorithm for Two-Dimensional Channel Flow Fluid Flow and Heat Transfer by Professor Jung-Yang San Mechanical Engineering Department National Chung Hsing University Two-dimensional, transient, incompressible
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
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 informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationThe effect of Entry Region on Thermal Field
The effect of Entry Region on Thermal Field Flow Fractionation Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract Achieving
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
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 informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationDaniel J. Jacob, Models of Atmospheric Transport and Chemistry, 2007.
1 0. CHEMICAL TRACER MODELS: AN INTRODUCTION Concentrations of chemicals in the atmosphere are affected by four general types of processes: transport, chemistry, emissions, and deposition. 3-D numerical
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 informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
More informationConvective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay
Convective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay Module No. # 01 Lecture No. # 32 Stefan Flow Model We are now familiar with
More informationPartial Differential Equations
Next: Using Matlab Up: Numerical Analysis for Chemical Previous: Ordinary Differential Equations Subsections Finite Difference: Elliptic Equations The Laplace Equations Solution Techniques Boundary Conditions
More informationComparison of some approximation schemes for convective terms for solving gas flow past a square in a micorchannel
Comparison of some approximation schemes for convective terms for solving gas flow past a square in a micorchannel Kiril S. Shterev and Sofiya Ivanovska Institute of Mechanics, Bulgarian Academy of Sciences,
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM. Prof. Peter Bermel February 6, 2017
ECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM Prof. Peter Bermel February 6, 2017 Outline Static Equilibrium Dynamic Equilibrium Thermal transport mechanisms Thermal transport modeling
More informationA Numerically Stable Method for Integration of the Multicomponent Species Diffusion Equations
Journal of Computational Physics 174, 460 472 (2001) doi:101006/jcph20016930, available online at http://wwwidealibrarycom on A Numerically Stable Method for Integration of the Multicomponent Species Diffusion
More informationHeat Transfer Modeling using ANSYS FLUENT
Lecture 2 - Conduction Heat Transfer 14.5 Release Heat Transfer Modeling using ANSYS FLUENT 2013 ANSYS, Inc. March 28, 2013 1 Release 14.5 Agenda Introduction Energy equation in solids Equation solved
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationAE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 1...in the phrase computational fluid dynamics the word computational is simply an adjective to fluid dynamics.... -John D. Anderson 2 1 Equations of Fluid
More informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
More informationAME 513. " Lecture 8 Premixed flames I: Propagation rates
AME 53 Principles of Combustion " Lecture 8 Premixed flames I: Propagation rates Outline" Rankine-Hugoniot relations Hugoniot curves Rayleigh lines Families of solutions Detonations Chapman-Jouget Others
More informationChemical reactors. H has thermal contribution, pressure contribution (often negligible) and reaction contribution ( source - like)
Chemical reactors - chemical transformation of reactants into products Classification: a) according to the type of equipment o batch stirred tanks small-scale production, mostly liquids o continuous stirred
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationCFD modeling of combustion
2018-10 CFD modeling of combustion Rixin Yu rixin.yu@energy.lth.se 1 Lecture 8 CFD modeling of combustion 8.a Basic combustion concepts 8.b Governing equations for reacting flow Reference books An introduction
More informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationChapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
More informationThe Bernoulli Equation
The Bernoulli Equation The most used and the most abused equation in fluid mechanics. Newton s Second Law: F = ma In general, most real flows are 3-D, unsteady (x, y, z, t; r,θ, z, t; etc) Let consider
More informationTwo-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of
More informationMultivariable Calculus
Multivariable Calculus In thermodynamics, we will frequently deal with functions of more than one variable e.g., P PT, V, n, U UT, V, n, U UT, P, n U = energy n = # moles etensive variable: depends on
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationInvestigation of ignition dynamics in a H2/air mixing layer with an embedded vortex
Paper # 070LT-0211 The 8th US National Meeting of the Combustion Institute, Park City, UT, May 19-22, 2013 Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex S.K. Menon
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationConvective Vaporization and Burning of Fuel Droplet Arrays
Convective Vaporization and Burning of Fuel Droplet Arrays Guang Wu and William A. Sirignano Mechanical & Aerospace Engineering University of California, Irvine May 19, 2009 Background and Task In many
More information7.2 Sublimation. The following assumptions are made in order to solve the problem: Sublimation Over a Flat Plate in a Parallel Flow
7..1 Sublimation Over a Flat Plate in a Parallel Flow The following assumptions are made in order to solve the problem: 1.. 3. The flat plate is very thin and so the thermal resistance along the flat plate
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationCE 329, Fall 2015 Second Mid-Term Exam
CE 39, Fall 15 Second Mid-erm Exam You may only use pencils, pens and erasers while taking this exam. You may NO use a calculator. You may not leave the room for any reason if you do, you must first turn
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationPROCESS SYSTEMS ENGINEERING Dr.-Ing. Richard Hanke-Rauschenbach
Otto-von-Guerice University Magdeburg PROCESS SYSTEMS ENGINEERING Dr.-Ing. Richard Hane-Rauschenbach Project wor No. 1, Winter term 2011/2012 Sample Solution Delivery of the project handout: Wednesday,
More information