Numerical studies on quasilinear and linear elliptic equations

We explore here the acceleration of convergence of iterative methods for the solution of a class of quasilinear and linear algebraic equations. The specific systems are the finite difference form of the Navier-Stokes equations and the energy equation for recirculating flows. The acceleration procedures considered are: the successive over relaxation scheme; several implicit methods; and a second-order procedure. A new implicit method—the alternating direction line iterative method—is proposed in this paper. The method combines the advantages of the line successive over relaxation and alternating direction implicit methods. The various methods are tested for their computational economy and accuracy on a typical recirculating flow situation. The numerical experiments show that the alternating direction line iterative method is the most economical method of solving the Navier-Stokes equations for all Reynolds numbers in the laminar regime. The usual ADI method is shown to be not so attractive for large Reynolds numbers because of the loss of diagonal dominance. This loss can however be restored by a suitable choice of the relaxation parameter, but at the cost of accuracy. The accuracy of the new procedure is comparable to that of the well-tested successive overrelaxation method and to the available results in the literature. The second-order procedure turns out to be the most efficient method for the solution of the linear energy equation.

A computer code for enclosed natural convection in the boundary layer regime

A computer code is developed for the numerical prediction of natural convection in rectangular two-dimensional cavities at high Rayleigh numbers. The governing equations are retained in the primitive variable form. The numerical method is based on finite differences and an ADI scheme. Convective terms may be approximated with either central or hybrid differencing for greater stability. A non-uniform grid distribution is possible for greater efficiency. The pressure is dealt with via a SIMPLE type algorithm and the use of a fast elliptic solver for the solenoidal velocity correction field significantly reduces computing times. Preliminary results indicate that the code is reasonably accurate, robust and fast compared with existing benchmarks and finite difference based codes, particularly at high Rayleigh numbers. Extension to three-dimensional problems and turbulence studies in similar geometries is readily possible and indicated.

Unsteady laminar boundary layers in a compressible stagnation flow

The unsteady laminar compressible boundary-layer flow in the immediate vicinity of a two-dimensional stagnation point due to an incident stream whose velocity varies arbitrarily with time is considered. The governing partial differential equations, involving both time and the independent similarity variable, are transformed into new co-ordinates with finite ranges by means of a transformation which maps an infinite interval into a finite one. The resulting equations are solved by converting them into a matrix equation through the application of implicit finite-difference formulae. Computations have been carried out for two particular unsteady free-stream velocity distributions: (1) a constantly accelerating stream and (2) a fluctuating stream. The results show that in the former case both the skin-friction and the heat-transfer parameter increase steadily with time after a certain instant, while in the latter they oscillate thus responding to the fluctuations in the free-stream velocity.

Nonsimilar incompressible laminar boundary layers with magnetic field.

The solution of the steady laminar incompressible nonsimilar magneto-hydrodynamic boundary layer flow and heat transfer problem with viscous dissipation for electrically conducting fluids over two-dimensional and axisymmetric bodies with pressure gradient and magnetic field has been presented. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The computations have been carried out for flow over a cylinder and a sphere. The results indicate that the magnetic field tends to delay or prevent separation. The heat transfer strongly depends on the viscous dissipation parameter. When the dissipation parameter is positive (i.e. when the temperature of the wall is greater than the freestream temperature) and exceeds a certain value, the hot wall ceases to be cooled by the stream of cooler air because the ‘heat cushion’ provided by the frictional heat prevents cooling whereas the effect of the magnetic field is to remove the ‘heat cushion’ so that the wall continues to be cooled. The results are found to be in good agreement with those of the local similarity and local nonsimilarity methods except near the point of separation, but they are in excellent agreement with those of the difference-differential technique even near the point of separation.

Nonsimilar incompressible laminar boundary-layer flows in micropolar fluids

The solution of the steady laminar incompressible nonsimilar boundary-layer problem for micropolar fluids over two-dimensional and axisymmetric bodies has been presented. The partial differential equations governing the flow have been transformed into new co-ordinates having finite range. The resulting equations have been solved numerically using implicit finite-difference scheme. The computations have been carried out for a cylinder and a sphere. The results indicate that the separation in micropolar fluids occurs at earlier streamwise locations as compared to Newtonian fluids. The skin friction and velocity profiles depend on the shape of the body and are almost insensitive to microrotation or coupling parameter, provided the coupling parameter is small. On the other hand, the microrotation profiles and microrotation gradient depend on the microrotation parameter and they are insensitive to the coupling parameter.

