The kinetic-theory for dilute solid liquid 2-phase flow

On the basis of a brief review of the continuum theory for macroscopic descriptions and the kinetic theory for microscopic descriptions in solid/liquid two-phase flows, some suggestions are presented, i.e. the solid phase may be described by the Boltzmann equation and the liquid phase still be described by conservation laws in the continuum theory. Among them the action force on the particles by the liquid fluid is a coupling factor which connects the phases. For dilute steady solid/liquid two-phase flows, the particle velocity distribution function can be derived by analogy with the procedures in the kinetic theory of gas molecules for the equilibrium state instead of being assumed, as previous investigators did. This done, more detailed information, such as the velocity probability density distribution, mean velocity distribution and fluctuating intensity etc. can be obtained directly from the particle velocity distribution function or from its integration. Experiments have been performed for dilute solid/liquid two-phase flow in a 4 x 6 cm2 sized circulating square pipe system by means of laser Doppler anemometry so that the theories can be examined. The comparisons show that the theories agree very well with all the measured data.

Kinetic theory for particle concentration distribution in two phase flow

The results of experiments in open channels and closed pipelines show two kinds of patterns for the vertical distribution of particle concentration (i.e., pattern I and pattern II). The former shows a pattern of maximum concentration at some location above the bottom and the downward decay of the concentration below the location. The latter always shows an increase of the particle concentration downward over the whole vertical, with the maximum value at the bottom. Many investigations were made on the pattern II, but few were made on pattern I. In this paper, a particle velocity distribution function is first obtained in the equilibrium state or in dilute steady state for the particle in two-phase flows, then a theoretical model for the particle concentration distribution is derived from the kinetic theory. More attention is paid to the predictions of the concentration distribution of pattern I and comparisons of the present model are made with the data measured by means of laser doppler anemometry (LDA). Very good agreements are obtained between the measured and calculated results.

Kinetic model of continuous-wave flow chemical lasers

In this paper an analysis of the kinetic theory of the continuous-wave flow chemical lasers(CWFCL) is presented with emphasis being laid on the effects of inhomogeneous broadeningon CWFCL's performance. The results obtained are applicable to the case where laser fre-quency is either coincident or incoincident with that of the eenter of the line shape. This rela-tion has been,compared with that of the rate model in common use. These two models are almostidentical as the broadening parameter η is larger than 1. The smaller the value of η, thegreater the difference between the results of these two models will be. For fixed η, the dif-ferences between fhe results of the two models increase with the increase of the frequencyshift parameter ξ. When η is about less than 0.2. the kinetic model can predict exactly the in-homogeneous broadening effects,while the rate model cannot.

A Molecular Dynamics Simulation: the Effect of Finite Size on the Thermal Conductivity in a Single Crystal Silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.

高温气体反应流中的物理力学问题

从气体动理论(kinetic theory of gases)的原理出发导出了非平衡条件下高温气体反应流基本守恒方程组.讨论了出现在方程组中的各种通量在非平衡条件下的意义,给出了从微观出发的输运系数表达式.介绍各种物理动理论与化学动理论过程的微观理论,讨论了在非平衡条件下现有理论的困难.

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Rarefied gas dynamics: fundamentals, simulations and micro flows

Dislocation nucleation governed softening and maximum strength in nano-twinned metals

In conventional metals, there is plenty of space for dislocations-line defects whose motion results in permanent material deformation-to multiply, so that the metal strengths are controlled by dislocation interactions with grain boundaries(1,2) and other obstacles(3,4). For nano-structured materials, in contrast, dislocation multiplication is severely confined by the nanometre-scale geometries so that continued plasticity can be expected to be source-controlled. Nano-grained polycrystalline materials were found to be strong but brittle(5-9), because both nucleation and motion of dislocations are effectively suppressed by the nanoscale crystallites. Here we report a dislocation-nucleation-controlled mechanism in nano-twinned metals(10,11) in which there are plenty of dislocation nucleation sites but dislocation motion is not confined. We show that dislocation nucleation governs the strength of such materials, resulting in their softening below a critical twin thickness. Large-scale molecular dynamics simulations and a kinetic theory of dislocation nucleation in nano-twinned metals show that there exists a transition in deformation mechanism, occurring at a critical twin-boundary spacing for which strength is maximized. At this point, the classical Hall-Petch type of strengthening due to dislocation pile-up and cutting through twin planes switches to a dislocation-nucleation-controlled softening mechanism with twin-boundary migration resulting from nucleation and motion of partial dislocations parallel to the twin planes. Most previous studies(12,13) did not consider a sufficient range of twin thickness and therefore missed this strength-softening regime. The simulations indicate that the critical twin-boundary spacing for the onset of softening in nano-twinned copper and the maximum strength depend on the grain size: the smaller the grain size, the smaller the critical twin-boundary spacing, and the higher the maximum strength of the material.

