High‐order upwind compact scheme and simulation of turbulent premixed V‐flame

A high‐order accurate upwind compact difference scheme with an optimal control coefficient is developed to track the flame front of a premixed V‐flame. In multi‐dimensional problems, dispersion effect appears in the form of anisotropy. By means of Fourier analysis of the operators, anisotropic effects of the upwind compact difference schemes are analysed. Based on a level set algorithm with the effect of exothermicity and baroclinicity, the flame front is tracked. The high‐order accurate upwind compact scheme is employed to approximate the level set equation. In order to suppress numerical oscillations, the group velocity control technique is used and the upwind compact difference scheme is combined with the random vortex method to simulate the turbulent premixed V‐flame. Distributions of velocities and flame brush thickness are obtained by this technique and found to be comparable with experimental measurement. Copyright © 2005 John Wiley & Sons, Ltd.     

Adaptive methods of lines for one-dimensional reaction-diffusion equations

Adaptive and non-adaptive finite difference methods are used to study one-dimensional reaction-diffusion equations whose solutions are characterized by the presence of steep, fast-moving flame fronts. Three non-adaptive techniques based on the methods of lines are described. The first technique uses a finite volume method and yields a system of non-linear, first-order, ordinary differential equations in time. The second technique uses time linearization, discretizes the time derivatives and yields a linear, second-order, ordinary differential equation in space, which is solved by means of three-point, fourth-order accurate, compact differences. The third technique takes advantage of the disparity in the time scales of the reaction and diffusion processes, splits the reaction--diffusion operator into a sequence of reaction and diffusion operators and solves the diffusion operator by means of either a finite volume method or a three-point, fourth-order accurate compact difference expression. The non-adaptive methods of lines presented in this paper may use equaliy or non-equally spaced fixed grids and require a large number of grid points to solve accurately one-dimensional problems characterized by the presence of steep, fast-moving fronts. Three adaptive methods for the solution of reaction-diffusion equations are considered. The first adaptive technique is static and uses a subequidistribution principle to determine the grid points, avoid mesh tangling and node overtaking and obtain smooth grids. The second adaptive technique is dynamic, uses an equidistribution principle with spatial and temporal smoothing and yields a system of first-order, non-linear, ordinary differential equations for the grid point motion. The third adaptive technique is hybrid, combines some features of static and dynamic methods, and uses a predictor-corrector strategy to predict the grid and solve for the dependent variables, respectively. The three adaptive techniques presented in this paper use physical co-ordinates and may employ finite volume or three-point, compact methods. The adaptive and non-adaptive finite difference methods presented in the paper are used to study a decomposition chemical reaction characterized by a scalar, one-dimensional reaction-diffusion equation, the propagation of a one-dimensional, confined, laminar flame in Cartesian co-ordinates and the Dwyer-Sanders model of one-dimensional flame propagation. It is shown that the adaptive moving method presented in this paper requires a smaller number of grid points than adaptive static, adaptive hybrid and non-adaptive methods. The adaptive hybrid method requires a smaller time step than adaptive static techniques, due to the lag between the grid prediction and the solution of the dependent variables. Non-adaptive methods of lines may yield temperature oscillations in front of and behind the flame front if Crank-Nicolson techniques are used to evaluate the time derivatives. Fourth-order accurate methods of lines in space yield larger temperature oscillations than second-order accurate methods of lines, and the magnitude of these oscillations decreases as the time step is decreased. It is also shown that three-point, fourth-order accurate discretizations of the spatial derivatives require the same number of grid points as second-order accurate, finite volume methods, in order to resolve accurately the structure of steep, fast-moving flame fronts.     

On some accurate finite-difference methods for laminar flame calculations

The ozone-decomposition flame has been studied by means of fourth- and second-order accurate schemes. The fourth-order methods include a method of lines, a time-linearization algorithm, and a majorant operatorsplitting technique. The second-order schemes include two time-linearization methods which use different temporal approximations. It is shown that the fourth-order techniques yield comparable results to those obtained with very accurate finite element and adaptive grid finite-difference algorithms. The results of the second-order methods are in good agreement with second-order explicit predictor-corrector methods but predict a lower flame speed than that obtained by means of fourth-order techniques. It is also shown that the temporal approximations are not as important as the spatial approximations in flame propagation problems characterized by the presence of several small time scales.     

A compact finite difference method for solving Burgers' equation

