Results of a CLM4 Land Surface Simulation over China Using a Multisource Integrated Land Cover Dataset

In this study, the high-accuracy multisource integrated Chinese land cover （MICLCover） dataset was used in version 4 of the Community Land Model （CLM4） to assess how the new land cover information affected land surface simulation over China. Compared to the default land cover dataset in CLM4, the MICL data indicated lower values for bare soil （14.6% reduction）, nee- dleleaf tree （3.6%）, and broadleaf tree （1.9%）; higher values for shrub cover （1.8% increase）, grassland （9.9%）, cropland （5.0%）, glaciers （0.5%）, lakes （1.6%）, and wetland （1.1%）; and unchanged for urban areas. Two comparative CLM4 simulations were conducted for the 33-yr period from 1972 to 2004, one using the MICL dataset and the other using the default dataset. The results revealed that the MICL data produced a 0.3% lower mean annual surface albedo over China than the original data. The largest contributor to the reduced value was semiarid regions （2.1% reduction）. The MICL-data albedo value agreed more closely with observations （MODIS broad- band black-sky albedo products） over arid and semiarid regions than for the original data to some extent. The simulated average sensible heat flux over China increased by only 0.1 W m 2 owing to the reduced values in arid and semiarid regions, as opposed to increases in humid and semihumid regions, while an increased latent heat flux of I W m-2 was reflected in almost identical changes over the whole region. In addition, the mean annual runoff simulated by CLM4 using MICL data decreased by 6.8 mm yr-1, primarily due to large simulated decreases in humid regions.

Determination of Slip Length in Couette Flow Based on an Analytical Simulation Incorporating Surface Interaction

An analytical simulation based on a new model incorporating surface interaction is conducted to study the slip phenomenon in the Couette flow at different scales. The velocity profile is calculated by taking account of the micro-force between molecules and macro-force from the viscous shearing effect, as they contribute to the achieve- ment of the slip length. The calculated results are compared with those obtained from the molecular dynamics simulation, showing an excellent agreement. Further, the effect of the shear rate on the slip is investigated. The results can well predict the fluid flow behaviors on a solid substrate, but has to be proved by experiment.

Simulation of Far-Field Properties of Coherent Vertical Cavity Surface Emitting Laser Array

Far-field properties dependent on array scale, separation, element width and emitted wavelength are system atically analyzed theoretically and experimentally. An array model based on the finite-difference method is established to simulate the far-field profile of the coherent arrays. Some important conclusions are obtained. To achieve a higher quality beam, it is necessary to decrease separation between elements, or to increase the element width. Higher brightness can be achieved in the array with larger scale. Emitted wavelength also has an influence on the far-field profile. These analyses can be extended to the future design of coherent vertical cavity surface emitting laser arrays.

Three-dimensional dynamic sea surface modeling based on ocean wave spectrum

In conventional marine seismic exploration data processing,the sea surface is usually treated as a horizontal free boundary.However,the sea surface is affected by wind and waves and there often exists dynamic small-range fluctuations.These dynamic fluctuations will change the energy propagation path and affect the final imaging results.In theoretical research,different sea surface conditions need to be described,so it is necessary to study the modeling method of dynamic undulating sea surface.Starting from the commonly used sea surface mathematical simulation methods,this paper mainly studies the realization process of simple harmonic wave and Gerstner wave sea surface simulation methods based on ocean wave spectrum,and compares their advantages and disadvantages.Aiming at the shortcomings of the simple harmonic method and Gerstner method in calculational speed and sea surface simulation effect,a method based on wave equation and using dynamic boundary conditions for sea surface simulation is proposed.The calculational speed of this method is much faster than the commonly used simple harmonic method and Gerstner wave method.In addition,this paper also compares the new method with the more commonly used higher-order spectral methods to show the characteristics of the improved wave equation method.

Influence of depth-thickness ratio of mining on the stability of a bedding slope with its sliding surface in concave deformation

In order to study the influence of depth-thickness ratio on bedding slope stability, whose sliding surface is flexural concave in shape under mining conditions, this paper aims to study the characteristics ofdeformarion and damage of bedding sliding with depth-thickness ratios of 200：1,150：1,120：1,100：1 and 50：1 by adopting numerical simulation analysis software combined with laboratory-made ＂under the influence of mining variable sliding surface slope similar simulation test bed＂, and to propose identification methods for slope stability under the infuence of mining. The results show that mining activities under the slope reduce slope stability. With a decrease in the mining depth ratio, the influence of mining on the slope increases gradually, and the damage to the slope gradually expands, the stability of the slope grad- ually reduces, fracture occurs on the slope toe and the central fissure gradually develops to the surface, and reaches slide threshold when the depth-thickness ratio is 50：1.

