A modified lattice Boltzmann model for shallow water flows over complex topography

A modified lattice Boltzmann model is proposed to describe shallow water flows over complex topography. In the proposed model, the quadratic depth term is excluded from the equilibrium distribution functions (EDFs), and the hydrostatic pressure term is combined with the bed slope term to be treated as a part of the sourcing term in the lattice Boltzmann equation (LBE). Therefore, it is unnecessary to match the coefficients of the quadratic depth term in the EDFs with those of the bed slope term in the sourcing terms in the LBE. This would bring more flexibility to the treatment of the sourcing terms in the LBE. In order to recover the shallow water equations (SWEs), the basic constraints are redefined, and under these constraints, the coefficients of the EDFs are derived afterwards. Several benchmark problems are used to validate the proposed model, including stationary case, steady flows over a two‐dimensional bump and tidal wave flows over irregular bed elevation. The computed results are in excellent agreement with the results of the other numerical methods and the analytical solutions, indicating that the proposed model is capable of simulating shallow water flows over complex bathymetry. It also proves that the proposed model has potential to produce competitive solutions to shallow water flows over complex bed topography. Copyright © 2014 John Wiley & Sons, Ltd.     

An evaluation of force terms in the lattice Boltzmann models in simulating shallow water flows over complex topography

Proper approximation of the force terms, especially the bed slope term, is of crucial importance to simulating shallow water flows in lattice Boltzmann (LB) models. However, there is little discussion on the schemes of adding force terms to LB models for shallow water equations (SWEs). In this study, we evaluate the performance of forcing schemes coupled with different LB models (LABSWE and MLBSWE) in simulating shallow water flows over complex topography and try to find out their intrinsic characteristics and applicability. Three cases are adopted for evaluation, including a stationary case, a one-dimensional tidal wave flow over an irregular bed, and a steady flow over a two-dimensional seamount. The simulating results are compared with analytical solutions or the results produced by the finite difference method. For LABSWE, all the forcing schemes, except for the weighting factor method, fail to produce accurate solutions for the test cases; this is probably due to the mismatch between the bed slope term in source terms and the quadratic depth term of the equilibrium distribution functions in these forcing schemes. For MLBSWE, all the forcing schemes are capable of simulating flows over the complex topography accurately; furthermore, those schemes taking into account the collision effect τ to eliminate the momentum induced by forces provide more accurate solutions with quicker convergence as the lattice size decreases. In this view, MLBSWE can bring more flexibility in treating the force terms and thus can be a better tool to simulate shallow water flows over complex topography in practical application.     

An unstructured finite volume model for dam‐break floods with wet/dry fronts over complex topography

A robust, well‐balanced, unstructured, Godunov‐type finite volume model has been developed in order to simulate two‐dimensional dam‐break floods over complex topography with wetting and drying. The model is based on the nonlinear shallow water equations in hyperbolic conservation form. The inviscid fluxes are calculated using the HLLC approximate Riemann solver and a second‐order spatial accuracy is achieved by implementing the MUSCL reconstruction technique. To prevent numerical oscillations near shocks, slope‐limiting techniques are used for controlling the total variation of the reconstructed field. The model utilizes an explicit two‐stage Runge–Kutta method for time stepping, whereas implicit treatments for friction source terms. The novelties of the model include the flux correction terms and the water depth reconstruction method both for partially and fully submerged cells, and the wet/dry front treatments. The proposed flux correction terms combined with the water depth reconstruction method are necessary to balance the bed slope terms and flux gradient in the hydrostatical steady flow condition. Especially, this well‐balanced property is also preserved in partially submerged cells. It is found that the developed wet/dry front treatments and implicit scheme for friction source terms are stable. The model is tested against benchmark problems, laboratory experimental data, and realistic application related to dam‐break flood wave propagation over arbitrary topography. Numerical results show that the model performs satisfactorily with respect to its effectiveness and robustness and thus has bright application prospects. Copyright © 2010 John Wiley & Sons, Ltd.     

Zero mass error using unsteady wetting–drying conditions in shallow flows over dry irregular topography

A wetting–drying condition (WDC) for unsteady shallow water flow in two dimensions leading to zero numerical error in mass conservation is presented in this work. Some applications are shown which demonstrate the effectiveness of the WDC in flood propagation and dam break flows over real geometries. The WDC has been incorporated into a cell centred finite volume method based on Roe's approximate Riemann solver across the edges of both structured and unstructured meshes. Previous wetting–drying condition based on steady‐state conditions lead to numerical errors in unsteady cases over configurations with strong variations on bed slope. A modification of the wetting–drying condition including the normal velocity to the cell edge enables to achieve zero numerical errors. The complete numerical technique is described in this work including source terms discretization as a complete and efficient 2D river flow simulation tool. Comparisons of experimental and numerical results are shown for some of the applications. Copyright © 2004 John Wiley & Sons, Ltd.     