In this paper, a high‐order accurate compact finite difference method using the Hopf–Cole transformation is introduced for solving one‐dimensional Burgers' equation numerically. The stability and convergence analyses for the proposed method are given, and this method is shown to be unconditionally stable. To demonstrate efficiency, numerical results obtained by the proposed scheme are compared with the exact solutions and the results obtained by some other methods. The proposed method is second‐ and fourth‐order accurate in time and space, respectively. Copyright © 2009 John Wiley & Sons, Ltd.     

Parameter variation methods for cell mapping

In the study of nonlinear dynamic systems, the influence of system parameters on the long term behaviour plays an important role. In this paper, parameter variation methods are presented which can be used when investigating a nonlinear dynamic system by means of simple or interpolated cell mapping. In the case of coexisting attractors, the proposed methods determine the evolution of the basin boundaries when a system parameter is varied. Application of the methods to a modified Duffing equation is performed. It is concluded that the proposed methods are very efficient and accurate.     

三维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解三维非定常变系数对流扩散方程的一种高精度全隐紧致差分格式,该格式在空间上具有四阶精度,时间具有二阶精度。为了克服传统迭代法在每一个时间步上迭代求解隐格式时收敛速度慢的缺点,采用多重网格加速技术,建立了适用于本文高精度全隐紧致格式的多重网格算法,从而大大加快了迭代收敛速度。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

An upwind compact mpact approach with group velocity control for compressible flow fields

In this paper we construct an upwind compact finite difference scheme with group velocity control for better simulation of compressible flow fields. Compared with traditional difference schemes, compact schemes have higher accuracy for the same stencil width. By means of the characteristic analysis of the operators, the group velocity of wave packets will be controlled to suppress the non‐physical oscillations in numerical solutions. In numerical simulation of the 3D compressible flow fields the third‐order accurate upwind compact operator is used to approximate the derivatives in the convection terms of the compressible N–S equations, the traditional finite difference scheme is used to approximate the viscous terms. Numerical solutions indicate that the method is satisfactory. Copyright © 2004 John Wiley & Sons, Ltd.     

KdV方程孤立子及其相互作用的Petrov-Galerkin有限元数值模拟

本文给出一种Petrov-Galerkin有限元方法,并用以求Kdv方程各种初值问题的数值解,包括孤立波进波解,多个孤立于的相互作用,孤立子与振荡尾波等,所得结果与分析解及其它数值结果作了比较,表明本方法精度高、稳定性好,几乎没有高频伪振荡,计算程序简洁、明瞭,经济实用。     

Governing equations of envelopes created by nearly bichromatic waves on deep water

In this paper, the author derives the modified Schrödinger equation that governs the envelope created by nearly bichromatic waves, which are defined by the waves whose energy is almost concentrated in two closely approached wavenumbers. The stability of the solution of the modified Schrödinger equation for nearly bichromatic waves on deep water is discussed and the fact that the Benjamin–Feir instability occurs in a condition is shown. Moreover, the solutions of the modified Schrödinger equation for nearly bichromatic waves on deep water are obtained and, in a special case, the solution becomes the standing wave solution is shown.     

A non‐diffusive,divergence‐free,finite volume‐based double projection method on non‐staggered grids

Second‐order accurate projection methods for simulating time‐dependent incompressible flows on cell‐centred grids substantially belong to the class either of exact or approximate projections. In the exact method, the continuity constraint can be satisfied to machine‐accuracy but the divergence and Laplacian operators show a four‐dimension nullspace therefore spurious oscillating solutions can be introduced. In the approximate method, the continuity constraint is relaxed, the continuity equation being satisfied up to the magnitude of the local truncation error, but the compact Laplacian operator has only the constant mode. An original formulation for allowing the discrete continuity equation to be satisfied to machine‐accuracy, while using a finite volume based projection method, is illustrated. The procedure exploits the Helmholtz–Hodge decomposition theorem for deriving an additional velocity field that enforces the discrete continuity without altering the vorticity field. This is accomplished by solving a second elliptic field for a scalar field obtained by prescribing that its additional discrete gradients ensure discrete continuity based on the previously adopted linear interpolation of the velocity. The resulting numerical scheme is applied to several flow problems and is proved to be accurate, stable and efficient. This paper has to be considered as the companion of: 'F. M. Denaro, A 3D second‐order accurate projection‐based finite volume code on non‐staggered, non‐uniform structured grids with continuity preserving properties: application to buoyancy‐driven flows. IJNMF 2006; 52 (4):393–432. Now, we illustrate the details and the rigorous theoretical framework. Copyright © 2006 John Wiley & Sons, Ltd.     

