Effect of liquid viscoelasticity on the interactions between cavitation bubbles and free surfaces

The influence of liquid viscoelasticity on the interaction between cavitation bubbles and free surfaces is of great practical significance in understanding bubble dynamics in biological systems. A series of millimeter cavitation bubbles were induced by laser near the free surfaces of the water and viscoelastic polyacrylamide (PAM) solutions with different concentrations. The effects of liquid viscoelasticity on the interactions of cavitation bubbles with free surfaces are analyzed from the perspectives of the evolution of free surface and bubble dynamics. The experimental results show that as the dimensionless standoff distance increases, the evolutions of free surface behaviors in all experimental fluids can be divided into six types of water mounds, i.e., breaking wrinkles, spraying water film, crown, swallowed water spike, hillock, and slight bulge. All the critical values of the dimensionless distance dividing different types decrease with increasing concentration. The evolutions of first four types of water mounds in PAM solutions differ from those in the water. Water droplets splashing in different directions are produced around the breaking wrinkles in the water. Meanwhile, the breaking wrinkles in PAM solution move with the “liquid filaments” towards the central axis. The water spike in the pattern of spraying water film in PAM solution is more stable than that in the water. As the solution concentration increases, the water skirt in the pattern of crown contracts earlier and faster, and the rate of increase in the height of the water skirt decreases. For swallowed water spike in PAM solution, the upper part of the newly formed water spike is not significantly thicker than the middle part, and thus the water waist structure does not form. Liquid viscoelasticity inhibits the bubble growth and collapse, and the bubble migration as well, especially in the second period. Shorter and thicker cavities are formed in PAM solutions with higher concentration, while slender and stable cavities formed in the water at the same dimensionless distance. The velocity and displacement of the tip of bullet jet both decrease as the solution concentration increases.

Investigation of the Interaction Between the Free Surface and A Semi/Shallowly Submerged Underwater Vehicle

The present study aims to investigate the interaction between the free surface and a semi/shallowly submerged underwater vehicle,especially when the submergence depth h is smaller than 0.75D(D:submarine maximum diameter).In this respect,the straight-ahead simulations of the generic SUBOFF underwater vehicle geometry are conducted with constant forward velocities using the Unsteady Reynolds-Averaged Navier-Stokes(URANS)equations with a Shear-Stress Transport(SST)k-ωturbulence model in commercial code Fluent,at submergence depths and Froude numbers ranging from h=0 to h=3.3D and from Fn=0.205 to Fn=0.512,respectively.The numerical models are verified against the existing experimental data.The analysis of the obtained results indicates that in the case of the semi and shallowly submerged underwater vehicle(UV),both the submergence depth and forward velocity have a great effect on the behaviors of hydrodynamic forces acting on the UV.The magnitude of maximum total resistance may reach almost five times the value of resistance exerted on the totally submerged hull.Both the forces acting on the UV and the generated waves when the submergence depth h is smaller than 0.75D are significantly different from those whenr h is larger than 0.75D.The conclusions can be used as reference for future research on near free surface motions of underwater vehicles and the design of small water-plane area twin hull.

STUDY OF NON-CONTACT MEASUREMENT OF FREE-FORM SURFACES

A non contact three dimensional measurement method is presented in this paper.This system consists of a laser triangulation probe,a probe head and a coordinate measuring machine (CMM).The measurement principle of the system is discussed,and a system calibration method employing a reference ball is proposed.The geometric model involving four frames is established to calculate the data points based on the reading of the laser probe and position information from the CMM.A measuring experiment for gesso free form surface using this system is carried out.

Elastic wave scattering from rough free surfaces of solids

Elastic wave scattering by a rough free surface of solids is analyzed. The analysis is based on the concept of scattering amplitude (SA) and perturbation approximation. The SA method is very convenient for rough surface scattering problems. By solving the boundary equations, the first and the second order solutions of approximate scattering amplitude are obtained. The general solutions are used for, as an example, the wave scattering by rough surfaces with Gaussian distribution. The mean field and variance are given. Finally, an experiment is designed to verify the theoretical predications.

A High-Order Conservative Semi-Lagrangian Solver for 3D Free Surface Flows with Sediment Transport on Voronoi Meshes

In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons.

Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties

Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.

SUPPORT SYSTEM FOR FREE SURFACE GRINDING BASED ON PATCH DIVISION

To obtain accurate forms and surfaces in free surface grinding, it is important to provide grinding conditions suitable for a curved surface. A grinding support system for the free surface （GSX-F） is proposed to help the operator grind a free surface with the high accuracy and the high productivity. To succeed in free surface grinding, the property of a ball type wheel must be known. Therefore, a basic study of free surface grinding with a ball type wheel is carried out based on the grinding center （GC）. Some working points for achieving sufficient accuracy in free surface grinding are discussed. GSX-F is constructed using the patch division method and is used to test grinding. Reasonable results are obtained.

