An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.

Appraisal of observational method for consolidation analysis

Two of the aims of laboratory one-dimensional consolidation tests are prediction of the end of primary settlement, and determination of the coefficient of consolidation of soils required for the time rate of consolidation analysis from time-compression data. Of the many methods documented in the literature to achieve these aims, Asaoka's method is a simple and useful tool, and yet the most neglected one since its inception in the geotechnical engineering literature more than three decades ago. This paper appraises Asaoka's method, originally proposed for the field prediction of ultimate settlement, from the perspective of laboratory consolidation analysis along with recent developments. It is shown through experimental illustrations that Asaoka's method is simpler than the conventional and popular methods, and makes a satisfactory prediction of both the end of primary compression and the coefficient of consolidation from laboratory one-dimensional consolidation test data.

NUMERICAL SIMULATIONS OF LIQUID JET BREAK UP IN A CROSSFLOW

Atomization is the process of disintegration of a liquid jet into ligaments and subsequently into smaller droplets. A liquid jet injected from a circular orifice into cross flow of air undergoes atomization primarily due to the interaction of the two phases rather than an intrinsic break up. Direct numerical simulation of this process resolving the finest droplets is computationally very expensive and impractical. In the present study, we resort to multiscale modelling to reduce the computational cost. The primary break up of the liquid jet is simulated using Gerris, an open source code, which employs Volume-of-Fluid (VOF) algorithm. The smallest droplets formed during primary atomization are modeled as Lagrangian particles. This one-way coupling approach is validated with the help of the simple test case of tracking a particle in a Taylor-Green vortex. The temporal evolution of the liquid jet forming the spray is captured and the flattening of the cylindrical liquid column prior to breakup is observed. The size distribution of the resultant droplets is presented at different distances downstream from the location of injection and their spatial evolution is analyzed.

Secondary breakup of a drop at moderate Weber numbers

We present volume of fluid based numerical simulations of secondary breakup of a drop with high density ratio (approx. 1000) and also perform experiments by injecting monodisperse water droplets in a continuous jet of air and capture the breakup regimes, namely, bag formation, bag-stamen, multibag and shear breakup, observed in the moderate Weber number range (20-120). We observe an interesting transition regime between bag and shear breakup for We = 80, in both simulations as well as experiments, where the formation of multiple lobes, is observed, instead of a single bag, which are connected to each other via thicker rim-like threads that hold them. We show that the transition from bag to shear breakup occurs owing to the rim dynamics which shows retraction under capillary forces at We = 80, whereas the rim is sheared away with flow at We = 120 thus resulting in a backward facing bag. The drop characteristics and timescales obtained in simulations are in good agreement with experiments. The drop size distribution after the breakup shows bimodal nature for the single-bag breakup mode and a unimodal nature following lognormal distribution for higher Weber numbers.

Cross-section fluctuations in open microwave billiards and quantum graphs: The counting-of-maxima method revisited

The fluctuations exhibited by the cross sections generated in a compound-nucleus reaction or, more generally, in a quantum-chaotic scattering process, when varying the excitation energy or another external parameter, are characterized by the width Gamma(corr) of the cross-section correlation function. Brink and Stephen Phys. Lett. 5, 77 (1963)] proposed a method for its determination by simply counting the number of maxima featured by the cross sections as a function of the parameter under consideration. They stated that the product of the average number of maxima per unit energy range and Gamma(corr) is constant in the Ercison region of strongly overlapping resonances. We use the analogy between the scattering formalism for compound-nucleus reactions and for microwave resonators to test this method experimentally with unprecedented accuracy using large data sets and propose an analytical description for the regions of isolated and overlapping resonances.

Application of enthalpy method for moving boundary problem in laser surface melting

A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization. The enthalpy method was applied to solve this two-phase axisymmetrical melting problem Computational results of temperature fields were obtained, which provide useful information to practical laser treatment processing. The validity of enthalpy method in solving such problems is presented.

Nonlinear Analysis of m Method for Single Pile Under Lateral Loading

通过对大量模型试验结果的分析,提出了简便合理的水平荷载作用下单桩的计算方法。该方法将常用m法中单一m值与位移建立指数关系,并由b_m及K两参数来描述。将这一关系引入常用单桩m法计算中,则计算结果可考虑单桩非线性响应的影响。

An improvement and proof of OGY method

OGY method is the most important method of controlling chaos. It stabilizes a hyperbolic periodic orbit by making small perturbations for a system parameter. This paper improves the method of choosing parameter, and gives a mathematics proof of it.

Preliminary Investigation of Fluid Flow in Meso-Scale Rectangular Tube

