A boundary element method for the simulation of non-spherical bubbles and their interactions near a free surface

The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4 π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demonstrated in our calculation results.     

Validation of a CIP-based tank for numerical simulation of free surface flows

A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the present paper.The tangent of hyperbola for interface capturing(THINC) scheme is applied for capturing complex free surfaces.The new model is capable of simulating a flow with violently varied free surface.A series of computations are conducted to assess the developed algorithm and its versatility.These tests include the collapse of water column with and without an obstacle,sloshing in a fixed tank,the generation of regular waves in a tank,the generation of extreme waves in a tank.Excellent agreements are obtained when numerical results are compared with available analytical,experimental,and other numerical results.     

Stability of aneurysm solutions in a fluid-filled elastic membrane tube

When a hyperelastic membrane tube is inflatedby an internal pressure, a localized bulge will form when thepressure reaches a critical value. As inflation continues thebulge will grow until it reaches a maximum size after whichit will then propagate in both directions to form a hat-likeprofile. The stability of such bulging solutions has recentlybeen studied by neglecting the inertia of the inflating fluidand it was shown that such bulging solutions are unstableunder pressure control. In this paper we extend this recentstudy by assuming that the inflation is by an inviscid fluidwhose inertia we take into account in the stability analysis.This reflects more closely the situation of aneurysm forma-tion in human arteries which motivates the current series ofstudies. It is shown that fluid inertia would significantly re-duce the growth rate of the unstable mode and thus it has astrong stabilizing effect.     

Wind tunnel simulation of wind loading on a solid structure of revolution

The wind tunnel simulations of wind loading on a solid structure of revolution with one smooth and five rough surfaces were conducted using wind tunnel tests. Timemean and fluctuating pressure distributions on the surface were obtained, and the relationships between the roughness Reynolds number and pressure distributions were analyzed and discussed. The results show that increasing the surface roughness can significantly affect the pressure distribution, and the roughness Reynolds numbers play an important role in the change of flow patterns. The three flow patterns of subcritical, critical and supercritical flows can be classified based on the changing patterns of both the mean and the fluctuating pressure distributions. The present study suggests that the wind tunnel results obtained in the supercritical pattern reflect more closely those of full-scale solid structure of revolution at the designed wind speed.     

Emission of dimers from a free surface of heated water

The emission rate of water dimers from a free surface and a wetted solid surface in various cases was calculated by a simplified Monte Carlo method with the use of the binding energy of water molecules. The binding energy of water molecules obtained numerically assuming equilibrium between the free surface of water and vapor in the temperature range of 298–438 K corresponds to the coordination number for liquid water equal to 4.956 and is close to the reference value. The calculation results show that as the water temperature increases, the free surface of water and the wetted solid surface become sources of free water dimers. At a temperature of 438 K, the proportion of dimers in the total flow of water molecules on its surface reaches 1%. It is found that in the film boiling mode, the emission rate of dimers decreases with decreasing saturation vapor. Two mechanisms of the emission are described.     

变深度浅水域中非定常船波

以Green—Naghdi(G—N)方程为基础，采用波动方程／有限元法计算船舶经过变深度浅水域时非定常波浪特性．把运动船舶对水面的扰动作为移动压强直接加在Green-Naghdi方程里，以描述运动船体和水面的相互作用．以Series60 CB＝0．6船为算例，给出自由面坡高，波浪阻力在船舶经过一个水下凸包时变化规律，并与浅水方程的结果进行了比较．计算结果表明，当船舶经过凸包时，波浪阻力先增加，后减少，并逐渐趋于正常．同时发现，当船速小于临界速度时(Fr＝√gh＜1．0)，G—N方程给出的船后尾波比浅水方程的结果明显，波浪阻力也比浅水方程的结果有所提高，频率散射必须考虑．当船速大于临界速度时(Fr＝√gh＞1．0)，G—N方程的计算结果与浅水方程差别不大，频率散射的影响可以忽略．     

Experimental and numerical studies of micro PEM fuel cell