Temperature Drop of Molten Metals in Open Channels

The temperature drop of molten metal flowing in open channels is numerically determined. Rectangular, trapezoidal and triangular geometries are considered. The overall heat transfer coefficients for the bottom, side walls and free surface of the channel have been taken from the literature. For each geometry, the volumetric flow rate, mean residence time and temperature drop as a function of the channel inclination angle were determined. The rectangular and trapezoidal geometries present the smallest temperature drops, while the triangular geometry presents the greatest temperature drop. The factors that most affect this drop are the value of the free surface area of the channel, and the average residence time of the molten metal in the channel.

A Mesh-Independent Brute-Force Approach for Traction-Free Corrections in Dislocation Problems

Simulation of dislocation dynamics opens the opportunity for researchers and scientists to observe in-depth many plastic deformation phenomena. In 2D or 3D media, modeling of physical boundary conditions accurately is one of the keys to the success of dislocation dynamics (DD) simulations. The scope of analytical solutions is restricted and applies to specific configurations only. But in dynamics simulations, the dislocations’ shape and orientation change over time thus limiting the use of analytical solutions. The authors of this article present a mesh-based generalized numerical approach based on the collocation point method. The method is applicable to any number of dislocations of any shape/orientation and to different computational domain shapes. Several verifications of the method are provided and successful implementation of the method in 3D DD simulations have been incorporated. Also, the effect of free surfaces on the Peach-Koehler force has been computed. Lastly, the effect of free surfaces on the flow stress of the material has been studied. The results clearly showed a higher force with increased closeness to the free surface and with increased dislocation segment length. The simulations’ results also show a softening effect on the flow stress results due to the effect of the free surfaces.

Numerical study of three-dimensional free surface dynamics

The dynamic problem of three-dimensional free surface is numerically studied in this paper. The ALE （Arbitrary Lagrange-Euler） kinematic description is introduced into the control equation system. The ALE description method is used to track free surface. Accurate formulations for calculating the normal vector on the free surface are presented. The discrete numerical equations by finite element method are developed by Galerkin weighted residual method. The boundary condition about free-surface tension is represented in the form of weak integration that can be computed by a differential geometry method derived in the present paper. The effect of contact angle is incorporated in the numerical algorithm. Furthermore, the numerical computations are performed and the comparison between computational and analytical results validated the effectiveness of the method. The results of this paper provide a fundamental understandings of the dynamics of liquid free surfaces, in which the surface tension and contact angle boundary conditions are taken into account. Finally, numerical simulation of largescale amplitude sloshing of liquid in a cylindrical container is performed and a numerical analysis of the effect of an annular ring-shaped rigid damping baffle on liquid sloshing oscillations in a cylindrical tank is also carried out.

Numerical and experimental investigation of rock breaking method under free surface by TBM disc cutter

To study the rock breaking method under the free surface induced by disc cutter,the rock breaking simulations were first conducted based on the discrete element method,and the dynamic process of rock breaking under the free surface was studied including stressed zone,crush zone,crack initiation and propagation.Then the crack propagation conditions,specific energy,etc.under different free surface distance(S)were also investigated combined with linear cutting experiments.The results show that the rock breaking process under the free surface induced by disc cutter is dominated by tension failure mode.There exists a critical S to promote crack propagation to free surface effectively.And this rock breaking method can improve the rock breaking force and breaking efficiency significantly when proper.

Prolonged simulation of near-free surface underwater explosion based on Eulerian finite element method

n the area of naval architecture and ocean engineering,the research about the underwater xplosion problem is of great significance.To achieve prolonged simulation of near-free surface underwater explosion,the underwater explosion transient numerical model is established in this paper based on compressible Eulerian finite element method(EFEM).Compared with Geers Hunter formula,EFEM is availably validated by simulating the free-field underwater xplosion case.Then,the bubble pulsation and flow field dynamic characteristics of the cases with different underwater explosive depth are compared in this work.Lastly,the height of the water hump and the pressure of flow flied are analyzed quantitatively through the simulation results.

A CIP-Based Numerical Simulation of Free Surface Flow Related to Freak Waves

An enhanced numerical model for simulating two-dimensional incompressible viscous flow with distorted free surface is reported. The numerical simulation is carried out through the CIP （Constrained Interpolation Profile）-based method, which is described in the paper. A more accurate interface capturing scheme, the VOF/WLIC scheme （VOF：Volume-of-Fluid;WLIC：weighed line interface calculation）, is adopted as the interface capturing method. To assess the developed algorithm and its versatility, a selection of test problems are examined, i.e. the square wave propagation, the Zalesak’s rigid body rotation, dam breaking problem with and without obstacles, wave sloshing in an excited wave tank and interaction between extreme waves and a floating body. Excellent agreements are obtained when numerical results are compared with available analytical, experimental, and other numerical results. These examples demonstrate that the use of the VOF/WLIC scheme in the free surface capturing makes better results and also the proposed CIP-based model is capable of predicting the freak wave-related phenomena.