Heat transfer to power-law fluids in thermal entrance region with viscous dissipation for constant heat flux conditions

An analytical solution of the heat transfer problem with viscous dissipation for non-Newtonian fluids with power-law model in the thermal entrance region of a circular pipe and two parallel plates under constant heat flux conditions is obtained using eigenvalue approach by suitably replacing one of the boundary conditions by total energy balance equation. Analytical expressions for the wall and the bulk temperatures and the local Nusselt number are presented. The results are in close agreement with those obtained by implicit finite-difference scheme. It is found that the role of viscous dissipation on heat transfer is completely different for heating and cooling conditions at the wall. The results for the case of cooling at the wall are of interest in the design of the oil pipe line.

Effect of soil variability on the bearing capacity of clay and in slope stability problems

Site-specific geotechnical data are always random and variable in space. In the present study, a procedure for quantifying the variability in geotechnical characterization and design parameters is discussed using the site-specific cone tip resistance data (qc) obtained from static cone penetration test (SCPT). The parameters for the spatial variability modeling of geotechnical parameters i.e. (i) existing trend function in the in situ qc data; (ii) second moment statistics i.e. analysis of mean, variance, and auto-correlation structure of the soil strength and stiffness parameters; and (iii) inputs from the spatial correlation analysis, are utilized in the numerical modeling procedures using the finite difference numerical code FLAC 5.0. The influence of consideration of spatially variable soil parameters on the reliability-based geotechnical deign is studied for the two cases i.e. (a) bearing capacity analysis of a shallow foundation resting on a clayey soil, and (b) analysis of stability and deformation pattern of a cohesive-frictional soil slope. The study highlights the procedure for conducting a site-specific study using field test data such as SCPT in geotechnical analysis and demonstrates that a few additional computations involving soil variability provide a better insight into the role of variability in designs.

Nonsimilar laminar incompressible boundary layers with vectored mass transfer

The effect of vectored mass transfer on the flow and heat transfer of the steady laminar incompressible nonsimilar boundary layer with viscous dissipation for two-dimensional and axisymmetric porous bodies with pressure gradient has been studied. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The computations have been carried out for a cylinder and a sphere. The skin friction is strongly influenced by the vectored mass transfer, and the heat transfer both by the vectored mass transfer and dissipation parameter. It is observed that the vectored suction tends to delay the separation whereas the effect of the vectored injection is just the reverse. Our results agree with those of the local nonsimilarity, difference-differential and asymptotic methods but not with those of the local similarity method.

Compressible magnetohydrodynamic boundary layer swirling nozzle and diffuser flows with a magnetic field

The flow, heat and mass transfer problem for boundary layer swirling flow of a laminar steady compressible electrically conducting gas with variable properties through a conical nozzle and a diffuser with an applied magnetic field has been studied. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme after they have been transformed into dimensionless form using the modified Lees transformation. The results indicate that the skin friction and heat transfer strongly depend on the magnetic field, mass transfer and variation of the density-viscosity product across the boundary layer. However, the effect of the variation of the density-viscosity product is more pronounced in the case of a nozzle than in the case of a diffuser. It has been found that large swirl is required to produce strong effect on the skin friction and heat transfer. Separationless flow along the entire length of the diffuser can be obtained by applying appropriate amount of suction. The results are found to be in good agreement with those of the local nonsimilarity method, but they differ quite significantly from those of the local similarity method.

Unsteady Incompressible Three-Dimensional Asymmetric Stagnation-Point Boundary Layers

The unsteady laminar incompressible boundary-layer flow near the three-dimensional asymmetric stagnation point has been studied under the assumptions that the free-stream velocity, wall temperature, and surface mass transfer vary arbitrarily with time. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. It is found that in contrast with the symmetric flow, the maximum heat transfer occurs away from the stagnation point due to the decrease in the boundary-layer thickness. The effect of the variation of the wall temperature with time on heat transfer is strong. The skin friction and heat transfer due to asymmetric flow only are comparatively less affected by the mass transfer as compared to those of symmetric flow.

Steady mixed convection stagnation-point flow of upper convected Maxwell fluids with magnetic field

