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.

适用于中等耦合等离子体的碰撞算子

完全电离等离子体中，当试探粒子分布函数fα是关于试探粒子速度vα的偶函数时，导出了一个新的动力学方程的碰撞算子．该碰撞算子同时包括了大角散射（库仑近碰撞）和小角散射（库仑远碰撞）的二体碰撞的贡献，因此，该碰撞算子同时适用于弱耦合（库仑对数ln∧≥10）和中等耦合（库仑对数2≤ln∧≤10）等离子体．而且经过修改的碰撞算子和Rosenbluth势有直接的联系，当试探粒子和场粒子满足条件mα＜mβ（如电子-离子碰撞或Lorentz气体模型）和｜vα｜〉｜vβ｜时，经约化的电子-离子碰撞算子同最初的Fokker

Nonequilibrium statistical mechanics of two-dimensional turbulence

The vorticity dynamics of two-dimensional turbulence are investigated analytically, applying the method of Qian (1983). The vorticity equation and its Fourier transform are presented; a set of modal parameters and a modal dynamic equation are derived; and the corresponding Liouville equation for the probability distribution in phase space is solved using a Langevin/Fokker-Planck approach to obtain integral equations for the enstrophy and for the dynamic damping coefficient eta. The equilibrium spectrum for inviscid flow is found to be a stationary solution of the enstrophy equation, and the inertial-range spectrum is determined by introducing a localization factor in the two integral equations and evaluating the localized versions numerically.

Finite element beam propagation method for analysis of plasmonic waveguide

The basic idea of the finite element beam propagation method (FE-BPM) is described. It is applied to calculate the fundamental mode of a channel plasmonic polariton (CPP) waveguide to confirm its validity. Both the field distribution and the effective index of the, fundamental mode are given by the method. The convergence speed shows the advantage and stability of this method. Then a plasmonic waveguide with a dielectric strip deposited on a metal substrate is investigated, and the group velocity is negative for the fundamental mode of this kind of waveguide. The numerical result shows that the power flow direction is reverse to that of phase velocity.

Parallel Fully-Implicit Computation of Magnetohydrodynamics Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETS,: package. The results of convergence tests arc plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

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.

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.

Variational approach to the closure problem of turbulence theory

A new method is proposed to solve the closure problem of turbulence theory and to drive the Kolmogorov law in an Eulerian framework. Instead of using complex Fourier components of velocity field as modal parameters, a complete set of independent real parameters and dynamic equations are worked out to describe the dynamic states of a turbulence. Classical statistical mechanics is used to study the statistical behavior of the turbulence. An approximate stationary solution of the Liouville equation is obtained by a perturbation method based on a Langevin-Fokker-Planck (LFP) model. The dynamic damping coefficient eta of the LFP model is treated as an optimum control parameter to minimize the error of the perturbation solution; this leads to a convergent integral equation for eta to replace the divergent response equation of Kraichnan's direct-interaction (DI) approximation, thereby solving the closure problem without appealing to a Lagrangian formulation. The Kolmogorov constant Ko is evaluated numerically, obtaining Ko = 1.2, which is compatible with the experimental data given by Gibson and Schwartz, (1963).

Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

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

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.

龙滩工程左岸边坡监测资料分析与稳定性评价研究

On the basis of the geological analysis and rock mass toppling deformation and failure mechanism analysis of Longtan engineering left bank slope, the synthetic space-time analysis and influence factors analysis on the surface monitoring data and deep rock mass monitoring data of B-zone of left bank slope are carried on. At the same time, based on the monitoring data analysis in conjunction with the predecessor's mechanics analysis results, the deformation state of B-zone of the left bank slope is discussed and its stability is synthetically evaluated. The detailed research contents and results are as following: According to the monitoring drill histogram analysis of Longtan engineering left bank slope, numerical simulation analysis and model experimentation analysis of bedded counter-inclined steep slope, a new type of toppling deformation and failure mode is proposed, that is "up-slope warping". Then the deformation and failure mode of bedded counter-inclined steep slope is summarized as "down-slope toppling" type, "up-slope warping" type and "complex fold" type. On the basis of synthetic space-time analysis to surface monitoring data and deep rock mass deformation monitoring data of B-zone of Longtan left bank slope;, we can get the conclusion that there exists potential instability rock mass over 520m altitude, especially over 560m altitude of slope B, and the rock mass of around strong-weathering line or creep rock mass breaking band controls the deformation of the whole slope. 1. According to the synthetic space-time analysis and influence factors analysis to the surface monitoring data of B-zone of Longtan left bank slope, a dynamical index, accumulative total acceleration index, which is used to analyze the influence factors of slope surface deformation, is raised. The principle and method of accumulative acceleration index are explained, and the index can be used for the influence factors analysis of the similar slope. 2. Summarize the results of geologic analysis, monitoring analysis and mechanics analysis, the following conclusion can be gotten: the stability of B-zone of the slope is basically good. However, on the condition of drainage and slope toe loading engineering, there is still some creep deformation in the rock mass over 520m altitude, especially over 560m altitude. So, better measures of the monitoring and timely maintenance of the drainage system are suggested in the paper.

Kramers Escape Rate in Nonlinear Diffusive Media

In this paper, we study nonlinear Kramers problem by investigating overdamped systems ruled by the one-dimensional nonlinear Fokker-Planck equation. We obtain an analytic expression for the Kramers escape rate under quasistationary conditions by employing

Analysis of super compact finite difference method and application to simulation of vortex-shock interaction

Turbulence and aeroacoustic noise high-order accurate schemes are required, and preferred, for solving complex flow fields with multi-scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth-order traditional and compact finite difference approximation. The comparison shows that the sixth-order accurate super compact method has higher resolving efficiency. The sixth-order super compact method, with a three-stage Runge-Kutta method for approximation of the compressible Navier-Stokes equations, is used to solve the complex flow structures induced by vortex-shock interactions. The basic nature of the near-field sound generated by interaction is studied.

Statistical and stochastic description of microvoid evolution observed in dynamic failure of ductile materials

A brief review is presented of statistical approaches on microdamage evolution. An experimental study of statistical microdamage evolution in two ductile materials under dynamic loading is carried out. The observation indicates that there are large differences in size and distribution of microvoids between these two materials. With this phenomenon in mind, kinetic equations governing the nucleation and growth of microvoids in nonlinear rate-dependent materials are combined with the balance law of void number to establish statistical differential equations that describe the evolution of microvoids' number density. The theoretical solution provides a reasonable explanation of the experimentally observed phenomenon. The effects of stochastic fluctuation which is influenced by the inhomogeneous microscopic structure of materials are subsequently examined (i.e. stochastic growth model). Based on the stochastic differential equation, a Fokker-Planck equation which governs the evolution of the transition probability is derived. The analytical solution for the transition probability is then obtained and the effects of stochastic fluctuation is discussed. The statistical and stochastic analyses may provide effective approaches to reveal the physics of damage evolution and dynamic failure process in ductile materials.