Numerical Simulation of Sloshing with Large Deforming Free Surface by MPS-LES Method

Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.

Model experiment of rock blasting with single borehole and double free-surface

In order to improve the quality of laneways and tunnel excavation by drilling and blasting and by making effective use of explosive energy, a model experiment of rock blasting with a single borehole and a double free-surface was performed with the objective of studying the effect of parameters such as charge structure, free-surface and rock compressive strength on rock blasting.The model experiments indicate that:1) the smaller the rock compressive strength and density, the more distinctive the cavity expanding action by blasting;2) the powder factor in an air-decoupling charge structure is larger than that in a coupling charge structure, i.e., the explosive energy in an air-decoupling charge structure transferred to the rocks is less than that in a coupling charge structure;3) a free-surface improves the utilizations of explosive energy;4) an air-decoupling charge structure helps to maintain the integrity and stability of wall rock in controlled perimeter blasting, such as in roadways and tunnel excavation by drilling and blasting.

Surface free energy of copper-zinc alloy for energy-saving of boiler

Cu-Zn,Cu-Zn-Sn,Cu-Zn-Ni alloys were melted by vacuum smelter.The effect factors to the surface free energy of the alloys such as chemical composition,crystal structure and surface crystal lattice distortion etc.were investigated by OCA30 automatic contact angle test instrument,metallography microscope and XRD instrument etc.Results suggests:adding alloy element to Cu may increase its surface free energy,and the more kinds of alloy elements are added,the more surface free energy increases;the alloy element Sn an increase the surface free energy of Cu-Zn alloy;Cu-Zn alloy with fir-tree crystal structure,great phase discrepancy and obvious composition aliquation has greater surface free energy;Cu-Zn alloy with compounds and serious surface crystal lattice distortion has greater surface free energy.

Three-dimensional numerical simulation of crown spike due to coupling effect between bubbles and free surface

The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble＇s effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.

Vertical 2D Modeling of Free Surface Flow with Hydrodynamic Pressure Using SIMPLE Arithmetic in σ Coordinates

A numerical model for shallow water flow has been developed based on the unsteady Reynolds-averaged Navier-Stokes equations with the hydrodynamic pressure instead of hydrostatic pressure assumption. The equations are transformed into the σ-coordinate system and the eddy viscosity is calculated with the standard k-ε turbulence model. The control volume method is used to discrete the equations, and the boundary conditions at the bed for shallow water models only include vertical diffusion terms expressed with wall functions. And the semi-implicit method for pressure linked equation arithmetic is adopted to solve the equations. The model is applied to the 2D vertical plane flow of a current over two steep-sided trenches for which experiment data are available for comparison and good agreement is obtained. And the model is used to predicting the flow in a channel with a steep-sided submerged breakwater at the bottom, and the streamline is drawn.

EXAMINATION OF THE SURFACE FREE ENERGY AND ACID-BASE PROPERTIES OF CELLULOSE BY THE COLUMN WICKING TECHNIQUE AND THE CRITICAL PACKING HEIGHT/DENSITY

The column wicking technique was applied to estimate the surface free energy of cellulose, the importance of which is to obtain a real effective capillary radius, Reff, initially from the plot of Washburn penetration distance versus time. Since the cellulose sample could not be packed with good reproducibility, therefore, Reff can not be obtained readily from the slope of the plot. A method was developed in this paper by uniting all apparent packing heights with a unique value to deduce a real effective capillary radius. Based on the defined critical packing height related to the critical packing density, the surface free energy and acid-base properties of cellulose Sigma C8002 were estimated.

Model for seawater fouling and effects of temperature,flow velocity and surface free energy on seawater fouling

A kinetic model was proposed to predict the seawater fouling process in the seawater heat exchangers.The new model adopted an expression combining depositional and removal behaviors for seawater fouling based on the Kern–Seaton model.The present model parameters include the integrated kinetic rate of deposition(k d)and the integrated kinetic rate of removal(k r),which have clear physical signi ficance.A seawater-fouling monitoring device was established to validate the model.The experimental data were well fitted to the model,and the parameters were obtained in different conditions.SEM and EDX analyses were performed after the experiments,and the results show that the main components of seawater fouling are magnesium hydroxide and aluminum hydroxide.The effects of surface temperature,flow velocity and surface free energy were assessed by the model and the experimental data.The results indicate that the seawater fouling becomes aggravated as the surface temperature increased in a certain range,and the seawater fouling resistance reduced as the flow velocity of seawater increased.Furthermore,the effect of the surface free energy of metals was analyzed,showing that the lower surface free energy mitigates the seawater fouling accumulation.