An Improved Ce/Se Scheme For Multi-Material Elastic-Plastic Flows And Its Applications

In this paper, a new definition of SE and CE, which is based on the hexahedron mesh and simpler than Chang's original CE/SE method (the space-time Conservation Element and Solution Element method), is proposed and an improved CE/SE scheme is constructed. Furthermore, the improved CE/SE scheme is extended in order to solve the elastic-plastic flow problems. The hybrid particle level set method is used for tracing the interfaces of materials. Proper boundary conditions are presented in interface tracking. Two high-velocity impact problems are simulated numerically and the computational results are carefully compared with the experimental data, as well as the results from other literature and LS-DYNA software. The comparisons show that the computational scheme developed currently is clear in physical concept, easy to be implemented and high accurate and efficient for the problems considered. (C) 2008 Elsevier Ltd. All rights reserved.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

A CE/SE Scheme for Flows in Porous Media and Its Applications

In this paper, a new computational scheme for solving flows in porous media was proposed. The scheme was based on an improved CE/SE method (the space-time Conservation Element and Solution Element method). We described porous flows by adopting DFB (Brinkman-Forchheimer extended Darcy) equation. The comparison between our computational results and Ghia's confirmed the high accuracy, resolution, and efficiency of our CE/SE scheme. The proposed first-order CE/SE scheme is a new reliable way for numerical simulations of flows in porous media. After investigation of effects of Darcy number on porous flow, it shows that Darcy number has dominant influence on porous flow for the Reynolds number and porosity considered.

Kinetic studies of gas flows

Our recent studies on kinetic behaviors of gas flows are reviewed in this paper. These flows have a wide range of background, but share a common feature that the flow Knudsen number is larger than 0.01. Thus kinetic approaches such as the direct simulation Monte Carlo method are required for their description. In the past few years, we studied several micro/nano-scale flows by developing novel particle simulation approach, and investigated the flows in low-pressure chambers and at high altitude. In addition, the microscopic behaviors of a couple of classical flow problems were analyzed, which shows the potential for kinetic approaches to reveal the microscopic mechanism of gas flows.

Exact solutions to collisionless external gas flows

This paper presents exact density, velocity and temperature solutions for two problems of collisionless gas flows around a flat plate or a spherical object. At any point off the object, the local velocity distribution function consists of two pieces of Maxwellian distributions: one for the free stream which is characterized by free stream density, temperature and average velocity, n0, T0, U0; and the other is for the wall and it is characterized by density at wall and wall temperature, nw,Tw. Directly integrating the distribution functions leads to complex but exact flowfield solutions. To validate these solutions, we perform numerical simulations with the direct simulation Monte Carlo (DSMC) method. In general, the analytical and numerical results are virtually identical. The evaluation of these analytical solutions only requires less than one minute while the DSMC simulations require several days.

两层液体Bénard-Marangoni对流的实验研究

采用粒子图像测速技术(PIV)对两层流体B(en)ard-Marangoni对流进行了实验研究.研究了各种不同厚度比下的临界对流模式,同时研究了在温差变大时向超临界对流模式转化过程.实验结果表明,界面张力对各种对流模式的形成和转变具有重要的作用.

A Simple Model For Predicting The Void Fraction Of Gas/Non-Newtonian Fluid Intermittent Flows In Upward Inclined Pipes

In this work, a simple correlation, which incorporates the mixture velocity, drift velocity, and the correction factor of Farooqi and Richardson, was proposed to predict the void fraction of gas/non-Newtonian intermittent flow in upward inclined pipes. The correlation was based on 352 data points covering a wide range of flow rates for different CMC solutions at diverse angles. A good agreement was obtained between the predicted and experimental results. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase intermittent flows.

Construction of habitat suitability models (HSMs) for benthic macroinvertebrate and their applications to instream environmental flows: A case study in Xiangxi River of Three Gorges Reservior region, China

