Effect of polymer chain structure and chain distance on luminescent efficiency

Three kinds of PPV-based copolymers were synthesized and characterized. Their luminescent properties were investigated and discussed by PL spectrum, and time-dependence luminescent spectrum in film and solution stare. The results show that in the range of our study, the PL intensity and lifetime of luminescent decay increase with the increasing length of flexible segments and the solution diluting, indicating the tendency of the increase of luminescent efficiency.

Refining near-native protein-protein docking decoys by local resampling and energy minimization

How to refine a near-native structure to make it closer to its native conformation is an unsolved problem in protein-structure and protein-protein complex-structure prediction. In this article, we first test several scoring functions for selecting locally resampled near-native protein-protein docking conformations and then propose a computationally efficient protocol for structure refinement via local resampling and energy minimization. The proposed method employs a statistical energy function based on a Distance-scaled Ideal-gas REference state (DFIRE) as an initial filter and an empirical energy function EMPIRE (EMpirical Protein-InteRaction Energy) for optimization and re-ranking. Significant improvement of final top-1 ranked structures over initial near-native structures is observed in the ZDOCK 2.3 decoy set for Benchmark 1.0 (74% whose global rmsd reduced by 0.5 angstrom or more and only 7% increased by 0.5 angstrom or more). Less significant improvement is observed for Benchmark 2.0 (38% versus 33%). Possible reasons are discussed.

磁层-电离层耦合现象的分析研究

The magnetosphere-ionosphere coupling is mainly manifested by the trans- porting processes of energy into the ionosphere , the energy is carried by solar wind and firstly accumulate at the magnetosphere, and the coupling processes also significantly include the interaction between the magnetosphere and ionosphere for mass and energy. At the quiet condition, energy is delivered by the large-scale convection of the geomagnetic field; the huge energy from solar wind bulk will be injected into and consumed at the near magnetosphere and ionosphere by the geomagnetic storm and substorm activities. Aurorae and FACs (Field-aligned currents) are the important phenomena in the coupling processes. In the present work, firstly, we analyze the activity characteristics of auroral precipitating particle, secondly, we study the distribution characters of large-scale field aligned currents (LS FACs) at storm-time using the observations from different satellites at different altitudes. Finally, we investigate the evolution of the geomagnetic field configuration at the nightside sector on the onset of the expansion phase in a substorm event, the substorm event happened at 0430UT to 0630UT on 8th Nov. 2004. The main results as follows： At the first, the data of the estimated power input (EPI) of auroral particles from NOAA/POES (Polar orbiting environmental satellite) for some 30 years have been analyzed. The variation tendencies of the EPI generally coincide with aa, AE and Dst indices. The annual variation of EPI shows equinox peaks and an asymmetric-activity with a higher peak in the winter-hemisphere than in the summer-hemisphere. The diurnal UT variations are different from north and south hemisphere: for north hemisphere, the peak appears at 1200UT, and the relative deviation is 22% to the daily average of the north hemisphere. For south hemisphere, the maximal deviation is 22% at 2000UT. So the diurnal variation of EPI is more dominant than the annual variation which maximal deviation is 3% to 12% for different seasons. Studies on correlations of the hourly average of EPI, Pa, with AE and Dst indices show a correlation coefficient r=0.74 of Pa and AE, and r=-0.55 of Pa and Dst. The hourly EPIs for north and south polar regions, NPa and SPa, show a north-south asymmetry with a higher correlation of SPa and AE (or Dst). Time delays of EPI with respect to magnetic indices are examined, the maximum correlation coefficient of Pa with AE (r=0.78) occurs when the time delay =0, suggesting a synchronous activity of auroral electrojet and auroral precipitating particles, while =1-2h, the correlation coefficient of Pa with Dst is maximum (r=0.57), suggesting that the activity of auroral particle precipitating may influence the ring current on some extent. Sencondly, we use the high-resolution magnetic field vector data of the CHAMP satellite to investigate the distribution of large-scale FACs during the great magnetic storm on 7th to 8th Nov. 2004. The results show that, whether in the northern or southern hemisphere, the number and density of large-scale FACs during the main-phase are more and bigger than these during the recover-phase, and the number of large-scale FACs in morning sector obviously is more than that in afternoon sector. In terms of the magnetic indices, we find that large-scale FACs in morning sector significantly affected by the substorm activities, while in afternoon sector the large-scale FACs mainly indicate the fluctuations of the ring-current in storm time. Accordingly to the former studies, similarly, we find that in the morning sector, the scale of the large-scale FACs move to the high-latitude region, and in the afternoon sector, large-scale FACs distinctly expand to the low-latitude region. During the time periods that the NOAA/POES auroral precipitating particle power data temporally correspond to the large-scale FACs, the more the power of auroral particle is, the more and bigger the number and density of FACs are. At the same time, we use the magnetic field vector data of POLAR obtain a good form of region 1, region 2, and three pieces of cusp FACs during a single transit at 1930UT-2006UT on 07th. And the characteristics of simultaneous electric field and energy particles observations on Polar are coincide with the five FACs pieces. Finally, by means of the observation of Cluster 4 and Goes 10、 Goes 12, we analyze the evolution process of the change of the magnetic field configuration at night sector at the expansion phase of a substorm event which happened during 0430UT to 0630UT on 8th Nov. 2004, we find that the times of the beginning of the polarizations of magnetic field are observed from Goes 10 to Goes 12 then to Cluster 4. So, at the synchronous orbit ( 6.6 RE) to 10RE distance scale of the neutral sheet, the current disruption spread tailward. Simultaneously, the strengthen of the FACs deduced from these satellites’ magnetic field observations are almost consistent with the times of polarizations, as well as the high energy particles injection and the electric field dominant variation. The onset times determined by the magnetic field polarizations from these satellites are all ahead of the onset time that confirmed from the auroral electrojet indices. So, these characters of different observations can be used as the criterions to determine the onset time for the substorms of such type as we studied.

