Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.

Response amplitude and hydrodynamic force for a buoy over a convex

A buoy as an offshore structure is often placed over a convex such as a caisson or a submerged island. The hydrodynamic fluid/solid interaction becomes more complex due to the convex compared with that on the flat. Both the buoy and the convex are idealized as vertical cylinders. Linear potential theory is used to investigate the response amplitude and the hydrodynamic force for a buoy over a convex due to diffraction and radiation in water of finite depth. These are derived from the total velocity potential. A set of theoretical added mass, damping coefficient, and exciting force expressions have been proposed. Analytical results of the response amplitude and hydrodynamic force are given. Finally, the numerical results show that the effect of the convex on the response amplitude and hydrodynamic force for the buoy is ignored if the size of the convex is relatively smaller.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

复杂地质条件下高精度地震勘探技术研究与应用-以泌阳凹陷为例

In order to carry out high-precision three-dimensional "integration" for the characteristics of the secondary seismic exploration for Biyang Depression, in the implementation process, through a combination of scientific research and production, summed up high-precision seismic acquisition, processing and interpretation technologies suitable for the eastern part of the old liberated areas, achieved the following results: 1. high-precision complex three-dimensional seismic exploration technology series suitable for shallow depression Biyang block group. To highlight the shallow seismic signal, apply goal-based observing system design, trail from the small panel to receive and protect the shallow treatment of a range of technologies; to explain the use of three-dimensional visualization and coherent combination of full-body three-dimensional fine interpretation identification of the 50-100 m below the unconformity surface and its formation of about 10 meters of the distribution of small faults and improve the small block and stratigraphic unconformity traps recognition. 2. high-precision series of three-dimensional seismic exploration technology suitable for deep depression Biyang low signal to noise ratio of information. Binding model using forward and lighting technology, wide-angle observation system covering the design, multiple suppression and raise the energy of deep seismic reflection processing and interpretation of detailed, comprehensive reservoir description, such as research and technology, identified a number of different types of traps. 3. high-precision seismic exploration technology series for the southern Biyang Depression high steep three-dimensional structure. The use of new technology of seismic wave scattering theory and high-precision velocity model based on pre-stack time migration and depth migration imaging of seismic data and other high-precision processing technology, in order to identify the southern steep slope of the local structure prediction and analysis of sandstone bedrock surface patterns provide a wealth of information.

印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress,intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored.

Fast growing mode of stimulated Raman scattering in a pure three-wave process

Stimulated Raman scattering (SRS) of a relativistic laser in plasmas is studied in the framework of the standard equation set of a three-wave process. As far as every wave involved in the process is concerned, its evolution has two aspects: time-dependent amplitude and time-dependent frequency. These two aspects affect each other. Strict analysis and numerical experiment on the full three-wave equation set reveal that a fast growing mode of the instability, which could reach a balance or saturation point during a period far shorter than an estimation based on conventional analysis, could take place in a standard three-wave process without coupling with a fourth wave. This fast growing mode is found to stem from the constraint set by the background density on the amplitude of the driven Langmuir wave. The effect of various parameters on the development of the SRS instability is studied by numerical calculation of the history of the instability in different cases. (c) 2007 American Institute of Physics.

Wave propagation and the frequency domain Green's functions in viscoelastic Biot/squirt (BISQ) media

In this paper, we examine the characteristics of elastic wave propagation in viscoelastic porous media, which contain simultaneously both the Biot-flow and the squirt-flow mechanisms (BISQ). The frequency-domain Green's functions for viscoelastic BISQ media are then derived based on the classic potential function methods. Our numerical results show that S-waves are only affected by viscoelasticity, but not by squirt-flows. However, the phase velocity and attenuation of fast P-waves are seriously influenced by both viscoelasticity and squirt-flows; and there exist two peaks in the attenuation-frequency variations of fast P-waves. In the low-frequency range, the squirt-flow characteristic length, not viscoelasticity, affects the phase velocity of slow P-waves, whereas it is opposite in the high-frequency range. As to the contribution of potential functions of two types of compressional waves to the Green's function, the squirt-flow length has a small effect, and the effects of viscoelastic parameter are mainly in the higher frequency range. Crown Copyright (C) 2006 Published by Elsevier Ltd. All rights reserved.

Photon-photon scattering in a plasma channel

The Heisenberg-Euler correction due to photon-photon scattering, a still unverified quantum electrodynamics effect, on electromagnetic wave interaction inside a plasma channel is investigated theoretically. From a signal laser beam in the relativistic overdense plasma channel, photon-photon scattering can produce a detectable output beam of different frequency and polarization. (C) 2003 American Institute of Physics.

Imaging properties of coherent anti-Stokes Raman scattering microscope

The coherent anti-Stokes Raman scattering (CARS) microscope with the combination of confocal and CARS techniques is a remarkable alternative for imaging chemical or biological specimens that neither fluoresce nor tolerate labelling. CARS is a nonlinear optical process, the imaging properties of CARS microscopy will be very different from the conventional confocal microscope. In this paper, the intensity distribution and the polarization property of the optical field near the focus was calculated. By using the Green function, the precise analytic solution to the wave equation of a Hertzian dipole source was obtained. We found that the intensity distributions vary considerably with the different experimental configurations and the different specimen shapes. So the conventional description of microscope (e.g. the point spread function) will fail to describe the imaging properties of the CARS microscope.

The intensity distributions of collected signals in coherent anti-Stokes Raman scattering microscopy

Coherent anti-Stokes Raman scattering (CARS) microscopy with the combining of confocal and CARS techniques is a remarkable alternative for imaging chemical or biological specimens that neither fluoresce nor tolerate labeling. The CARS is a nonlinear optical process, the imaging properties of CARS microscopy will be very different from the conventional confocal microscopy. In this paper, we calculated the propagation of CARS signals by using the wave equation in medium and the slowly varying envelope approximation (SVEA), and find that the intensity angular distributions vary considerably with the different experimental configurations and the different specimen shapes. So the conventional description of microscopy (e.g.. the point spread function) will fail to descript the imaging properties of CARS microscopy. (c) 2004 Elsevier B.V. All rights reserved.

Attosecond X-ray pulse obtained from linear Thomson scattering

Linear Thomson scattering by a relativistic electron of a short pulse laser has been investigated by computer simulation. Under a laser field with a pulse of 33.3-fs full-width at half-maximum, and the initial energy of an electron of gamma(0) = 10, the motion of the electron is relativistic and generates an ultrashort radiation of 76-as with a photon wave length of 2.5-nm in the backward scattering. The radiation under a high relativistic energy electron has better characteristic than under a low relativistic energy electron in terms of the pulse width and the angular distribution. (c) 2005 Elsevier GrnbH. All rights reserved.

Scattering by two degenerate anisotropic modes in square-lattice dielectric photonic crystals

We systematically investigate the square-lattice dielectric photonic crystals that have been used to demonstrate flat slab imaging experimentally. A right-handed Bloch mode is found in the left-handed frequency region by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Using the multiple scattering theory, numerical simulations demonstrate that the left-handed mode and the right-handed mode are excited simultaneously by a point source and result in two kinds of transmitted waves. Impacted by the evanescent waves, superposition of these transmitted waves brings on complicated near field distributions such as the so-called imaging and its disappearance.