Miscibility and crystallization of poly(beta-hydroxybutyrate)/poly(vinyl acetate-co-vinyl alcohol) blends

Poly(vinyl acetate-co-vinyl alcohol) copolymers (P(VAc-co-VA)) were synthesized by hydrolysis-alcoholysis of PVAc. The miscibility, crystallization, and morphology of poly(P-hydroxybutyrate) (PHB) and P(VAc-co-VA) blends were studied by differential scanning calorimetry, optical microscopy (OM), and SAXS. It is found that the P(VAc-co-VA)s with vinyl alcohol content of 9, 15, and 22 mol % will form a miscible phase with the amorphous part of PHB in the solution-cast samples. The melting-quenched samples of PHB/P(VAc-co-VA) blends with different vinyl alcohol content show different phase behavior. PHB and P(VAc-co-VA9) with low vinyl alcohol content (9% mel) will form a miscible blend in the melt state. PHB and P(VAc-co-VA15) with 15 mol % vinyl alcohol will not form miscible blends while PHB/P(VAc-co-VA15) blend with 20/80 composition will form a partially miscible blend in the melt state. PHB and P(VAc-co-VA22) with 22 mol % vinyl alcohol are not miscible in the whole composition range. The single glass transition temperature of the blends within the whole composition range suggests that PHB and P(VAc-co-VA9) are totally miscible in the melt. The crystallization kinetics was studied from the whole crystallization and spherulite growth for the miscible blends. The equilibrium melting point of PHB in the PHB/P(VAc-co-VA9) blends, which was obtained from DSC results using the Hoffman-Weeks equation, decreases with the increase in P(VAc-co-VA9) content. The negative value of the interaction parameter determined from the equilibrium melting point depression supports the miscibility between the components. The kinetics of spherulitic crystallization of PHB in the blends was analyzed according to nucleation theory in the temperature range studied in this work. The best fit of the data to the kinetic theory is obtained by employing WLF parameters and the equilibrium melting points obtained by DSC. The addition of P(VAc-co-VA) did not affect the crystalline structure of PHB, as shown by the WAXD results. The long periods of blends obtained from SAXS increase with the increase in P(VAc-co-VA) content. It indicates that the amorphous P(VAc-co-VA) was rejected to interlamellar phase corporating with the amorphous part of PHB.

Kinetic study for hopping conduction through deoxyribonucleic acid molecules

Recent experiments indicated that disorder effect in deoxyribonucleic acid (DNA) may lead to a transition of the electronic hole transport mechanism from band resonant tunneling to thermally activated hopping. In this letter, based on Mott's variable-range hopping theory, we present a kinetic study for the hole transport properties of DNA molecules. Beyond the conventional argument in large-scale systems, our numerical study for finite-size DNA molecules reveals a number of unique features for: (i) the current-voltage characteristics, (ii) the temperature and length dependence, and (iii) the transition from conducting to insulating behaviors. (c) 2005 American Institute of Physics.

Kinetic study of MOCVD III-V quaternary antimonides

The kinetics of MOCVD GaInAsSb and AlGaAsSb was studied by the growth rate as a function of growth temperature and partial pressure of III and V MO species. The diffusion theory was used to explain the mass transport processes in MOCVD III-V quaternary antimonides. On the basis of the discussion about their growth kinetics and epilayer properties, the good quality multi-epilayers of these two quaternary antimonides and their photodetectors and arrays with wavelength of 1.8 similar to 2.3 mu m and detectivities of D* > 10(9) cm Hz(1/2) W-1 were obtained.

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion.Including the configurational dependence will challenge the transition state theory of protein folding.

Primary and secondary crystallization kinetic analysis of nylon 1212

The isothermal and non-isothermal melt-crystallization kinetics of nylon 1212 were investigated by differential scanning calorimetry. Primary and secondary crystallization behaviors were analysed based on different approaches. The results obtained suggested that primary crystallization under isothermal conditions involves three-dimensional spherulite growth initiated by athermal nucleation, while under non-isothermal conditions, the mechanism of primary crystallization is more complex. Secondary crystallization displays a lower-dimensional crystal growth, both in the isothermal and non-isothermal processes. The crystallite morphology of nylon 1212, isothermally crystallized at various temperatures, was observed by polarized optical microscopy. The activation energies of crystallization under isothermal and non-isothermal conditions were also calculated based on different approaches.