Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

Wave loading on floating platforms by internal solitary waves

Morison's equation is used for estimating internal solitary wave-induced forces exerted on SPAR and semi-submersible platforms. And the results we got have also been compared to ocean surface wave loading. It is shown that Morison's equation is an appropriate approach to estimate internal wave loading even for SPAR and semi-submersible platforms, and the internal solitary wave load on floating platforms is comparable to surface wave counterpart. Moreover, the effects of the layers with different thickness on internal solitary wave force are investigated.

Characteristics of flow fields induced by interfacial waves in two-layer fluid

When designing deep ocean structures, it is necessary to estimate the effects of internal waves on the platform and auxiliary parts such as tension leg, riser and mooring lines. Up to now, only a few studies are concerned with the internal wave velocity fields. By using the most representative two-layer model, we have analyzed the behavior of velocity field induced by interfacial wave in the present paper. We find that there may exist velocity shear of fluid particles in the upper and lower layers so that any structures in the ocean are subjected to shear force nearby the interface. In the meantime, the magnitude of velocity for long internal wave appears spatially uniform in the respective layer although they still decay exponentially. Finally, the temporal variation for Stokes and solitary waves are shown to be of periodical and pulse type.

An algorithm for coarse particle sedimentation simulation by Stokesian dynamics

Coarse Particle sedimentation is studied by using an algorithm with no adjustable parameters based on stokesian dynamics. Only inter-particle interactions of hydrodynamic force and gravity are considered. The sedimentation of a simple cubic array of spheres is used to verify the computational results. The scaling and parallelism with OpenMP of the method are presented. Random suspension sedimentation is investigated with Mont Carlo simulation. The computational results are shown in good agreement with experimental fitting at the lower computational cost of O(N In N).

Effects of rainfall infiltration on the stability of soil slopes

Slope failure due to rainfall is a common geotechnical problem. The mechanics of rainfall induced landslides involves the interaction of a number of complex hydrologic and geotechnical factors. This study attempts to identify the influence of some of these factors on the stability of soil slope including rainfall intensity, hydraulic conductivity and the strength parameters of soil.

Comparisons of static, quasi-static and dynamic 3D porous media scale network models for two-phase immiscible flow in porous media

A dynamic 3D pore-scale network model is formulated for investigating the effect of interfacial tension and oil-water viscosity during chemical flooding. The model takes into account both viscous and capillary forces in analyzing the impact of chemical properties on flow behavior or displacement configuration, while the static model with conventional invasion percolation algorithm incorporates the capillary pressure only. From comparisons of simulation results from these models. it indicates that the static pore scale network model can be used successfully when the capillary number is low. With the capillary increases due to the enhancement of water viscosity or decrease of interfacial tension, only the quasi-static and dynamic model can give insight into the displacement mechanisms.

Lateral motion and departure of vapor bubbles in nucleate pool boiling on thin wires in microgravity

A space experiment on bubble behavior and heat transfer in subcooled pool boiling phenomenon has been performed utilizing the temperature-controlled pool boiling (TCPB) device both in normal gravity in the laboratory and in microgravity aboard the 22(nd) Chinese recoverable satellite. The fluid is R113 at 0.1 MPa and subcooled by 26 degrees C nominally. A thin platinum wire of 60 mu m in diameter and 30mm in length is simultaneously used as heater and thermometer. Only the lateral motion and the departure of discrete vapor bubbles in nucleate pool boiling are reported and analyzed in the present paper. A scale analysis on the Marangoni convection surrounding a bubble in the process of subcooled nucleate pool boiling leads to formulas of the characteristic velocity of the lateral motion and its observability. The predictions consist with the experimental observations. Considering the Marangoni effect, a new qualitative model is proposed to reveal the mechanism underlying the bubble departure processes and a quantitative agreement can also be acquired.

Effects of temperature on the threshold of phosphorus for algal blooms

Algal blooms, worsening marine ecosystems and causing great economic loss, have been paid much attention to for a long time. Such environmental factors as light penetration, water temperature, and nutrient concentration are crucial in blooms processes. Among them, only nutrients can be controlled. Therefore, the threshold of nutrients for algal blooms is of great concern. To begin with, a dynamic eutrophication model has been constructed to simulate the algal growth and phosphorus cycling. The model encapsulates the essential biological processes of algal growth and decay, and phosphorus regeneration due to algal decay. The nutrient limitation is based upon commonly used Monod's kinetics. The effects of temperature and phosphorus limitation are particularly addressed. Then, we have endeavored to elucidate the threshold of phosphorus at different temperature for algal blooms. Based on the numerical simulation, the isoquant contours of change rate of alga as shown in the figure are obtained, which obviously demonstrate the threshold of nutrient at an arbitrary reasonable temperature. The larger the change rate is, the more rapidly the alga grows. If the phosphorus concentration at a given temperature remains larger than the threshold the algal biomass may increase monotonically, leading to the algal blooming. With the rising of temperature, the threshold is apparently reduced, which may explain why likely red tide disasters occur in a fine summer day. So, high temperature and sufficient phosphorus supply are the major factors which result in algal growth and blowout of red tide.

Piezoelectricity of ZnO films prepared by sol-gel method

