Search results

The purpose of this paper is to provide a brief introduction of the finite difference based parabolic equation (PE) modeling to the advanced engineering students and academic researchers.

A three-dimensional parabolic equation (3DPE) model is developed from the ground up for modeling wave propagation in the tunnel via a rectangular waveguide structure. A discussion of vector wave equations from Maxwell’s equations followed by the paraxial approximations and finite difference implementation is presented for the beginners. The obtained simulation results are compared with the analytical solution.

It is shown that the alternating direction implicit finite difference method (FDM) is more efficient in terms of accuracy, computational time and memory than the explicit FDM. The reader interested in maximum details of individual contributions such as the latest achievements in PE modeling until 2021, basic PE derivation, PE formulation’s approximations, finite difference discretization and implementation of 3DPE, can learn from this paper.

For the purpose of this paper, a simple 3DPE formulation is presented. For simplicity, a rectangular waveguide structure is discretized with the finite difference approach as a design problem. Future work could use the PE based FDM to study the possibility of utilization of meteorological techniques, including the effects of backward traveling waves as well as making comparisons with the experimental data.

The proposed work is directly applicable to typical problems in the field of tunnel propagation modeling for both national commercial and military applications.

The purpose of this paper is to both determine the effects of the nonlinearity on the wave dynamics and assess the temporal and spatial accuracy of five finite difference methods for the solution of the inviscid generalized regularized long-wave (GRLW) equation subject to initial Gaussian conditions.

Two implicit second- and fourth-order accurate finite difference methods and three Runge-Kutta procedures are introduced. The methods employ a new dependent variable which contains the wave amplitude and its second-order spatial derivative. Numerical experiments are reported for several temporal and spatial step sizes in order to assess their accuracy and the preservation of the first two invariants of the inviscid GRLW equation as functions of the spatial and temporal orders of accuracy, and thus determine the conditions under which grid-independent results are obtained.

It has been found that the steepening of the wave increase as the nonlinearity exponent is increased and that the accuracy of the fourth-order Runge-Kutta method is comparable to that of a second-order implicit procedure for time steps smaller than 100th, and that only the fourth-order compact method is almost grid-independent if the time step is on the order of 1,000th and more than 5,000 grid points are used, because of the initial steepening of the initial profile, wave breakup and solitary wave propagation.

This is the first study where an accuracy assessment of wave breakup of the inviscid GRLW equation subject to initial Gaussian conditions is reported.

Detailed results of numerical calculations of transient, 2D incompressible flow around and in the wake of a square prism at Re = 100, 200 and 500 are presented. An implicit finite‐difference operator‐splitting method, a version of the known SIMPLEC‐like method on a staggered grid, is described. Appropriate theoretical results are presented. The method has second‐order accuracy in space, conserving mass, momentum and kinetic energy. A new modification of the multigrid method is employed to solve the elliptic pressure problem. Calculations are performed on a sequence of spatial grids with up to 401 × 321 grid points, at sequentially halved time steps to ensure grid‐independent results. Three types of flow are shown to exist at Re = 500: a steady‐state unstable flow and two which are transient, fully periodic and asymmetric about the centre line but mirror symmetric to each other. Discrete frequency spectra of drag and lift coefficients are presented.

This paper aims to improve the computational efficiency of higher-order accurate Noye–Hayman [NH (9,9)] implicit finite difference scheme for the solution of electromagnetic scattering problems in tunnel environments.

The proposed method consists of two major steps: First, the higher-order NH (9,9) scheme is numerically discretized using the finite-difference method. The second step is to use an algorithm based on hierarchical interpolative factorization (HIF) to accelerate the solution of this scheme.

It is observed that the simulation results obtained from the numerical tests illustrate very high accuracy of the NH (9,9) method in typical tunnel environments. HIF algorithm makes the NH (9,9) method computationally efficient for two-dimensional (2D) or three-dimensional (3D) problems. The proposed method could help in reducing the computational cost of the NH (9,9) method very close to O(n) usual O(n3) for a full matrix.

For simplicity, in this study, perfect electric conductor boundary conditions are considered. Future research may also include the utilization of meteorological techniques, including the effects of backward traveling waves, and make comparisons with the experimental data.

This study is directly applicable to typical problems in the field of tunnel propagation modeling for both national commercial and military applications.

The problem of laminar natural convection from a vertical circular cone maintained at either a uniform surface temperature or a uniform surface heat flux, and placed in a thermally stratified medium is considered. The governing non‐similarity boundary layer equation for uniform surface temperature are analyzed by using two distinct solution methodologies; namely, (i) a finite difference method and (ii) a local non‐similarity method. For uniform surface heat flux case, the solutions of the governing non‐similarity boundary layer equations are obtained by using three distinct solution methodologies, namely, (i) a finite difference method, (ii) a series solution method and (iii) an asymptotic solution method. The solutions are presented in terms of local skin‐friction and local Nusselt number for different values of Prandtl number and are displayed graphically. Effects of variations in the Prandtl number and stratification parameter on the velocity and temperature profiles are also shown graphically. Solutions obtained by finite difference method are compared with the other methods and found to be in excellent agreement.

