Characteristics of Micronozzle Gas Flows

The fluid characteristics of gas flows in the micronozzle whose throat height is 20 μm were investigated by the direct simulation Monte Carlo (DSMC) method. In a series of cases, the dependence of mass flux on the pressure difference was gained, and the DSMC's results show good agreement with the experimental data. The comparison of mass flux and the Mach number contours between the DSMC and Navier-Stokes equations adding slip boundary also reveals quantitatively that the continuum model will be invalid gradually even when the average Knudsen number is smaller than 0.01. As one focus of the present paper, the phenomenon of the multiple expansion-compression waves that comes from the nozzle's divergent part was analyzed in detailed.

MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Surface waves in a vertically excited circular cylindrical container

The nonlinear free surface amplitude equation, which has been derived from the inviscid fluid by solving the potential equation of water waves with a singular perturbation theory in a vertically oscillating rigid circular cylinder, is investigated successively in the fourth-order Runge-Kutta approach with an equivalent time-step. Computational results include the evolution of the amplitude with time, the characteristics of phase plane determined by the real and imaginary parts of the amplitude, the single-mode selection rules of the surface waves in different forced frequencies, contours of free surface displacement and corresponding three-dimensional evolution of surface waves, etc. In addition, the comparison of the surface wave modes is made between theoretical calculations and experimental measurements, and the results are reasonable although there are some differences in the forced frequency.

Effects of temperature on the threshold of phosphorus for algal blooms

Algal blooms, worsening marine ecosystems and causing great economic loss, have been paid much attention to for a long time. Such environmental factors as light penetration, water temperature, and nutrient concentration are crucial in blooms processes. Among them, only nutrients can be controlled. Therefore, the threshold of nutrients for algal blooms is of great concern. To begin with, a dynamic eutrophication model has been constructed to simulate the algal growth and phosphorus cycling. The model encapsulates the essential biological processes of algal growth and decay, and phosphorus regeneration due to algal decay. The nutrient limitation is based upon commonly used Monod's kinetics. The effects of temperature and phosphorus limitation are particularly addressed. Then, we have endeavored to elucidate the threshold of phosphorus at different temperature for algal blooms. Based on the numerical simulation, the isoquant contours of change rate of alga as shown in the figure are obtained, which obviously demonstrate the threshold of nutrient at an arbitrary reasonable temperature. The larger the change rate is, the more rapidly the alga grows. If the phosphorus concentration at a given temperature remains larger than the threshold the algal biomass may increase monotonically, leading to the algal blooming. With the rising of temperature, the threshold is apparently reduced, which may explain why likely red tide disasters occur in a fine summer day. So, high temperature and sufficient phosphorus supply are the major factors which result in algal growth and blowout of red tide.

Space-time correlations of fluctuating velocities in turbulent shear flows

Space-time correlations or Eulerian two-point two-time correlations of fluctuating velocities are analytically and numerically investigated in turbulent shear flows. An elliptic model for the space-time correlations in the inertial range is developed from the similarity assumptions on the isocorrelation contours: they share a uniform preference direction and a constant aspect ratio. The similarity assumptions are justified using the Kolmogorov similarity hypotheses and verified using the direct numerical simulation DNS of turbulent channel flows. The model relates the space-time correlations to the space correlations via the convection and sweeping characteristic velocities. The analytical expressions for the convection and sweeping velocities are derived from the Navier-Stokes equations for homogeneous turbulent shear flows, where the convection velocity is represented by the mean velocity and the sweeping velocity is the sum of the random sweeping velocity and the shearinduced velocity. This suggests that unlike Taylor’s model where the convection velocity is dominating and Kraichnan and Tennekes’ model where the random sweeping velocity is dominating, the decorrelation time scales of the space-time correlations in turbulent shear flows are determined by the convection velocity, the random sweeping velocity, and the shear-induced velocity. This model predicts a universal form of the spacetime correlations with the two characteristic velocities. The DNS of turbulent channel flows supports the prediction: the correlation functions exhibit a fair good collapse, when plotted against the normalized space and time separations defined by the elliptic model.

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.

Numerical Study on Critical Wedge Angle of Cellular Detonation Reflections

The critical wedge angle (CWA) for the transition from regular reflection (RR) to Mach reflection (MR) of a cellular detonation wave is studied numerically by an improved space-time conservation element and solution element method together with a two-step chemical reaction model. The accuracy of that numerical way is verified by simulating cellular detonation reflections at a 19.3∘ wedge. The planar and cellular detonation reflections over 45∘–55∘ wedges are also simulated. When the cellular detonation wave is over a 50∘ wedge, numerical results show a new phenomenon that RR and MR occur alternately. The transition process between RR and MR is investigated with the local pressure contours. Numerical analysis shows that the cellular structure is the essential reason for the new phenomenon and the CWA of detonation reflection is not a certain angle but an angle range.

透镜自重变形引起波像差的有限元分析

在用大口径、长焦距平行光管模拟激光远场特性时，其光组透镜在重力作用下的变形不能忽略，为了分析对出射光束质量的影响，采用有限元分析软件“ANSYS”建立了平行光管光组中声400mm平凸透镜的有限元模型，给出一种分析透镜轴向变形引起的波像差的方法，在不同工况下，计算了平凸透镜在重力作用下轴向变形的峰谷值和均方根值，对轴向变形量均方根值最小工况画出了透镜表面变形的等值线图，计算了声350mm通光口径内的波像差峰谷值和均方根值，对平行光管光组的波像差做出估计，验证了设计的合理性。