A single micro proton exchange membrane fuel cell (PEMFC) has been produced using Micro-electromechanical systems (MEMS) technology with the active area of 2.5 cm 2 and channel depth of about 500 μ m.A theoretical analysis is performed in this study for a novel MEMS-based design of a micro PEMFC.The model consists of the conservation equations of mass,momentum,species and electric current in a fully integrated finite-volume solver using the CFD-ACE+ commercial code.The polarization curves of simulation are well correlated with experimental data.Three-dimensional simulations are carried out to treat prediction and analysis of micro PEMFC temperature,current density and water distributions in two different fuel flow rates (15 cm 3 /min and 40 cm 3 /min).Simulation results show that temperature distribution within the micro PEMFC is affected by water distribution in the membrane and indicate that low and uniform temperature distribution in the membrane at low fuel flow rates leads to increased membrane water distribution and obtains superior micro PEMFC current density distribution under 0.4 V operating voltage.Model predictions are well within those known for experimental mechanism phenomena.     

Flutter and resonances of a flag near a free surface

We investigate the effects of a nearby free surface on the stability of a flexible plate in axial flow. Confinement by rigid boundaries is known to affect flag flutter thresholds and fluttering dynamics significantly, and this work considers the effects of a more general confinement involving a deformable free surface. To this end, a local linear stability is proposed for a flag in axial uniform flow and parallel to a free surface, using one-dimensional beam and potential flow models to revisit this classical fluid–structure interaction problem. The physical behaviour of the confining free surface is characterized by the Froude number, corresponding to the ratio of the incoming flow velocity to that of the gravity waves. After presenting the simplified limit of infinite span (i.e. two-dimensional problem), the results are generalized to include finite-span and lateral confinement effects. In both cases, three unstable regimes are identified for varying Froude number. Rigidly-confined flutter is observed for low Froude number, i.e. when the free surface behaves as a rigid wall, and is equivalent to the classical problem of the confined flag. When the flow and wave velocities are comparable, a new instability is observed before the onset of flutter (i.e. at lower reduced flow speed) and results from the resonance of a structural bending wave and one of the fundamental modes of surface gravity waves. Finally, for large Froude number (low effect of gravity), flutter is observed with significant but passive deformation of the free surface in response of the flag’s displacement.     

A new model for wrinkling and collapse analysis of membrane inflated beam

A new model is proposed to accurately predict the wrinkling and collapse loads of a membrane inflated beam. In this model, the pressure effects are considered and a modified factor is introduced to obtain an accurate prediction. The former is achieved by modifying the pressure-related structural parameters based on elastic small strain considerations, and the modified factor is determined by our test data. Compared with previous models and our test data, the present model, named as shell-membrane model, can accurately predict the wrinkling and collapse loads of membrane inflated beams.     

强约束球形装药反应裂纹传播和反应烈度表征实验

炸药燃烧的高温高压气体产物可以进入基体裂纹中引发炸药表面热传导燃烧,形成所谓的对流燃烧。在一定约束条件下,不断上升的气体压力反过来又使炸药基体产生更多的裂纹,为对流燃烧提供更多的通道和燃烧表面积,快速生成大量产物气体导致高烈度反应现象的产生。本文中设计了一种新型强约束球形装药中心点火实验,针对一种HMX为基的PBX炸药,对高烈度反应条件下燃烧裂纹传播和反应增长过程进行了观测,实验中采用测得的反应压力和壳体速度历程对反应烈度进行了量化表征。在带窗口结构中,早期炸药中的燃烧裂纹不可见;中期燃烧裂纹扩展到药球表面时,先形成4条沿经线方向近似对称的主裂纹,随后环向贯通并扩展到整个药球表面;最后的剧烈反应造成强烈发光。上述反应演化经历低压增长阶段约为100 μs,之后伴随着壳体变形膨胀产生剧烈的反应,此时产物压力在约10 μs时间内超过1 GPa,并形成约20%相对于裸炸药爆轰的超压输出。在全钢结构中,20 mm厚的壳体膨胀速度最大可达到500 m/s,此时壳体完全破裂。     

Some aspects of sheet-metal forming using impact loading through water columns

This paper reports on work carried out on the free forming of 50.8-mm-diam low-carbon-steel sheet specimens using a water column as the energy-transmilling medium. A ‘simple’ theory is presented that allows an estimation of the pressure variation, within the water column, with time to be made. The energy is supplied from a falling weight that impacts a punch testing on the top surface of the water column. The large mass of the falling weight relative to the mass of the water column allows an energy balance to be made between the kinetic energy of the weight and the strain energy absorbed by the column. Experimental results are presented that show the theroy to be capable of providing a reasonable estimation of pressure variation.     