On the local stability analysis of the approximate harmonic balance solutions

Periodic response of nonlinear oscillators is usually determined by approximate methods. In the "steady state" type methods, first an approximate solution for the steady state periodic response is determined, and then the local stability of this solution is determined by analyzing the equation of motion linearized about this predicted "solution". An exact stability analysis of this linear variational equation can provide erroneous stability type information about the approximate solutions. It is shown that a consistent stability type information about these solutions can be obtained only when the linearized variational equation is analyzed by approximate methods, and the level of accuracy of this analysis is consistent with that of the approximate solutions. It is demonstrated that these consistent stability results do not imply that the approximate solution is qualitatively correct. It is also shown that the difference between an approximate and the next higher order stability analysis can be used to "guess" the role of higher harmonics in the periodic response. This trial and error procedure can be used to ensure the qualitatively correct and numerically accurate nature of the approximate solutions and the corresponding stability analysis.     

Effect of Temperature on Slip Coefficients of Molecular Gases

Results obtained by accurate analytical methods applied to the problem of molecular-gas slip over a rigid spherical surface are reported. The Boltzmann equation is modified to take into account rotational degrees of freedom in the BGK model is used as a master kinetic equation. The calculated slip coefficients are shown to depend on the Prandtl number and on the gas temperature. Slip coefficients for several molecular gases are plotted as functions of temperature. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 58–65, January–February, 2006.     

Symmetry analysis of modified 2D Burgers vortex equation for unsteady case

In this paper, a symmetry analysis of the modified 2D Burgers vortex equation with a flow parameter is presented. A general form of classical and non-classical symmetries of the equation is derived. These are fundamental tools for obtaining exact solutions to the equation. In several physical cases of the parameter, the specific classical and non-classical symmetries of the equation are then obtained. In addition to rediscovering the existing solutions given by different methods, some new exact solutions are obtained with the symmetry method, showing that the symmetry method is powerful and more general for solving partial differential equations(PDEs).     

New analytic solutions for static problems of rectangular thin plates point-supported at three corners

This paper deals with the bending of rectangular thin plates point-supported at three corners using an analytic symplectic superposition method. The problems are of fundamental importance in both civil and mechanical engineering, but there were no accurate analytic solutions reported in the literature. This is attributed to the difficulty in seeking the solutions that satisfy the governing fourth-order partial differential equation with the free boundary conditions at all the edges as well as the support conditions at the corners. In the following, the Hamiltonian system-based equation for plate bending is formulated, and two types of fundamental problems are analytically solved by the symplectic method. The analytic solutions of the plates point-supported at three corners are then obtained by superposition, where the constants are obtained by a set of linear equations. The solution procedure presented in this paper offers a rigorous way to yield analytic solutions of similar problems. Some numerical results, validated by the finite element method, are shown to provide useful benchmarks for comparison and validation of other solution methods.     

Development and application of an adaptive finite element method to reaction-diffusion equations

An adaptive finite element method is developed and applied to study the ozone decomposition laminar flame. The method uses a semidiscrete, linear Galerkin approximation in which the size of the elements is controlled by an integral which minimizes the changes in mesh spacing. The sizes and locations of the elements are controlled by the location and magnitude of the largest temperature gradient. The numerical results obtained with this adaptive finite element method are compared with those obtained using fixed-node finite-difference schemes and an adaptive finite-difference method. It is shown that the adaptive finite element method developed here using 36 elements can yield as accurate flame speeds as fourth-order accurate, fixed-node, finite-difference methods when 272 collocation points are employed in the calculations.     

Forced dissipative Boussinesq equation for solitary waves excited by unstable topography

In the paper, the effects of topographic forcing and dissipation on solitary Rossby waves are studied. Special attention is given to solitary Rossby waves excited by unstable topography. Based on the perturbation analysis, it is shown that the nonlinear evolution equation for the wave amplitude satisfies a forced dissipative Boussinesq equation. By using the modified Jacobi elliptic function expansion method and the pseudo-spectral method, the solutions of homogeneous and inhomogeneous dissipative Boussinesq equation are obtained, respectively. With the help of these solutions, the evolutional character of Rossby waves under the influence of dissipation and unstable topography is discussed.     

Development of level set method with good area preservation to predict interface in two‐phase flows

