Fast monostatic GPR modeling

We propose the exploding-reflector method to simulate a monostatic survey with a single simulation. The exploding reflector, used in seismic modeling, is adapted for ground-penetrating radar (GPR) modeling by using the analogy between acoustic and electromagnetic waves. The method can be used with ray tracing to obtain the location of the interfaces and estimate the properties of the medium on the basis of the traveltimes and reflection amplitudes. In particular, these can provide a better estimation of the conductivity and geometrical details. The modeling methodology is complemented with the use of the plane-wave method. The technique is illustrated with GPR data from an excavated tomb of the nineteenth century.

Charge disproportionation in CaCu3Fe4O12

First-principles calculations using the APW+lo method, as implemented in the WIEN2K code, have been used to investigate the structural, electronic, and magnetic properties of the perovskite CaCu3Fe4O12, including the high-temperature Im-3 and low-temperature Pn-3 phase. The high-temperature phase presents a homogeneous valence and an orbital degenerate half-metallic behavior, which is consistent with the previous theoretical result. Instead orbital ordering, charge ordering, or disproportionation on Fe sites occur in the low-temperature phase, leading to the insulating character. More importantly, the charge disproportionation is of 2d(5)L -> d(5)L(2) + d(5) type (where L denotes an oxygen hole or a ligand hole), and the origin for the phenomenon is discussed in detail.

利用台阵面波资料研究华北克拉通北缘燕山褶皱带及其邻区的地壳与上地幔结构

The continent of eastern China, especially the North China Craton (NCC), has endured intensive tectonic renovation during Mesozoic and Cenozoic, with the presence of widespread magmatism, high heat flow and development of large sedimentary basins and mountain ranges. The cratonic lithosphere of the region has been destroyed remarkably, which is characterized by not only a significant reduction in thickness but also complex modifications in physical and chemical properties of the lithosphere. As for the tectonic regime controlling the evolution of the NCC, various models have been put forward, including the impingement of mantle plumes (“mushroom cloud” model), the collision of south China block and north China block, the subduction of the Pacific plate, etc. Lithosphere delamination and thermal erosion were proposed as the two end-member mechanisms of the lithospheric thinning. However, given the paucity of the data, deep structural evidence is currently still scarce for distinguishing and testifying these models. To better understand the deep structure of the NCC, from 2000 to the present, temporary seismic array observations have been conducted in the NCC by the Seismological Laboratory of the Institute of the Geology and Geophysics, Chinese Academy of Sciences under the North China Interior Structure Project (NCISP). Many arrays extend from the North China Craton and the off-craton regions, and traverse a lot of main tectonic boundaries. A total of more than 300 broadband seismic stations have been deployed along several profiles that traversed the major tectonic units within the craton’s interior, at the boundary areas and in the neighboring off-craton regions. These stations recorded abundant high-quality data, which provides an unprecedented opportunity for us to unravel the deep structural features of the NCC using seismological methods. Among all the seismological methods, the surface wave method appears to be an efficient and widely adopted technique in studying the crustal and upper mantle structures. In particular, it can provide the absolute values of S-wave velocity that are difficult to obtain with other methods. Benefiting from the deployment of dense seismic arrays, progresses have been made in improving the spatial resolution of surface wave imaging, which makes it possible to resolve the fine-scale velocity structures of the crust and upper mantle based on surface wave analysis. Meanwhile, the differences in the S-wave velocities derived from Rayleigh and Love wave data can provide information on the radial anisotropy beneath the seismic arrays. In this thesis, using the NCISP-III broadband data and based on phase velocity dispersion analysis and inversion of fundamental mode Rayleigh and Love waves, I investigated the lateral variations in the S-wave velocity structure of the crust and uppermost mantle beneath the Yanshan Belt and adjacent regions at the northeastern boundary of the NCC. Based on the constructed structural images, I discussed possible deep processes of the craton destruction in the study region.

A New Version of Smoothed Point Method to Simulate Linear Elastic Wave Propagation

By using the kernel function of the smoothed particle hydrodynamics (SPH) and modification of statistical volumes of the boundary points and their kernel functions, a new version of smoothed point method is established for simulating elastic waves in solid. With the simplicity of SPH kept, the method is easy to handle stress boundary conditions, especially for the transmitting boundary condition. A result improving by de-convolution is also proposed to achieve high accuracy under a relatively large smooth length. A numerical example is given and compared favorably with the analytical solution.

Improved Analysis Method for Wave-Induced Pipeline Stability on Sandy Seabed

