Lattice hydrodynamic model for traffic flow on curved road with passing

In order to investigate the effect of passing upon traffic flow on curved road, in this paper, an extended one-dimensional lattice hydrodynamic model for traffic flow on curved road with passing is proposed. The stability condition is obtained by the use of linear stability analysis. The result of stability analysis shows that passing behavior plays an important role in influencing the stability of traffic flow as well as radian of curved road. The nonlinear wave equations including Burgers, Korteweg-de Vries and modified Korteweg-de Vries equations are derived to describe the nonlinear traffic behavior in different regions, respectively. The analytical results show that reducing the coefficient of passing may enhance the stability of traffic flow. Jamming transition occurs between uniform flow and kink jam when the coefficient of passing is less than the critical value. When the coefficient of passing is larger than the critical value, jamming transition occurs from uniform flow to irregular wave through chaotic phase with decreasing sensitivity parameter. In addition, compared with other segments of curved road, traffic flow with passing easily becomes unstable and complicated at the entrance and exit of curved road, especially at the entrance of curved road. The numerical simulations are given to illustrate and clarify the analytical results.     

局部狭窄血管中血液振荡流的速度和切应力

本文建立一种分析局部缓慢狭窄血管中血液振荡流的数学模型，给出了血液的轴向流速，径向流速和切应力的包含压力梯度项的解析表达式，并讨论了血管内由局部狭窄引起的压力梯度沿轴向变化的规律。文章以局部余弦狭窄为例进行数值计算，详细讨论上游均匀管段压力梯度的定常部分和不同次谐波对狭窄管段内流速和切应力的影响。数值结果表明，与均匀管情况相比，在狭窄段内，血液振荡流轴向流速无论平均值还是脉动幅值均明显增大，且径向流速不再为零。但径向流速仍远小于轴向流速。同时，切应力也不再仅由轴向流速梯度提供，径向流速梯度也将产生切应力，但是在计算管壁切向上的切应力时，径向流速梯度的贡献仍相当大。与均匀管管壁切应力沿流运方向保持恒定不同。狭窄管管壁切应力（平均值和脉动值）将随着狭窄高度的增大而增大，在狭窄最大高度处达到最大，因而沿流动方向产生了较大的切应力梯度。     

Modeling Variable Density Flow and Solute Transport in Porous Medium: 2. Re‐Evaluation of the Salt Dome Flow Problem

Case 5, Level 1 of the international HYDROCOIN groundwater flow modeling project is an example of idealized flow over a salt dome. The groundwater flow is strongly coupled to solute transport since density variations in this example are large (20%).Several independent teams simulated this problem using different models. Results obtained by different codes can be contradictory. We develop a new numerical model based on the mixed hybrid finite elements approximation for flow, which provides a good approximation of the velocity, and the discontinuous finite elements approximation to solve the advection equation, which gives a good approximation of concentration even when the dispersion tensor is very small. We use the new numerical model to simulate the salt dome flow problem.In this paper we study the effect of molecular diffusion and we compare linear and nonlinear dispersion equations. We show the importance of the discretization of the boundary condition on the extent of recirculation and the final salt distribution. We study also the salt dome flow problem with a more realistic dispersion (very small dispersion tensor). Our results are different to prior works with regard to the magnitude of recirculation and the final concentration distribution. In all cases, we obtain recirculation in the lower part of the domain, even for only dispersive fluxes at the boundary. When the dispersion tensor becomes very small, the magnitude of recirculation is small. Swept forward displacement could be reproduced by using finite difference method to compute the dispersive fluxes instead of mixed hybrid finite elements.     

Multi-mode solitons in a long-short range traffic lattice model with time delay

