Microoptical research in Shanghai Institute of Optics and Fine Mechanics

This paper summarized the recent research results of Changhe Zhou's group of Information Optics Lab in Shanghai Institute of Optics and Fine Mechanics (SIOM). The first is about the Talbot self-imaging research. We have found the symmetry rule, the regular-rearranged neighboring phase difference rule and the prime-number decamping rule, which is briefly summarized in a recent educational publication of Optics and Photonics News, pp.46-50, November 2004. The second is about four novel microoptical gratings designed and fabricated in SIOM. The third is about the design and fabrication of novel supperresolution phase plates for beam shaping and possible use in optical storage. The fourth is to develop novel femtosecond laser information processing techniques by incorporating microoptical elements, for example, use of a pair of reflective Dammann gratings for splitting the femtosecond laser pulses. The most attractive feature of this approach is that the conventional beam splitter is avoided. The conventional beam splitter would introduce the unequal dispersion due to the broadband spectrum of ultrashort laser pulses, which will affect the splitting result. We implemented the Dammann splitting apparatus by using two-layered reflective Dammann gratings, which generates the almost same array without angular dispersion. We believe that our device is highly interesting for splitting femtosecond laser pulses.

High-accuracy same-beam VLBI observations using Shanghai and Urumqi telescopes

The same-beam VLBI observations of Rstar and Vstar, which were two small satellites of Japanese lunar mission, SELENE, were successfully performed by using Shanghai and Urumqi 25-m telescopes. When the separation angle between Rstar and Vstar was less than 0.1 deg, the differential phase delay of the X-band signals between Rstar and Vstar on Shanghai-Urumqi baseline was obtained with a very small error of 0.15 mm rms, which was reduced by 1-2 order compared with the former VLBI results. When the separation angle was less than 0.56 deg, the differential phase delay of the S-band signals was also obtained with a very small error of several mm rms. The orbit determination for Rstar and Vstar was performed, and the accuracy was improved to a level of several meters by using VLBI and Doppler data. The high-accuracy same-beam differential VLBI technique is very useful in orbit determination for a spacecraft, and will be used in orbit determination for Mars missions of China Yinghuo-1 and Russia Phobos-grunt.

Research progress on high-enthalpy and hypersonic flows

The research progress on high-enthalpy and hypersorlic flows having been achieved in the Institute of Mechanics, Chinese Academy of Sciences, is reported in this paper. The paper consists of three main parts: The first part is on the techniques to develop advanced hypersonic test facilities, in which the detonation-driven shock-reflected tunnel and the detonation-driven shock-expanded tube are introduced. The shock tunnel can be used for generating hypersonic flows of a Mach number ranging from 10 to 20, and the expansion tube is applicable to simulate the flows with a speed of 7 similar to 10km/s. The second part is dedicated to the shock tunnel nozzle flow diagnosis to examine properties of the hypersonic flows thus created. The third part is on experiments and numerical simulations. The experiments include measuring the aerodynamic pitching moment and heat transfer in hypersonic flows, and the numerical work reports nozzle flow simulations and flow non-equilibrium effects on the possible experiments that may be carried out on the above-mentioned hypersonic test facilities.

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

Analytical calculations of Eulerian and Lagrangian time correlations in turbulent shear flows

The Taylor series expansion method is used to analytically calculate the Eulerian and Lagrangian time correlations in turbulent shear flows. The short-time behaviors of those correlation functions can be obtained from the series expansions. Especially, the propagation velocity and sweeping velocity in the elliptic model of space-time correlation are analytically calculated and further simplified using the sweeping hypothesis and straining hypothesis. These two characteristic velocities mainly determine the space-time correlations.

DNS and LES of turbulent channel flow with hydrophobic surface

Hydrophobic surface benefits for drag reduction. Min and Kim[1] do the first Direct Numerical Simulation on drag reduction in turbulent channel flow. And Fukagata and Kasagi[2] make some theoretical analysis based on Dean[3]'s formula and some observations in the DNS results. Using their theory, they conclude that drag reduction is possible in large Reynolds number. Both Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed in our research. How the LES behaving in the turbulent channel flow with hydrophobic surface is examined. Original Smagorinsky model and its Dynamical model are used in LES. The slip velocities predicted by LES using Dynamical model are in good agreement with DNS as shown in the Figure. Although the percentage of drag reduction predicted by LES shows some discrepancies, it is in the error limit for industrial flow. First order and second order moments of LES are also examined and compared with DNS's results. The first-order moments is calculated well by LES. But there are some discrepancies of second-order moments between LES and DNS. [GRAPHICS]

Receptivity to free-stream disturbance waves for blunt cone axial symmetry hypersonic boundary layer

Based on high-order compact upwind scheme, a high-order shock-fitting finite difference scheme is studied to simulate the generation of boundary layer disturbance waves due to free-stream waves. Both steady and unsteady flow solutions of the receptivity problem are obtained by resolving the full Navier-Stokes equations. The interactions of bow-shock and free-stream disturbance are researched. Direct numerical simulation (DNS) of receptivity to free-stream disturbances for blunt cone hypersonic boundary layers is performed.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

Wave loading on floating platforms by internal solitary waves

Morison's equation is used for estimating internal solitary wave-induced forces exerted on SPAR and semi-submersible platforms. And the results we got have also been compared to ocean surface wave loading. It is shown that Morison's equation is an appropriate approach to estimate internal wave loading even for SPAR and semi-submersible platforms, and the internal solitary wave load on floating platforms is comparable to surface wave counterpart. Moreover, the effects of the layers with different thickness on internal solitary wave force are investigated.