Implementation of high‐order compact finite‐difference method to parabolized Navier–Stokes schemes

The numerical solution to the parabolized Navier–Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth‐order compact finite‐difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite‐difference algorithm of Beam and Warming type with a high‐order compact accuracy. A shock‐fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously computed with the fourth‐order accuracy. The computations are carried out for a benchmark case: hypersonic axisymmetric flow over a blunt cone at Mach 8. A sensitivity study is performed for the basic flowfield, including profiles and their derivatives obtained from the fourth‐order compact PNS and IPNS solutions, and the effects of grid size and numerical dissipation term used are discussed. The present results for the flowfield variables and also their derivatives are compared with those of other basic flow models to demonstrate the accuracy and efficiency of the proposed method. The present work represents the first known application of a high‐order compact finite‐difference method to the PNS schemes, which are computationally more efficient than Navier–Stokes solutions. Copyright © 2008 John Wiley & Sons, Ltd.     

The nonlinear response of Cattaneo-type thermal loading of a laser pulse on a medium using the generalized thermoelastic model

The nonlinear thermoelastic responses of an elastic medium exposed to laser generated shortpulse heating are investigated in this article. The thermal wave propagation of generalized thermoelastic medium under the impact of thermal loading with energy dissipation is the focus of this research. To model the thermal boundary condition(in the form of thermal conduction),generalized Cattaneo model(GCM) is employed. In the reference configuration, a nonlinear coupled Lord-Shulman-type generalized thermoelasticity formulation using finite strain theory(FST) is developed and the temperature dependency of the thermal conductivity is considered to derive the equations. In order to solve the time-dependent and nonlinear equations, Newmark's numerical time integration technique and an updated finite element algorithm is applied and to ensure achieving accurate continuity of the results, the Hermitian elements are used instead of Lagrangian's. The numerical responses for different factors such as input heat flux and nonlinear terms are expressed graphically and their impacts on the system's reaction are discussed in detail.The results of the study are presented for Green–Lindsay model and the findings are compared with Lord-Shulman model especially with regards to heat wave propagation. It is shown that the nature of the laser's thermal shock and its geometry are particularly determinative in the final stage of deformation. The research also concluded that employing FST leads to achieving more accuracy in terms of elastic deformations; however, the thermally nonlinear analysis does not change the results markedly. For this reason, the nonlinear theory of deformation is required in laser related reviews, while it is reasonable to ignore the temperature changes compared to the reference temperature in deriving governing equations.     

Finite volume scheme for the solution of fluid flow problems on unstructured non‐staggered grids

A recently developed non‐staggered methodology which uses the principle of applying fourth‐order dissipation to the governing pressure‐correction equation is developed so it can be applied to unstructured grids. A finite volume methodology is used for discretization. The fourth‐order dissipation term is found using second‐order gradient operators. This makes it straightforward to incorporate the dissipation term on unstructured grids. The new methodology is compared with solutions from a standard finite volume second‐order flow solver and is also tested for a standard laminar driven‐lid flow problem with grids systems that do not have a uniform structure. Finally, we demonstrate how the new methodology can be used to predict flow over a wavy boundary. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical simulation of shock wave propagation in microchannels using continuum and kinetic approaches

Numerical simulations of shock wave propagation in microchannels and microtubes (viscous shock tube problem) have been performed using three different approaches: the Navier–Stokes equations with the velocity slip and temperature jump boundary conditions, the statistical Direct Simulation Monte Carlo method for the Boltzmann equation, and the model kinetic Bhatnagar–Gross–Krook equation with the Shakhov equilibrium distribution function. Effects of flow rarefaction and dissipation are investigated and the results obtained with different approaches are compared. A parametric study of the problem for different Knudsen numbers and initial shock strengths is carried out using the Navier–Stokes computations.      

A method for solving the factorized vorticity-stream function equations by finite elements