We consider a new form of solutions of a special lattice model for traffic system. By analyzing nearest sites’ interactions, time delay, and bumpy effects, we deduce the bifurcation lines and surfaces for stable and unstable regions and show how they vary as parameters vary. It shows that keeping other conditions unchanged, as the incoming flow increases, the traffic flow becomes unstable, opposite to when outgoing flow increases, it becomes stable. Besides, considering delayed optimal flow, multiple sites effect or artificial parameters can also help stabilize the traffic road condition. Moreover, by putting it into the framework of mKdV equations, we obtain the kink–antikink solitons involving all parameters, which show the feature of the traffic congestion. The result is original, and our model in differential or difference form can be reduced into the previous ones by choosing appropriate parameters. Since the optimal velocity function we considered involves finitely or infinitely many sites, the density waves can be in multi-mode and high dimension forms and can also be quasi-periodic, we show a new feature of the traffic lattice system.

     

An eigenfunction solution of entry flow in a semi-infinite pipe at low Reynolds numbers

The developing flow in a semi-infinite pipe is solved for by an eigenfunction method. The results are applicable to zero and low Reynolds number flows. The eigenfunction method yields a solution which accurately predicts overshoots in the axial velocity for small axial positions. Results are presented for both uniform and non-uniform inlet conditions. The latter condition has important implications for numerical simulations of the developing flow problem.This research was performed while Mr Benson was a graduate student at Clarkson College, Potsdam, NY, USA.     

金属切屑塑性流动的稳定性

对金属正交切削过程中切屑形成机制和材料塑性流动行为进行实验研究和理论分析. 通过对4 种常用金属材料正交切削过程的实验研究和切屑形貌的微观观察,确定了连续切屑转变成锯齿切屑的临界速度. 结果表明该临界速度与材料性能相关. 在实验观察基础上,提出描述材料正交切削过程的二维分析模型. 该模型假设切屑形成区为包括主剪切区和次剪切区的一个平行四边形. 载荷有主剪切区中的剪应力和次剪切区中的正压力;通过量纲分析得到描述材料正交切削过程的无量纲主控参数和无量纲形式的基本控制方程;应用线性稳定性分析方法建立平面应变状态下评价材料塑性流动稳定性的普遍准则;求得切屑形成区内材料塑性变形的速度和应力近似解. 讨论切屑形成、形貌转变以及相关的塑性失稳机制. 分析结果表明, 表征材料惯性与阻尼之比的无量纲参数— 雷诺数可以作为主控参数描述金属切削过程以及切屑材料塑性流动的稳定性.      

THE INLET EFFECT ON THE DRAG FACTOR OF A SPHERE IN A TUBE

The creeping motion Ground a sphere situated axisymmetrically near the entrance of asemi-infinite circular cylindrical tube is analyzed using infinite series solutions for thevelocity components. pressure and the stream function. Truncating the infinite series. thecorresponding coefficients in the series are determined by a collocation technique. The dragfactor and the stress distribution on the surface of the sphere are calculated for the sphere inmotion in quiescent fluid and for the flow with uniform velocity at the entrance past a rigidlyheld sphere. The results indicate that a sphere near the entrance which has a uniformentrance velocity profile will suffer larger drag than that in infinite tube.Theconvergence of the collocation technique is tested by numerical calculation. It is shown thatthe technique has good convergence properties.     

Numerical study of single bubble motion in liquid metal exposed to a longitudinal magnetic field

The paper presents numerical simulations modeling the ascent of an argon bubble in liquid metal with and without an external magnetic field. The governing equations for the fluid and the electric potential are discretized on a uniform Cartesian grid and the bubble is represented with a highly efficient immersed boundary method. The simulations performed were conducted matching experiments under the same conditions so that sound validation is possible. The three-dimensional trajectory of the bubble is analyzed quantitatively and related to the flow structures in the wake. Indeed, the substantial impact of the magnetic field in the bubble trajectory results from its influence on the wake. Quantitative data describing the selective damping of vortex structures are provided and discussed. As a result of applying a longitudinal field, the time-averaged bubble rise velocity increases for large bubbles, it reaches a maximum and then decreases when further increasing the magnetic interaction parameter. For small bubbles, the time-averaged bubble rise velocity decreases when increasing the magnetic field. The bubble Strouhal number as a dimensionless frequency is reduced with the application of a magnetic field for all bubbles considered and the zig–zag trajectory of the bubble becomes more rectilinear.     

