On the computation of space-time correlations by large-eddy simulation

The effect of subgrid-scale (SGS) modeling on velocity (space-) time correlations is investigated in decaying isotropic turbulence. The performance of several SGS models is evaluated, which shows superiority of the dynamic Smagorinsky model used in conjunction with the multiscale large-eddy simulation (LES) procedure. Compared to the results of direct numerical simulation, LES is shown to underpredict the (un-normalized) correlation magnitude and slightly overpredict the decorrelation time scales. This can lead to inaccurate solutions in applications such as aeroacoustics. The underprediction of correlation functions is particularly severe for higher wavenumber modes which are swept by the most energetic modes. The classic sweeping hypothesis for stationary turbulence is generalized for decaying turbulence and used to analyze the observed discrepancies. Based on this analysis, the time correlations are determined by the wavenumber energy spectra and the sweeping velocity, which is the square root of the total energy. Hence, an accurate prediction of the instantaneous energy spectra is most critical to the accurate computation of time correlations. (C) 2004 American Institute of Physics.

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.

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.

地质体科学计算可视化研究及其应用

Describing visually space-time properties of geological phenomena consists of one of the most important parts in geology research. Such visual images are of usually helpful for analyzing geological phenomena and for discovering the regulations behind geological phenomena. This report studies mainly three application problems of scientific visualization in geology: (Dvisualizing geological body A new geometric modeling technique with trimmed surface patches has been eveloped to visualize geological body. Constructional surfaces are represented as trimmed surfaces and a constructional solid is represented by the upper and lower surface composed of trimmed surface patches from constructional surfaces. The technique can completely and definitely represent the structure of geological body. It has been applied in visualization for the coal deposit in Huolinhe, the aquifer thermal energy storage in Tianjin and the structure of meteorite impact in Cangshan et al. (2)visualizing geological space field Efficient visualization methods have been discussed. Marching-Cube algorithm used has been improved and is used to extract iso~surface from 3D data set, iso-line from 2D data set and iso-point from ID data set. The improved method has been used to visualize distribution and evolution of the abnormal pressures in Zhungaer Basin. (3)visualizing porous space a novel way was proposed to define distance from any point to a convex set. Thus a convex set skeleton-based implicit surface modeling technique is developed and used to construct a simplified porous space model. A Buoyancy Percolation numerical simulation platform has been developed to simulate the process of migration of oil in the porous media saturated with water.

An Improved Ce/Se Scheme And Its Application To Detonation Propagation

An improved two-dimensional space-time conservation element and solution element ( CE/ SE) method with second-order accuracy is proposed, examined and extended to simulate the detonation propagations using detailed chemical reaction models. The numerical results of planar and cellular detonation are compared with corresponding results by the Chapman-Jouguet theory and experiments, and prove that the method is a new reliable way for numerical simulations of detonation propagation.

Wavelike patterns in one-dimensional coupled map lattices

We investigate the existence of wavelike solution for the logistic coupled map lattices for which the spatiotemporal periodic patterns can be predicted by a simple two-dimensional mapping. The existence of such wavelike solutions is proved by the implicit function theorem with constraints. We also examine the stabilities of these wave solutions under perturbations of uniform small deformation type. We show that in some specific cases these perturbations are completely general. The technique used in this paper is also applicable to investigate other space-time regular patterns.

Experimental Study and Numerical Simulation of Cellular Structures and Mach Reflection of Gaseous Detonation Wave

In this paper the Deflagration to Detonation Transition (DDT) process of gaseous H-2-O-2 mixture and Mach reflection of gaseous detonation wave on a wedge have been conducted experimentally. The cellular pattern of DDT process and Mach reflection were obtained from experiments with wedge angle theta = 10(0) similar to 40(0) and initial pressure of gaseous mixture 16kPa similar to 26.7kPa. The 2-D numerical simulations of DDT process and Mach reflection of detonation wave were performed by using the simplified ZND model and improved space-time conservation element and solution element (CE/SE) method. The numerical cellular structures were compared with the cellular patterns of soot track. Compared results were shown that it is satisfactory. The characteristic comparisons on Mach reflection of air shock wave and detonation wave were carried also out and their differences were given.

An Improved Ce/Se Scheme For Multi-Material Elastic-Plastic Flows And Its Applications

In this paper, a new definition of SE and CE, which is based on the hexahedron mesh and simpler than Chang's original CE/SE method (the space-time Conservation Element and Solution Element method), is proposed and an improved CE/SE scheme is constructed. Furthermore, the improved CE/SE scheme is extended in order to solve the elastic-plastic flow problems. The hybrid particle level set method is used for tracing the interfaces of materials. Proper boundary conditions are presented in interface tracking. Two high-velocity impact problems are simulated numerically and the computational results are carefully compared with the experimental data, as well as the results from other literature and LS-DYNA software. The comparisons show that the computational scheme developed currently is clear in physical concept, easy to be implemented and high accurate and efficient for the problems considered. (C) 2008 Elsevier Ltd. All rights reserved.

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

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.

对数值模拟气相爆轰中二阶段化学反应模型再研究

采用改进的高精度时-空守恒元解元算法(the space-time conservation element and solution element method,CE/SE method)和考虑组分的二阶段化学反应模型(Sichel的二步模型)对气相爆轰问题的数值模拟进行了分析.分析发现采用Sichel的二步模型得到的数值结果虽然比早期二阶段化学反应模型(旧二步模型)更接近实验值,但是仍然不能得到爆轰过程准确气体动力学参数.为此通过修改组分的质量分数分布形式对Sichel的二步模型进行了改造,然后采用新的二步模型对平面爆轰波进行了数值模拟.数值结果表明采用新的二步模型计算得到气体动力学参数更接近于实验值和基元反应模型的计算值,在计算精度上有较大提高.

A CE/SE Scheme for Flows in Porous Media and Its Applications

In this paper, a new computational scheme for solving flows in porous media was proposed. The scheme was based on an improved CE/SE method (the space-time Conservation Element and Solution Element method). We described porous flows by adopting DFB (Brinkman-Forchheimer extended Darcy) equation. The comparison between our computational results and Ghia's confirmed the high accuracy, resolution, and efficiency of our CE/SE scheme. The proposed first-order CE/SE scheme is a new reliable way for numerical simulations of flows in porous media. After investigation of effects of Darcy number on porous flow, it shows that Darcy number has dominant influence on porous flow for the Reynolds number and porosity considered.

Self-compression of femtosecond pulses in argon with a power close to the self-focusing threshold

Self-compression of femtosecond pulses in noble gases with an input power close to the self-focusing threshold has been investigated experimentally and theoretically. It is demonstrated that either multiphoton ionization (MPI) or space time focusing and self-steepening effects can induce pulse shortening, but they predominate at different beam intensities during the propagation. The latter effects play a key role in the final pulse self-compression. By choosing an appropriate focusing parameter, action distance of the space time focusing and self-steepening effects can be lengthened, which can promote a shock pulse structure with a duration as short as two optical cycles. It is also found that, for our calculation cases in which an input pulse power is close to the self-focusing threshold, either group velocity dispersion (GVD) or multiphoton absorption (MPA) has a negligible influence on pulse characteristics in the propagation process.

飞秒光谱全息技术及飞秒整形技术的应用

飞秒光谱全息是飞秒脉冲整形技术中非常重要的一种方法，它可以实现飞秒脉冲信号的记录、再现和处理。主要介绍飞秒光谱全息技术及飞秒脉冲时空变换整形在飞秒化学领域的最新应用。