A Generalization of the Arnold-Wendland Lemma to a Modified Collocation Method for Boundary Integral Equations in ℝ3

We prove quasioptimal and optimal order estimates in various Sobolev norms for the approximation of linear strongly elliptic periodic pseudodifferential equations in two independent variables by a modified method of nodal collocation by odd degree polynomial splines. In the one-dimensional case, our method coincides with the method of nodal collocation when odd degree polynomial splines are employed for the trial functions. The convergence analysis is based on an equivalence which we establish between our method and a nonstandard Galerkin method for an operator closely related to the given operator. This equivalence is realized through a crucial intermediate result (which we now term the Arnold-Wendland lemma) to connect the solution of central finite difference equations and that of certain nonstandard Galerkin equations. The results of this paper are genuine two-dimensional generalizations of the results obtained by ARNOLD and WENDLAND in [2] for the one-dimensional equations.     

Modeling Sediment Deposition Patterns Arising From Suddenly Released Fixed-Volume Turbulent Suspensions

Models presented in several recent papers [1–3] dealing with particle transport by, and deposition from, bottom gravity currents produced by the sudden release of dilute, well‐mixed fixed‐volume suspensions have been relatively successful in duplicating the experimentally observed long‐time, distal, areal density of the deposit on a rigid horizontal bottom. These models, however, fail in their ability to capture the experimentally observed proximal pattern of the areal density with its pronounced dip in the region initially occupied by the well‐mixed suspension and its equally pronounced local maximum at roughly the one‐third point of the total reach of the deposit. The central feature of the models employed in [1–3] is that the particles are always assumed to be vertically well‐mixed by fluid turbulence and to settle out through the bottom viscous sublayer with the Stokes settling velocity for a fluid at rest with no re‐entrainment of particles from the floor of the tank. Because this process is assumed from the outset in the models of [1–3], the numerical simulations for a fixed‐volume release will not take into account the actual experimental conditions that prevail at the time of release of a well‐mixed fixed‐volume suspension. That is, owing to the vigorous stirring that produces the well‐mixed suspension, the release volume will initially possess greater turbulent energy than does an unstirred release volume, which may only acquire turbulent energy as a result of its motion after release through various instability mechanisms. The eddy motion in the imposed fluid turbulence reduces the particle settling rates from the values that would be observed in an unstirred release volume possessing zero initial turbulent energy. We here develop a model for particle bearing gravity flows initiated by the sudden release of a fixed‐volume suspension that takes into account the initial turbulent energy of mixing in the release volume by means of a modified settling velocity that, over a time scale characteristic of turbulent energy decay, approaches the full Stokes settling velocity. Thereafter, in the flow regime, we assume that the turbulence persists and, in accord with current understanding concerning the mechanics of dense underflows, that this turbulence is most intense in the wall region at the bottom of the flow and relatively coarse and on the verge of collapse (see [22]) at the top of the flow where the density contrast is compositionally maintained. We capture this behavior by specifying a “shape function” that is based upon experimental observations and provides for vertical structure in the volume fraction of particles present in the flow. The assumption of vertically well‐mixed particle suspensions employed in [1–5] corresponds to a constant shape function equal to unity. Combining these two refinements concerning the settling velocity and vertical structure of the volume fraction of particles into the conservation law for particles and coupling this with the fluid equations for a two‐layer system, we find that our results for areal density of deposits from sudden releases of fixed‐volume suspensions are in excellent qualitative agreement with the experimentally determined areal densities of deposit as reported in [1, 3, 6]. In particular, our model does what none of the other models do in that it captures and explains the proximal depression in the areal density of deposit.     

Prediction of concentration and particle size distribution in the flow of multi-sized particulate slurry through rectangular duct

The analytical model proposed by Karabelas [AICHe J. 23 (4) (1977) 426] has been modified and upgraded to predict the concentration profile and particle size distribution across the cross-section of rectangular duct for the flow of multi-sized particulate slurry. The predictions have been compared with the experimental data reported by Kaushal [Prediction of particle distribution in the flow of multi-sized particulate slurries through closed ducts and open channels, 1995]. The limitations of the original model have been identified and modifications incorporated to take into account the effect of solid concentration on settling rate and particle diffusivity. The predictions by the modified model are in good agreement with the experimental results.     

半无穷长圆管内的低雷诺数入口流