A two‐step conservative level set method is proposed in this study to simulate the gas/water two‐phase flow. For the sake of accuracy, the spatial derivative terms in the equations of motion for an incompressible fluid flow are approximated by the coupled compact scheme. For accurately predicting the modified level set function, the dispersion‐relation‐preserving advection scheme is developed to preserve the theoretical dispersion relation for the first‐order derivative terms shown in the pure advection equation cast in conservative form. For the purpose of retaining its long‐time accurate Casimir functionals and Hamiltonian in the transport equation for the level set function, the time derivative term is discretized by the sixth‐order accurate symplectic Runge–Kutta scheme. To resolve contact discontinuity oscillations near interface, nonlinear compression flux term and artificial damping term are properly added to the second‐step equation of the modified level set method. For the verification of the proposed dispersion‐relation‐preserving scheme applied in non‐staggered grids for solving the incompressible flow equations, three benchmark problems have been chosen in this study. The conservative level set method with area‐preserving property proposed for capturing the interface in incompressible fluid flows is also verified by solving the dam‐break, Rayleigh–Taylor instability, bubble rising in water, and droplet falling in water problems. Good agreements with the referenced solutions are demonstrated in all the investigated problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Oscillations in high‐order finite difference solutions of stiff problems on non‐uniform grids

This work investigates the mitigation and elimination of scheme‐related oscillations generated in compact and classical fourth‐order finite difference solutions of stiff problems, represented here by the Burgers and Reynolds equations. The regions where severe gradients are anticipated are refined by the use of subdomains where the grid is distributed according to a geometric progression. It is observed that, for multi‐domain solutions, both the classical and compact fourth‐order finite difference schemes can exhibit spurious oscillations. When present, the oscillations are initially generated around the interface between the uniform and non‐uniform grid subdomains. Based on a thorough study of the grid distribution effects, it is shown that the numerical oscillations are caused by inadequate geometric progression ratios within the non‐uniformly discretized subdomains. Indeed, accurate solutions are obtainable if and only if the grid ratios in the non‐uniform subdomains are greater than a critical threshold ratio. It is concluded that high‐order classical and compact schemes can be used with confidence to efficiently solve one‐ or two‐dimensional problems whose solutions exhibit sharp gradients in very thin regions, provided that the numerically generated oscillations are eliminated by an appropriate choice of grid distribution within the non‐uniformly discretized subdomains. Copyright © 1999 John Wiley & Sons, Ltd.     

Two‐objective optimization strategies using the adjoint method and game theory in inverse natural convection problems

This paper considers the problem of estimating the strengths of two time‐varying heat sources simultaneously, from measurements of the temperature inside the square domain in a porous medium, when prior knowledge of the source functions is not available. This problem is an inverse natural convection problem. In order to circumvent this problem, we define several optimization criteria (objective functionals) that measure discrepancies between model and measured data, where objective functionals depend on two heat sources and use multi‐criteria optimization to identify Nash equilibria, which are solutions to the non‐cooperative game according to game theory. Two non‐cooperative game strategies are considered: competitive (Nash) game and hierarchical (modified Stackelberg) game. The methodology that we employ relies on a combination of mixed finite element space approximations, finite difference time discretizations, adjoint equation and sensitivity equation techniques, and nonlinear conjugate gradient algorithms for the solutions of estimating two heat sources. Applying the Sobolev gradient for the noise removal is investigated. The performance of the present technique of inverse analysis is evaluated, by means of several numerical experiments, and is found to be very accurate as well as efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite element methods for one-dimensional combustion problems

Three adaptive finite element methods based on equidistribution, elliptic grid generation and hybrid techniques are used to study a system of reaction–diffusion equations. It is shown that these techniques must employ sub-equidistributing meshes in order to avoid ill-conditioned matrices and ensure the convergence of the Newton method. It is also shown that elliptic grid generation methods require much longer computer times than hybrid and static rezoning procedures. The paper also includes characteristic, Petrov–Galerkin and flux-corrected transport algorithms which are used to study a linear convection–reaction–diffusion equation that has an analytical solution. The flux-corrected transport technique yields monotonic solutions in good agreement with the analytical solution, whereas the Petrov–Galerkin method with quadratic upstream-weighted functions results in very diffused temperature profiles. The characteristic finite element method which uses a Lagrangian–Eulerian formulation overpredicts the flame front location and exhibits overshoots and undershoots near the temperature discontinuity. These overshoots and undershoots are due to the interpolation of the results of the Lagrangian operator onto the fixed Eulerian grid used to solve the reaction–diffusion operator, and indicate that characteristic finite element methods are not able to eliminate numerical diffusion entirely.