DENTING AND FAILURE OF LIQUID-FILLED TUBES UNDER LATERAL IMPACT

In this paper, numerical simulation of impact cases of liquid-filled tube impacted by missiles is conducted with a commercial finite element code LS-DYNA, and the results obtained are compared with the experimental data to verify the validity of the numerical simulation model adopted. With the verified numerical method, the processes of dynamic response of a blunt indenter impacting an empty or liquid-filled three-span continuous tubular beam are studied when the parameter such as the indenter’s mass, liquid’s density or impact velocity is varied and the other conditions are kept the same. The simulation results indicate that the critical perforation energy and the deformation of the wall of the pipe are significantly influenced by the presence of the liquid and the pressure. The liquid filling the tube provides a ‘foundation’ pressure to resist and localize the deformation, which may affect the perforation process and lead to a reduction of the ballistic limit. The simulation results also indicate that the increase of the fluid density filled in the tube will decrease the ballistic limit, but the fluid density must be in some scope. The relationship between the ballistic limit velocity of the tube and the mass of the impact missile is nonlinear in the Cartesian coordinate while it becomes linear through logarithmic transformation.     

The effect of inflow conditions on the transition to turbulence in large eddy simulations of spatially developing mixing layers

Large Eddy Simulations (LES) of spatially developing turbulent mixing layers have been performed for flows of uniform density and Reynolds numbers of up to 50,000 based on the visual thickness of the layer and the velocity difference across it. On a fine LES grid, a validation simulation performed with a hyperbolic tangent inflow profile produces flow statistics that compare extremely well with reference Direct Numerical Simulation (DNS) data. An inflow profile derived from laminar Blasius profiles produces a flow that is significantly different to the reference DNS, particularly with respect to the initial development of the flow. When compared with experimental data, however, it is the boundary layer-type inflow simulation produces the better prediction of the flow statistics, including the mean transition location. It is found that the boundary layer inflow condition is more unstable than the hyperbolic tangent inlet profile. A suitably designed coarse LES grid produces good predictions of the mean transition location with boundary layer inflow conditions at a low computational cost. The results suggest that hyperbolic tangent functions may produce unreliable DNS data when used as the initial condition for studies of the transition in the mixing layer flow.     

油田注水开发中水嘴节流与地层渗流耦合研究

配水器的选取会对油田注水效果产生明显影响。为精确地描述注水井注入过程中井底压力和流量的动态响应,本文在传统的内边界井处理模型的基础上建立了考虑水嘴压力损失影响的井处理模型。耦合求解井底压力、流量及地层流动,得到了注入过程中考虑水嘴压力损失影响的井底压力和流量响应情况。以此来分析不同水嘴、渗透率对井底压力、注入量随时间变...     

Turbulent travel speeds in nonlinear vehicle dynamics

Quarter car models of vehicles rolling on wavy roads lead to limit cycles of travel speed and acceleration with period doublings and bifurcation effects for appropriate driving force parameters. In case of narrow-banded road excitations, speed jumps occur, additionally. This has the consequence that the driving speed becomes turbulent. Bifurcation and jump effects vanish with growing vehicle damping. The same happens for increasing bandwidth of road excitations when, e.g., on flat highways there are no big road waves but only small noisy slope processes generated by rough road surfaces. The paper derives a new stability condition in mean. Numerical time integrations are stabilized by means of polar coordinates. Equivalently, Fourier series expansions are introduced in the angle domain. Phase portraits of travel speed and acceleration show new period-doublings of limit cycles when speed gets stuck before resonance. The paper extends these investigations to the stochastic case that road surfaces are random generated by filtered white noise. By means of Gaussian closure, a nonlinear mean speed equation is derived which includes the extreme cases of wavy roads and road noise.

     