We determine the effects of micro‐inertia density and the vortex viscosity on laminar free convection boundary layer flow of a thermomicropolar fluid past a vertical plate with exponentially varying surface temperature as well as surface heat flux. The governing nonsimilarity boundary layer equations are analyzed using: first, a series solution for small ξ (a scaled streamwise distribution of micro‐inertia density), second, an asymptotic solution for large ξ and, third, a full numerical solution implicit finite difference method together with Keller‐box scheme. Results are expressed in terms of local skin friction and local Nusselt number. The effects of varying the vortex viscosity parameter, Δ, surface temperature and the surface heat flux gradient n and m respectively against ξ for fluids having Prandtl number equals 0.72 and 7.0 are determined.

A mixed implicit semi Lagrangian finite difference‐finite volume method for numerical simulation of 2D air motion inside cylinders is derived and discussed. A conformal mapping from a physical (moving) domain to a computational (fixed) one is resorted in order to deal with a grid independent of time, making the numerical code very efficient. The numerical method is mass and energy conservative, unconditionally stable and at each timestep requires the solution of two well structured five‐band systems of linear equations. Its accuracy is first order in time and second one in space where the solution is smooth, while due to FCT space accuracy drops to the first order where the solution is steep. Stability of the method is proved both by a classical Von Newmann analysis and analysis of the matrices involved in the systems of linear equations. All these elements make the numerical method particularly fast. Numerical experiments are performed that show the influence of the maximum Courant number (with respect to the fluid speed) on the performance of the numerical method; moreover, comparison of simulations with a major existing code for engines is worked out.

A steady, two‐dimensional natural convection flow of a viscous, incompressible fluid having temperature‐dependent viscosity and thermal conductivity about a truncated cone is considered. We use suitable transformations to obtain the equations governing the flow in convenient form and integrate them by using an implicit finite difference method. Perturbation solutions are employed to obtain the solution in the regimes near and far away from the point of truncation. The results are obtained in terms of the local skin friction and the local Nusselt number. Perturbation solutions are compared with the finite difference solutions and found to be in excellent agreement. The dimensionless velocity, viscosity and thermal conductivity distributions are also displayed graphically, showing the effects of various values of the pertinent parameter for smaller values of Prandtl number.

The purpose of this paper is to study numerically the influence of a magnetic field on the transient free convective boundary layer flow of a nanofluid over a moving semi-infinite vertical cylinder with heat transfer

The problem is governed by the coupled non-linear partial differential equations with appropriate boundary conditions. The fluid is a water-based nanofluid containing nanoparticles of copper. The Brinkman model for dynamic viscosity and Maxwell–Garnett model for thermal conductivity are used. The governing boundary layer equations are written according to The Tiwari–Das nanofluid model. A robust, well-tested, implicit finite difference method of Crank–Nicolson type, which is unconditionally stable and convergent, is used to find the numerical solutions of the problem. The velocity and temperature profiles are studied for significant physical parameters such as the magnetic parameter, nanoparticles volume fraction and the thermal Grashof number Gr. The local skin-friction coefficient and the Nusselt number are also analysed and presented graphically.

The present computations have shown that an increase in the values of either magnetic parameter M or nanoparticle volume fraction decreases the local skin-friction coefficient, whereas the opposite effect is observed for thermal Grashof number Gr. The local Nusselt number increases with a rise in Gr and ϕ values. But an increase in M reduces the local Nusselt number.

This paper is relatively original and presents numerical investigation of transient two-dimensional laminar boundary layer free convective flow of a nanofluid over a moving semi-infinite vertical cylinder in the presence of an applied magnetic field. The present study is of immediate application to all those processes which are highly affected by heat enhancement concept and a magnetic field. Further the present study is relevant to nanofluid materials processing, chemical engineering coating operations exploiting nanomaterials and others.

The purpose of this paper is to focus on the numerical modelling of transient natural convection flow of an incompressible viscous nanofluid past an impulsively started semi-infinite vertical plate with variable surface temperature.

The problem is governed by the coupled non-linear partial differential equations with appropriate boundary conditions. A robust, well-tested, Crank-Nicolson type of implicit finite-difference method, which is unconditionally stable and convergent, is used to solve the governing non-linear set of partial differential equations.

The local and average values of the skin-friction coefficient (viscous drag) and the average Nusselt number (the rate of heat transfer) decreased, while the local Nusselt number increased for all nanofluids, namely, aluminium oxide-water, copper-water, titanium oxide-water and silver-water with an increase in the temperature exponent m. Selecting aluminium oxide as the dispersing nanoparticles leads to the maximum average Nusselt number (the rate of heat transfer), while choosing silver as the dispersing nanoparticles leads to the minimum local Nusselt number compared to the other nanofluids for all values of the temperature exponent m. Also, choosing silver as the dispersing nanoparticles leads to the minimum skin-friction coefficient (viscous drag), while selecting aluminium oxide as the dispersing nanoparticles leads to the maximum skin-friction coefficient (viscous drag) for all values of the temperature exponent m.

The Brinkman model for dynamic viscosity and Maxwell-Garnett model for thermal conductivity are employed. The governing boundary layer equations are written according to The Tiwari-Das nanofluid model. A range of nanofluids containing nanoparticles of aluminium oxide, copper, titanium oxide and silver with nanoparticle volume fraction range less than or equal to 0.04 are considered.

The present simulations are relevant to nanomaterials thermal flow processing in the chemical engineering and metallurgy industries. This study also provides an important benchmark for further simulations of nanofluid dynamic transport phenomena of relevance to materials processing, with alternative computational algorithms (e.g. finite element methods).

This paper is relatively original and illustrates the influence of variable surface temperature on transient natural convection flow of a viscous incompressible nanofluid and heat transfer from an impulsively started semi-infinite vertical plate.