地震P-S波时差耦合作用的斜坡崩滑效应研究

China locates between the circum-Pacific and the Mediterranean-Himalayan seismic belt. The seismic activities in our country are very frequent and so are the collapses and slides of slope triggered by earthquakes. Many collapses and slides of slope take place mainly in the west of China with many earthquakes and mountains, especially in Sichuan and Yunnan Provinces. When a strong earthquake happening, the damage especially in mountains area caused by geological hazards it triggered such as rock collapses, landslides and debris flows is heavier than that it caused directly. A conclusion which the number of lives lost caused by geological hazards triggered by a strong earthquake in mountains area often accounts for a half even more of the total one induced by the strong earthquake can be made by consulting the statistical loss of several representative earthquakes. As a result, geological hazards such as collapses and slides of slope triggered by strong earthquakes attract wide attention for their great costs. Based on field geological investigation, engineering geological exploration and material data analysis, chief conclusions have been drawn after systematic research on formation mechanism, key inducing factors, dynamic characteristics of geological hazards such as collapses and slides of slope triggered by strong earthquakes by means of engineering geomechanics comprehensive analysis, finite difference numerical simulation test, in-lab dynamic triaxial shear test of rock, discrete element numerical simulation. Based on research on a great number of collapses and landslides triggered by Wenchuan and Xiaonanhai Earthquake, two-set methods, i.e. the method for original topography recovering based on factors such as lithology and elevation comparing and the method for reconstructing collapsing and sliding process of slope based on characteristics of seism tectonic zone, structural fissure, diameter spatial distribution of slope debris mass, propagation direction and mechanical property of seismic wave, have been gotten. What is more, types, formation mechanism and dynamic characteristics of collapses and slides of slope induced by strong earthquakes are discussed comprehensively. Firstly, collapsed and slided accumulative mass is in a state of heavily even more broken. Secondly, dynamic process of slope collapsing and sliding consists of almost four stages, i.e. broken, thrown, crushed and river blocked. Thirdly, classified according to failure forms, there are usually four types which are made up of collapsing, land sliding, land sliding-debris flowing and vibrating liquefaction. Finally, as for key inducing factors in slope collapsing and sliding, they often include characteristics of seism tectonic belts, structure and construction of rock mass, terrain and physiognomy, weathering degree of rock mass and mechanical functions of seismic waves. Based on microscopic study on initial fracturing of slope caused by seismic effect, combined with two change trends which include ratio of vertical vs. horizontal peak ground acceleration corresponding to epicentral distance and enlarging effect of peak ground acceleration along slope, key inducing factor of initial slope fracturing in various area with different epicentral distance is obtained. In near-field area, i.e. epicentral distance being less than 30 km, tensile strength of rock mass is a key intrinsic factor inducing initial fracturing of slope undergoing seismic effect whereas shear strength of rock mass is the one when epicentral distance is more than 30 km. In the latter circumstance, research by means of finite difference numerical simulation test and in-lab dynamic triaxial shear test of rock shows that initial fracture begins always in the place of slope shoulder. The fact that fracture strain and shear strength which are proportional to buried depth of rock mass in the place of slope shoulder are less than other place and peak ground acceleration is enlarged in the place causes prior failure at slope shoulder. Key extrinsic factors inducing dynamic fracture of slope at different distances to epicenter have been obtained through discrete element numerical simulation on the total process of collapsing and sliding of slope triggered by Wenchuan Earthquake. Research shows that combined action of P and S seismic waves is the key factor inducing collapsing and sliding of slope at a distance less than 64 km to initial epicenter along earthquake-triggering structure. What is more, vertical tensile action of P seismic wave plays a leading role near epicenter, whereas vertical shear action of S seismic wave plays a leading role gradually with epicentral distance increasing in this range. On the other hand, single action of P seismic wave becomes the key factor inducing collapsing and sliding of slope at a distance between 64 km and 216 km to initial epicenter. Horizontal tensile action of P seismic wave becomes the key factor gradually from combined action between vertical and horizontal tensile action of P seismic wave with epicentral distance increasing in this distance range. In addition, initial failure triggered by strong earthquakes begins almost in the place of slope shoulder. However, initial failure beginning from toe of slope relates probably with gradient and rock occurrence. Finally, starting time of initial failure in slope increases usually with epicentral distance. It is perhaps that the starting time increasing is a result of attenuating of seismic wave from epicenter along earthquake-triggering structure. It is of great theoretical and practical significance for us to construct towns and infrastructure in fragile geological environment along seism tectonic belts and conduct risk management on earthquake-triggered geological hazards by referring to above conclusions.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