Based on a long-term ecological monitoring, the present study chose the most dominant benthic macroinvertebrate (Baetis spp.) as target organisms in Xiangxi River, built the habitat suitability models (HSMs) for water depth, current velocity and substrate, respectively, which is the first aquatic organisms model for habitat suitability in the Chinese Mainland with a long-term consecutive in situ measurement. In order to protect the biointegrity and function of the river ecosystem, the theory system of instream environmental flow should be categorized into three hierarchies, namely minimum required instream flow (hydrological level), minimum instream environmental flow (biospecies level), and optimum instream environmental flow (ecosystem level). These three hierarchies of instream environmental flow models were then constructed with the hydrological and weighted usable area (WUA) method. The results show that the minimum required instream flow of Xiangxi River calculated by the Tennant method (10% of the mean annual flow) was 0.615 m(3) s(-1); the minimum instream environmental flow accounted for 19.22% of the mean annual flow (namely 1.182 m(3) s(-1)), which was the damaged river channel. ow in the dry season; and 42.91% of the mean annual flow (namely 2.639 m(3) s(-1)) should be viewed as the optimum instream environmental flow in order to protect the health of the river ecosystem, maintain the instream biodiversity, and reduce the impact of small hydropower stations nearby the Xiangxi River. We recommend that the hydrological and biological methods can help establish better instream environmental. ow models and design best management practices for use in the small hydropower station project. (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.

Structural and optical properties of quaternary AlInGaN epilayers grown by MOCVD with various TMGa flows

AlInGaN quaternary epilayers have been grown with various TMGa flows by metalorganic chemical vapor deposition to investigate the influence of growth rate on the structural and optical properties. Triple-axis X-ray diffraction measurements show AlInGaN epilayers have good crystalline quality. Photolummescence (PL) measurements show that the emission intensity of AlInGaN epilayers is twenty times stronger than that of AlGaN epilayer with comparable Al content. V-shaped pits are observed at the surface of AlInGaN epilayers by atomic force microscopy (AFM) and transmission electron microscopy (TEM). High growth rate leads to increased density and size of V-shaped pits, but crystalline quality is not degraded. (C) 2003 Elsevier B.V. All rights reserved.

Experimental investigation on TBAB clathrate hydrate slurry flows in a horizontal tube: Forced convective heat transfer behaviors

Tetra-n-butyl-ammonium bromide (TBAB) clathrate hydrate slurry (CHS) is one kind of secondary refrigerants, which is promising to be applied into air-conditioning or latent-heat transportation systems as a thermal storage or cold carrying medium for energy saving. It is a solid-liquid two phase mixture which is easy to produce and has high latent heat and good fluidity. In this paper, the heat transfer characteristics of TBAB slurry were investigated in a horizontal stainless steel tube under different solid mass fractions and flow velocities with constant heat flux. One velocity region of weakened heat transfer was found. Moreover, TBAB CHS was treated as a kind of Bingham fluids, and the influences of the solid particles, flow velocity and types of flow on the forced convective heat transfer coefficients of TBAB CHS were investigated. At last, criterial correlations of Nusselt number for laminar and turbulent flows in the form of power function were summarized, and the error with experimental results was within 20%.

Micromixing enhanced by pulsating flows

We consider the micromixing enhancement by pulsating flows. Dimensionless governing equations and boundary conditions were developed for T-type micromixers with two inlet pulsating flows. The problem involves a set of parameters. Three key dimensionless parameters are identified: the Reynolds number, the Strouhal number, and the disturbance amplitude. Suitable Strouhal number or disturbance amplitude causes symmetrical meniscus-shape mixing interfaces, separating the whole mixing channel into a set of segments. Thus uniform exit species concentration can be reached. Too large or too small Strouhal number or disturbance amplitude yields the meniscus-shape mixing interfaces deviating from the centerline of the mixing channel, deteriorating the mixing performance. The optimized disturbance amplitude is increased with increases in Strouhal numbers. Low Reynolds number needs larger disturbance amplitude.

Statistical analysis of microscale gas flows and stock price changes

First, recent studies on the information preservation (IP) method, a particle approach for low-speed micro-scale gas flows, are reviewed. The IP method was validated for benchmark issues such as Couette, Poiseuille and Rayleigh flows, compared well with measured data for typical internal flows through micro-channels and external flows past micro flat plates, and combined with the Navier-Stokes equations to be a hybrid scheme for subsonic, rarefied gas flows. Second, the focus is moved to the microscopic characteristic of China stock market, particularly the price correlation between stock deals. A very interesting phenomenon was found that showed a reverse transition behaviour between two neighbouring price changes. This behaviour significantly differs from the transition rules for atomic and molecular energy levels, and it is very helpful to understand the essential difference between stock markets and nature.

Mechanism Responsible for the Complete Suppression of Karman Vortex in Flows Past a Wavy Square-Section Cylinder

The Karman vortex shedding is totally suppressed in flows past a wavy square-section cylinder at a Reynolds number of 100 and the wave steepness of 0.025. Such a phenomenon is illuminated by the numerical simulations. In the present study, the mechanism responsible for it is mainly attributed to the vertical vorticity. The geometric disturbance on the rear surface leads to the appearance of spanwise flow near the base. The specific vertical vorticity is generated on the rear surface and convecting into the near wake. The wake flow is recirculated with the appearance of the pair of recirculating cells. The interaction between the upper and lower shear layers is weakened by such cells, so that the vortex rolls could not be formed and the near wake flow becomes stable.