Solute Transport Analysis Through Heterogeneous Media in Nonuniform in the Average Flow by a Stochastic Approach

The stochastic approach has been shown to be an excellent tool for the characterisation and analysis of velocity fields and transport processes through heterogeneous porous formations. The main results (linear theory) have been obtained for problems with simplified flow conditions, usually in the assumption of uniform in the average flow, but a great effort is spent to reach theoretical results for more complex situations.This paper deals with 2D heterogeneous aquifers subject to uniform recharge; the stochastic approach is adopted to characterise, as ensemble behaviour, the velocity field and transport processes of a nonreactive solute. The impact of transmissivity conditioning on solute particles trajectories is analysed and an application is carried out. The analytical formulations, obtained by a first order analysis, are compared to the one resulting from constant in the average hydraulic gradient, and their reliability is investigated with numerical tests performed by a Monte Carlo method.The result of this study is that strong non-stationarities are present in the flow and transport process. A detailed analysis shows that the theoretical results cannot be extended to cases with high heterogeneity level, unlike the uniform in the average flow fields.     

Finite volume solution of the Navier–Stokes equations in velocity–vorticity formulation

For the incompressible Navier–Stokes equations, vorticity‐based formulations have many attractive features over primitive‐variable velocity–pressure formulations. However, some features interfere with the use of the numerical methods based on the vorticity formulations, one of them being the lack of a boundary conditions on vorticity. In this paper, a novel approach is presented to solve the velocity–vorticity integro‐differential formulations. The general numerical method is based on standard finite volume scheme. The velocities needed at the vertexes of each control volume are calculated by a so‐called generalized Biot–Savart formula combined with a fast summation algorithm, which makes the velocity boundary conditions implicitly satisfied by maintaining the kinematic compatibility of the velocity and vorticity fields. The well‐known fractional step approaches are used to solve the vorticity transport equation. The paper describes in detail how we accurately impose no normal‐flow and no tangential‐flow boundary conditions. We impose a no‐flux boundary condition on solid objects by the introduction of a proper amount of vorticity at wall. The diffusion term in the transport equation is treated implicitly using a conservative finite update. The diffusive fluxes of vorticity into flow domain from solid boundaries are determined by an iterative process in order to satisfy the no tangential‐flow boundary condition. As application examples, the impulsively started flows through a flat plate and a circular cylinder are computed using the method. The present results are compared with the analytical solution and other numerical results and show good agreement. Copyright © 2005 John Wiley & Sons, Ltd.     

Fluid inertia in large amplitude oscillatory shear

Homogeneous shearing is required in sliding plate flow experiments with one plate fixed and the other oscillating. However, when fluid inertia becomes significant, the velocity gradient and the stress will not be uniform. MacDonald et al. (1969) and Schrag (1977) investigated this effect for a linear viscoelastic fluid. However, linear viscoelasticity does not describe the behavior of melts in large amplitude oscillatory shear (LAOS). Jeyaseelan et al. (1993) have shown that the Berkeley kinetic network model does accurately describe the LAOS behavior of polymer melts. In this work, the Berkeley model is solved for LAOS in sliding plate flow with fluid inertia, by numerical integration of spatially discretized forms of the governing equations. Nonlinear viscoelasticity is predicted to aggravate the effects of fluid inertia in LAOS and experiments confirm this. Specifically, fluid inertia amplifies the first harmonic and produces no even harmonics. Operating limits are presented graphically for minimizing inertial effects in LAOS experiments. Received: 2 January 1998 Accepted: 27 April 1998     

Influence of the inlet velocity profiles on the prediction of velocity distribution inside an electrostatic precipitator