人工源频率域大地电磁法正反演研究

The CSAMT method is playing an important role in the exploration of geothermal and the pre-exploration in tunnel construction project recently. In order to instruct the interpretation technique for the field data, the forward method from ID to 3D and inversion method in ID and 2D are developed in this paper for the artificial source magnetotelluric in frequency domain. In general, the artificial source data are inverted only after the near field is corrected on the basis of the assumption of half-homogeneous space; however, this method is not suitable for the complex structure because the assumption is not valid any more. Recently the new idea about inversion scheme without near field correction is published in order to avoid the near field correction error. We try to discuss different inversion scheme in ID and 2D using the data without near field correction.The numerical integration method is used to do the forward modeling in ID CSAMT method o The infinite line source is used in the 2D finite-element forward modeling, where the near-field effect is occurred as in the CSAMT method because of using artificial source. The pseudo-delta function is used to modeling the source distribution, which reduces the singularity when solving the finite-element equations. The effect on the exploration area is discussed when anomalous body exists under the source or between the source and exploration area; A series of digital test show the 2D finite element method are correct, the results of modeling has important significant for CSAMT data interpretation. For 3D finite-element forward modeling, the finite-element equation is derived by Galerkin method and the divergence condition is add forcedly to the forward equation, the forward modeling result of the half homogeneous space model is correct.The new inversion idea without near field correction is followed to develop new inversion methods in ID and 2D in the paper. All of the inversion schemes use the data without near field correction, which avoid introducing errors caused by near field correction. The modified grid parameter method and the layer-by-layer inversion method are joined in the ID inversion scheme. The RRI method with artificial source are developed and finite-element inversion method are used in 2D inversion scheme. The inversion results using digital data and the field data are accordant to the model and the known geology data separately, which means the inversion without near field correction is accessible. The feasibility to invert the data only in exploration area is discussed when the anomalous body exists between the source and the exploration area.

