Numerical experiments on one-dimensional model of turbulence

The initial-value problem of a forced Burgers equation is numerically solved by the Fourier expansion method. It is found that its solutions finally reach a steady state of 'laminar flow' which has no randomness and is stable to disturbances. Hence, strictly speaking, the so-called Burgers turbulence is not a turbulence. A new one-dimensional model is proposed to simulate the Navier-Stokes turbulence. A series of numerical experiments on this one-dimensional turbulence is made and is successful in obtaining Kolmogorov's (1941) k exp(-5/3) inertial-range spectrum. The (one-dimensional) Kolmogorov constant ranges from 0.5 to 0.65.

A transition layer model and its application to resonant tunneling in heterostructures

A transition layer model is proposed and used to calculate resonant tunneling in a double-barrier quantum well system. Compared with the ideal step of the potential at the interface, the studied system has transition layers that are composed by many thin rectangular barriers with a random height. It is found that these transition layers can improve the peak-to-valley ratio of the tunneling current and change the negative differential conductance.

Linear Stability Analysis of Convection in Two-Layer System with an Evaporating Vapor-Liquid Interface

Classical theories have successfully provided an explanation for convection in a liquid layer heated from below without evaporation. However, these theories are inadequate to account for the convective instabilities in an evaporating liquid layer, especially in the case when it is cooled from below. In the present paper, we study the onset of Marangoni convection in a liquid layer being overlain by a vapor layer.A new two-sided model is put forward instead of the one-sided model in previous studies. Marangoni-Bénard instabilities in evaporating liquid thin layers are investigated with a linear instability analysis. We define a new evaporation Biot number, which is different from that in previous studies and discuss the influences of reference evaporating velocity and evaporation Biot number on the vapor-liquid system. At the end, we explain why the instability occurs even when an evaporating liquid layer is cooled from below.

Physical mechanism and numerical simulations of surface layer temperature inversion in tropical ocean

The one-dimensional Kraus-Turner mixed layer model improved by Liu is developed to consider the effect of salinity and the equations of temperature and salinity under the mixed layer. On this basis, the processes of growth and death of surface layer temperature inversion is numerically simulated under different environmental parameters. At the same time, the physical mechanism is preliminarily discussed combining the observations at the station of TOGA-COARE 0 degrees N, 156 degrees E. The results indicate that temperature inversion sensitively depends on the mixed layer depth, sea surface wind speed and solar shortwave radiation, etc., and appropriately meteorological and hydrological conditions often lead to the similarly periodical occurrence of this inversion phenomenon.

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

Marangoni-Benard instability with the exchange of evaporation at liquid-vapour interface

A new two-sided model rather than the one-sided model in previous works is put forward. The linear instability analysis is performed on the Marangoni-Benard convection in the two-layer system with an evaporation interface. We define a new evaporation Biot number which is different from that in the one-sided model, and obtain the curves of critical Marangoni number versus wavenumber. The influence of evaporation velocity and Biot number on the system is discussed and a new phenomenon uninterpreted before is now explained from our numerical results.

Nonequilibrium statistical mechanics of one-dimensional turbulence

The method of statistical mechanics is applied to the study of the one-dimensional model of turbulence proposed in an earlier paper. The closure problem is solved by the variational approach which has been developed for the three-dimensional case, yielding two integral equations for two unknown functions. By solving the two integral equations, the Kolmogorov k−5/3 law is derived and the (one-dimensional) Kolmogorov constant Ko is evaluated, obtaining Ko=0.55, which is in good agreement with the result of numerical experiments on one-dimensional turbulence.

Influence of Evaporation on Marangoni-Benard Instability

Evaporative convection and instability give rise to both scientific and technological interests. Practically, a number of the industrial applications such as thin-film evaporators, boiling technologies and heat pipes concern with the evaporation process of which through the vapor-liquid interface the heat and mass transfer occur. From a physical viewpoint, one of interesting questions is the mechanisms of convection instability in thin-liquid layers induced by the coupling of evaporation phenomenon and Marangoni effect at the mass exchanged interface. Classical theories, including Rayleigh’s and Pearson’s, have only successfully explained convection in a liquid layer heated from below without evaporation. However these theories are unable to explain the convection in an evaporating thin layer, especially liquid layer is cooled from below. In present paper, a new two-sided model is put forward rather than the one-sided model in previous works. In previous works, the vapor is treated as passive gas and dynamics of vapor has been ignored. In this case, the vapor liquid system can be described by one-sided model. In our two-sided model, the dynamics of vapor should be considered. Linear instability analysis of the Marangoni-Bénard convection in the two-layer system with an evaporation interface is performed. We define a new evaporating Biot number which is different from the Biot number in one-sided model and obtain the curves of critical Marangoni number versus wave number. In our theoretical results, the Biot number and the evaporating velocity play a major role in the stability of the vapor-liquid system.

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.

Mechanisms underlying two kinds of surface effects on elastic constants

