Modified method of characteristics for solving population balance equations

This paper presents a new numerical method for solving the population balance equation using the modified method of characteristics. Aggregation and break‐up are neglected but the density function variations in the three‐dimensional space and its dependence on the external fields are accounted for. The method is an interpretation of the Lagrangian approach. Based on a pre‐specified grid, it follows the particles backward in time as opposed to forward in the case of traditional method of characteristics. Unlike the direct marching method, the inverse marching method uses a fixed grid thus, making it compatible with other numerical schemes (e.g. finite‐volume, finite elements) that may be used to solve other coupled equations such as the mass, momentum, and energy conservation equations. The numerical solutions are compared with the exact analytical solutions for simple one‐dimensional flow cases. Very good agreement between the numerical and the theoretical solutions has been obtained confirming the validity of the numerical procedure and the associated computer program. Copyright © 2003 John Wiley & Sons, Ltd.     

A local domain‐free discretization method to simulate three‐dimensional compressible inviscid flows

In this paper, the domain‐free discretization method (DFD) is extended to simulate the three‐dimensional compressible inviscid flows governed by Euler equations. The discretization strategy of DFD is that the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior‐dependent points are updated at each time step by extrapolation along the wall normal direction in conjunction with the wall boundary conditions and the simplified momentum equation in the vicinity of the wall. Spatial discretization is achieved with the help of the finite element Galerkin approximation. The concept of ‘osculating plane’ is adopted, with which the local DFD can be easily implemented for the three‐dimensional case. Geometry‐adaptive tetrahedral mesh is employed for three‐dimensional calculations. Finally, we validate the DFD method for three‐dimensional compressible inviscid flow simulations by computing transonic flows over the ONERA M6 wing. Comparison with the reference experimental data and numerical results on boundary‐conforming grid was displayed and the results show that the present DFD results compare very well with the reference data. Copyright © 2009 John Wiley & Sons, Ltd.     

An unstructured finite‐volume method for coupled models of suspended sediment and bed load transport in shallow‐water flows

The aim of this work is to develop a well‐balanced finite‐volume method for the accurate numerical solution of the equations governing suspended sediment and bed load transport in two‐dimensional shallow‐water flows. The modelling system consists of three coupled model components: (i) the shallow‐water equations for the hydrodynamical model; (ii) a transport equation for the dispersion of suspended sediments; and (iii) an Exner equation for the morphodynamics. These coupled models form a hyperbolic system of conservation laws with source terms. The proposed finite‐volume method consists of a predictor stage for the discretization of gradient terms and a corrector stage for the treatment of source terms. The gradient fluxes are discretized using a modified Roe's scheme using the sign of the Jacobian matrix in the coupled system. A well‐balanced discretization is used for the treatment of source terms. In this paper, we also employ an adaptive procedure in the finite‐volume method by monitoring the concentration of suspended sediments in the computational domain during its transport process. The method uses unstructured meshes and incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep sediment concentrations and bed load gradients that may form in the approximate solutions. Details are given on the implementation of the method, and numerical results are presented for two idealized test cases, which demonstrate the accuracy and robustness of the method and its applicability in predicting dam‐break flows over erodible sediment beds. The method is also applied to a sediment transport problem in the Nador lagoon.Copyright © 2013 John Wiley & Sons, Ltd.     

A fast and robust fictitious domain method for modelling viscous flows in complex mixers: The example of propellant make‐down

This work deals with the development of a fast three‐dimensional numerical strategy for the simulation of viscous fluid flow in complex mixing systems. The proposed method is based on a distributed Lagrange multiplier fictitious domain method and the use of the low‐cost MINI finite element. Contrary to the previous fictitious domain method developed by our group a few years ago, the underlying partial differential equations are solved here in a coupled manner using a consistent penalty technique. The method is discussed in detail and its precision is assessed by means of experimental data in the case of an agitated vessel. A comparison made with our existing fictitious domain method and its decoupled Uzawa‐based solver clearly shows the advantages of resorting to the MINI finite element and fully coupled solution strategy. The new technique is then applied to the simulation of the flow of a Newtonian viscous fluid in a three‐blade planetary mixer in the context of the production of solid propellants. Copyright © 2008 John Wiley & Sons, Ltd.     

Velocity–pressure coupling in finite difference formulations for the Navier–Stokes equations