ZnO piezoelectric thin films were prepared on crystal substrate Si(111) by sol-gel technology, then characterized by scanning electron microscopy, X-ray diffraction and atomic force microscopy (AFM). The ZnO films characterized by X-ray diffraction are highly oriented in (002) direction with the growing of the film thickness. The morphologies, roughness and grain size of ZnO film investigated by AFM show that roughness and grain size of ZnO piezoelectric films decrease with the increase of the film thickness. The roughness dimension is 2.188-0.914 nm. The piezoelectric coefficient d(33) was investigated with a piezo-response force microscope (PFM). The results show that the piezoelectric coefficient increases with the increase of thickness and (002) orientation. When the force reference is close to surface roughness of the films, the piezoelectric coefficient measured is inaccurate and fluctuates in a large range, but when the force reference is big, the piezoelectric coefficient d(33) changes little and ultimately keeps constant at a low frequency.

Preparation and optical properties of CdTe/CdO center dot nH(2)O core/shell nano-composites in aqueous solution

The deposition of CdO center dot nH(2)O On CdTe nanoparticles was studied in an aqueous phase. The CdTe nanocrystals (NCs) were prepared in aqueous solution through the reaction between Cd2+ and NaHTe in the presence of thioglycolic acid as a stabilizer. The molar ratio of the Cd2+ to Te2- in the precursory solution played an important role in the photoluminescence of the ultimate CdTe NCs. The strongest photoluminescence was obtained under 4.0 of [Cd2+]/[Te2-] at pH similar to 8.2. With the optimum dosage of Cd(II) hydrous oxide deposited on the CdTe NCs, the photoluminescence was enhanced greatly. The photoluminescence of these nanocomposites was kept constant in the pH range of 8.0-10.0, but dramatically decreased with an obvious blue-shifted peak while the pH was below 8.0. In addition, the photochemical oxidation of CdTe NCs with cadmium hydrous oxide deposition was markedly inhibited.

Space experimental studies of microgravity fluid science in China

Microgravity fluid physics is an important part of microgravity sciences, which consists of simple fluids of many new systems, gas-liquid two-phase flow and heat transfer, and complex fluid mechanics. In addition to the importance of itself in sciences and applications, microgravity fluid physics closely relates to microgravity combustion, space biotechnology and space materials science, and promotes the developments of interdisciplinary fields. Many space microgravity experiments have been per- formed on board the recoverable satellites and space ships of China and pushed the rapid development of microgravity sciences in China. In the present paper, space experimental studies and the main re- sults of the microgravity fluid science in China in the last 10 years or so are introduced briefly.

Measurements Of In Situ Stress And Mining-Induced Stress In Beiminghe Iron Mine Of China

In order to obtain the distribution rules of in situ stress and mining-induced stress of Beiminghe Iron Mine, the stress relief method by overcoring was used to measure the in situ stress, and the MC type bore-hole stress gauge was adopted to measure the mining-induced stress. In the in situ stress measuring, the technique of improved hollow inclusion cells was adopted, which can realize complete temperature compensation. Based on the measuring results, the distribution model of in situ stress was established and analyzed. The in situ stress measuring result shows that the maximum horizontal stress is 1.75-2.45 times of vertical stress and almost 1.83 times of the minimum horizontal stress in this mineral field. And the mining-induced stress measuring result shows that, according to the magnitude of front abutment pressure the stress region can be separated into stress-relaxed area, stress-concentrated area and initial stress area. At the -50 m mining level of this mine, the range of stress-relaxed area is 0-3 m before mining face; the range of stress-concentrated area is 3-55 m before mining face, and the maximum mining-induced stress is 16.5-17.5 MPa, which is 15-20 m from the mining face. The coefficient of stress concentration is 1.85.

MODELING EVOLUTION OF INTERNAL WAVES IN SOUTH CHINA SEA

Internal waves are an important factor in the design of drill operations and production in deep water, because the waves have very large amplitude and may induce large horizontal velocity. How the internal waves occur and propagate over benthal terrain is of great concern for ocean engineers. In the present paper, we have formulated a mathematical model of internal wave propagation in a two-layer deep water, which involves the effects of friction, dissipation and shoaling, and is capable of manifesting the variation of the amplitude and the velocity pattern. After calibration by field data measured at the Continental Slope in the Northern South China Sea, we have applied the model to the South China Sea, investigating the westward propagation of internal waves from the Luzon Strait, where internal waves originate due to the interaction of benthal ridge and tides. We find that the internal wave induced velocity profile is obviously characterized by the opposite flow below and above the pycnocline, which results in a strong shear, threatening safety of ocean structures, such as mooring system of oil platform, risers, etc. When internal waves propagate westwards, the amplitude attenuates due to the effects of friction and dissipation. The preliminary results show that the amplitude is likely to become half of its initial value at Luzon Strait when the internal waves propagate about 400 kilometers westwards.

PIXE and ICP-AES analysis of early glass unearthed from Xinjiang (China)

Early glasses (about 1066 BC-220 AD) unearthed from Xinjiang of China were chemically characterized by using PIXE and ICP-AES. It was found that these glasses were basically attributed to PbO-BaO-SiO2 system, K2O-SiO2 system, Na2O-CaO-SiO2 system and Na2O-CaO-PbO-SiO2 system. The results from the cluster analysis showed that some glasses had basically similar recipe and technology. The PbO-BaO-SiO2 glass and the K2O-SiO2 glass were thought to come from the central area and the south of ancient China, respectively. The part of the Na2O-CaO-SiO2 glass (including the Na2O-CaO-PbO-SiO2 glass) might be imported from Mesopotamia, while the other part might be locally produced. (c) 2005 Elsevier B.V. All rights reserved.