1969年Lew及Fung[1]计算了圆管内的低雷诸数入口流.1982年Dagan等人[2]得到了有限长圆柱形孔道内蠕动流的级数解.[1]中所得的数值解实质上代表有限长圆管内的低雷诺数入口流,因为一般解中的富氏积分已用富氏级数代替.本文直接计算富氏积分,更精确地求出了真正的半无穷长圆管内Stokes入口流的速度分布,压力分布以及流函数,与此对应的入口段长度为圆管半径的1.2倍,接近于Lew及Fung得到的结果1.3倍.此外,本文还研究了配置法的收敛性,证明了此法在入口流问题中具有很好的收敛性,因此可以在其他类似的问题中采用.     

一维Allen-Cahn方程紧差分格式的离散最大化原则和能量稳定性研究

本文主要研究相场模拟中的Allen-Cahn模型,考虑一维Allen-Cahn方程紧差分方法的数值逼近.建立具有O(∫τ2+h4)精度的全离散紧差分格式,证明在合理的步长比和时间步长的约束下,其数值解满足离散最大化原则,在此基础上,研究了全离散格式的能量稳定性.最后给出数值算例.     

A decoupled finite particle method for modeling incompressible flows with free surfaces

Smoothed particle hydrodynamics (SPH) is a meshfree Lagrangian particle method, and it has been applied to different areas in engineering and sciences. One concern of the conventional SPH is its low accuracy due to particle inconsistency, which hinders the further methodology development. The finite particle method (FPM) restores the particle consistency in the conventional SPH and thus significantly improves the computational accuracy. However, as pointwise corrective matrix inversion is necessary, FPM may encounter instability problems for highly disordered particle distribution. In this paper, through Taylor series analyses with integration approximation and assuming diagonal dominance of the resultant corrective matrix, a new meshfree particle approximation method, decoupled FPM (DFPM), is developed. DFPM is a corrective SPH method, and is flexible, cost-effective and easy in coding with better computational accuracy. It is very attractive for modeling problems with extremely disordered particle distribution as no matrix inversion is required. One- and two-dimensional numerical tests with different kernel functions, smoothing lengths and particle distributions are conducted. It is demonstrated that DFPM has much better accuracy than conventional SPH, while particle distribution and the selection of smoothing function and smoothing length have little influence on DFPM simulation results. DFPM is further applied to model incompressible flows including Poiseuille flow, Couette flow, shear cavity and liquid sloshing. It is shown that DFPM is as accurate as FPM while as flexible as SPH, and it is very attractive in modeling incompressible flows with possible free surfaces.     

ON THE CENTRAL RELAXING SCHEME Ⅱ: SYSTEMS OF HYPERBOLIC CONSERVATION LAWS

1. IntroductionWe are interested in construction of the central reltalng sChemes for system of noIilinearhyperbolic conservation lawswith initial data U(0, x) = Uo(x), x = (x1 ? ...! xd), based on the local relaJxation approkimationof Eq.(1.1) [2, 3, 6, 8, 9, 12].To i11ustrate the basic idea of the relaalng schemes, for the sake of simplicity in the presentation, we restrict our attention to onedimensional scalar conservaioll lawsFirst, introduce a linear hyperbollc system with a stiff sourc…     

不同大小粒子之间相互作用的直接数值模拟

应用改进的拉格朗日乘子/虚拟区域算法对不同大小的两个圆形粒子在二维方槽中的沉降过程和相互作用进行了直接数值模拟,并进行了实验验证．结果表明不同大小的两个粒子在沉降过程中的相互作用可以描述为追尾、接触、旋转和分离4个过程,只有当两个粒子尺度差异很小时,才会重复进行DKT过程．在两个粒子相互作用的过程中,小粒子的运动受到大粒子的影响更剧烈一些,而相反大粒子运动包括运动轨迹和速度所受到的影响则相对较小．     

两平行圆板间径向层流进口段效应分析

本文首先将B.B.方法[1]推广至二平行圆板间的径向扩散流动,由边界层运动方程式同时推导出动量积分方程式和能量积分方程式,而后再用Picard逐次逼近法[2]解能量积分方程式,求得进口段通道长随边界层厚度而改变的二级近似显函数表达式.从而为进口段效应诸系数的直接解析分析提供了可能.特别是当圆板外径小于进口段长度时,更加突出地表现了本方法的优越性.由于采用了能量积分方程式,则压力损失系数的各项才得以从理论上独立地推导出来.本文所提供的压力损失系数计算值,在进口修正雷诺数Re<100时,和文献[3]比较与实验值更为接近.因此在该范围内本文的结果既可靠又简便.     

