The Millimeter Wave Tomography Application for the Subsurface Imaging

In the present paper new results of modeling and experimental investigation on millimeter wave subsurface tomography are submitted. Tomographic algorithm is employed for imaging of subsurface objects in the case when conductivity of probed medium is not equal to zero. The possibilities and restrictions of this algorithm for image processing are shown both as a result of modeling and as a result of experiments. A new tomography setup allowing obtaining images of different inhomogeneities in dielectric media is considered.     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

An efficient sampling method for stochastic inverse problems

A general framework is developed to treat inverse problems with parameters that are random fields. It involves a sampling method that exploits the sensitivity derivatives of the control variable with respect to the random parameters. As the sensitivity derivatives are computed only at the mean values of the relevant parameters, the related extra cost of the present method is a fraction of the total cost of the Monte Carlo method. The effectiveness of the method is demonstrated on an example problem governed by the Burgers equation with random viscosity. It is specifically shown that this method is two orders of magnitude more efficient compared to the conventional Monte Carlo method. In other words, for a given number of samples, the present method yields two orders of magnitude higher accuracy than its conventional counterpart.     

On the covariance of the level sizes in random recursive trees

In this paper we study the covariance structure of the number of nodes k and l steps away from the root in random recursive trees. We give an analytic expression valid for all k, l and tree sizes N. The fraction of nodes k steps away from the root is a random probability distribution in k. The expression for the covariances allows us to show that the total variation distance between this (random) probability distribution and its mean converges in probability to zero. © 2002 Wiley Periodicals, Inc. Random Struct. Alg., 20: 519–539, 2002     

Possible triggering of solar activity to big earthquakes (Ms ≥ 8) in faults with near west-east strike in China

A big earthquake (Ms = 8.1) occurred at western Kunlun Mts gap on November 14, 2001. It is the biggest earthquake in China in the past 50 years and the only one with a magnitude bigger than 8 in China after Dangxiong Earthquake (M = 8) in Tibet in 1951. The big earthquake has drawn wide attention of scientists in seismological, geophysical, geological and astro-geodynamic circles. Scientists are interested in the occurrence regularity of the very big earthquakes and try to find out the fac…     

Fastiterative methods for least squares estimations

Least squares estimations have been used extensively in many applications, e.g. system identification and signal prediction. When the stochastic process is stationary, the least squares estimators can be found by solving a Toeplitz or near-Toeplitz matrix system depending on the knowledge of the data statistics. In this paper, we employ the preconditioned conjugate gradient method with circulant preconditioners to solve such systems. Our proposed circulant preconditioners are derived from the spectral property of the given stationary process. In the case where the spectral density functions() of the process is known, we prove that ifs() is a positive continuous function, then the spectrum of the preconditioned system will be clustered around 1 and the method converges superlinearly. However, if the statistics of the process is unknown, then we prove that with probability 1, the spectrum of the preconditioned system is still clustered around 1 provided that large data samples are taken. For finite impulse response (FIR) system identification problems, our numerical results show that annth order least squares estimator can usually be obtained inO(n logn) operations whenO(n) data samples are used. Finally, we remark that our algorithm can be modified to suit the applications of recursive least squares computations with the proper use of sliding window method arising in signal processing applications.Research supported in part by HKRGC grant no. 221600070, ONR contract no. N00014-90-J-1695 and DOE grant no. DE-FG03-87ER25037.     

The covariance matrix of the Potts model: A random cluster analysis