The steady MHD mixed convection flow of a viscoelastic fluid in the vicinity of two-dimensional stagnation point with magnetic field has been investigated under the assumption that the fluid obeys the upper-convected Maxwell (UCM) model. Boundary layer theory is used to simplify the equations of motion. induced magnetic field and energy which results in three coupled non-linear ordinary differential equations which are well-posed. These equations have been solved by using finite difference method. The results indicate the reduction in the surface velocity gradient, surface heat transfer and displacement thickness with the increase in the elasticity number. These trends are opposite to those reported in the literature for a second-grade fluid. The surface velocity gradient and heat transfer are enhanced by the magnetic and buoyancy parameters. The surface heat transfer increases with the Prandtl number, but the surface velocity gradient decreases.

Influence of spatial variability of permeability property on steady state seepage flow and slope stability analysis

In recent years, spatial variability modeling of soil parameters using random field theory has gained distinct importance in geotechnical analysis. In the present Study, commercially available finite difference numerical code FLAC 5.0 is used for modeling the permeability parameter as spatially correlated log-normally distributed random variable and its influence on the steady state seepage flow and on the slope stability analysis are studied. Considering the case of a 5.0 m high cohesive-frictional soil slope of 30 degrees, a range of coefficients of variation (CoV%) from 60 to 90% in the permeability Values, and taking different values of correlation distance in the range of 0.5-15 m, parametric studies, using Monte Carlo simulations, are performed to study the following three aspects, i.e., (i) effect ostochastic soil permeability on the statistics of seepage flow in comparison to the analytic (Dupuit's) solution available for the uniformly constant permeability property; (ii) strain and deformation pattern, and (iii) stability of the given slope assessed in terms of factor of safety (FS). The results obtained in this study are useful to understand the role of permeability variations in slope stability analysis under different slope conditions and material properties. (C) 2009 Elsevier B.V. All rights reserved.

Conservation law with the flux function discontinuous in the space variable- II - Convex-concave type fluxes and generalized entropy solutions

We deal with a single conservation law with discontinuous convex-concave type fluxes which arise while considering sign changing flux coefficients. The main difficulty is that a weak solution may not exist as the Rankine-Hugoniot condition at the interface may not be satisfied for certain choice of the initial data. We develop the concept of generalized entropy solutions for such equations by replacing the Rankine-Hugoniot condition by a generalized Rankine-Hugoniot condition. The uniqueness of solutions is shown by proving that the generalized entropy solutions form a contractive semi-group in L-1. Existence follows by showing that a Godunov type finite difference scheme converges to the generalized entropy solution. The scheme is based on solutions of the associated Riemann problem and is neither consistent nor conservative. The analysis developed here enables to treat the cases of fluxes having at most one extrema in the domain of definition completely. Numerical results reporting the performance of the scheme are presented. (C) 2006 Elsevier B.V. All rights reserved.

Simulation of Dust Loaded V-I Characteristics of a Commercial Thermal Power Plant Precipitator

A new approach based on finite difference method, is proposed for the simulation of electrical conditions in a dc energized wire-duct electrostatic precipitator with and without dust loading. Simulated voltage-curren characteristics with and without dust loading were compared with the measured characteristics for analyzing the performance of a precipitator. The simple finite difference method gives sufficiently accurate results with reduced mesh size. The results for dust free simulation were validated with published experimental data. Further measurements were conducted at a thermal power plant in India and the results compares well with the measured ones.

Unsteady MHD mixed convection flow over an impulsively stretched permeable vertical surface in a quiescent fluid

The unsteady mixed convection flow of an incompressible laminar electrically conducting fluid over an impulsively stretched permeable vertical surface in an unbounded quiescent fluid in the presence of a transverse magnetic field has been investigated. At the same time, the surface temperature is suddenly increased from the surrounding fluid temperature or a constant heat flux is suddenly imposed on the surface. The problem is formulated in such a way that for small time it is governed by Rayleigh type of equation and for large time by Crane type of equation. The non-linear coupled parabolic partial differential equations governing the unsteady mixed convection flow under boundary layer approximations have been solved analytically by using the homotopy analysis method as well as numerically by an implicit finite difference scheme. The local skin friction coefficient and the local Nusselt number are found to decrease rapidly with time in a small time interval and they tend to steady-state values for t* >= 5. They also increase with the buoyancy force and suction, but decrease with injection rate. The local skin friction coefficient increases with the magnetic field, but the local Nusselt number decreases. There is a smooth transition from the unsteady state to the steady state. (C) 2010 Elsevier Ltd. All rights reserved.