Three-dimensional Monte Carlo simulation of bulk fin field effect transistor

In this paper, we investigate the performance of the bulk fin field effect transistor （FinFET） through a three- dimensional （3D） full band Monte Carlo simulator with quantum correction. Several scattering mechanisms, such as the acoustic and optical phonon scattering, the ionized impurity scattering, the impact ionization scattering and the surface roughness scattering are considered in our simulator. The effects of the substrate bias and the surface roughness scattering near the Si/SiO2 interface on the performance of bulk FinFET are mainly discussed in our work. Our results show that the on-current of bulk FinFET is sensitive to the surface roughness and that we can reduce the substrate leakage current by modulating the substrate bias voltage.

Design optimization of transonic compressor stage using CFD and response surface model

In order to shorten the design period, the paper describes a new optimization strategy for computationally expensive design optimization of turbomachinery, combined with design of experiment (DOE), response surface models (RSM), genetic algorithm (GA) and a 3-D Navier-Stokes solver(Numeca Fine). Data points for response evaluations were selected by improved distributed hypercube sampling (IHS) and the 3-D Navier-Stokes analysis was carried out at these sample points. The quadratic response surface model was used to approximate the relationships between the design variables and flow parameters. To maximize the adiabatic efficiency, the genetic algorithm was applied to the response surface model to perform global optimization to achieve the optimum design of NASA Stage 35. An optimum leading edge line was found, which produced a new 3-D rotor blade combined with sweep and lean, and a new stator one with skew. It is concluded that the proposed strategy can provide a reliable method for design optimization of turbomachinery blades at reasonable computing cost.

Generation of linear and nonlinear waves in numerical wave tank using clustering technique-volume of fluid method

A two-dimensional （2D） numerical model is developed for the wave sim- ulation and propagation in a wave flume. The fluid flow is assumed to be viscous and incompressible, and the Navier-Stokes and continuity equations are used as the governing equations. The standard k-e model is used to model the turbulent flow. The Navier- Stokes equations are discretized using the staggered grid finite difference method and solved by the simplified marker and cell （SMAC） method. Waves are generated and propagated using a piston type wave maker. An open boundary condition is used at the end of the numerical flume. Some standard tests, such as the lid-driven cavity, the constant unidirectional velocity field, the shearing flow, and the dam-break on the dry bed, are performed to valid the model. To demonstrate the capability and accuracy of the present method, the results of generated waves are compared with available wave theories. Finally, the clustering technique （CT） is used for the mesh generation, and the best condition is suggested.

Optimizing the Qusai-static Folding and Deploying of Thin-Walled Tube Flexure Hinges with Double Slots

The thin-walled tube flexure（TWTF） hinges have important potential application value in the deployment mechanisms of satellite and solar array, but the optimal design of the TWTF hinges haven＇t been completely solved, which restricts their applications. An optimal design method for the qusai-static folding and deploying of TWTF hinges with double slots is presented based on the response surface theory. Firstly, the full factorial method is employed to design of the experiments. Then, the finite element models of the TWTF hinges with double slots are constructed to simulate the qusai-static folding and deploying non-linear analysis. What＇s more, the mathematical model of the TWTF flexure hinge quasi-static folding and deploying properties are derived by the response surface method. Considering of small mass and high stability, the peak moment of quasi-static folding and deploying as well as the lightless are set as the objectives to get the optimal performances. The relative errors of the objectives between the optimal design results and the FE analysis results are less than 7%, which demonstrates the precision of the surrogate models. Lastly, the parameter study shows that both the slots length and the slots width both have significant effects to the peak moment of quasi-static folding and deploying of TWTF hinges with double slots. However, the maximum Mises stress of quasi-static folding is more sensitive to the slots length than the slots width. The proposed research can be applied to optimize other thin-walled flexure hinges under quasi-static folding and deploying, which is of great importance to design of flexure hinges with high stability and low stress.

A simulation study on the relation between muscle motor unit numbers and the non-Gaussianity/non-linearity levels of surface electromyography