A new finite element method for solving the time-dependent incompressible Navier-Stokes equations with general boundary conditions is presented. The two second-order partial differential equations for the vorticity and the stream function are factorized, apart from the non-linear advection term, by eliminating the coupling due to the double specification on the stream function at (a part of) the boundary. This is achieved by reducing the no-slip boundary conditions to projection integral conditions for the vorticity field and by evaluating the relevant quantities involved according to an extension of the method of Glowinski and Pironneau for the biharmonic problem. Time integration schemes and iterative algorithms are introduced which require the solution only of banded linear systems of symmetric type. The proposed finite element formulation is compared with its finite difference equivalent by means of a few numerical examples. The results obtained using 4-noded bilinear elements provide an illustration of the superiority of the finite element based spatial discretization.     

Postprocessing algorithm for particle-tracking velocimetry based on ellipsoidal equations

 Accurate postprocessing methods are required in order to analyze the detailed flow structures from the scattered data of particle-tracking velocimetry (PTV). In particular, vorticity distributions and stream functions are not reasonably obtained by conventional methods. This paper proposes a new postprocessing algorithm based on ellipsoidal differential equations; this method utilizes data as discrete boundary conditions. The results obtained by the proposed algorithm fully satisfy the equation of continuity and simultaneously correspond well with the raw data. The performance of the algorithm is examined by applying it to two-dimensional vortex flows and isotropic turbulent flows. The results reveal that the present algorithm has the highest accuracy among several conventional methods for detecting vorticity and streamlines. Received: 29 January 2001 / Accepted: 7 August 2001     

Comment on the paper “Transient MHD free convective flow past an infinite vertical plate embedded in a porous medium with viscous dissipation,Siva Reddy Sheri,R. Srinivasa Raju,Meccanica, DOI 10.1007/s11012-015-0285-y”

A numerical investigation of transient magnetohydrodynamic free convection flow past an infinite vertical plate embedded in a porous medium with viscous dissipation is presented in the above paper. The governing differential equations are transformed into a set of non-linear coupled partial differential equations and are solved numerically using the finite element method. Numerical results for the velocity, temperature and concentration profiles within the boundary layer are presented and discussed.     

A fractional-step method for the incompressible Navier–Stokes equations related to a predictor–multicorrector algorithm

An implicit fractional-step method for the numerical solution of the time-dependent incompressible Navier–Stokes equations in primitive variables is studied in this paper. The method, which is first-order-accurate in the time step, is shown to converge to an exact solution of the equations. By adequately splitting the viscous term, it allows the enforcement of full Dirichlet boundary conditions on the velocity in all substeps of the scheme, unlike standard projection methods. The consideration of this method was actually motivated by the study of a well-known predictor–multicorrector algorithm, when this is applied to the incompressible Navier–Stokes equations. A new derivation of the algorithm in a general setting is provided, showing in what sense it can also be understood as a fractional-step method; this justifies, in particular, why the original boundary conditions of the problem can be enforced in this algorithm. Two different finite element interpolations are considered for the space discretization, and numerical results obtained with them for standard benchmark cases are presented. © 1998 John Wiley & Sons, Ltd.     

Wave induced thermal boundary layers in a compressible fluid: analysis and numerical simulations

The response of a semi-infinite compressible fluid to a step-wise change in temperature of its boundary is investigated analytically and numerically. Numerical results of the boundary layer structure are compared with Clarke’s analytical solution for a gas with thermal conductivity proportional to temperature. To avoid unwanted numerical dissipation in the numerical analysis, the space-time conservation element and solution element (CESE) method has been adopted to solve the unsteady 1-D Navier-Stokes equations. Good agreement between analytical and numerical results has been found for the development of the thermal boundary layer on a long time scale. Weak shock waves and expansion waves induced by the thermal boundary layer due to its compressibility, are observed in the numerical simulation. Finally, the numerical method has been applied to the reflection of a non-linear expansion wave and to a shock wave from an isothermal wall, thereby illustrating the effect of the boundary layer on the external flow field.     

Free Convection Boundary Layer Flow Past a Vertical Surface in a Porous Medium with Temperature-Dependent Properties

