Finite volume method for unsteady flow


Unsteady flows are characterized as flows in which the properties of the fluid are time dependent. It gets reflected in the governing equations as the time derivative of the properties are absent. For Studying Finite-volume method for unsteady flow there is some governing equations [1]>

Governing Equation Edit

The conservation equation for the transport of a scalar in unsteady flow has the general form as [2]


  is density and   is conservative form of all fluid flow,
  is the Diffusion coefficient and   is the Source term.   is Net rate of flow of   out of fluid element(convection),
  is Rate of increase of   due to diffusion,
  is Rate of increase of   due to sources.

  is Rate of increase of   of fluid element(transient),

The first term of the equation reflects the unsteadiness of the flow and is absent in case of steady flows. The finite volume integration of the governing equation is carried out over a control volume and also over a finite time step ∆t.


The control volume integration of the steady part of the equation is similar to the steady state governing equation's integration. We need to focus on the integration of the unsteady component of the equation. To get a feel of the integration technique, we refer to the one-dimensional unsteady heat conduction equation.[3]




Now, holding the assumption of the temperature at the node being prevalent in the entire control volume, the left side of the equation can be written as [4]


By using a first order backward differencing scheme, we can write the right hand side of the equation as


Now to evaluate the right hand side of the equation we use a weighting parameter   between 0 and 1, and we write the integration of  


Now, the exact form of the final discretised equation depends on the value of  . As the variance of   is 0<   <1, the scheme to be used to calculate   depends on the value of the  

Different Schemes Edit

1. Explicit Scheme in the explicit scheme the source term is linearised as  . We substitute   to get the explicit discretisation i.e.:[5]


where  . One thing worth noting is that the right side contains values at the old time step and hence the left side can be calculated by forward matching in time. The scheme is based on backward differencing and its Taylor series truncation error is first order with respect to time. All coefficients need to be positive. For constant k and uniform grid spacing,   this condition may be written as


This inequality sets a stringent condition on the maximum time step that can be used and represents a serious limitation on the scheme. It becomes very expensive to improve the spatial accuracy because the maximum possible time step needs to be reduced as the square of   [6]

2. Crank-Nicolson scheme : the Crank-Nicolson method results from setting  . The discretised unsteady heat conduction equation becomes



Since more than one unknown value of T at the new time level is present in equation the method is implicit and simultaneous equations for all node points need to be solved at each time step. Although schemes with   including the Crank-Nicolson scheme, are unconditionally stable for all values of the time step it is more important to ensure that all coefficients are positive for physically realistic and bounded results. This is the case if the coefficient of   satisfies the following condition


which leads to


the Crank-Nicolson is based on central differencing and hence is second order accurate in time. The overall accuracy of a computation depends also on the spatial differencing practice, so the Crank-Nicolson scheme is normally used in conjunction with spatial central differencing

3. Fully implicit Scheme when the value of Ѳ is set to 1 we get the fully implicit scheme. The discretised equation is: [7]



Both sides of the equation contain temperatures at the new time step, and a system of algebraic equations must be solved at each time level. The time marching procedure starts with a given initial field of temperatures  . The system of equations is solved after selecting time step  . Next the solution   is assigned to   and the procedure is repeated to progress the solution by a further time step. It can be seen that all coefficients are positive, which makes the implicit scheme unconditionally stable for any size of time step. Since the accuracy of the scheme is only first-order in time, small time steps are needed to ensure the accuracy of results. The implicit method is recommended for general purpose transient calculations because of its robustness and unconditional stability

References Edit

  1. ^ Retrieved November 10, 2013. {{cite web}}: Missing or empty |title= (help)[dead link]
  2. ^ An Introduction to Computational Fluid Dynamics H. K. Versteeg and W Malalasekra Chapter 8 page 168
  3. ^ An Introduction to Computational Fluid Dynamics H. K. Versteeg and W Malalasekera Chapter 8 page 169
  4. ^ Kim, Dongjoo; Choi, Haecheon (2000-08-10). "A Second-Order Time-Accurate Finite Volume Method for Unsteady Incompressible Flow on Hybrid Unstructured Grids". Journal of Computational Physics. 162 (2): 411–428. Bibcode:2000JCoPh.162..411K. doi:10.1006/jcph.2000.6546.
  5. ^ An Introduction to Computational Fluid Dynamics H. K. Versteeg and W Malalasekera Chapter 8 page 171
  6. ^ topic 7
  7. ^ topic 7