两相流中柱状固粒对流体湍动特性影响的研究

对含柱状固粒的两相流场,建立了包含柱状固粒对流场影响的流体脉动速度方程,在求解脉动速度方程的基础上,经平均得到流体的湍流强度和雷诺应力.将该方法用于槽流湍流场的求解,并与单相流实验结果进行了比较.计算中变化柱状固粒的参数,给出了固粒的体积分数、长径比、松驰时间对流场湍动特性的影响,说明粒子对流场的湍动特性起着抑制作用,其抑制的程度与粒子的体积分数、长径比成正比,与粒子的松弛时间成反比.     

Numerical modeling on inlet aperture effects on flow pattern in primary settling tanks

Inlets should be designed to dissipate the kinetic energy or velocity head of the mixed liquor and to prevent short-circuiting, mitigate the effects of density currents, and minimize blanket disturbances. Flow in primary settling tank is simulated by means of computational fluid dynamics. The fluid is assumed incompressible and non-buoyant. A two-dimensional computational and one phase fluid dynamics model was built to simulate the flow properties in the settling tank including the velocity profiles, the flow separation area and kinetic energy. In this study, the RNG turbulent model was solved with the Navier–Stokes equations. In order to evaluate hydraulic influences on the velocity profile, separation length and kinetic energy, three different of opening positions and two and three aperture in inlets were simulated. The flow model uses to apply a fixed-grid of cells that are all rectangular faces; the fluid moves through the grid and free surfaces are tracked with the volume-of-fluid (VOF) technique. Effects of numbers and locations of inlet apertures on the flow field are presented and the results show the positions of inlet apertures are affected on the flow pattern in the settling basin and increasing the numbers of slots can reduce kinetic energy in the inlet zone and produce uniform flow.     

二相平面渗流问题的不变流管解法及其理论分析

本文总结和改进了工程上广泛应用的求解二相平面渗流问题的不变流管近似方法.对其核心部份即给定压差时一维变截面流管中的二相驱替问题作了全面的考察,证明了解的存在唯一性,给出了准确解、数值解及其收敛性和稳定性分析.     

Approximate analytical solution for supercritical flow in rectangular curved channels

In this study, we present an asymptotical mathematical model and an analytical solution for a supercritical flow in curved rectangular open channels. An original approach is proposed for solving the free-surface configuration and features of the flow in the presence of cross shock waves. The two-dimensional steady depth-averaged shallow water equations are transformed into an equivalent one-dimensional (1D) unsteady flow problem and a first order approximation is then obtained using small perturbation theory. Furthermore, the 1D asymptotic model is solved analytically by Laplace integral transformation and the two-dimensional flow field solution is reconstructed according to the translating planes. The free-surface profile along the outer chute wall and downstream channel was compared with the available experimental data, and the results indicated the satisfactory agreement of the maximum flow depth, peak positions, and wavelength. The proposed approach provides accurate predictions of the flow features and it facilitates the safe design of curved channel transitions.     

The Jin–Xin relaxation approximation of scalar conservation laws in several dimensions with large initial perturbation

This paper is concerned with nonlinear stability of strong planar rarefaction waves for the Jin–Xin relaxation approximation of scalar conservation laws in several dimensions. For such a problem, local stability of weak or strong planar rarefaction waves have been obtained in Luo (1997) [20] and Zhao (2000) [43] respectively. For the global stability results, to the best of our knowledge, the only result available now is on the one-dimensional case, cf. Zhao (2000) [43], which is based on the maximum principle established in Natalini (1996) [30]. The main purpose of this paper is try to deduce some nonlinear stability results with large initial perturbation. Our analysis is based on the elementary energy method and the continuation argument.     

Nonstationary mass transfer of a reacting spherical particle in laminar stream of a viscous fluid : PMM vol. 38, n 6, 1974, pp. 1041–1047

The process of the formation of a stationary mass transfer mode for a moving reacting particle is examined. An analytic expression valid for a nonstationary distribution of the concentration of matter in a steady stream of viscous fluid, flowing past a spherical particle, was obtained for the case when at a certain instant a chemical reaction of the first order begins at the surface of the sphere. The problem is solved for small finite Reynolds and Péclet numbers. The solution of the corresponding stationary problem has been obtained in [1]. Paper [2] examined a nonstationary heat transfer of a fluid spherical drop in an inviscid flow with spasmodic change of initial temperature at high Péclet numbers. Paper [3] contains an analysis of the problem of a nonstationary heat transfer of a rigid spherical particle for small Reynolds and Péclet numbers at spasmodic change of temperature of the particle surface. The results obtained in [3] can be used to describe the mass transfer for a moving reacting particle only in the case of a diffusion mode of the chemical reaction.     