ABSTRACT Recently, people are confused with two opposite variations of elastic modulus with decreasing size of nano scale sample: elastic modulus either decreases or increases with decreas- ing sample size. In this paper, based on intermolecular potentials and a one dimensional model, we provide a unified understanding of the two opposite size effects. Firstly, we analyzed the mi- crostructural variation near the surface of an fcc nanofilm based on the Lennard-Jones potential. It is found that the atomic lattice near the surface becomes looser in comparison with the bulk, indicating that atoms in the bulk are located at the balance of repulsive forces, resulting in the decrease of the elastic moduli with the decreasing thickness of the film accordingly. In addition, the decrease in moduli should be attributed to both the looser surface layer and smaller coor- dination number of surface atoms. Furthermore, it is found that both looser and tighter lattice near the surface can appear for a general pair potential and the governing mechanism should be attributed to the surplus of the nearest force to all other long range interactions in the pair po- tential. Surprisingly, the surplus can be simply expressed by a sum of the long range interactions and the sum being positive or negative determines the looser or tighter lattice near surface re- spectively. To justify this concept, we examined ZnO in terms of Buckingham potential with long range Coulomb interactions. It is found that compared to its bulk lattice, the ZnO lattice near the surface becomes tighter, indicating the atoms in the bulk located at the balance of attractive forces, owing to the long range Coulomb interaction. Correspondingly, the elastic modulus of one- dimensional ZnO chain increases with decreasing size. Finally, a kind of many-body potential for Cu was examined. In this case, the surface layer becomes tighter than the bulk and the modulus increases with deceasing size, owing to the long range repulsive pair interaction, as well as the cohesive many-body interaction caused by the electron redistribution.

Instabilities of a liquid film flowing down an inclined porous plane

The problem of a film flowing down an inclined porous layer is considered. The fully developed basic flow is driven by gravitation. A careful linear instability analysis is carried out. We use Darcy's law to describe the porous layer and solve the coupling equations of the fluid and the porous medium rather than the decoupled equations of the one-sided model used in previous works. The eigenvalue problem is solved by means of a Chebyshev collocation method. We compare the instability of the two-sided model with the results of the one-sided model. The result reveals a porous mode instability which is completely neglected in previous works. For a falling film on an inclined porous plane there are three instability modes, i.e., the surface mode, the shear mode, and the porous mode. We also study the influences of the depth ratio d, the Darcy number delta, and the Beavers-Joseph coefficient alpha(BJ) on the instability of the system.

Study of turbulent heat transfer of aviation kerosene flows in a curved pipe at supercritical pressure

The now and heat transfer characteristics of China No. 3 aviation kerosene in a heated curved tube under supercritical pressure are numerically investigated by a finite volume method. A two-layer turbulence model, consisting of the RNG k-epsilon two-equation model and the Wolfstein one-equation model, is used for the simulation of turbulence. A 10-species kerosene surrogate model and the NIST Supertrapp software are applied to obtain the thermophysical and transport properties of the kerosene at various temperature under a supercritical pressure of 4 MPa. The large variation of thermophysical properties of the kerosene at the supercritical pressure make the flow and heat transfer more complicated, especially under the effects of buoyancy and centrifugal force. The centrifugal force enhances the heat transfer, but also increases the friction factors. The rise of the velocity caused by the variation of the density does not enhance the effects of the centrifugal force when the curvature ratios are less than 0.05. On the contrary, the variation of the density increases the effects of the buoyancy. (C) 2010 Elsevier Ltd. All rights reserved.

Instabilities and transient behaviors of a liquid film flowing down a porous inclined plane

The nonmodal linear stability of a falling film over a porous inclined plane has been investigated. The base flow is driven by gravity. We use Darcy's law to describe the flow in the porous medium. A simplified one-sided model is used to describe the fluid flow. In this model, the influence of the porous layer on the flow in the film can be identified by a parameter beta. The instabilities of a falling film have traditionally been investigated by linearizing the governing equations and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In the present paper, we have studied the linear stability of three-dimensional disturbances using the nonmodal stability theory. Particular attentions are paid to the transient behavior rather than the long time behavior of eigenmodes predicted by traditional normal mode analysis. The transient behaviors of the response to external excitations and the response to initial conditions are studied by examining the pseudospectral structures and the energy growth function G(t) Before we study the nonmodal stability of the system, we extend the results of long-wave analysis in previous works by examining the linear stabilities for streamwise and spanwise disturbances. Results show that the critical conditions of both the surface mode and the shear mode instabilities are dependent on beta for streamwise disturbances. However, the spanwise disturbances have no unstable eigenvalue. 2010 American Institute of Physics. [doi:10.1063/1.3455503]

中太平洋海区富钴结壳地球化学特征及成因机制

本文利用ICP-AES、EPMA、X-ray衍射等测试技术以及聚类分析和因子分析等多种数理统计方法，系统地对中太平洋海区富钴结壳的元素地球化学特征、矿物组成和微观组构进行了研究，并探讨了其成因机制，获得以下主要认识： 1中太平洋海区富钴结壳类型多样，均为水成成因，其矿物相主要由锰矿物相、铁矿物相和非金属矿物相组成；富钴结壳壳层发育多种原生和次生构造类型。 2磷酸盐化作用不仅强烈改变富钴结壳元素初始含量，而且造成富钴结壳某些元素间的相关性发生改变，这些敏感型元素对可用于指示富钴结壳是否发生磷酸盐化。在不同水深段内，未磷酸盐化型富钴结壳的主要元素随经向、纬向的变化趋势相似，表明其受水体化学障、表面生产力和物质来源等环境参数控制；而随水深的变化则具有区域一致性，表明水体化学具有区域成层性。 3未磷酸盐化富钴结壳稀土元素含量和轻重稀土分馏程度随水深发生规律变化，这种变化不仅与它们在海洋中的含量和行为有关，也与海洋背景颗粒的吸附有关；铈(Ce)在富钴结壳中基本上呈4价，且动力学因素控制了其富集过程，因此Ce异常不能用于指示富钴结壳形成环境的氧化程度。 4基于富钴结壳微层呈锯齿状且同一微层生长速率不同，提出了富钴结壳在各种基底表面生长以及后继发育过程受固液界面双电层控制的发育模式。在富钴结壳整个发育过程中，经历了从贫氧环境向富氧环境的转变，但微环境则呈富氧-低氧过程的交替。 关键词： 富钴结壳；中太平洋；元素地球化学；界面双电层