Perturbational Finite Volume Method For The Solution of 2-D Navier-Stokes Equations on Unstructured and Structured Colocated Meshes

Based on the first-order upwind and second-order central type of finite volume( UFV and CFV) scheme, upwind and central type of perturbation finite volume ( UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from 0.3 to 0. 8 and convergence perform excellent with Reynolds number variation from 102 to 104.

Steady-state crack growth and fracture work based on the theory of mechanism-based strain gradient plasticity

Mode I steady-state crack growth is analyzed under plane strain conditions in small scale yielding. The elastic-plastic solid is characterized by the mechanism-based strain gradient (MSG) plasticity theory [J. Mech. Phys. Solids 47 (1999) 1239, J. Mech. Phys. Solids 48 (2000) 99]. The distributions of the normal separation stress and the effective stress along the plane ahead of the crack tip are computed using a special finite element method based on the steady-state fundamental relations and the MSG flow theory. The results show that during the steady-state crack growth, the normal separation stress on the plane ahead of the crack tip can achieve considerably high value within the MSG strain gradient sensitive zone. The results also show that the crack tip fields are insensitive to the cell size parameter in the MSG theory. Moreover, in the present research, the steady-state fracture toughness is computed by adopting the embedded process zone (EPZ) model. The results display that the steady-state fracture toughness strongly depends on the separation strength parameter of the EPZ model and the length scale parameter in the MSG theory. Furthermore, in order for the results of steady crack growth to be comparable, an approximate relation between the length scale parameters in the MSG theory and in the Fleck-Hutchinson strain gradient plasticity theory is obtained.

Interface Layer Effect On The Stress Distribution Of A Wafer-Bonded Bilayer Structure

The interface layer plays an important role in stress transfer in composite structures. However, many interface layer properties such as the modulus, thickness, and uniformity are difficult to determine. The model developed in this article links the influence of the interface layer on the normal stress distribution along the layer thickness with the layer surface morphology before bonding. By doing so, a new method of determining the interfacial parameter(s) is suggested. The effects of the layer thickness and the surface roughness before bonding on the normal stress distribution and its depth profile are also discussed. For ideal interface case with no interfacial shear stress, the normal stress distribution pattern can only be monotonically decreased from the interface. Due to the presence of interfacial shear stress, the normal stress distribution is much more complex, and varies dramatically with changes in the properties of the interface layer, or the dimensions of the bonding layers. The consequence of this dramatic stress field change, such as the shift of the maximum stress from the interface is also addressed. The size-dependent stress distribution in the thickness direction due to the interface layer effect is presented. When the interfacial shear stress is reduced to zero, the model presented in this article is also demonstrated to have the same normal stress distribution as obtained by the previous model, which does not consider the interface layer effect.

Fracture statistics of brittle materials: Weibull or normal distribution

The fit of fracture strength data of brittle materials (Si3N4, SiC, and ZnO) to the Weibull and normal distributions is compared in terms of the Akaike information criterion. For Si3N4, the Weibull distribution fits the data better than the normal distribution, but for ZnO the result is just the opposite. In the case of SiC, the difference is not large enough to make a clear distinction between the two distributions. There is not sufficient evidence to show that the Weibull distribution is always preferred to other distributions, and the uncritical use of the Weibull distribution for strength data is questioned.

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

The plastic instability behaviour of laser-textured steel sheet

As the production of a new technique that can offer both good formability and high image clarity for texturing metal sheet, laser-textured sheet has attracted the attention of many manufacturers and users. Among the many subjects to be studied, plastic instability behaviour of the laser-textured sheet is one of most important to understand its ability in extending material ductility and to appropriately control this technique. Experimental investigations are carried out in this paper to study the macroscopic behaviour and microstructural mechanism of the laser-textured sheet, and comparison is made with the normal sheet taken from the same coil of metal sheet. It is demonstrated that, the difference in the behaviour of plastic instability obviously shows tendency to delay strain localization and the onset of thickness necking. Shear banding and internal void damage are spread to a much wider region in the sheet being laser-textured. The prestrained microcraters enforced on the surface of the textured sheet act as hardening spots, which are likely to share out deformation and inhibit the increasing rate of voiding, and eventually favouring the ductility of the material used.

Computation of turbulent-flow in a square duct - aspects of the secondary flow and its origin

A numerical study of turbulent flow in a straight duct of square cross-section is made. An order-of-magnitude analysis of the 3-D, time-averaged Navier-Stokes equations resulted in a parabolic form of the Navier-Stokes equations. The governing equations, expressed in terms of a new vector-potential formulation, are expanded as a multi-deck structure with each deck characterized by its dominant physical forces. The resulting equations are solved using a finite-element approach with a bicubic element representation on each cross-sectional plane. The numerical integration along the streamwise direction is carried out with finite-difference approximations until a fully-developed state is reached. The computed results agree well with other numerical studies and compare very favorably with the available experimental data. One important outcome of the current investigation is the interpretation analytically that the driving force of the secondary flow in a square duct comes mainly from the second-order terms of the difference in the gradients of the normal and transverse Reynolds stresses in the axial vorticity equation.

High-speed streams from coronal holes and the accelerating mechanism of the solar-wind

In this paper, the general Mach number equation is derived, and the influence of typical energy forms in the solar wind is analysed in detail. It shows that the accelerating process of the solar wind is influenced critically by the form of heating in the corona, and that the transonic mechanism is mainly the result of the adjustment of the variation of the crosssection of flowing tubes and the heat source term.The accelerating mechanism for both the high-speed stream from the coronal hole and the normal solar wind is similar. But, the temperature is low in the lower level of the coronal hole and more heat energy supply in the outside is required, hence the high speed of the solar wind; while the case with the ordinary coronal region is just the opposite, and the velocity of the solar wind is therefore lower. The accelerating process for various typical parameters is calculated, and it is found that the high-speed stream may reach 800 km/sec.

Generalization of the JTZ model to open plane wakes

The JTZ model [C. Jung, T. T¶el and E. Ziemniak, Chaos 3, (1993) 555], as a theoretical model of a plane wake behind a circular cylinder in a narrow channel at a moderate Reynolds number, has previously been employed to analyze phenomena of chaotic scattering. It is ex- tended here to describe an open plane wake without the con¯ned nar- row channel by incorporating a double row of shedding vortices into the intermediate and far wake. The extended JTZ model is found in qualitative agreement with both direct numerical simulations and ex- perimental results in describing streamlines and vorticity contours. To further validate its applications to particle transport processes, the in- teraction between small spherical particles and vortices in an extended JTZ model °ow is studied. It is shown that the particle size has signif- icant in°uences on the features of particle trajectories, which have two characteristic patterns: one is rotating around the vortex centers and the other accumulating in the exterior of vortices. Numerical results based on the extended JTZ model are found in qualitative agreement with experimental ones in the normal range of particle sizes.

Simulated pathogenic conformational switch regions matched well with the biochemical findings

Pathogenic conformational conversion is a general causation of many disease, such as transmissible spon- giform encephalopathy (TSE) caused by misfolding of prion, sickle cell anemia, and etc. In such structural changes, misfolding occurs in regions important for the stability of native structure firstly. This destabi- lizes the normal conformation and leads to subsequent errors in folding pathway. Sites involved in the first stage can be deemed switch regions of the protein, and are vital for conformational conversion. Namely it could be a switch of disease at residue level. Here we report an algorithm that can identify such sites computationally with an accuracy of 93%, by calculating the probability of the native structure of a short segment jumping to a mistake one. Knowledge of such switch sites could be used to target clinical therapy, study physiological and pathologic mechanism of protein, and etc.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

Effect of bismuth oxide on the thermal stability and Judd-Ofelt parameters of Er 3(+)/Yb3(+) co-doped aluminophosphate glasses

The Er3+/Yb3+ co-doped glasses with compositions of xBi(2)O(3)-(65-x)P2O5-4Yb(2)O(3)-11Al(2)O(3)-5BaO-15Na(2)O (where x = 0, 2.5, 5, 7.5 and 10 mol%) were prepared using the normal melt quench technique. The optical absorption spectra of the glasses were recorded in the wavelength range 300-1700 nm. The effect of Bi2O3 content on the thermal stability and absorption spectra of glasses was investigated. In addition, the Judd-Ofelt parameters and oscillator strengths were calculated by employing Judd-Ofelt theory. It was observed that the positions of the fundamental absorption edge and cut-off wavelength shifted towards red as the content of Bi2O3 increased. However, there were no red shifts found both in the peak wavelength and in the center of mass wavelength of all absorption bands with Bi2O3 content increasing. The results of Judd-Ofelt theory analysis showed that Judd-Ofelt parameters Omega(t), (t = 2, 4, 6) changed sharply when Bi2O3 concentration exceeded 5 mol%. The variation trends of experimental oscillator strength were similar with those of Judd-Ofelt parameters as function of Bi2O3 concentrations. Moreover, differential scanning calorimetry experiments showed that the increases of Bi2O3 content weakened the network structure and then lowered the thermal stability of the glasses. The spontaneous emission probability A(rad), branching ratio beta and the radiative lifetime tau(rad) were also calculated and analyzed. The stimulated emission cross-section of Er3+ was calculated according to the McCumber theory. It was found that the stimulated emission cross-section of Er3+ was monotonically increases with Bi2O3 content increasing. (C) 2006 Elsevier B.V. All rights reserved.

The computer simulation of primate (Macaca mulatta) spatial visual delayed response and cerebral prefrontal control function

Cerebral prefrontal function is one of the important aspects in neurobiology. Based on the experimental results of neuroanatomy, neurophysiology, behavioral sciences, and the principles of cybernetics and information theory after constructed a simple model simulating prefrontal control function, this paper simulated the behavior of Macaca mulatta completing delayed tasks both before and after its cerebral prefrontal cortex being damaged. The results indicated that there is an obvious difference in the capacity of completing delayed response tasks for the normal monkeys and those of prefrontal cortex cut away. The results are agreement with experiments. The authors suggest that the factors of affecting complete delayed response tasks might be in information keeping and extracting of memory including information storing, keeping and extracting procedures rather than in information storing process.

The Influence of Age and Sex on the Cell Counts of Peripheral Blood Leukocyte Subpopulations in Chinese Rhesus Macaques

Non-human primates such as Chinese rhesus macaques are the favorable models for preclinical study of potential therapeutic drugs, vaccines and mechanisms of human diseases. Little is known about the normal levels of leukocyte subpopulations of Chinese rhe

A STUDY ON THE EFFECT OF LESIONS OF AREA-7 OF THE PARIETAL CORTEX ON THE SHORT-TERM VISUAL SPATIAL MEMORY OF RHESUS-MONKEYS (MACACA-MULATTA)

This research is focused on the contribution of area 7 to the short-term visual spatial memory. Three rhesus monkeys (Macaca mulatta) were trained in the direct delayed response task in which 5 delay intervals were used in each session. When each monkey reached the criterion of 90% correct responses in 5 successive sessions, two monkeys underwent a surgery while the other one received a sham operation as a control. In the first stage of the surgery, bilateral areas 7a, 7b and 7ip of the parietal cortex of two monkeys were precisely lesioned. After 7 days of recuperation, the monkeys were required to do the same task. The average percentage of correct responses in the lesioned animals decreased from 94.7% to 89.3% and 93.3% to 82.0% respectively (no significance, P > 0.05, n = 2). In addition, the monkeys' complex movements were mildly impaired. The lesioned monkeys were found to have difficulty picking up food from the wells. In the second stage, bilateral area 7m was lesioned. In the 5 postoperative sessions, the average percentage of correct responses in one monkey, with a relatively precise 7m lesion, decreased from 94.7% to 92.2% (no significance, P > 0.05), while the other monkey, with widely spread necrosis of lateral parietal cortex, showed an. obvious decline in performance, but still over the chance level. After 240 trials this monkey reattained the normal criterion. The results of this research suggest that the lesions of area 7 of the parietal cortex did not significantly affect the short-term visual spatial memory, which has been shown to be sensitive to lesions of the prefrontal cortex; they also support the notion of dissociation of spatial functions in the prefrontal and parietal cortices.