A new numerical procedure for solving the two‐dimensional, steady, incompressible, viscous flow equations on a staggered Cartesian grid is presented in this paper. The proposed methodology is finite difference based, but essentially takes advantage of the best features of two well‐established numerical formulations, the finite difference and finite volume methods. Some weaknesses of the finite difference approach are removed by exploiting the strengths of the finite volume method. In particular, the issue of velocity–pressure coupling is dealt with in the proposed finite difference formulation by developing a pressure correction equation using the SIMPLE approach commonly used in finite volume formulations. However, since this is purely a finite difference formulation, numerical approximation of fluxes is not required. Results presented in this paper are based on first‐ and second‐order upwind schemes for the convective terms. This new formulation is validated against experimental and other numerical data for well‐known benchmark problems, namely developing laminar flow in a straight duct, flow over a backward‐facing step, and lid‐driven cavity flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Application of a preconditioned high‐order accurate artificial compressibility‐based incompressible flow solver in wide range of Reynolds numbers

In the present study, the preconditioned incompressible Navier‐Stokes equations with the artificial compressibility method formulated in the generalized curvilinear coordinates are numerically solved by using a high‐order compact finite‐difference scheme for accurately and efficiently computing the incompressible flows in a wide range of Reynolds numbers. A fourth‐order compact finite‐difference scheme is utilized to accurately discretize the spatial derivative terms of the governing equations, and the time integration is carried out based on the dual time‐stepping method. The capability of the proposed solution methodology for the computations of the steady and unsteady incompressible viscous flows from very low to high Reynolds numbers is investigated through the simulation of different 2‐dimensional benchmark problems, and the results obtained are compared with the existing analytical, numerical, and experimental data. A sensitivity analysis is also performed to evaluate the effects of the size of the computational domain and other numerical parameters on the accuracy and performance of the solution algorithm. The present solution procedure is also extended to 3 dimensions and applied for computing the incompressible flow over a sphere. Indications are that the application of the preconditioning in the solution algorithm together with the high‐order discretization method in the generalized curvilinear coordinates provides an accurate and robust solution method for simulating the incompressible flows over practical geometries in a wide range of Reynolds numbers including the creeping flows.     

Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

The two‐dimensional streamline upwind scheme for the convection–reaction equation

This paper is concerned with the development of the finite element method in simulating scalar transport, governed by the convection–reaction (CR) equation. A feature of the proposed finite element model is its ability to provide nodally exact solutions in the one‐dimensional case. Details of the derivation of the upwind scheme on quadratic elements are given. Extension of the one‐dimensional nodally exact scheme to the two‐dimensional model equation involves the use of a streamline upwind operator. As the modified equations show in the four types of element, physically relevant discretization error terms are added to the flow direction and help stabilize the discrete system. The proposed method is referred to as the streamline upwind Petrov–Galerkin finite element model. This model has been validated against test problems that are amenable to analytical solutions. In addition to a fundamental study of the scheme, numerical results that demonstrate the validity of the method are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

A kernel gradient free (KGF) SPH method

The finite particle method (FPM) is a modified SPH method with high order accuracy while retaining the advantages of SPH in modeling problems with free surfaces, moving interfaces, and large deformations. In both SPH and FPM, kernel gradient is necessary in kernel and particle approximation of a field function and its derivatives. In this paper, a new FPM is presented, which only involves kernel function itself in kernel and particle approximation. The kernel gradient is not necessary in the whole computation, and this approach is thus referred to as a kernel gradient free (KGF) SPH method. This is helpful when a kernel function is not differentiable or the resultant kernel gradient is not sufficiently smooth, and thus it is more general in selecting a kernel function. Moreover, different from the original FPM with an asymmetric corrective matrix, in the new FPM, the resultant corrective matrix is symmetric, and this is advantageous in particle approximations. A series of numerical examples have been conducted to show the efficiencies of KGF‐SPH including one‐dimensional mathematical tests of polynomial functions with equal or variable smoothing length and two‐dimensional incompressible fluid flow of shear cavity. It is found that KGF‐SPH is comparable with FPM in accuracy and is flexible as SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

An unstructured finite volume technique for the 3D Poisson equation on arbitrary geometry using a σ‐coordinate system

We present a solver for a three‐dimensional Poisson equation issued from the Navier–Stokes equations applied to model rivers, estuaries, and coastal flows. The three‐dimensional physical domain is composed of an arbitrary domain in the horizontal direction and is bounded by an irregular free surface and bottom in the vertical direction. The equations are transformed vertically to the σ‐coordinate system to obtain an accurate representation of top and bottom topographies. The method is based on a second‐order finite volume technique on prisms consisting of triangular grids in the horizontal direction. The algorithm is accompanied by an analysis of different linear system solvers in order to achieve fast solutions. Numerical experiments are conducted to test the numerical accuracy and the computational efficiency of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.     