Cryptography using multiple one-dimensional chaotic maps

Recently, Pareek et al. [Phys. Lett. A 309 (2003) 75] have developed a symmetric key block cipher algorithm using a one-dimensional chaotic map. In this paper, we propose a symmetric key block cipher algorithm in which multiple one-dimensional chaotic maps are used instead of a one-dimensional chaotic map. However, we also use an external secret key of variable length (maximum 128-bits) as used by Pareek et al. In the present cryptosystem, plaintext is divided into groups of variable length (i.e. number of blocks in each group is different) and these are encrypted sequentially by using randomly chosen chaotic map from a set of chaotic maps. For block-by-block encryption of variable length group, number of iterations and initial condition for the chaotic maps depend on the randomly chosen session key and encryption of previous block of plaintext, respectively. The whole process of encryption/decryption is governed by two dynamic tables, which are updated time to time during the encryption/decryption process. Simulation results show that the proposed cryptosystem requires less time to encrypt the plaintext as compared to the existing chaotic cryptosystems and further produces the ciphertext having flat distribution of same size as the plaintext.     

Particle transport in a bottom-feed separation vessel

A two-dimensional, axisymmetric numerical model of particle separation in a bottom-feed separation vessel is presented. The model includes six separate particle classes and assumes that the settling velocity of each particle class is sufficiently small when compared to the high inflow turbulence levels that the effect of the particles on turbulence can be neglected. Low particle settling velocities coupled with low particle volume fractions allows application of a drift-flux multi-phase model. The comparison between numerical results and measured plant data is in good agreement for overflow of all particle classes. Results of simulations show that bottom feeding results in a negatively buoyant, particle-laden jet being formed in the core of the vessel. The fraction of large particles that is carried out through the overflow is found to be critically dependent on the inlet velocity. The most effective way to reduce carry-over of large particles at the same time as maintaining through-put is to increase the diameter of the inlet feed pipe.     

Simulation of suspension gravity currents with different initial aspect ratio and layout of turbidity fence

A three-dimensional (3D) numerical model, using large eddy simulation (LES), is developed for simulating the motion of suspension gravity currents. The suitable values of model parameters are determined using the existing experimental data of a two-dimensional (2D) suspension (a mixture composed of water and glass bead particles) cloud. The simulated gravity current with different initial aspect ratio (length/breadth) of the suspension is compared with the reported data of 3D laboratory experiments to investigate the effect of initial aspect ratio on the flow characteristics and the diffusion of turbidity under the presence of a turbidity fence. The comparison of simulated results of such main flow characteristics as front height, front propagation velocity and particle deposition with the experimental data reveals that the model is capable of simulating the complex behavior of the 3D suspension gravity currents to a reasonably good accuracy under complex conditions.     

One-dimensional viscoelastic fluid model where viscosity and normal stress coefficients depend on the shear rate

We study the unsteady motion of a viscoelastic fluid modeled by a second-order fluid where normal stress coefficients and viscosity depend on the shear rate by using a power-law model. To study this problem, we use the one-dimensional nine-director Cosserat theory approach which reduces the exact three-dimensional equations to a system depending only on time and on a single spatial variable. Integrating the equation of conservation of linear momentum over the tube cross-section, with the velocity field approximated by the Cosserat theory, we obtain a one-dimensional system. The velocity field approximation satisfies both the incompressibility condition and the kinematic boundary condition exactly. From this one-dimensional system we obtain the relationship between average pressure and volume flow rate over a finite section of the tube with constant and variable radius. Also, we obtain the correspondent equation for the wall shear stress which enters directly in the formulation as a dependent variable. Attention is focused on some numerical simulation of unsteady/steady flows for average pressure, wall shear stress and on the analysis of perturbed flows.     

带激波的一维非定常两相流动的数值模拟

本文将处理带激波的单相气体非定常流动问题的两种高分辨数值方法(随机取样法和二阶GRP有限差分法)推广应用于气固悬浮体(亦称含灰气体)两相情况,计算了含灰气体激波管中两相激波特性、波后流场结构及气固两相流动参数随时间的变化.数值结果表明:这两种方法均能给出带有尖锐间断阵面的两相激波松弛结构.二阶GRP方法在计算精度和机时耗用等方面优于随机取样法.