Recent research has demonstrated that surface electromyography （sEMG） signals have non-Gaussianity and non-linearity properties. It is known that more muscle motor units are recruited and firing rates （FRs） increase as exertion increases. A hy- pothesis was proposed that the Gaussianity test （Sg） and linearity test （St） levels of sEMG signals are associated with the num- ber of active motor units （nMUs） and the FR. The hypothesis has only been preliminarily discussed in experimental studies. We used a simulation sEMG model involving spatial （active MUs） and temporal （three FRs） information to test the hypothesis. Higher-order statistics （HOS） from the bi-frequency domain were used to perform Sg and St. Multivariate covariance analysis and a correlation test were employed to determine the nMUs-Sg relationship and the nMUs-St relationship. Results showed that nMUs, the FR, and the interaction of nMUs and the FR all influenced the Sg and St values. The nMUs negatively correlated to both the Sg and St values. That is, at the three FRs, sEMG signals tended to a more Gaussian and linear distribution as exertion and nMUs increased. The study limited experiment factors to the sEMG non-Gaussianity and non-linearity levels. The study quantitatively described nMUs and the FR of muscle that are not directly available from experiments. Our finding has guiding significance for muscle capability assessment and prosthetic control.

NUMERICAL SIMULATION OF HIGH TURBULENT FREE SURFACE FLOW UNDER THE CURVILINEAR COORDINATES

This paper presents a new numerical method to simulate the high velocity turbulent flow with free surface by solving two-dimensional incompressible unsteady Navier-Stokes Eqs. , together with the k-ε turbulence model. In order to treat the non-rectangular boundary (or curvilinear boundary), orthogonal boundary-fitted grid is used and the Navier-Stokes Eqs. and k-ε turbulence model are rewritten and discreted in orthogonal curvilinear coordinates. Meanwhile, gas-liquid two-field model theory is introduced to treat the free-surface problem.

NUMERICAL SIMULATION OF THE SURFACE RADIATION FIELD IN A HILLY COUNTRY

Based on the simulations of the topographic parameters(mean orientation,slope,and terrain screen angle)in a hilly land,discussion is made term by term of the technique for modelling of all components is done in the surface radiation balance over a rugged terrain,thus presenting a computer model for each com- ponent.In terms of experimental data,the components is calculated for the mesh grids of 100m×100 m each over 3.0×3.5 km^2 in the Zhaogongling,southern Dabie Mountains and a map is prepared showing the distribution of these components for January and July.Results show that the hilly-land surface radiation field matches the terrain element field quite well,which reveals the determining function of the latter,with orientational effect dominant in winter and terrain screen effect most significant in summer.The simulation technique presented is in principle applicable to the calculation for a radiation field in any kind of topogra- phy,thus providing a means for further exploration of hilly-land climate resource.

Truncated separation method for characterizing and reconstructing bi-Gaussian stratified surfaces

Existing ISO segmented and continuous separation methods for differentiating the two components contained within a bi-Gaussian stratified surface were developed based on the fit of the probability material ratio curve.In the present study,because of the significant effect of the plateau component on tribological behavior such as asperity contact,wear and friction,a truncated separation method is proposed based on the truncation of the upper Gaussian component defined by zero skewness.The three separation methods are applied to real worn surfaces.Surface-separation and surface-reconstruction results show that the truncated method accurately captures the upper component identically to the ISO and continuous ones.The identification of the lower component characteristics requires performing a curve fit procedure on the data left after truncation.However,the truncated method fails in identifying the upper component when the material ratio of the transition is less than 9％.

Analysis of the shape of heavy droplets on flat and spherical surface

In this study,a theoretical model was established for predicting the equilibrium shape of the droplet on flat and spherical surfaces.The theoretical equilibrium shape of heavy droplets could be obtained once contact angle and volume of droplets were given.It showed that the predictions of the theoretical flat model were in good agreement with the shape obtained by Surface Evolver when the contact angle is below 120 and the droplet size is on the order of capillary length.This available range will decrease and increase when the heavy droplet is on convex and concave spherical surface,respectively,in contrast to that on flat surface.The available range will decrease more for higher curvature of convex spherical surfaces.

Water-trapping and drag-reduction effects of fish Ctenopharyngodon idellus scales and their simulations