Kinetic and Size Control of Polystyrene and Polyacrylic Octadecyl Ester Lattices Via Polymerization in O/W Microemulsions

The kinetic studies of the acrylic octadecyl ester and styrene polymerization in microemulsion systems, (1) cetyl pyridine bromide (CPDB)/t-butanol/styrene/water; (2) CPDB/t-butanol/toluene + acrylic octadecyl ester (1:1, w/v)/ water; (3) cetyl pyridine bromide/styrene/formamide, were made by using dynamic laser light scattering techniques (DLS). The mechanisms of nucleation of latex particles were discussed. The most possible nucleation location of the styrene and acrylic octadecyl ester microlatex particles in aqueous microemulsion system is in aqueous phase via homogeneous nucleation. Meanwhile, parts of microlatex particles are possibly produced via swollen micelles (microemulsions) and monomer droplets nucleation. On the other hand, the most possible nucleation location of the styrene microlatex particles in nonaqueous microemulsion system is inside monomer droplets. The relationship between the amount of monomer and the size of microlatex was also investigated. It has been found that the size of microlatex particles could be controlled by changing the amount of monomer. (C) 2002 Elsevier Science B.V. All rights reserved.

配液结晶法制备溶菌酶蛋白质晶体的生长机理研究

采用配液结晶法制取了溶菌酶蛋白质晶体,使用动态光散射测量了溶液中聚集体的颗粒度几率分布;使用Zeiss显微镜测定了溶菌酶(110)晶面的生长速度.实验表明:随着蛋白质和NaCl浓度的增加,溶液中聚集体的颗粒尺寸也相应增加.随着反应时间的增加,溶菌酶分子在溶液中的聚集反应,逐渐达到平衡;在蛋白质和NaCl浓度较高时,溶菌酶晶体的(110)面生长较快,而在蛋白质和NaCl浓度较低时,该晶面生长较慢.基于二维成核生长机理,从晶体生长动力学理论方程出发,计算了二维成核的形成能a=4.01×10-8J?cm-2.

NaCl对液-液扩散法生长溶菌酶晶体的影响

用动态光散射法研究了不同浓度NaCl对液—液扩散法生长溶菌酶晶体的影响，并测量了晶体生长前后体系的Zeta电势．结果表明，NaCl浓度较高时，在溶菌酶溶液—凝胶界面处会发生液液分层现象，溶液中一直存在较大的聚集体，生长出的晶体质量较差．而在合适的NaCl浓度下，随着溶液Zeta电势降低，溶液中溶菌酶的大的聚集体发生解聚集，生长出的晶体质量较高．

Study of Growth Mechanism of Lysozyme Crystal by Batch Crystallization Method

The lysozyme crystals were made by batch crystallization method and the distribution of aggregate in solution were measured by dynamic light scattering. The results showed that the dimension of aggregate increased with the increase of the concentration of lysozyme and NaCl, lysozyme molecules aggregated gradually in solution and finally arrived at balance each other. The higher the concentrations of lysozyme and NaCl were, the faster the growth rate of (I 10) face was. The growth rates of lysozyme crystal were obtained by a Zeiss microscope, and the effective surface energy (a) of growing steps were calculated about 4.01 X 10(-8) J.cm(-2) according to the model of multiple two-dimensional nucleation mechanism.

光散射技术在蛋白质晶体生长研究中的应用和进展

光散射技术广泛应用于生物大分子的晶体生长研究中，它包括静态光散射和动态光散射两种。利用静态光散射可以测定蛋白质溶液渗透的第二维里系数；利用动态光散射可以测定蛋白质溶液的平动扩散系数，获得溶液中蛋白质粒子的流体力学半径及分布情况，分离蛋白质结晶的成核与生长过程，研究大分子的聚集行为和晶体生长的动力学。借助光散射技术可以实现蛋白质晶体生长过程的动态控制。近些年光散射仪器向着小型化、轻便化的方向发展，光散射技术不断得到改进，日益完善，不仅用于地面实验，也应用于空间领域蛋白质晶体生长的研究中。