Shoreline tracking and implicit source terms for a well balanced inundation model

The HyFlux2 model has been developed to simulate severe inundation scenario due to dam break, flash flood and tsunami‐wave run‐up. The model solves the conservative form of the two‐dimensional shallow water equations using the finite volume method. The interface flux is computed by a Flux Vector Splitting method for shallow water equations based on a Godunov‐type approach. A second‐order scheme is applied to the water surface level and velocity, providing results with high accuracy and assuring the balance between fluxes and sources also for complex bathymetry and topography. Physical models are included to deal with bottom steps and shorelines. The second‐order scheme together with the shoreline‐tracking method and the implicit source term treatment makes the model well balanced in respect to mass and momentum conservation laws, providing reliable and robust results. The developed model is validated in this paper with a 2D numerical test case and with the Okushiri tsunami run up problem. It is shown that the HyFlux2 model is able to model inundation problems, with a satisfactory prediction of the major flow characteristics such as water depth, water velocity, flood extent, and flood‐wave arrival time. The results provided by the model are of great importance for the risk assessment and management. Copyright © 2009 John Wiley & Sons, Ltd.     

A robust 2D shallow water model for solving flow over complex topography using homogenous flux method

A robust Godunov‐type numerical scheme solver is proposed for solving 2D SWEs and is applied to simulate flow over complex topography with wetting and drying. In reality, the topography is usually complex and irregular; therefore, to avoid the numerical errors generated by such features, a Homogenous Flux Method is used to handle the bed slope term in the SWEs. The method treats the bed slope term as a flux to be incorporated into the flux gradient and so maintains the balance between the two in a Godunov‐type shock‐capturing scheme. The main advantages of the method are: first, it is simple and easy to implement; second, numerical experiments demonstrate that it can handle discontinuous or vertical bed topography without any special treatment and third, it is applicable to both steady and unsteady flows. It is demonstrated how the approach set out here can be applied to the nonlinear hyperbolic system of the SWEs. The two‐dimensional hyperbolic system is then solved by use of a second‐order total‐variation‐diminishing version of the weighted average flux method in conjunction with a Harten‐Lax‐van Leer‐Contract approximate Riemann solver incorporating the new flux gradient term. Several benchmark tests are presented to validate the model and the approach is verified against experimental measurements from the European Union Concerted Action on Dam Break Modelling project. These show very good agreement. Finally, the method is applied to a volcano‐induced outburst flood over an initially dry channel with complex irregular topography to demonstrate the technique's capability in simulating a real flood. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical and experimental study on the motion characteristics of single bubble in a complex channel

This paper is an extended study from previous work. In this study, the focus is paid to the dynamics of bubble rising and deformation in a complex channel, while the previous work is in straight channel. For this purpose, a three-dimensional lattice Boltzmann method (LBM) is employed to simulate the dynamics behaviour of a bubble rising in a complex channel consisting of three half-round throats. To validate the numerical method, a visual experiment was carried out by means of a high-speed digital camera and computer image processing technology. The behaviour of the rising bubble through glycerine solution in a complex channel was recorded. Some physical parameters such as rising velocities, trajectory and shapes of the bubble were calculated and processed based on the experimental data. In the same conditions, the trajectory, shapes and rising velocities of the bubble were simulated during its rising process by the proposed LBM. The numerical results are in good agreement with the experimental results. It demonstrates that LBM used in this work is feasible for simulating two-phase flow in such a complex channel.     

An immersed boundary method wall model for high‐Reynolds‐number channel flow over complex topography

High‐Reynolds‐number channel flows regularly encounter topographies composed of multiple length scales and that protrude into the boundary layer. Physically, the presence of immersed obstacles leads to increased velocity gradients, turbulence production, and manifestation of wakes. Considerable challenges are associated with numerically describing the presence of obstacles in channel flows. Common approaches include generation of a computational mesh that is uniquely designed for the flow and obstacle, the immersed boundary method, and terrain‐following coordinates. There are challenges and limitations associated with each of these techniques. Specification of boundary conditions representing the perimeter of solid obstacles is a primary challenge of the immersed boundary method. In this document, a simplistic canopy stress‐like wall model is used to impose boundary conditions. The model isolates aerodynamically relevant local frontal areas through evaluation of the gradient of the topographic height field. The gradient of the height field describes both the surface‐normal direction and the frontal area, making it ideal for detecting areas on which the flow impinges. The model is tested in numerical simulations of turbulent half‐channel flow over topographies with different obstacles affixed–right prisms, rectangular prisms, ellipsoidal mounds, and sinusoids. In all cases, the performance is strong relative to datasets presented in the literature. Results are finally presented for numerical simulation of flow over complex synthetic fractal‐like topography and a synthetic city. These results show interesting trends in how the turbulent multiscale flow field responds to multiscale topography. Copyright © 2012 John Wiley & Sons, Ltd.     