In the last decades, surface drag reduction has been re-emphasized because of its practical values in engineering applications,including vehicles, aircrafts, ships, and fuel pipelines. The bionic study of drag reduction has been attracting scholars' attentions. Here, it was determined that the delicate microstructures on the scales of the fish Ctenopharyngodon idellus exhibit remarkable drag-reduction effect. In addition, the underlying drag-reduction mechanism was carefully investigated. First,exceptional morphologies and structures of the scales were observed and measured using a scanning electron microscope and3-dimensional(3D) microscope. Then, based on the acquired data, optimized 3D models were created. Next, the mechanism of the water-trapping effect of these structures was analyzed through numerical simulations and theoretical calculations. It was determined that there are many microcrescent units with certain distributions on its surface. In fact, these crescents are effective in generating the "water-trapping" effect and forming a fluid-lubrication film, thus reducing the skin friction drag effectively.Contrasting to a smooth surface, the dynamics finite-element analysis indicated that the maximum drag-reduction rate of a bionic surface is 3.014% at 0.66 m/s flow rate. This study can be used as a reference for an in-depth analysis on the bionic drag reduction of boats, underwater vehicles, and so forth.

Simulation of the Effects of the Preceding SST Anomalies over the Tropical Eastern Pacific on Precipitation to the South of the Yangtze River in June

Numerical experiments are performed to simulate the response of the atmospheric circulation and pre-cipitation over East China in June to the sea surface temperature（SST）anomalies over the tropical eastern Pacific（TEP）from preceding September to June by using an atmospheric general circulation model （AGCM）.We constructed composite positive/negative SST anomalies（P-SSTAs/N-SSTAs）based on the observational SST anomalies over the TEP from September 1997 to June 1998.The results show that：（1） the response of the precipitation in the Yangtze River basin and its southern area（YRBS）to El Nino with different durations varies with the maximum amplitude of the precipitation anomalies appearing when the imposed duration is from November to next June,and the minimum appearing when the SST anomalies is only imposed in June.The anomalies of the precipitation are reduced when the duration of the forcing SST anomalies over the TEP is shortened and the positive SST anomalies in the preceding autumn tend to cause significantly more rainfall in the YRBS.This is in agreement with previous diagnostic analysis results.（2）The simulated precipitation anomalies over the YRBS are always obviously positive under strong or weak positive SST anomalies over the TEP.The intensity of the precipitation anomalies increases with increasing intensity of the SST anomalies in the experiments.The simulation results are consistent with the observations during the warm SST events,suggesting reasonable modeling results.（3）When negative SST anomalies in the TEP are put into the model,the results are different from those of the diagnostic analysis of La Nina events.Negative precipitation anomalies in YRBS could be reproduced only when the negative SST anomalies are strong enough.

Influence of heating parameters on properties of the Al-Si coating applied to hot stamping

The Al-Si coating of ultra-high strength steel has been applied to hot stamping more and more widely, owing to solving the problem of oxidation and decarburization. However, the evolution of Al-Si coating during the heating process was rarely studied in the previous study. The tests about the influence of heating parameters, such as heating temperature, heating rates and dwell time, on properties of the Al-Si coating were carried out on the Gleeble-3500 thermal simulator. The properties of the Al-Si coating, for instance, volume fraction of FeAl intermetallics, α-Fe layer as well as porosity and 3D surface topography, were explored in the study. Results showed that more and more Kirkendall voids and cracks appeared in the Al-Si coating when the heating temperature exceeded 600°C. The heating rates almost had no influence on properties of the Al-Si coating when the temperature was equal to or lower than 500°C. The volume fraction of FeAl intermetallics in the coating with dwell time from 3 s to 8 min at 930°C was0, 6.19%, 17.03% and 20.65%, separately. The volume fraction of the α-Fe layer in the coating changed from zero to 31.52%with the prolonged dwell time. The porosity of the coating ranged from 0.51% to 4.98% with the extension of dwell time. The unsmooth degree of the surface of the coating rose gradually with the increasing of heating rates and the extension of dwell time.The 3D surface topography of the coating was determined by the comprehensive effect of atoms diffusion, new formed phases,surface tension and the degree of oxidation of the coating surface. Experiments indicated that rapid heating was not suitable for the coating when the temperature exceeded 500°C. Experiments also demonstrated that enough dwell time was essential to obtain the superior properties of the coating.