The existing Det Norske Veritas DNV Recommended Practice RP E305 for pipeline on-bottom stability is mainly based on the pipe–soil interaction model reported by Wagner et al. in 1987, and the wake model reported by Lambrakos et al. in 1987, to calculate the soil resistance and the hydrodynamic forces upon pipeline, respectively. Unlike the methods in the DNV Practice, in this paper, an improved analysis method is proposed for the on-bottom stability of a submarine pipeline, which is based on the relationships between Um/ gD 0.5 and Ws / D2 for various restraint conditions obtained by the hydrodynamic loading experiments, taking into account the coupling effects between wave, pipeline, and sandy seabed. The analysis procedure is illustrated with a detailed flow chart. A comparison is made between the submerged weights of pipeline predicted with the DNV Practice and those with the new method. The proposed analysis method may provide a helpful tool for the engineering practice of pipeline on-bottom stability design.

Optimization of tried wave-function in quantum monte-carlo method

A method for optimizing tried wave functions in quantum Monte Carlo method has been found and used to calculate the energies of molecules, such as H-2, Li-2, H-3+, H-3 and H-4. Good results were obtained.

A shock tube method in the kinetic study of nonequilibrium ionization-recombination by rarefaction wave cooling

A shock tube method is introduced to study the ionization–recombination kinetics of high temperature gas, in which a test gas is heated and ionized by a reflected shock wave and subsequently quenched by a strong rarefaction wave reflected on the end wall of the driver section as the main cooling wave associated with a rarefaction wave incident back into region 5 when the reflected shock wave interacts with the contact surface. As the quenching rate of the strong rarefaction wave reaches 106 K/s, a nonequilibrium ionization-recombination process occurs, during which the ion recombination with electrons dominates.

Wave-front analysis method of circular aperture sampling for collimation testing

A new type of wave-front analysis method for the collimation testing of laser beams is proposed. A concept of wave-front height is defined, and, on this basis, the wave-front analysis method of circular aperture sampling is introduced. The wave-front height of the tested noncollimated wave can be estimated from the distance between two identical fiducial diffraction planes of the sampled wave, and then the divergence is determined. The design is detailed, and the experiment is demonstrated. The principle and experiment results of the method are presented. Owing to the simplicity of the method and its low cost, it is a promising method for checking the collimation of a laser beam with a large divergence. © 2005 Optical Society of America.

Calculation of Valence Subband Structures for Strained Quantum-Wells by Plane Wave Expansion Method Within 6 * 6 Luttinger-Kohn Model

The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6 * 6 Luttinger-Kohn model. The effect of the number and period of plane-waves used for expansion on the stability of energy eigenvalues is examined. For practical calculation, it should choose the period large sufficiently to ensure the envelope functions vanish at the boundary and the number of plane waves large enough to ensure the energy eigenvalues keep unchanged within a prescribed range.

An augmented method for free boundary problems with moving contact lines

An augmented immersed interface method (IIM) is proposed for simulating one-phase moving contact line problems in which a liquid drop spreads or recoils on a solid substrate. While the present two-dimensional mathematical model is a free boundary problem, in our new numerical method, the fluid domain enclosed by the free boundary is embedded into a rectangular one so that the problem can be solved by a regular Cartesian grid method. We introduce an augmented variable along the free boundary so that the stress balancing boundary condition is satisfied. A hybrid time discretization is used in the projection method for better stability. The resultant Helmholtz/Poisson equations with interfaces then are solved by the IIM in an efficient way. Several numerical tests including an accuracy check, and the spreading and recoiling processes of a liquid drop are presented in detail. (C) 2010 Elsevier Ltd. All rights reserved.

A New High Order Accurate Shock Capture Method with Wave Booster

A new kind of shock capturing method is developed. Before applying the high order accurate traditional scheme which is called as base scheme in this paper the fluid parameters are preconditioned in order to control the group velocity. The newly constructed scheme is high order accurate, simple, has high resolution of the shock, and less computer time consumed.

A parametric method of retrieving ocean wave spectra from synthetic aperture radar images

A parametric method that extracts the ocean wave directional spectra from synthetic aperture radar (SAR) image is presented. The 180 degrees ambiguity of SAR image and the loss of information beyond the azimuthal cutoff can be overcome with this method. The ocean wave spectra can be obtained from SAR image directly by using iteration inversion mapping method with forward nonlinear mapping. Some numerical experiments have been made by using ERS-1 satellite SAR imagette data. The ocean wave direction retrieved from SAR imagette data is in agreement with the wind direction from the scatterometer data.