simulation and interaction of fluid dynamics

In the fluid simulation, the fluids and their surroundings may greatly change properties such as shape and temperature simultaneously, and different surroundings would characterize different interactions, which would change the shape and motion of the fluids in different ways. On the other hand, interactions among fluid mixtures of different kinds would generate more comprehensive behavior. To investigate the interaction behavior in physically based simulation of fluids, it is of importance to build physically correct models to represent the varying interactions between fluids and the environments, as well as interactions among the mixtures. In this paper, we will make a simple review of the interactions, and focus on those most interesting to us, and model them with various physical solutions. In particular, more detail will be given on the simulation of miscible and immiscible binary mixtures. In some of the methods, it is advantageous to be taken with the graphics processing unit (GPU) to achieve real-time computation for middle-scale simulation.

Gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, Northern South China Sea

Based on the comprehensive interpretation and study of the Neogene fracture system and diapiric structure, it can be concluded that the diapiric structures, high-angle fractures and vertical fissure system are the main gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, northern South China Sea. The Neogene fractures widely developed in the study area may be classed into two groups: NW (NNW)-trending and NE (NNE)-trending. The first group was active in the Late Miocene, while the second one was active since the Pliocene. The NE (NNE)-trending fractures were characterized by lower activity strength and larger scale, and cut through the sediment layers deposited since the Pliocene. Within the top sediment layers, the high-angle fracture and vertical fissure system was developed. The diapiric structures display various types such as a turtle-back-like arch, weak piercing, gas chimney, and fracture (or crack, fissure). On the seismic profile, some diapiric structures show the vertical chimney pathway whose top is narrow and the bottom is wide, where some ones extend horizontally into pocket or flower-shaped structures and formed the seismic reflection chaotic zones. Within the overlying sediment layers of the diapiric structure, the tree branch, flower-shaped high-angle fractures and vertical fissures were developed and became the pathway and migration system of the gas-bearing fluid influx. In the study area, the diapiric structures indicate a high temperature/over pressure system ever developed. Closely associated and abundant bright-spots show the methane-bearing fluid influx migrated vertically or horizontally through the diapiric structures, high-angle fractures and vertical fissures. In the place where the temperature and pressure conditions were favor for the formation of gas hydrate, the hydrate reservoir deposition sub-system was developed.

Uniform mems chip temperatures in the nucleate boiling heat transfer region by selecting suitable, medium boiling number range