复杂地形网格生成研究

在研究羽流扩散的过程中,根据实验中提供的复杂地形数据,采用一系列的格式转化,将其变为流体力学计算软件Fluent和其建模软件Gambit能够读取的数据,并生成相应的计算网格,为计算复杂地形上的羽流扩散创造了条件。     

An experimental study on boundary layer transition detection over a pitching supercritical airfoil using hot-film sensors

In the present work, experimental tests are conducted to study boundary layer transition over a supercritical airfoil undergoing pitch oscillations using hot-film sensors. Tests have been undertaken at an incompressible flow. Three reduced frequencies of oscillations and two mean angles of attack are studied and the influences of those parameters on transition location are discussed. Different algorithms are examined on the hot-film signals to detect the transition point. Results show the formation of a laminar separation bubble near the leading edge and at relatively higher angles of attack which leads to the transition of the boundary layer. However, at lower angles of attack, the amplification of the peaks in voltage signal indicate the emergence of the vortical structures within the boundary layer, introducing a different transition mechanism. Moreover, an increase in reduced frequency leads to a delay in transition onset, postponing it to a higher angle of attack, which widens the hysteresis between the upstroke and downstroke motions. Rising the reduced frequency yields in weakening or omission of vortical disturbances ensuing the removal of spikes in the signals. Of the other important results observed, is faster movement of the relaminarization point in the higher mean angle of attack. Finally, a time–frequency analysis of the hot-film signals is performed to investigate evolution of spectral features of the transition due to the pitching motion. An asymmetry is clearly observed in frequency pattern of the signals far from the bubble zone towards the trailing edge; this may reflect the difference between the transition and relaminarization physics. Also, various ranges of frequency were obtained for different transition mechanisms.     

Biomechanical research of joint III. An experimental biomechanical research on the Femur's articular surface of knee of pongidae

In this paper, moiré contour fringes is applied to study the femur's articular surface of the knee of Pongidae. The preliminary division of the femur's articular surface of knee into three types is proposed. The moiré contour fringes ofthe medial condyle is taken as a mark according to the references. Owing to the fact that the moiré contour fringes obtained from experiments after the 2nd order of fringe basically follow a certain rule, an investigation is made on the distribution of the angle α which is defined as the angle of the major axis of the 2nd order's near-oval shaped moiré contour fringe on the medial condyle with the horizontal axis. Preliminary distribution graphs are given in the paper.     

爆炸冲击波在长直坑道和复杂坑道内传播规律的对比研究

利用非线性动力分析软件LS-DYNA,对炸药在坑道入口处爆炸情况下,爆炸流场在长直坑道、＂L＂形坑道、＂S＂形坑道内的变化过程进行了数值仿真,得到了爆炸冲击波在不同形状坑道中的传播规律。研究结果表明：带弯道设计的坑道可以有效地削弱冲击波强度,提高坑道的抗破坏能力;弯道外侧峰值压力比同位置长直坑道中的峰值压力大,因此在进行带弯道的坑道设计时,应相应地加强弯道外侧坑道壁的强度。该方法适合于解决爆炸冲击波在复杂弯曲坑道内传播的问题,为今后的坑道防护设计提供了参考和依据。     

高密度B炸药的燃烧转爆轰实验研究

采用电探针及压力传感器测试技术对密度为1.597 g/cm3的固体B炸药（TNT/RDX=40/60）的燃烧转爆轰性能进行研究。实验结果表明,在较强的约束条件下（45号钢管,内径20 mm,外径64 mm,长500 mm）,B炸药形成了DDT现象,诱导爆轰距离为295~310 mm。     

纯钽动态本构关系的实验研究

利用材料试验机和分离式霍布金森压杆实验装置测量了纯钽在较宽温度(-100～550℃)和应变率(10-5～103s-1)范围内的应力 应变曲线,该曲线表明钽的流动应力对温度和应变率都很敏感。利用所测得的应力 应变曲线拟合了钽在不同应变率和不同温度下的Johnson Cook与Zerilli Armstrong本构方程,Zeril li Armstrong本构方程与实验数据基本吻合,Johnson Cook本构方程与实验数据吻合得不太好。对Johnson Cook本构方程中的应变率强化项作了修正,修正后的Johnson Cook本构方程与实验数据吻合比修正前的好。     

An experimental study on the thermal performance of ground heat exchanger

A knowledge of ground thermal properties is most important for the proper design of large GHE (ground heat exchanger) systems. Thermal response tests have so far been used primarily not only for in situ determination of design data for GHE systems, but also for the evaluation of grout material, heat exchanger types and groundwater effects. The main purpose has been to determine in situ values of effective ground thermal conductivity, including the effect of groundwater flow and natural convection in boreholes.     