Steady/unsteady aerodynamic analysis of wings at subsonic,sonic and supersonic Mach numbers using a 3D panel method

This paper treats the kernel function of an integral equation that relates a known or prescribed upwash distribution to an unknown lift distribution for a finite wing. The pressure kernel functions of the singular integral equation are summarized for all speed range in the Laplace transform domain. The sonic kernel function has been reduced to a form, which can be conveniently evaluated as a finite limit from both the subsonic and supersonic sides when the Mach number tends to one. Several examples are solved including rectangular wings, swept wings, a supersonic transport wing and a harmonically oscillating wing. Present results are given with other numerical data, showing continuous results through the unit Mach number. Computed results are in good agreement with other numerical results. Copyright © 2003 John Wiley & Sons, Ltd.     

An improved KGF‐SPH with a novel discrete scheme of Laplacian operator for viscous incompressible fluid flows

The kernel gradient free (KGF) smoothed particle hydrodynamics (SPH) method is a modified finite particle method (FPM) which has higher order accuracy than the conventional SPH method. In KGF‐SPH, no kernel gradient is required in the whole computation, and this leads to good flexibility in the selection of smoothing functions and it is also associated with a symmetric corrective matrix. When modeling viscous incompressible flows with SPH, FPM or KGF‐SPH, it is usual to approximate the Laplacian term with nested approximation on velocity, and this may introduce numerical errors from the nested approximation, and also cause difficulties in dealing with boundary conditions. In this paper, an improved KGF‐SPH method is presented for modeling viscous, incompressible fluid flows with a novel discrete scheme of Laplacian operator. The improved KGF‐SPH method avoids nested approximation of first order derivatives, and keeps the good feature of ‘kernel gradient free’. The two‐dimensional incompressible fluid flow of shear cavity, both in Euler frame and Lagrangian frame, are simulated by SPH, FPM, the original KGF‐SPH and improved KGF‐SPH. The numerical results show that the improved KGF‐SPH with the novel discrete scheme of Laplacian operator are more accurate than SPH, and more stable than FPM and the original KGF‐SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

Finite element framework for describing dynamic wetting phenomena

The finite element simulation of dynamic wetting phenomena, requiring the computation of flow in a domain confined by intersecting a liquid–fluid free surface and a liquid–solid interface, with the three‐phase contact line moving across the solid, is considered. For this class of flows, different finite element method (FEM) implementations have been used in the literature, and in some cases, these produced apparently contradictory results. In the present paper, a robust framework for the FEM simulation of dynamic wetting flows is developed, which, by consistently adhering to the FEM methodology, leaves no room for ad hoc ‘optional’ variations in the numerical handling of these flows. The developed approach makes it possible to conduct a convergence study, assess the spatial resolution required to achieve a preset accuracy and provide the corresponding benchmark calculations. This analysis allows one to identify numerical artefacts, which had previously been interpreted as physical effects, and demonstrates that suppressing numerical errors using a ‘strong’ implementation of a boundary condition creates bigger and less detectable errors elsewhere in the computational domain. We provide practical recommendations on the spatial resolution required by a numerical scheme for a given set of non‐dimensional similarity parameters and give a user‐friendly step‐by‐step guide specifying the entire implementation, which allows the reader to easily reproduce all presented results including the benchmark calculations. It is also shown how the developed framework accommodates generalizations of the mathematical model accounting for additional physical effects, such as gradients in surface tensions. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling of aggregation kernels for fluidized beds using discrete particle model simulations

Aggregation is one of the many important processes in chemical and process engineering. Several researchers have attempted to understand this complex process in fluidized beds using the macro-model of population balance equations （PBEs）. The aggregation kernel is an effective parameter in PBEs, and is defined as the product of the aggregation efficiency and collision frequency functions. Attempts to derive this kernel have taken different approaches, including theoretical, experimental, and empirical techniques. The present paper calculates the aggregation kernel using micro-model computer simulations, i.e., a discrete particle model. We simulate the micro-model without aggregation for various initial conditions, and observe that the collision frequency function is in good agreement with the shear kernel. We then simulate the micro-model with aggregation and calculate the aggregation efficiency rate.     