The not only lower but also uniform MEMS chip temperatures can he reached by selecting suitable boiling number range that ensures the nucleate boiling heat transfer. In this article, boiling heat transfer experiments in 10 silicon triangular microchannels with the hydraulic diameter of 55.4 mu m were performed using acetone as the working fluid, having the inlet liquid temperatures of 24-40 degrees C, mass fluxes of 96-360 kg/m(2)s, heat fluxes of 140-420 kW/m(2), and exit vapor mass qualities of 0.28-0.70. The above data range correspond to the boiling number from 1.574 x 10(-3) to 3.219 x 10(-3) and ensure the perfect nucleate boiling heat transfer region, providing a very uniform chip temperature distribution in both streamline and transverse directions. The boiling heat transfer coefficients determined by the infrared radiator image system were found to he dependent on the heat Axes only, not dependent on the mass Axes and the vapor mass qualities covering the above data range. The high-speed flow visualization shows that the periodic flow patterns take place inside the microchannel in the time scale of milliseconds, consisting of liquid refilling stage, bubble nucleation, growth and coalescence stage, and transient liquid film evaporation stage in a full cycle. The paired or triplet bubble nucleation sites can occur in the microchannel corners anywhere along the flow direction, accounting for the nucleate boiling heat transfer mode. The periodic boiling process is similar to a series of bubble nucleation, growth, and departure followed by the liquid refilling in a single cavity for the pool boiling situation. The chip temperature difference across the whole two-phase area is found to he small in a couple of degrees, providing a better thermal management scheme for the high heat flux electronic components. Chen's [11 widely accepted correlation for macrochannels and Bao et al.'s [21 correlation obtained in a copper capillary tube with the inside diameter of 1.95 mm using R11 and HCFC123 as working fluids can predict the present experimental data with accepted accuracy. Other correlations fail to predict the correct heat transfer coefficient trends. New heat transfer correlations are also recommended.

Investigation of Interfacial Tensions for Carbon Dioxide Aqueous Solutions by Perturbed-Chain Statistical Associating Fluid Theory Combined with Density-Gradient Theory

The perturbed-chain statistical associating fluid theory and density-gradient theory are used to construct an equation of state (EOS) applicable for the phase behaviors of carbon dioxide aqueous solutions. With the molecular parameters and influence parameters respectively regressed from bulk properties and surface tensions of pure fluids as input, both the bulk and interfacial properties of carbon dioxide aqueous solutions are satisfactorily correlated by adjusting the binary interaction parameter (k(ij)). Our results show that the constructed EOS is able to describe the interfacial properties of carbon dioxide aqueous solutions in a wide temperature range, and illustrate the influences of temperature, pressure, and densities in each phase on the interfacial properties.

Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.

Are the available boiling heat transfer coefficients suitable for silicon microchannel heat sinks?

The typical MEMS fabrication of micro evaporators ensures the perfect smooth wall surface that is lack of nucleation sites, significantly decreasing the heat transfer coefficients compared with miniature evaporators fabricated using copper or stainless steel. In the present paper, we performed the boiling heat transfer experiment in silicon triangular microchannel heat sink over a wide parameter range for 102 runs. Acetone was used as the working fluid. The measured boiling heat transfer coefficients versus the local vapor mass qualities are compared with the classical Chen’s correlation and other correlations for macro and miniature capillary tubes. It is found that most of these correlations significantly over-predict the measured heat transfer coefficients. New correlations are given. There are many reasons for such deviations. The major reason is coming from the perfect smooth silicon surface that lowers the heat transfer performances. New theory is recommended for the silicon microchannel heat sink that should be different from metallic capillary tubes.

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.

A NEW TYPE OF SILICON-ON-INSULATOR WITH A PERFECT SURFACE SILICON LAYER

A novel silicon structure consisting of a silicon-on-defect layer (SODL), with enhanced surface Hall mobility in the surface layer on a buried defect layer (DL), has been discovered [J. Li, Nucl. Instr. and Meth. B59/60 (1991) 1053]. SODL material was formed by using proton implantation and subsequent two-step annealing. The implantation was carried out with a Varian 350D ion implanter. Based on the discovery, a standard measurement method (current-voltage curve method) was adopted to measure the true resistivity value of the DL in order to replace the spreading resistivity measurement by which the true resistivity in seriously defective silicon cannot be obtained. By adopting the current-voltage current method, the true resistivity value of the DL is measured to be 4.2 x 10(9) OMEGA cm. The SODL material was proved to be a silicon-on-insulator substrate.

Quasi-3D Cytoskeletal Dynamics of Osteocytes under Fluid Flow

Osteocytes respond to dynamic fluid shear loading by activating various biochemical pathways, mediating a dynamic process of bone formation and resorption. Whole-cell deformation and regional deformation of the cytoskeleton may be able to directly regulate this process. Attempts to image cellular deformation by conventional microscopy techniques have been hindered by low temporal or spatial resolution. In this study, we developed a quasi-three-dimensional microscopy technique that enabled us to simultaneously visualize an osteocyte's traditional bottom-view profile and a side-view profile at high temporal resolution. Quantitative analysis of the plasma membrane and either the intracellular actin or microtubule (MT) cytoskeletal networks provided characterization of their deformations over time. Although no volumetric dilatation of the whole cell was observed under flow, both the actin and MT networks experienced primarily tensile strains in all measured strain components. Regional heterogeneity in the strain field of normal strains was observed in the actin networks, especially in the leading edge to flow, but not in the MT networks. In contrast, side-view shear strains exhibited similar subcellular distribution patterns in both networks. Disruption of MT networks caused actin normal strains to decrease, whereas actin disruption had little effect on the MT network strains, highlighting the networks' mechanical interactions in osteocytes.