An experimental study on turbulent coherent structures near a sheared air-water interface

The turbulence structures near a sheared air-water interface were experimentally investigated with the hydrogen bubble visualization technique. Surface shear was imposed by an airflow over the water flow which was kept free from surface waves. Results show that the wind shear has the main influence on coherent structures under air-water interfaces. Low- and high- speed streaks form in the region close to the interface as a result of the imposed shear stress. When a certain airflow velocity is reached, “turbulent spots” appear randomly at low-speed streaks with some characteristics of hairpin vortices. At even higher shear rates, the flow near the interface is dominated primarily by intermittent bursting events. The coherent structures observed near sheared air-water interfaces show qualitative similarities with those occurring in near-wall turbulence. However, a few distinctive phenomena were also observed, including the fluctuating thickness of the instantaneous boundary layer and vertical vortices in bursting processes, which appear to be associated with the characteristics of air-water interfaces. The project supported by the National Natural Science Foundation of China (Grant No.19672070)     

高温下轻质泡沫铝动态力学性能实验

对传统的分离式Hopkinson压杆装置加以改进,设计了一种长杆直接撞击Hopkinson杆的实验方案,检测出低波阻抗材料在高温动态加载下的应力均匀性。对轻质泡沫铝材料的实验表明,在同一撞击速度下,温度越高,试件两端的应力均匀性越差,增加温度与提高撞击速度均会导致泡沫铝材料冲击端与支撑端的应力不均匀性。根据高温下应力均匀性的实验结果,确定高温下试件均匀变形对应的冲击速度,再通过传统的分离式Hopkinson压杆实验得出泡沫铝在高温动态下的力学性能。

     

An experimental study on the static and impact torsional buckling of cylindrical shells

In this paper, experimental studies are carried out on the buckling of circular cylindrical thin shells under impact torque. Experiments of impact buckling are made on a Hopkinson torsional bar. The torsional bar gives a step torque on the shells. Through an analysis of the strain-time curve obtained in experiment, the dynamic buckling critical torqueM _er and buckling waves numbern of the shell with different geometric data and some qualitative results are obtained. The buckling behavior of circular cylindrical thin shells under static and impact torque is compared.     

An experimental study on the interactions between surface waves and floating viscoelastic covers

A possible mathematical ice model for the wave interactions in polar seas was developed based on the assumption that an ice cover behaved as a Voigt viscoelastic material. The dispersion relation was found to depend on the rheological properties of the cover. In the present study, an experimental approach was developed that can enable the verification of the theoretical predictions in the laboratory. The approach utilized the blended mixture of white oil and Polydimethylsiloxane (PDMS) material with various mass percentages of a curing agent, to create a floating layer with a range of targeted viscoelastic properties. Due to the large coverage required for wave flume experiments, special curing procedures were also established for the preparation of PDMS material. The rheological results showed that the mechanical behavior of the floating cover was close to a Voigt material. Experiments were conducted to analyze the wave interactions with the floating viscoelastic cover. The measured data showed an obvious change of wavelength when waves propagated along the cover region. It is observed that the change in wavelength can be linked quantitatively to the viscoelastic properties based on the numerical predictions by Wang and Shen (2010). Some differences were however noted for less viscous covers under longer wave periods. A direct comparison of the PDMS covers with a polypropylene (PP) cover was also performed for verification. Only wave lengthening was observed under the PP cover. With a shear modulus more than three orders of magnitude greater than that of PDMS, the theoretical wavelength for the PP cover from Wang and Shen (2010) is very close to that of the thin elastic plate theory from Fox and Squire (1990). Comparison between these two theoretical results and the measured data again deviated with longer wave periods. In both PDMS and PP cases, edge effects and pitching motion of the covers were present at various degrees. In addition, the materials were not strictly a Voigt type. The small deviation from the idealized rheological behavior could also contribute to differences between theoretical and experimental results.     

氢气/丙烷/空气预混气体爆轰性能的实验研究

通过采用压力传感器和烟灰板两种测试设备,开展了常温常压下氢气/丙烷和空气混合气体爆轰性能的实验研究。实验过程中观察到自持爆轰波,爆轰速度比值在0.99~1之间,爆轰压力比值在0.8~1.2之间。爆轰胞格尺寸在10~50 mm范围内,建立了爆轰胞格尺寸和化学诱导长度的关系式。随着丙烷不断添加,爆轰速度减小,而爆轰压力和胞格尺寸增加。这种变化趋势起初较快,而后变缓。因为起初氢气摩尔分数较大,混合气体趋向于氢气/空气的爆轰性能;而后因丙烷摩尔质量较大,丙烷逐渐起主要作用,混合气体表现出丙烷/空气的爆轰性能。