The influence of the velocity profile at the inlet boundary on the simulation of air velocity distribution inside an electrostatic precipitator is presented in this study. Measurements and simulations were performed in a duct and an electrostatic precipitator (ESP). A four-hole cobra probe was used for the measurement of velocity distribution. The flow simulation was performed by using the computational fluid dynamics (CFD) code FLUENT. Numerical calculations for the air flow were carried out by solving the Reynolds-averaged Navier–Stokes equations coupled with the realizable k-ε turbulence model equations. Simulations were performed with two different velocity profiles at the inlet boundary – one with a uniform (ideal) velocity profile and the other with a non-uniform (real) velocity profile to demonstrate the effect of velocity inlet boundary condition on the flow simulation results inside an ESP. The real velocity profile was obtained from the velocity measured at different points of the inlet boundary whereas the ideal velocity profile was obtained by calculating the mean value of the measured data. Simulation with the real velocity profile at the inlet boundary was found to predict better the velocity distribution inside the ESP suggesting that an experimentally measured velocity profile could be used as velocity inlet boundary condition for an accurate numerical simulation of the ESP.     

Casson流体在旋转圆盘上的流动

本文研究Casson流体在旋转圆盘上的涂层流动特性,得到基本流动的速度分布,并且用Runge-Kutta法进行数值计算,得到薄膜厚度随时间和流态参数的变化规律,还用能量法检验了流动稳定性。     

Large eddy simulation of polydispersed inertial particles using two-way coupled PDF-PBE

This paper presents large eddy simulation (LES) of polydispersed turbulent recirculating flows using a two-way coupled probability density function of the population balance equation (PDF-PBE). A stochastic Monte Carlo method is adopted to solve the PDF-PBE on an ensemble of notional Lagrangian particles and the method of Stokes binning, that was recently developed by Salehi et al. (2017) is employed to explicitly treat effects of inertia. The PDF-PBE is applied to the experiment of Boŕee et al. (2001) which studied dispersion of polysized inertial particles in a bluff body configuration. The particle mass loading is 22% where the dispersed elements affect the carrier phase velocity. The simulations are performed using both the standard Smagorinsky and the Wall-Adapting Local Eddy-viscosity (WALE) subgird turbulent models. It is found that the subgrid model has a significant impact on the results, particularly on the carrier velocity. The WALE model shows a better agreement with the measurements. Different boundary conditions are tested for injection of notional particles. It is demonstrated both dispersed and carrier phase velocities are initially sensitive to the particle boundary condition whereas the difference between tested conditions becomes marginal further downstream. The best results are obtained for the particle boundary condition that accounts for effects of inertia at the inlet. Finally, the sensitivity of the PDF-PBE simulations to the number of Stokes bins is studied. It is found that eight Stokes bins are enough to accurately model the polysized particles dispersion in a bluff body configuration. The results indicate that particle dispersion is notably sensitive to the number of Stokes bins but the particle velocity predictions are much less sensitive with small variation in velocity results over the range of Stokes bins that were tested.     

基于离散单元法和人工神经网络的近壁颗粒动力学特征研究

颗粒与壁面的相互作用往往对颗粒流动具有显著影响. 为研究颗粒与壁面作用机理, 对滚筒内颗粒流动过程进行离散单元法(DEM)数值模拟. 基于模拟结果统计分析靠近壁面处颗粒的运动特征, 结果表明, 小摩擦系数时颗粒平动和旋转速度均近似满足正态分布, 但由于壁面影响, 摩擦系数增大时颗粒沿滚筒轴向的旋转速度偏离正态分布, 颗粒动力学理论推导壁面边界条件时应考虑速度正态分布的修正及速度脉动的各向异性. 采用人工神经网络(ANN)构建了颗粒无因次旋转温度、滑移速度和平动温度之间的函数模型, 进而可以在常规双流模型壁面边界条件中考虑颗粒旋转的影响. 基于DEM模拟及结果分析可以为壁面边界条件的理论构造和半经验修正提供基础数据和封闭模型.