Numerical Investigation On Evolution Of Cylindrical Cellular Detonation

Cylindrical cellular detonation is numerically investigated by solving two-dimensional reactive Euler equations with a finite volume method on a two-dimensional self-adaptive unstructured mesh. The one-step reversible chemical reaction model is applied to simplify the control parameters of chemical reaction. Numerical results demonstrate the evolution of cellular cell splitting of cylindrical cellular detonation explored in experimentas. Split of cellular structures shows different features in the near-field and far-field from the initiation zone. Variation of the local curvature is a key factor in the behavior of cell split of cylindrical cellular detonation in propagation. Numerical results show that split of cellular structures comes from the self-organization of transverse waves corresponding to the development of small disturbances along the detonation front related to detonation instability.

A Novel Inverse Method For Determining The Refractive Indices Of Medium And Dispersed Particles Simultaneously By Turbidity Measurement

The refractive indices of particles and dispersion medium are important parameters in many colloidal experiments using optical techniques, such as turbidity and light scattering measurements. These data are in general wavelength-dependent and may not be available at some wavelengths fitting to the experimental requirement. in this Study we present a novel approach to inversely determine the refractive indices of particles and dispersion medium by examining the consistency of measured extinction cross sections of particles with their theoretical values using a series of trial values of the refractive indices. The colloidal suspension of polystyrene particles dispersed in water was used as an example to demonstrate how this approach works and the data obtained via such a method are compared with those reported in literature, showing a good agreement between both. Furthermore, the factors that affect the accuracy of measurements are discussed. We also present some data of the refractive indices of polystyrene over a range of wavelengths smaller than 400 nm that have been not reported in the available literature. (C) 2008 Elsevier Inc. All rights reserved.

Short-surface waves radiated by a partially-immersed 3-dimensional oscillating body

The short-surface waves generated by a 3-D arbitrarily oscillating body floating onwater are discussed. In the far-field off the body, the phase and the amplitude functions ofthe radiated waves are determined by the ray method. An undetermined constant is includ-ed in the amplitude function. From the result of Ref. [1], the near-field boundary layersolution near the body waterline is obtained. The amplitude of this solution depends on thewhole wall shape of the body and the slope at the body waterline on the cross-sections per-pendicular to the waterline. By matching the far-field solution with the near-field bound-ary layer solution, the undetermined constant in the amplitude function of the far-fieldradiated waves is determined. For the special case of a half-submerged sphere which per-forms vertical oscillating motion, the result obtained in this paper is in agreement withthat of Ref. [ 2 ].

Vortex motion in the early stages of unsteady flow around a circular cylinder

In this paper, processes in the early stages of vortex motion and the development of flow structure behind an impulsively-started circular cylinder at high Reynolds number are investigated by combining the discrete vortex model with boundary layer theory, considering the separation of incoming flow boundary layer and rear shear layer in the recirculating flow region. The development of flow structure and vortex motion, particularly the formation and development of secondary vortex and a pair of secondary vortices and their effect on the flow field are calculated. The results clearly show that the flow structure and vortices motion went through a series of complicated processes before the symmetric main vortices change into asymmetric: development of main vortices induces secondary vortices; growth of the secondary vortices causes the main vortex sheets to break off and causes the symmetric main vortices to become “free” vortices, while a pair of secondary vortices is formed; then the vortex sheets, after breaking off, gradually extend downstream and the structure of a pair of secondary vortices becomes relaxed. These features of vortex motion look very much like the observed features in some available flow field visualizations. The action of the secondary vortices causes the main vortex sheets to break off and converts the main vortices into free vortices. This should be the immediate cause leading to the instability of the motion of the symmetric main vortices. The flow field structure such as the separation position of boundary layer and rear shear layer, the unsteady pressure distributions and the drag coefficient are calculated. Comparison with other results or experiments is also made.