Absolute left-handed behaviors in a triangular elliptical-rod photonic crystal

In this paper we theoretically study the left-handed behaviors in a two-dimensional triangular photonic crystal made of elliptical rods in air. An absolute left-handed region is found in the second photonic band by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Typical left-handed behaviors such as negative refraction, flat superlensing and plano-concave lensing are demonstrated by the finite-difference time-domain simulations. These behaviors are also compared with the quasi-negative refraction and the resulted focusing effects in a square-lattice two-dimensional photonic crystal. (c) 2005 Optical Society of America

Scattering by two degenerate anisotropic modes in square-lattice dielectric photonic crystals

We systematically investigate the square-lattice dielectric photonic crystals that have been used to demonstrate flat slab imaging experimentally. A right-handed Bloch mode is found in the left-handed frequency region by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Using the multiple scattering theory, numerical simulations demonstrate that the left-handed mode and the right-handed mode are excited simultaneously by a point source and result in two kinds of transmitted waves. Impacted by the evanescent waves, superposition of these transmitted waves brings on complicated near field distributions such as the so-called imaging and its disappearance.

Negative refraction in a honeycomb-lattice photonic crystal

In this paper we theoretically investigate a photonic crystal with dielectric rods in a honeycomb lattice. Two left-handed frequency regions are found in the second and third photonic band by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Subwavelength imaging by the photonic crystal flat lens are systematically studied by numerical simulations using the multiple scattering method. Different from the photonic crystals with noncircular dielectric rods in air, this structure is almost isotropic at the optimal frequency for superlensing. As a comparison, flat slab focusing is also demonstrated at other frequencies in the two left-handed regions. (c) 2006 Elsevier Ltd. All rights reserved.

基于滤波方法的DEM尺度变换方法研究

针对水土流失评价中地形因子的尺度效应问题,引入数字图像处理中的滤波分析方法,对5 m分辨率DEM进行了滤波处理,并通过地形综合特征分析,坡度、曲率制图统计分析,派生等高线与河流提取统计分析等方法,对滤波处理的综合效果进行了分析论证。结果表明,利用滤波方法可以实现对高分辨率DEM的综合,利用比例尺因子为5或者10的处理,其结果与10 m分辨率DEM或1∶5万地形图在反映地形起伏方面相似。

ANUDEM和TIN两种建立DEM方法的对比研究

针对流域水文和土壤侵蚀定量模拟分析需要,对基于TIN和Hutchinson方法建立的DEM表现地形形态和起伏的能力进行了对比研究。结果表明,基于TIN建立的DEM始终存在一些平顶现象,一些较小的侵蚀沟被忽略,其上提取的河流不完全连续,多处出现多重线条河流,因而不能如实地反映地形起伏的细部特征。而基于ANUDEM建立的DEM,其派生等高线的形状与输入等高线吻合较好,较好地表现了地形的形态和起伏,对地形和坡度的反映更加连续和光滑,其上提取的河流信息基本与地形图上的河流一致。

Short-term morphological changes after sand-digging during a sand sculpture activity

The morphology of the beach backshore and foreshore at Huiquan Bay, Qingdao, China, is characterized by a single intertidal sandbar system with a spring tide range of 4.59 m. The beach was measured with a laser total station of Leica TPS402. Contours of the beach were generated using data collected in March and November 2005. The survey method provided 2 mm measuring accuracy and 4-10 m horizontal spacing. The net accretion volume of the foreshore was about 11, 215 m(3) from March to November. After sand sculpture activity, the axis of the sand trough migrated onshore from about 3.5 m to 17.5 m on the foreshore beach in November. At the same time, the axis of the sandbar crest migrated onshore no more than 42.25 m on the northwest foreshore; and it migrated offshore no more than 23.75 m on the southeast foreshore. On the northwest and southeast foreshore beach, two strips of erosion areas with a thickness of 0-0.2 m appeared on the sandbar crest. Accretion occurred at the bottom of the sand trough with a thickness of similar to 0.2-0.6 m. The sandbar height decreased after sand sculpture activity, and it was no more than 0.7 m in March and 0.6 m in November. Human activities, such as sand digging on the sandbar crest during sand sculpture activity, also can disturb the beach morphology of intertidal bar systems. This phenomenon also was validated by comparison of beach morphology, the results of a color artificial tracer experiment and a sediment transportation trend prediction.

Impact of Brunswick river mouth training walls on adjacent beaches, Brunswick heads, New South Wales, Australia

River training walls have been built at scores of locations along the NSW coast and their impacts on shoreline change are still not fully understood. In this study, the Brunswick River entrance and adjacent beaches are selected for examination of the impact of the construction of major training walls. Thirteen sets of aerial photographs taken between 1947 and 1994 are used in a CIS approach to accurately determine tire shoreline Position, beach contours and sand volumes, and their changes in both time and space, and then to assess the contribution of both tire structures and natural hydrodynamic conditions to large scale (years-decades and kilometres) beach changes. The impact of the training walls can be divided into four stages: natural conditions prior to their construction (pre 1959), major downdrift erosion and updrift accretion during and. following the construction of the walls in 1959 similar to 1962 and 1966. diminishing impact of the walls between 1966 and 1987, and finally no apparent impact between 1987 similar to 1994. The impact extends horizontally about 8 km updrift and 17 km downdrift, and temporally up to 25 years..