MD simulation of a copper rod under thermal shock

In this paper, thermoelastic problem of onedimensional copper rod under thermal shock is simulated using molecular dynamics method by adopting embedded atom method potential. The rod is on axis x, the left outermost surface of which is traction free and the right outermost surface is fixed. Free boundary condition is imposed on the outermost surfaces in direction y and z. The left and right ends of the rod are subjected to hot and cold baths, respectively. Temperature, displacement and stress distributions are obtained along the rod at different moments, which are shown to be limited in the mobile region, indicating that the heat propagation speed is limited rather than infinite. This is consistent with the prediction given by generalized thermoelastic theory. From simulation results we find that the speed of heat conduction is the same as the speed of thermal stress wave. In the present paper, the simulations are conducted using the large-scale atomic/molecular massively parallel simulator and completed visualization software.     

Experimental study on correlation between turbulence and sound in a subsonic wind tunnel

In this paper, the effects of turbulence on sound generation and velocity fluctuations due to pressure waves in a large subsonic wind tunnel are studied. A trip strip located at different positions in the contraction part or at one position in the diffuser of a large wind tunnel is used to investigate the aforementioned phenomenon, and the results indicate that the trip strip has significant effects on sound reduction. The lowest turbulence intensity and sound are obtained from a trip strip with a diameter of 0.91 mm located either at X/L = 0.79 or at X/L = 0.115 in the wide portion of the contraction. Furthermore, the effect of monopole, dipole and quadrupole sources of aerodynamic noise at different velocities is investigated, and it is demonstrated that the contribution of the monopole is dominant, while the shares due to the dipole and quadrupole remain less important. In addition, it is found that the sound waves have a modest impact on the measured longitudinal turbulence and are generated essentially by eddies.     

Wave scattering by flexible porous vertical membrane barrier in a two-layer fluid

The scattering of water waves by a flexible porous membrane barrier in a two-layer fluid having a free surface is analysed in two dimensions. The membrane barrier is extended over the entire water depth in a two-layer fluid, each fluid being of finite depth. In the present analysis, linear wave theory and small amplitude membrane response are assumed. The porous membrane barrier is tensioned and pinned at both the free surface and the seabed. The associated mixed boundary value problem is reduced to a linear system of equations by utilizing a general orthogonality relation along with least-squares approximation method. Because of the flow discontinuity at the interface, the eigenfunctions involved have a discontinuity at the interface and the orthogonality relation used is a generalization of the classical one corresponding to a single-layer fluid. The reflection and transmission coefficients for the surface and internal modes, the free surface and interface elevations and the nondimensional membrane deflection are computed for various physical parameters like the nondimensional tension parameter, porous-effect parameter, fluid density ratio, ratio of water depths of the two fluids to analyse the efficiency of a porous membrane as a wave barrier in the two-layer fluid.     

Liquid crystal model of vesicle deformation in alternating electric fields

Considering the effects of osmotic pressure, elastic bending, Maxwell pressure, surface tension, as well as flexo-electric and dielectric properties of phospholipid membrane, the shape equation for sphere vesicle in alternation (AC) electric field is derived based on the liquid crystal model by minimizing the free energy due to coupled mechanical and AC electrical fields. Besides the effect of elastic bending, the influence of osmotic pressure and surface tension on the frequency dependent behavior of vesicle membrane in AC electric field is also discussed. Our theoretical results for membrane deformation are consistent with corresponding experiments. The present model provides the possibility to further disclose the frequency-depended behavior of biological cells in the coupled AC electric and different mechanical fields.     

Coupled bending-torsion vibration of a homogeneous beam with a single delamination subjected to axial loads and static end moments