Numerical solutions for the free convection heat transfer in a viscous fluid at a permeable surface embedded in a saturated porous medium, in the presence of viscous dissipation with temperature-dependent variable fluid properties, are obtained. The governing equations for the problem are derived using the Darcy model and the Boussinesq approximation (with nonlinear density temperature variation in the buoyancy force term). The coupled non-linearities arising from the temperature-dependent density, viscosity, thermal conductivity, and viscous dissipation are included. The partial differential equations of the model are reduced to ordinary differential equations by a similarity transformation and the resulting coupled, nonlinear ordinary differential equations are solved numerically by a second order finite difference scheme for several sets of values of the parameters. Also, asymptotic results are obtained for large values of | f _w|. Moreover, the numerical results for the velocity, the temperature, and the wall-temperature gradient are presented through graphs and tables, and are discussed. It is observed that by increasing the fluid variable viscosity parameter, one could reduce the velocity and thermal boundary layer thickness. However, quite the opposite is true with the non-linear density temperature variation parameter.     

Flux‐corrected transport technique for open channel flow

In modeling flow in open channels, the traditional finite difference/finite volume schemes become inefficient and warrant special numerical treatment in the presence of shocks and discontinuities. The numerical oscillations that arise by making use of a second‐ and higher‐order schemes require some additional smoothing mechanism. A characteristic feature of high‐resolution schemes lies in smooth capturing of the shock fronts. This paper provides a general formulation for a flux‐corrected transport algorithm to the one‐dimensional open channel flow equations. The preliminary results presented show that the present algorithm is an efficient, conservative and robust tool that can be easily coded. To demonstrate the robustness of the present formulation, results are compared with other published numerical results, experimental data and analytical solutions when available. In particular, a comprehensive study on the effect of the source term, dry bed, variable width channel, steep sloping channel and flow with mixed flow conditions (as in a hydraulic jump) has been carried out to test the efficacy of the present algorithm. Copyright © 1999 John Wiley & Sons, Ltd.     

A 3D incompressible Navier–Stokes velocity–vorticity weak form finite element algorithm

The velocity–vorticity formulation is selected to develop a time‐accurate CFD finite element algorithm for the incompressible Navier–Stokes equations in three dimensions.The finite element implementation uses equal order trilinear finite elements on a non‐staggered hexahedral mesh. A second order vorticity kinematic boundary condition is derived for the no slip wall boundary condition which also enforces the incompressibility constraint. A biconjugate gradient stabilized (BiCGSTAB) sparse iterative solver is utilized to solve the fully coupled system of equations as a Newton algorithm. The solver yields an efficient parallel solution algorithm on distributed‐memory machines, such as the IBM SP2. Three dimensional laminar flow solutions for a square channel, a lid‐driven cavity, and a thermal cavity are established and compared with available benchmark solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

THE APPLICATION OF A FINITE DIFFERENCE/FINITE VOLUME ALGORITHM TO SOLVE MAXWELL'S EQUATIONS IN THE TIME DOMAIN

Abstract

A finite volume/finite difference method based on Ni's multigrid formulation is introduced for the solution of Maxwell's equations. The scheme is presented for the cases of transverse magnetic scattering from two-dimensional circular and square cylinders, as well as from NACA 0012 airfoil. The codes are validated against the traditional Method of Moments, which is analogous to a panel method in CFD. The circular cylinder scattering is compared to the analytical series solution for better understanding how the roles of numerical dispersion and dissipation errors affect the solution. The reflecting boundary conditions are modeled by the idea of inducing fields inside the conductor and a method of modeling the singularities that arise at a sharp corner is presented. Absorbing boundary conditions are modeled by integrating along the characteristic compatibility equations in the direction of the outgoing wave.     

NND格式在非结构网格中的推广

在张涵信提出的无波动、无自由参数的差分格式（ＮＮＤ格式）的基础上，构造了适用于非结构网格的二阶精度ＮＮＤ有限体积格式，解决了现有非结构网格方法中为抑制激波附近的波动而必须引入含自由参数的人工粘性项的困难，并采用网格自适应技术以提高效率．通过对二维平板激波反射和前台阶在管道内的流动问题的计算，表明本方法可有效地用于Ｅｕｌｅｒ方程的求解．     

Contact Interaction of the Faces of a Rectangular Crack under Normally Incident Tension–Compression Waves