We consider the covariance matrix,G ^mm =q ²<(_x,m);(_y,m)>, of thed-dimensionalq-states Potts model, rewriting it in the random cluster representation of Fortuin and Kasteleyn. In any of theq ordered phases, we identify the eigenvalues of this matrix both in terms of representations of the unbroken symmetry group of the model and in terms of random cluster connectivities and covariances, thereby attributing algebraic significance to these stochastic geometric quantities. We also show that the correlation length corresponding to the decay rate of one of the eigenvalues is the same as the inverse decay rate of the diameter of finite clusers. For dimensiond=2, we show that this correlation length and the correlation length of the two-point function with free boundary conditions at the corresponding dual temperature are equal up to a factor of two. For systems with first-order transitions, this relation helps to resolve certain inconsistencies between recent exact and numerical work on correlation lengths at the self-dual point _o. For systems with second order transitions, this relation implies the equality of the correlation length exponents from above and below threshold, as well as an amplitude ratio of two. In the course of proving the above results, we establish several properties of independent interest, including left continuity of the inverse correlation length with free boundary conditions and upper semicontinuity of the decay rate for finite clusters in all dimensions, and left continuity of the two-dimensional free boundary condition percolation probability at _o. We also introduce DLR equations for the random cluster model and use them to establish ergodicity of the free measure. In order to prove these results, we introduce a new class of events which we call decoupling events and two inequalities for these events. The first is similar to the FKG inequality, but holds for events which are neither increasing nor decreasing; the second is similar to the van den Berg-Kesten inequality in standard percolation. Both inequalities hold for an arbitrary FKG measure.     

Dimer asymmetry in superoxide dismutase studied by molecular dynamics simulation

Summary Molecular dynamics (MD) simulations of 100 ps have been carried out to study the active-site behaviour of the Cu,Zn superoxide dismutase dimer (SOD) in water. The active site of each subunit was monitored during the whole simulation by calculating the distances between functional residues and the catalytic copper. The results indicate that charge orientation is maintained at each active site but the solvent accessibility varies. Analysis of the MD simulation, carried out by using the atomic displacement covariance matrix, has shown a different intra-subunit correlation pattern for the two monomers and the presence of inter-subunit correlations. The MD simulation presented here indicates an asymmetry in the two active sites and different dynamic behaviour of the two SOD subunits.     

Using symbolic computation software packages in production of multidimensional finite volume‐based large eddy simulation codes

The numerical simulation of turbulence is one of the most challenging tasks in the field of the modern computational science. At present, the most advanced approach is the large eddy simulation (LES) technique wherein a formal separation between resolved (large) and unresolved (small) scales of the motion is in effect by means of a filtering operation applied onto the governing equations. However, LES requires very sophisticated numerical discretizations in terms of both accuracy and efficiency. Often, the modelling of the unresolved subgrid scale terms adds further computational complexities. This paper illustrates the suitability in using software packages for symbolic computation (in the present case, Maple© for helping in the production of subroutines for a new multidimensional, high‐order accurate finite volume‐based LES code. Specifically, it will be detailed how producing, rapidly and efficiently, the routines for computing convective, diffusive as well as subgrid scale modelling fluxes. It is particularly detailed how exploiting the package for differential calculus and linear algebra for the analytical integration of the flux polynomials over the finite volume faces. The structure of the LES code is illustrated, and an accuracy analysis of the local truncation errors is performed comparing the third‐order accurate multidimensional upwind and the classical second‐order centred reconstruction in the wavenumbers space. Then, some numerical results for the turbulent plane channel and some brief points concerning the parallelization issue are addressed. Copyright © 2012 John Wiley & Sons, Ltd.     

Toward low‐noise synthetic turbulent inflow conditions for aeroacoustic calculations

The accuracy of boundary conditions for computational aeroacoustics is a well‐known challenge, due in part to the necessity of truncating the flow domain and replacing the analytical boundary conditions at infinity with numerical boundary conditions. In particular, the inflow boundary condition involving turbulent velocity or scalar fields is likely to introduce spurious waves into the domain, therefore degrading the flow behavior and deteriorating the physical acoustic waves. In this work, a method to generate low‐noise, divergence‐free, synthetic turbulence for inflow boundary conditions is proposed. It relies on the classical view of turbulence as a superposition of random eddies convected with the mean flow. Within the proposed model, the vector potential and the requirement that the individual eddies must satisfy the linearized momentum equations about the mean flow are used. The model is tested using isolated eddies convected through the inflow boundary and an experimental benchmark data for spatially decaying isotropic turbulence. Copyright © 2013 John Wiley & Sons, Ltd.