Depth‐averaged two‐dimensional curvilinear explicit finite analytic model for open‐channel flows

A depth‐averaged two‐dimensional model has been developed in the curvilinear co‐ordinate system for free‐surface flow problems. The non‐linear convective terms of the momentum equations are discretized based on the explicit–finite–analytic method with second‐order accuracy in space and first‐order accuracy in time. The other terms of the momentum equations, as well as the mass conservation equation, are discretized by the finite difference method. The discretized governing equations are solved in turn, and iteration in each time step is adopted to guarantee the numerical convergence. The new model has been applied to various flow situations, even for the cases with the presence of sub‐critical and supercritical flows simultaneously or sequentially. Comparisons between the numerical results and the experimental data show that the proposed model is robust with satisfactory accuracy. Copyright © 2000 John Wiley & Sons, Ltd.     

Local Heaviside‐weighted LRPIM meshless method and its application to two‐dimensional potential flows

In this paper, the local radial point interpolation meshless method (LRPIM) is used for the analysis of two‐dimensional potential flows, based on a local‐weighted residual method with the Heaviside step function as the weighting function over a local subdomain. Trial functions are constructed using radial basis functions. The present method is a truly meshless method based only on a number of randomly located nodes. Integration over the subdomains requires only a simple integration cell to obtain the solution. No element matrix assembly is required and no special treatment is needed to impose the essential boundary conditions. The novelty of the paper is the use of a local Heaviside weight function in the LRPIM, which does not need local domain integration and integrations only on the boundary of the local domains are needed. Effects of the sizes of local subdomain and interpolation domain on the performance of the present method are investigated. The behavior of shape parameters of multiquadrics has been systematically studied. Two numerical tests in groundwater and fluid flows are presented and compared with closed‐form solutions and finite element method. The results show that the use of a local Heaviside weight function in the LRPIM is highly accurate and possesses no numerical difficulties. Copyright © 2008 John Wiley & Sons, Ltd.     

??????????????SPH?????????

运用von Neumann稳定性分析方法对采用3次B-样条核函数的一维标准SPH进行了稳 定性分析,在此基础上对扰动波长的影响进行了讨论,并用一维弹性杆的SPH 计算对分析结论进行了数值验证. 理论分析和数值计算表明, von Neumann分析所得的稳定性随扰动波长变化的规律与实际计算结果完全相符, SPH稳定性与扰 动波长密切相关. 对于所分析的SPH方法,在扰动波长最小处稳定性最差,而在其它扰 动波长下,SPH的计算则可以在拉伸状态下稳定.     

A particle finite element method applied to long wave run‐up

This paper presents a Lagrangian–Eulerian finite element formulation for solving fluid dynamics problems with moving boundaries and employs the method to long wave run‐up. The method is based on a set of Lagrangian particles which serve as moving nodes for the finite element mesh. Nodes at the moving shoreline are identified by the alpha shape concept which utilizes the distance from neighbouring nodes in different directions. An efficient triangulation technique is then used for the mesh generation at each time step. In order to validate the numerical method the code has been compared with analytical solutions and a preexisting finite difference model. The main focus of our investigation is to assess the numerical method through simulations of three‐dimensional dam break and long wave run‐up on curved beaches. Particularly the method is put to test for cases where different shoreline segments connect and produce a computational domain surrounding dry regions. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical investigation of premixed combustion in a porous burner with integrated heat exchanger

In this paper, we perform a numerical analysis of a two-dimensional axisymmetric problem arising in premixed combustion in a porous burner with integrated heat exchanger. The physical domain consists of two zones, porous and heat exchanger zones. Two dimensional Navier–Stokes equations, gas and solid energy equations, and chemical species transport equations are solved and heat release is described by a multistep kinetics mechanism. The solid matrix is modeled as a gray medium, and the finite volume method is used to solve the radiative transfer equation to calculate the local radiation source/sink in the solid phase energy equation. Special attention is given to model heat transfer between the hot gas and the heat exchanger tube. Thus, the corresponding terms are added to the energy equations of the flow and the solid matrix. Gas and solid temperature profiles and species mole fractions on the burner centerline, predicted 2D temperature fields, species concentrations and streamlines are presented. Calculated results for temperature profiles are compared to experimental data. It is shown that there is good agreement between the numerical solutions and the experimental data and it is concluded that the developed numerical program is an excellent tool to investigate combustion in porous burner.     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.