Delaminations in structures may significantly reduce the stiffness and strength of the structures and may affect their vibration characteristics. As structural components, beams have been used for various purposes, in many of which beams are often subjected to axial loads and static end moments. In the present study, an analytical solution is developed to study the coupled bending-torsion vibration of a homogeneous beam with a single delamination subjected to axial loads and static end moments. Euler–Bernoulli beam theory and the "free mode" assumption in delamination vibration are adopted. This is the first study of the influences of static end moments upon the effects of delaminations on natural frequencies, critical buckling loads and critical moments for lateral instability. The results show that the effects of delamination on reducing natural frequencies, critical buckling load and critical moment for lateral instability are aggravated by the presence of static end moment. In turn, the effects of static end moments on vibration and instability characteristics are affected by the presence of delamination. The analytical results of this study can serve as a benchmark for finite element method and other numerical solutions.     

Boundary layer flow over a moving surface in a nanofluid with suction or injection

An analysis is performed to study the heat transfer characteristics of steady two-dimensional boundary layer flow past a moving permeable flat plate in a nanofluid. The effects of uniform suction and injection on the flow field and heat transfer characteristics are numerically studied by using an implicit finite difference method. It is found that dual solutions exist when the plate and the free stream move in the opposite directions. The results indicate that suction delays the boundary layer separation, while injection accelerates it.     

Dynamical behavior of a mobile system with two degrees of freedom near the resonance

The design of mobile robots that can move without wheels or legs is an important engineering and technological problem.Self-propelling mechanisms consisting of a body that has contact with a rough surface and moveable internal masses are considered.Mathematical models of such systems are presented in this paper.First,a model of a vibration driven robot that moves along a rough horizontal plane with isotropic dry friction is studied.It is shown that by changing the off-resonance frequency detuning in sign,one can control the direction of motion of the system.In addition,a locomotion system which moves in an environment with anisotropic viscous friction is considered.For all models,the method of averaging to obtain an algebraic equation for the steady-state"average"velocity of the system is used. Prototypes were constructed to compare the theoretical results with experimental ones.     

Mechanics of the implosion of cylindrical shells in a confining tube

A fundamental experimental investigation, with corresponding computational simulations, was conducted to understand the physical mechanisms of implosions of cylindrical shells occurring within a tubular confining space which has a limited potential energy reservoir. In particular, attention was focused on studying the generation of pressure waves from the implosion, the interaction of the pressure waves with the confining tube walls and end caps, and the collapse mechanisms of the implodable volume. Experiments were conducted with three implodable volume geometries which had similar critical collapse pressures. The implodable volumes were aluminum 6061-T6 cylindrical tubing and were placed concentrically within the confining tube. Pressure histories recorded along the length of the confining tube during the experiments were utilized to analytically evaluate the deformation of the implodable volume using fluid–structure coupled deformation models. Computational simulations were conducted using a coupled Eulerian–Lagrangian scheme to explicitly model the implosion process of the tubes along with the resulting compressible fluid flow. The numerical model developed in this study is shown to have high correlation with the experimental results and will serve as a predictive tool for the simulation of the implosion of different cylindrical geometries as well as various tube-in-tube implosion configurations. The experimental results show that the limited hydrostatic potential energy available in a confined environment, as compared to a free field, significantly influences the implosion process. The wall velocities of the implodable volume during the collapse, as well as the extent of the collapse progression, are largely affected by the sudden decrease in the available hydrostatic potential energy. This energy is shown to be partially transformed into elasto-plastic strain energy absorbed in the deformation of the implodable volume, as well as the kinetic energy of the water during the implosion process. Experiments also show that the extent of the collapse progression of an implodable volume can potentially be inhibited within a closed environment, which can lead to the arresting of an implosion event prior to completion for larger implodable volumes. The pressure waves generated during collapse comprise of waves emitted due to the impact of the implodable volume walls, the arrest of rushing water and contact propagation along the walls. These processes later evolve into water hammer type axial wave behavior.     

The SPH approach to the process of container filling based on non-linear constitutive models

In this work,the transient free surface of container filling with non-linear constitutive equation’s fluids is numerically investigated by the smoothed particle hydrodynamics(SPH) method.Specifically,the filling process of a square container is considered for non-linear polymer fluids based on the Cross model.The validity of the presented SPH is first verified by solving the Newtonian fluid and OldroydB fluid jet.Various phenomena in the filling process are shown,including the jet buckling,jet thinning,splashing or spluttering,steady filling.Moreover,a new phenomenon of vortex whirling is more evidently observed for the Cross model fluid compared with the Newtonian fluid case.