A three-dimensional problem on the contact interaction between the faces of a rectangular crack under a normally incident harmonic tension–compression wave is considered. The problem is solved by using the method of boundary integral equations and an iterative algorithm. The contact forces and the discontinuity in the displacement of the crack faces are studied. The results obtained are compared with those for a finite plane crack.     

Analysis of bonded anisotropic wedges with interface crack under anti-plane shear loading

The antiplane stress analysis of two anisotropic finite wedges with arbitrary radii and apex angles that are bonded together along a common edge is investigated. The wedge radial boundaries can be subjected to displacement-displacement boundary condi- tions, and the circular boundary of the wedge is free from any traction. The new finite complex transforms are employed to solve the problem. These finite complex transforms have complex analogies to both kinds of standard finite Mellin transforms. The traction free condition on the crack faces is expressed as a singular integral equation by using the exact analytical method. The explicit terms for the strength of singularity are extracted, showing the dependence of the order of the stress singularity on the wedge angle, material constants, and boundary conditions. A numerical method is used for solving the resul- tant singular integral equations. The displacement boundary condition may be a general term of the Taylor series expansion for the displacement prescribed on the radial edge of the wedge. Thus, the analysis of every kind of displacement boundary conditions can be obtained by the achieved results from the foregoing general displacement boundary condition. The obtained stress intensity factors （SIFs） at the crack tips are plotted and compared with those obtained by the finite element analysis （FEA）.     

A moving boundary hyperbolic problem for a stress impact in a bar of rate-type material

In this paper we present some results obtained by studying the mathematical model describing a moving boundary hyperbolic problem related to a time dependent stress impact in a bar of Maxwell-like material. Due to the impact a shock front propagates with a finite speed. Here our interest is to underline the influence of the dissipative term on the propagation of the shock front.

In the framework of the similarity analysis we are able to reduce the moving boundary hyperbolic problem to a free boundary value problem for an ordinary differential system. It is then possible, by applying two numerical transformation methods, to solve the free boundary value problem numerically. The influence of the dissipative term is evident: the free boundary (that defines the shock front propagation) is an increasing function of the dissipative coefficient.     

An investigation of artificial dissipation for the cell-vertex finite volume method

This paper is concerned with an investigation of artificial dissipation models that are used with the cell-vertex finite volume approximation of the compressible Euler and Navier-Stokes equations. Based on the observation that first and second-order upwind schemes can be written as a central discretization plus an appropriately scaled dissipative flux, a matrix scaling of second and fourth-differences is implemented in an artificial dissipation model that also uses a procedure to limit the applicability of shock-capturing dissipation. Analysis of the model and the associated boundary conditions is given to determine the effect on accuracy. Numerical results are given for transonic Euler flow past a NACA0012 aerofoil profile which demonstrate the improved shock-capturing capability of the model. Results for laminar subsonic viscous flow over a flat plate show that the matrix-dissipation model reduces the amount of spurious artificial dissipation within boundary layers.     

On the performance of adaptive mesh refinement computation

A finite difference scheme for the unsteady Euler equations using an adaptive mesh refinement (AMR) algorithm was applied to the time-dependent flowfield of shock diffraction problems. The effectiveness of the AMR algorithm was evaluated against a uniform mesh algorithm. Computational results showed that to obtain solutions with equivalent resolution, the AMR algorithm requires much less processing time, when compared with a uniform mesh algorithm.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Analysis of shock wave reflection from fixed and moving boundaries using a stabilized particle method

In the present paper, the efficiency of an enhanced formulation of the stabilized corrective smoothed particle method （CSPM） for simulation of shock wave propagation and reflection from fixed and moving solid boundaries in compressible fluids is investigated. The Lagrangian nature and its accuracy for imposing the boundary conditions are the two main reasons for adoption of CSPM. The governing equations are further modified for imposition of moving solid boundary conditions. In addition to the traditional artificial viscosity, which can remove numerically induced abnormal jumps in the field values, a velocity field smoothing technique is introduced as an efficient method for stabilizing the solution. The method has been implemented for one- and two-dimensional shock wave propagation and reflection from fixed and moving boundaries and the results have been compared with other available solutions. The method has also been adopted for simulation of shock wave propagation and reflection from infinite and finite solid boundaries.