The fluid flow associated with micro and meso scale devices is currently of interest. Experiments were performed to study the fluid flow in meso-scale channels. A straight flow tube was fabricated with 1.0x4.0mm^2 in rectangular cross section and 200mm in length, which was made of quartz for flow visualization and PIV measurements. Reynolds numbers were ranged from 311 to over 3105. The corresponding pressure drop was from 0.65KPa to over 16.58KPa between the inlet and outlet of the tube. The micro PIV was developed to measure the velocity distribution in the tube. A set of microscope object lens was mounted ahead of CCD camera to obtain optimized optical magnification on the CCD chip. The velocity distributions near the outlet of the tube were measured to obtain full-developed flow. A CW laser beam was focused directly on the test section by a cylinder lens to form a small light sheet. Thus, high power density of light was formed on the view region. It is very important to the experiment while the velocity of the flow reaches to a few meters per second within millimeter scale. In this case, it is necessary to reduce exposure time to microseconds for PIV measurements. In the present paper, the experimental results are compared with the classical theories.

Progress In Modeling Of Fluid Flows In Crystal Growth Processes

Modeling of fluid flows in crystal growth processes has become an important research area in theoretical and applied mechanics. Most crystal growth processes involve fluid flows, such as flows in the melt, solution or vapor. Theoretical modeling has played an important role in developing technologies used for growing semiconductor crystals for high performance electronic and optoelectronic devices. The application of devices requires large diameter crystals with a high degree of crystallographic perfection, low defect density and uniform dopant distribution. In this article, the flow models developed in modeling of the crystal growth processes such as Czochralski, ammonothermal and physical vapor transport methods are reviewed. In the Czochralski growth modeling, the flow models for thermocapillary flow, turbulent flow and MHD flow have been developed. In the ammonothermal growth modeling, the buoyancy and porous media flow models have been developed based on a single-domain and continuum approach for the composite fluid-porous layer systems. In the physical vapor transport growth modeling, the Stefan flow model has been proposed based on the flow-kinetics theory for the vapor growth. In addition, perspectives for future studies on crystal growth modeling are proposed. (c) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Large Scale Numerical Simulations For Multi-Phase Fluid Dynamics With Moving Interfaces

This short communication presents our recent studies to implement numerical simulations for multi-phase flows on top-ranked supercomputer systems with distributed memory architecture. The numerical model is designed so as to make full use of the capacity of the hardware. Satisfactory scalability in terms of both the parallel speed-up rate and the size of the problem has been obtained on two high rank systems with massively parallel processors, the Earth Simulator (Earth simulator research center, Yokohama Kanagawa, Japan) and the TSUBAME (Tokyo Institute of Technology, Tokyo, Japan) supercomputers.

Experimental study of thermocapillary convection on a liquid bridge consisting of fluid with low Prandtl number in a floating half-zone

A device of mercury liquid bridge of floating half-zone is designed to experimentally explore thermocapillary convection and its instability of a low Prandtl number liquid. Noncontacted diagnostic techniques were developed to monitor surface flow and surface deformation. The surface flow and the influence of a growing surface film (or skin) on the flow were observed. It is shown that the film is a key factor in changing the behavior associated with the thermocapillary convection. The experiment indicates that the critical Marangoni number should be much higher than that expected by the numerical simulation. The condition and process of surface film growth are discussed. The surface oscillation of the mercury bridge wrapped with ''dirt-film'' was observed, and the characteristics and the frequency associated with this oscillation are given.

Experimental Study on Surface Deformation and Surface Wave of Fluid Convection with Free Surface

An optical diagnostic system consisting of Michelson interferometer with image processor has been developed for study of the kinetics of thermal capillary convection and buoyancy convection. This optical interferometer has been used to observe and measure surface deformation and surface wave of capillary convection and buoyancy convection in a rectangular cavity with different temperature’s sidewalls. Fourier transformation is used to image processing. The quantitative results of surface deformation and surface wave have been calculated from the interference fringe pattern. With the increasing of temperature gradient, the liquid surface slant gradually. It’s deformation has been calculated, which is related directly with temperature gradient. This is one of the characters introducing convection. Another interesting phenomenon is the inclining direction, which is different when the liquid layer is thin or thick. When the liquid layer is thin, convection is mainly controlled by thermocapillary effect. However, When the liquid layer is thick, convection is mainly controlled by buoyancy effect. Surface deformation in the present experiment are more and more declining in this process. The present experiment proved that surface deformation appears before the appearance of surface wave on fluid convection, it is related with temperature gradient, and the height of liquid layer, and lies on capillary convection and buoyancy convection. The present experiment also demonstrates that the amplitude of surface wave of thermocapillary-buoyancy convection is much smaller than surface deformation, the wave is covered by deformation.