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.

Effect of substitution of chromium for nickel on structure and electrochemical characteristic of Ti0.26Zr0.07V0.24Mn0.1Ni0.33 multi-phase hydrogen storage alloy

Ti-Zr-V-Mn-Ni-based multi-component alloys demonstrate high discharge capacity in KOH electrolyte. However, the drastic decrease in their discharge capacities makes them unsuitable for use as negative electrode material in the Ni/MH battery. In present work, Ni is partially replaced by Cr in the Ti-Zr-V-Mn-Ni-based alloys to improve their cycle life. The effects of Cr substitution on microstructures and the electrochemical characteristics of the alloys are investigated. It is found that Cr substitution is very effective to improve the cyclic durability of the alloys although the discharge capacity decreases with changing x from 0.05 to 0.20. Some kinetic performances have been also investigated using electrochemical impedance spectroscopy (EIS) and potentiostatic discharge technique.

Application of Electrical Resistance Tomography System to Monitor Gas/Liquid Two-Phase Flow in a Horizontal Pipe

The Electrical Resistance Tomography (ERT) technique possesses great potential in monitoring widely exiting industrial two/multi-phase flow. For vertical pipe flow and inclined pipe flow, some application studies with exciting results have been reported, but there is rarely a paper regarding the application of ERT to horizontal gas/liquid pipe flow. This paper addresses this issue and proposes a smart method, Liquid Level Detection method, to conventional ERT system. The enhanced ERT system using the new method can monitor horizontal pipe flow effectively and its application is no longer restricted by the flow conditions. Some experimental results from monitoring an air/water slug pipe flow are presented.

Numerical Investigation on Laser Transformation Hardening with Different Temporal Pulse Shapes

A 3-D numerical model for pulsed laser transformation hardening (LTH) is developed using the finite element method. In this model, laser spatial and temporal intensity distribution, temperature-dependent thermophysical properties of material, and multi-phase transformations are considered. The influence of laser temporal pulse shape on connectivity of hardened zone, maximum surface temperature of material and hardening depth is numerically investigated at different pulse energy levels. Results indicate that these hardening parameters are strongly dependent on the temporal pulse shape. For the rectangular temporal pulse shape, the temperature field obtained from this model is in excellent agreement with analytical solution, and the predicted hardening depth is favorably compared with experimental one. It should be pointed out that appropriate temporal pulse shape should be selected according to pulse energy level in order to achieve desirable hardening quality under certain laser spatial intensity distribution.

Numerical Simulation of Pulsed Laser Transformation Hardening

A novel pulsed laser surface processing technology is introduced, which can make use of the spatial and temporal profile of laser pulse to obtain ideal hardening parameters. The intensity distribution of laser pulse is spatially and temporally controlled by using laser shape transformation technology. A 3D numerical model including multi-phase transformations is established to explore material microstructure evolution induced by temperature field evolution. The influences of laser spatial-temporal profiles on hardening parameters are investigated. Different from the continuous laser processing technology, results indicate that spatial and temporal profiles are important factors in determining processing quality during pulsed laser processing method.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

超图多级划分算法框架及对划分结果的多阶段优化

超图划分应用于大规模矩阵计算、大规模集成电路等领域.详细地阐述了超图多级划分的算法框架,并提出对划分结果进行优化的一种手段,通过进行多阶段的V循环优化,在可以接受的运行时间内得到对超图的一个较优的划分.

超图多级划分算法框架及对划分结果的多阶段优化

超图划分应用于大规模矩阵计算、大规模集成电路等领域.详细地阐述了超图多级划分的算法框架,并提出对划分结果进行优化的一种手段,通过进行多阶段的循环优化,在可以接受的运行时间内得到对超图的一个较优的划分.

Boundary Characterization Experiment Series Overview

Ocean acoustic propagation and reverberation in continental shelf regions is often controlled by the seabed and sea surface boundaries. A series of three multi-national and multi-disciplinary experiments was conducted between 2000-2002 to identify and measure key ocean boundary characteristics. The frequency range of interest was nominally 500-5000 Hz with the main focus on the seabed, which is generally considered as the boundary of greatest importance and least understood. Two of the experiments were conducted in the Mediterranean in the Strait of Sicily and one experiment in the North Atlantic with sites on the outer New Jersey Shelf (STRATAFORM area) and on the Scotian Shelf. Measurements included seabed reflection, seabed, surface, and biologic scattering, propagation, reverberation, and ambient noise along with supporting oceanographic, geologic, and geophysical data. This paper is primarily intended to provide an overview of the experiments and the strategies that linked the various measurements together, with detailed experiment results contained in various papers in this volume and other sources

Computer Simulation Tb(lll) Speciation in Human Blood Plasma

A multi phase model of human blood plasma was developed and the Tb(Ⅲ) speciation in this system was studied. The results show that the speciation of Tb(Ⅲ) depends on the concentration of Tb(Ⅲ). When the concentration of Tb(Ⅲ) is below 4.000×10 -8 mol/L, most of Tb(Ⅲ) exists as soluble species while the concentration of Tb(Ⅲ) is in between 4.000 ×10 -8 mol/L and 1.667×10 -2 mol/L, precipitates(TbPO 4 and Tb 2 (CO 3 ) 3 ) are the dominant species of Tb(Ⅲ). Among soluble Tb(Ⅲ) ...

Computer simulation of insoluble Pr(III) speciation in human interstitial fluid

A multi-phase model of Pr(III) speciation in human interstitial fluid was constructed and insoluble Pr(III) speciation was studied. When the total concentration of Pr(III) is below 8.401E-10 mol/L, soluble Pr(III) species are main species. With rising the total concentration of Pr(III), Pr(III) is firstly bound to phosphate to form precipitate of PrPO4, then bound to carbonate and another precipitate of Pr-2(CO3)(3) was obtained. When the total concentration is between 1.583E-9 mol/L and 4.000E-3 mol/L, the insoluble species are predominant Pr(III) species.

Tb(III) speciation in human blood plasma by computer simulation

A multi-phase model was developed and Tb(III) speciation in human blood plasma was studied. At a concentration below 3.744x 10(-4) mol/L (or at the concentration), Tb(III) is mostly bound to phosphate to form precipitate of TbPO4. As the concentration of Tb(III) increases, phosphate is exceeded and another kind of precipitate of Tb-2(CO3)(3) appears. Among soluble Tb(III) species, Tb(III) mainly distribute in [Tb (Tf)] at low concentration and in [Tb (HSAA, [Tb-2 (Tf)], [Th (IgG)], [Tb (Lactate)](2+), [Tb (CitArgH)] and free Tb(III) at high concentration.

复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.