Extremely large-scale array(XL-array)communications can significantly improve the transmission rate,spectral efficiency,and spatial resolution,and has great potential in next-generation mobile communication networks.A...Extremely large-scale array(XL-array)communications can significantly improve the transmission rate,spectral efficiency,and spatial resolution,and has great potential in next-generation mobile communication networks.A crucial problem in XLarray communications is to determine the boundary of applicable regions of the plane wave model(PWM)and spherical wave model(SWM).In this paper,we propose new PWM/SWM demarcations for XL-arrays from the viewpoint of channel gain and rank.Four sets of results are derived for four different array setups.First,an equi-power line is derived for a point-touniform linear array(ULA)scenario,where an inflection point is found at±π6 central incident angles.Second,an equi-power surface is derived for a point-touniform planar array(UPA)scenario,and it is proved that cos2(ϕ)cos2(φ)=12 is a dividing curve,where ϕ andφdenote the elevation and azimuth angles,respectively.Third,an accurate and explicit expression of the equi-rank surface is obtained for a ULA-to-ULA scenario.Finally,an approximated expression of the equirank surface is obtained for a ULA-to-UPA scenario.With the obtained closed-form expressions,the equirank surface for any antenna structure and any angle can be well estimated.Furthermore,the effect of scatterers is also investigated,from which some insights are drawn.展开更多
A nuclear explosion in the rock mass medium can produce strong shock waves,seismic shocks,and other destructive effects,which can cause extreme damage to the underground protection infrastructures.With the increase in...A nuclear explosion in the rock mass medium can produce strong shock waves,seismic shocks,and other destructive effects,which can cause extreme damage to the underground protection infrastructures.With the increase in nuclear explosion power,underground protection engineering enabled by explosion-proof impact theory and technology ushered in a new challenge.This paper proposes to simulate nuclear explosion tests with on-site chemical explosion tests in the form of multi-hole explosions.First,the mechanism of using multi-hole simultaneous blasting to simulate a nuclear explosion to generate approximate plane waves was analyzed.The plane pressure curve at the vault of the underground protective tunnel under the action of the multi-hole simultaneous blasting was then obtained using the impact test in the rock mass at the site.According to the peak pressure at the vault plane,it was divided into three regions:the stress superposition region,the superposition region after surface reflection,and the approximate plane stress wave zone.A numerical simulation approach was developed using PFC and FLAC to study the peak particle velocity in the surrounding rock of the underground protective cave under the action of multi-hole blasting.The time-history curves of pressure and peak pressure partition obtained by the on-site multi-hole simultaneous blasting test and numerical simulation were compared and analyzed,to verify the correctness and rationality of the formation of an approximate plane wave in the simulated nuclear explosion.This comparison and analysis also provided a theoretical foundation and some research ideas for the ensuing study on the impact of a nuclear explosion.展开更多
To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algor...To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is：（1） Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.（2） Define the topology relationships and related concepts of triangular unit ray tracing.（3） The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.（4） By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point＇s traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.展开更多
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 pla...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.展开更多
The Pre-Stack Depth Migration （PSDM） method based on wavefield continuation is the most reliable method for imaging complex structure in the subsurface, although there are large computational costs and poorly adapti...The Pre-Stack Depth Migration （PSDM） method based on wavefield continuation is the most reliable method for imaging complex structure in the subsurface, although there are large computational costs and poorly adaptive geometry. Plane wave shot migration is another method to perform exact wave equation prestack imaging with high computational efficiency and without the migration aperture problem. Moreover, wavefield energy can be compensated at the target zone by controlled illumination. In this paper, plane wave shot PSDM was implemented by the control of the plane down-going wavefield and selection of number and range of the raypaths in order to optimize the imaging effect. In addition, controlled illumination techniques are applied to enhance the imaging precision of interesting areas at different depths. Numerical calculation indicates that plane wave shot imaging is a rapid and efficient method with less computational cost and easy parallel computation compared to the single-square-root operator imaging for common shot gathers and double- square-root operator imaging for common midpoint gathers.展开更多
The super-cell plane wave expansion method is employed to calculate band structures for the design of a siliconbased one-dimensional phononic crystal plate with large absolute forbidden bands. In this method, a low im...The super-cell plane wave expansion method is employed to calculate band structures for the design of a siliconbased one-dimensional phononic crystal plate with large absolute forbidden bands. In this method, a low impedance medium is introduced to replace the free stress boundary, which largely reduces the computational complexity. The dependence of band gaps on structural parameters is investigated in detail. To prove the validity of the super-cell plane wave expansion, the transmitted power spectra of the Lamb wave are calculated by using a finite element method. With the detailed computation, the band-gap of a one-dimensional plate can be designed as required with appropriate structural parameters, which provides a guide to the fabrication of a Lamb wave phononic crystal.展开更多
A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the l...A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.展开更多
The electronic density of states and band structures of doped and un-doped anatase TiO2 were studied by the Linearized Augmented Plane Wave method based on the density functional theory. The calculation shows that the...The electronic density of states and band structures of doped and un-doped anatase TiO2 were studied by the Linearized Augmented Plane Wave method based on the density functional theory. The calculation shows that the band structures of TiO2 crystals doped with transition metal atoms become narrower. Interesting, an excursion towards high energy level with increasing atomic number in the same element period could be observed after doping with transition metal atoms.展开更多
The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV wave...The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV waves in the discrete model are first obtained by means of separating the characteristic equation of the motion equation, and then used to analyse the properties of P-and SV-homogeneous, inhomogeneous waves and other types of motion in the model. The dispersion characters, cut-off frequencies of P and SV waves, the polarization drift and appendent anisotropic property of wave motion caused by the discretization are finally discussed.展开更多
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 d...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.展开更多
A theoretical analysis is made, using plane wave expansion, on how the width of the first three band gaps is influenced by filling ratio, dielectric constant ratio, and periodic width in one-dimensional photonic cryst...A theoretical analysis is made, using plane wave expansion, on how the width of the first three band gaps is influenced by filling ratio, dielectric constant ratio, and periodic width in one-dimensional photonic crystals （PhCs）. From simulation and analysis, there are one, two, and three peak points on the first, second and third band gaps respectively with the changes of filling ratio un- der fixed dielectric constant ratio. When filling ratio is fixed, the bandwidth of the first band gap consistently increases with dielectric constant ratio. However, no similar trend is observed in the second and the third band gaps. Because of scaling properties, varying periodic width does not alter the relative bandwidth.展开更多
Based on the vectorial structure of an electromagnetic wave, the analytical and concise expressions for the TE and TM terms of a veetorial plane wave diffracted by a circular aperture are derived in the far-field. The...Based on the vectorial structure of an electromagnetic wave, the analytical and concise expressions for the TE and TM terms of a veetorial plane wave diffracted by a circular aperture are derived in the far-field. The expressions of the energy flux distributions of the TE term, the TM term and the diffracted plane wave are also presented. The ratios of the power of the TE and TM terms to that of the diffracted plane wave are examined in the far-field. In addition, the far-field divergence angles of the TE term, the TM term and the diffracted plane wave, which are related to the energy flux distribution, are investigated. The different energy flux distributions of the TE and TM terms result in the discrepancy of their divergence angles. The influences of the linearly polarized angle and the radius of the circular aperture on the far-field divergence angles of the TE term, the TM term and the diffracted plane wave are discussed in detail. This research may promote the recognition of the optical propagation through a circular aperture.展开更多
Contrary to the superposition principle, it is well known that photorefraction exists in the vacuum with the presence of a strong static field, a laser field, or a rotational magnetic field. Different from the classic...Contrary to the superposition principle, it is well known that photorefraction exists in the vacuum with the presence of a strong static field, a laser field, or a rotational magnetic field. Different from the classical optical crystals, the refractive index also depends on the phase of the strong electromagnetic field. We obtain the phase and direction dependence of the refractive index of a probe wave incident in the strong field of a circular-polarized plane wave by solving the Maxwell equations corrected by the effective Lagrangian. It may provide a valuable theoretical basis to calculate the polarization evolution of waves in the strong electromagnetic circumstances of pulsar or neutron stars.展开更多
The model of equations of micropolar generalized magneto-thermoelasticity is introduced within the context of the theory of two temperatures generalized thermoelasticity and we consider a problem of an isotropic homog...The model of equations of micropolar generalized magneto-thermoelasticity is introduced within the context of the theory of two temperatures generalized thermoelasticity and we consider a problem of an isotropic homogeneous micropolar medium taking into account the heat effects and allowing the magnetic field effects. A plane wave analysis is employed to obtain the exact formulas of the two temperatures （conductive and mechanical）, displacement components, micro-rotation components, stresses, couple stresses, induced electric current, electric field and magnetic field. Arbitrary application is chosen to enable us to get the complete solution. The considered variables are presented graphically and discussions are made for the results.展开更多
The sea-surface height （SSH） signatures of internal tides extracted from the TOPEX/Poseidon （T/P） altimeter data along satellite tracks are fitted with superposition of several plane waves which have different wav...The sea-surface height （SSH） signatures of internal tides extracted from the TOPEX/Poseidon （T/P） altimeter data along satellite tracks are fitted with superposition of several plane waves which have different wavenumber vectors. The key problem of plane wave fitting with iterative method is how to determine the initial value of wavenumber of each plane wave. The previous solving method is to analyze the internal tidal SSH signatures along each track with wavenumber spectrum. But it is found that the problem cannot be solved completely with the wavenumber spectrum analysis method only. The method based on the combination of wavenumber spectrum analysis method and the exhaustive method is proposed to determine the initial values of wavenumbers for iteration. Numerical results indicate that the proposed method is not only reasonable and feasible but also better than the previous method. The proposed method is an improvement of the previous one, which is beneficial to improving the precision of plane wave fitting of the T/P internal tidal SSH signatures and deepening the understanding of the internal tides in ocean.展开更多
For the last few decades,the parameter estimation of electromagnetic plane waves i.e.,far field sources,impinging on antenna array geometries has attracted a lot of researchers due to their use in radar,sonar and unde...For the last few decades,the parameter estimation of electromagnetic plane waves i.e.,far field sources,impinging on antenna array geometries has attracted a lot of researchers due to their use in radar,sonar and under water acoustic environments.In this work,nature inspired heuristics based on the flower pollination algorithm(FPA)is designed for the estimation problem of amplitude and direction of arrival of far field sources impinging on uniform linear array(ULA).Using the approximation in mean squared error sense,a fitness function of the problem is developed and the strength of the FPA is utilized for optimization of the cost function representing scenarios for various number of sources non-coherent located in the far field.The worth of the proposed FPA based nature inspired computing heuristic is established through assessment studies on fitness,histograms,cumulative distribution function and box plots analysis.The other worthy perks of the proposed scheme include simplicity of concept,ease in the implementation,extendibility and wide range of applicability to solve complex optimization problems.These salient features make the proposed approach as an attractive alternative to be exploited for solving different parameter estimation problems arising in nonlinear systems,power signal modelling,image processing and fault diagnosis.展开更多
We construct the two-flux colliding plane wave solutions in higher-dimensional gravity theory with dilaton,and two complementary fluxes. Two kinds of solutions have been obtained: Bell-Szekeres (BS) type and homogeneo...We construct the two-flux colliding plane wave solutions in higher-dimensional gravity theory with dilaton,and two complementary fluxes. Two kinds of solutions have been obtained: Bell-Szekeres (BS) type and homogeneous type. After imposing the junction condition, we find that only the BS type solution is physically well-defined. Furthermore, we show that the future curvature singularity is always developed for our solutions.展开更多
The plane wave numerical technique is recast from Ampere’s and Faraday’s laws for materials that are characterized with a bianisotropic form of the constitutive relations. The populating expressions are provided for...The plane wave numerical technique is recast from Ampere’s and Faraday’s laws for materials that are characterized with a bianisotropic form of the constitutive relations. The populating expressions are provided for the eigenvalue matrix system that can be directly solved for the angular frequencies and field profiles when bianisotropy is included. To demonstrate the computation process and expected state diagrams and field profiles, numerical computation examples are provided for a bianisotropic Bragg Array with central defect. It is shown that the location of the magnetoelectric tensor elements has a significant effect on the eigenstates of an equivalent isotropic (anisotropic) structure. One form of the magnetoelectric tensor (diagonal elements only) leads to the observation of merging states and the formation of exceptional points. The numerical approach presented can be implemented as an add-on to the familiar plane wave numerical technique.展开更多
Plane wave generators(PWGs) are used to accelerate flyer plates to high velocities with their generated plane waves,which are widely used in the test of dynamic properties of materials.The traditional PWG is composed ...Plane wave generators(PWGs) are used to accelerate flyer plates to high velocities with their generated plane waves,which are widely used in the test of dynamic properties of materials.The traditional PWG is composed of two explosives with different detonation velocities.It is difficult to implement the related fabrication processes and control the generated waves due to its complicated structures.A simple plane wave generator is presented in this paper,which is composed of two identical cylindrical high explosive(HE) charges and an air-metal barrier.A theoretical model was established based on two different paths of the propagation of detonation waves,based on which the size of air-metal barrier was calculated for a given charge.The corresponding numerical simulations were also carried out by AUTODYN-2D based on the calculated results,which were used to compare with the theoretical calculations.A detonation wave with a flatness of 0.039 μs within the range of 70-percent diameter of the main charge was obtained through the simulations.展开更多
The governing equations of a transversely isotropic dissipative medium are solved analytically to obtain the speeds of plane waves. The appropriate solutions satisfy the required boundary conditions at the stress-free...The governing equations of a transversely isotropic dissipative medium are solved analytically to obtain the speeds of plane waves. The appropriate solutions satisfy the required boundary conditions at the stress-free surface to obtain the expressions of the reflection coefficients of reflected quasi-P (qP) and quasi-SV (qSV) waves in closed form for the incidence of qP and qSV waves. A particular model is chosen for numerical computation of these reflection coefficients for a certain range of the angle of incidence. The numerical values of these reflection coefficients are shown graphically against the angle of incidence for different values of initial stress parameter. The impact of initial stress parameter on the reflection coefficients is observed significantly.展开更多
基金supported in part by the National Natural Science Foundation of China (NSFC) under Grants 62271310 and 62125108in part by the Fundamental Research Funds for the Central Universities of Chinain part by the NSFC under Grant 62431014
文摘Extremely large-scale array(XL-array)communications can significantly improve the transmission rate,spectral efficiency,and spatial resolution,and has great potential in next-generation mobile communication networks.A crucial problem in XLarray communications is to determine the boundary of applicable regions of the plane wave model(PWM)and spherical wave model(SWM).In this paper,we propose new PWM/SWM demarcations for XL-arrays from the viewpoint of channel gain and rank.Four sets of results are derived for four different array setups.First,an equi-power line is derived for a point-touniform linear array(ULA)scenario,where an inflection point is found at±π6 central incident angles.Second,an equi-power surface is derived for a point-touniform planar array(UPA)scenario,and it is proved that cos2(ϕ)cos2(φ)=12 is a dividing curve,where ϕ andφdenote the elevation and azimuth angles,respectively.Third,an accurate and explicit expression of the equi-rank surface is obtained for a ULA-to-ULA scenario.Finally,an approximated expression of the equirank surface is obtained for a ULA-to-UPA scenario.With the obtained closed-form expressions,the equirank surface for any antenna structure and any angle can be well estimated.Furthermore,the effect of scatterers is also investigated,from which some insights are drawn.
基金supported by the General Program of the National Natural Science Foundation of China(Grant No.52074295)the Special Fund for Basic Scientific Research Business Expenses of Central Universities(Grant No.2022YJSSB06)supported by State Key Laboratory for Geomechanics and Deep Underground Engineering,China University of Mining and technology,Beijing,China(Grant No.SKLGDUEK202217).
文摘A nuclear explosion in the rock mass medium can produce strong shock waves,seismic shocks,and other destructive effects,which can cause extreme damage to the underground protection infrastructures.With the increase in nuclear explosion power,underground protection engineering enabled by explosion-proof impact theory and technology ushered in a new challenge.This paper proposes to simulate nuclear explosion tests with on-site chemical explosion tests in the form of multi-hole explosions.First,the mechanism of using multi-hole simultaneous blasting to simulate a nuclear explosion to generate approximate plane waves was analyzed.The plane pressure curve at the vault of the underground protective tunnel under the action of the multi-hole simultaneous blasting was then obtained using the impact test in the rock mass at the site.According to the peak pressure at the vault plane,it was divided into three regions:the stress superposition region,the superposition region after surface reflection,and the approximate plane stress wave zone.A numerical simulation approach was developed using PFC and FLAC to study the peak particle velocity in the surrounding rock of the underground protective cave under the action of multi-hole blasting.The time-history curves of pressure and peak pressure partition obtained by the on-site multi-hole simultaneous blasting test and numerical simulation were compared and analyzed,to verify the correctness and rationality of the formation of an approximate plane wave in the simulated nuclear explosion.This comparison and analysis also provided a theoretical foundation and some research ideas for the ensuing study on the impact of a nuclear explosion.
基金the National Natural Science Foundation of China(Grant Nos.50974081,50774051, 51034003)the Foundation for the Author of National Excellent Doctoral Dissertation of PR China(Grant No.200958)+1 种基金the Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT0843)the Open Research Fund Program of the State Key Laboratory of Coal Resources and Safe Mining(China University of Mining and Technology, Beijing)(No.2007-07)
文摘To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is：（1） Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.（2） Define the topology relationships and related concepts of triangular unit ray tracing.（3） The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.（4） By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point＇s traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.
文摘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.
基金This project is sporspored by Fund item:the National Development and Innovation Committee Program (2005) 2372the National High-tech R&D Program (863 Program) 2006AA06Z241 of ChinaYouth Innovation Fund of CNPC (Program:Prestack Imaging Integral Study for Complex near Surface)
文摘The Pre-Stack Depth Migration （PSDM） method based on wavefield continuation is the most reliable method for imaging complex structure in the subsurface, although there are large computational costs and poorly adaptive geometry. Plane wave shot migration is another method to perform exact wave equation prestack imaging with high computational efficiency and without the migration aperture problem. Moreover, wavefield energy can be compensated at the target zone by controlled illumination. In this paper, plane wave shot PSDM was implemented by the control of the plane down-going wavefield and selection of number and range of the raypaths in order to optimize the imaging effect. In addition, controlled illumination techniques are applied to enhance the imaging precision of interesting areas at different depths. Numerical calculation indicates that plane wave shot imaging is a rapid and efficient method with less computational cost and easy parallel computation compared to the single-square-root operator imaging for common shot gathers and double- square-root operator imaging for common midpoint gathers.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 10874086 and 10834009)the National Basic Research Program of China (Grant No. 2010CB327803)
文摘The super-cell plane wave expansion method is employed to calculate band structures for the design of a siliconbased one-dimensional phononic crystal plate with large absolute forbidden bands. In this method, a low impedance medium is introduced to replace the free stress boundary, which largely reduces the computational complexity. The dependence of band gaps on structural parameters is investigated in detail. To prove the validity of the super-cell plane wave expansion, the transmitted power spectra of the Lamb wave are calculated by using a finite element method. With the detailed computation, the band-gap of a one-dimensional plate can be designed as required with appropriate structural parameters, which provides a guide to the fabrication of a Lamb wave phononic crystal.
基金the National Natural Science Foundation of China(50478014)the National 973 Program(2007CB714200)the Beijing Natural Science Foundation(8061003).
文摘A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.
文摘The electronic density of states and band structures of doped and un-doped anatase TiO2 were studied by the Linearized Augmented Plane Wave method based on the density functional theory. The calculation shows that the band structures of TiO2 crystals doped with transition metal atoms become narrower. Interesting, an excursion towards high energy level with increasing atomic number in the same element period could be observed after doping with transition metal atoms.
基金The project sponsored by the Earthquake Science Foundation under Contract No. 90141
文摘The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV waves in the discrete model are first obtained by means of separating the characteristic equation of the motion equation, and then used to analyse the properties of P-and SV-homogeneous, inhomogeneous waves and other types of motion in the model. The dispersion characters, cut-off frequencies of P and SV waves, the polarization drift and appendent anisotropic property of wave motion caused by the discretization are finally discussed.
基金This work was supported by the National Natural Science Foundation of China(No.19772064)by the project of CAS KJ 951-1-20
文摘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.
基金Supported by the National Natural Science Foundation of China (61036006)
文摘A theoretical analysis is made, using plane wave expansion, on how the width of the first three band gaps is influenced by filling ratio, dielectric constant ratio, and periodic width in one-dimensional photonic crystals （PhCs）. From simulation and analysis, there are one, two, and three peak points on the first, second and third band gaps respectively with the changes of filling ratio un- der fixed dielectric constant ratio. When filling ratio is fixed, the bandwidth of the first band gap consistently increases with dielectric constant ratio. However, no similar trend is observed in the second and the third band gaps. Because of scaling properties, varying periodic width does not alter the relative bandwidth.
基金Project supported by the National Natural Science Foundation of China (Grant No.10974179)the Natural Science Foundation of Zhejiang Province of China (Grant No.Y1090073)
文摘Based on the vectorial structure of an electromagnetic wave, the analytical and concise expressions for the TE and TM terms of a veetorial plane wave diffracted by a circular aperture are derived in the far-field. The expressions of the energy flux distributions of the TE term, the TM term and the diffracted plane wave are also presented. The ratios of the power of the TE and TM terms to that of the diffracted plane wave are examined in the far-field. In addition, the far-field divergence angles of the TE term, the TM term and the diffracted plane wave, which are related to the energy flux distribution, are investigated. The different energy flux distributions of the TE and TM terms result in the discrepancy of their divergence angles. The influences of the linearly polarized angle and the radius of the circular aperture on the far-field divergence angles of the TE term, the TM term and the diffracted plane wave are discussed in detail. This research may promote the recognition of the optical propagation through a circular aperture.
基金supported by the National Basic Research Program of China(Grant No.2011CB808104)the National Natural Science Foundation of China(Grant No.11105233)
文摘Contrary to the superposition principle, it is well known that photorefraction exists in the vacuum with the presence of a strong static field, a laser field, or a rotational magnetic field. Different from the classical optical crystals, the refractive index also depends on the phase of the strong electromagnetic field. We obtain the phase and direction dependence of the refractive index of a probe wave incident in the strong field of a circular-polarized plane wave by solving the Maxwell equations corrected by the effective Lagrangian. It may provide a valuable theoretical basis to calculate the polarization evolution of waves in the strong electromagnetic circumstances of pulsar or neutron stars.
文摘The model of equations of micropolar generalized magneto-thermoelasticity is introduced within the context of the theory of two temperatures generalized thermoelasticity and we consider a problem of an isotropic homogeneous micropolar medium taking into account the heat effects and allowing the magnetic field effects. A plane wave analysis is employed to obtain the exact formulas of the two temperatures （conductive and mechanical）, displacement components, micro-rotation components, stresses, couple stresses, induced electric current, electric field and magnetic field. Arbitrary application is chosen to enable us to get the complete solution. The considered variables are presented graphically and discussions are made for the results.
基金The National Natural Science Foundation of China under contract No. 41076006the State Ministry of Science and Technology of China under contract No. 2008AA09A402the Ministry of Education’s "111" Project of China under contract No. B07036
文摘The sea-surface height （SSH） signatures of internal tides extracted from the TOPEX/Poseidon （T/P） altimeter data along satellite tracks are fitted with superposition of several plane waves which have different wavenumber vectors. The key problem of plane wave fitting with iterative method is how to determine the initial value of wavenumber of each plane wave. The previous solving method is to analyze the internal tidal SSH signatures along each track with wavenumber spectrum. But it is found that the problem cannot be solved completely with the wavenumber spectrum analysis method only. The method based on the combination of wavenumber spectrum analysis method and the exhaustive method is proposed to determine the initial values of wavenumbers for iteration. Numerical results indicate that the proposed method is not only reasonable and feasible but also better than the previous method. The proposed method is an improvement of the previous one, which is beneficial to improving the precision of plane wave fitting of the T/P internal tidal SSH signatures and deepening the understanding of the internal tides in ocean.
基金the Deanship of Scientific Research at Majmaah University for supporting this work under Project Number No.R-2021-27.
文摘For the last few decades,the parameter estimation of electromagnetic plane waves i.e.,far field sources,impinging on antenna array geometries has attracted a lot of researchers due to their use in radar,sonar and under water acoustic environments.In this work,nature inspired heuristics based on the flower pollination algorithm(FPA)is designed for the estimation problem of amplitude and direction of arrival of far field sources impinging on uniform linear array(ULA).Using the approximation in mean squared error sense,a fitness function of the problem is developed and the strength of the FPA is utilized for optimization of the cost function representing scenarios for various number of sources non-coherent located in the far field.The worth of the proposed FPA based nature inspired computing heuristic is established through assessment studies on fitness,histograms,cumulative distribution function and box plots analysis.The other worthy perks of the proposed scheme include simplicity of concept,ease in the implementation,extendibility and wide range of applicability to solve complex optimization problems.These salient features make the proposed approach as an attractive alternative to be exploited for solving different parameter estimation problems arising in nonlinear systems,power signal modelling,image processing and fault diagnosis.
文摘We construct the two-flux colliding plane wave solutions in higher-dimensional gravity theory with dilaton,and two complementary fluxes. Two kinds of solutions have been obtained: Bell-Szekeres (BS) type and homogeneous type. After imposing the junction condition, we find that only the BS type solution is physically well-defined. Furthermore, we show that the future curvature singularity is always developed for our solutions.
文摘The plane wave numerical technique is recast from Ampere’s and Faraday’s laws for materials that are characterized with a bianisotropic form of the constitutive relations. The populating expressions are provided for the eigenvalue matrix system that can be directly solved for the angular frequencies and field profiles when bianisotropy is included. To demonstrate the computation process and expected state diagrams and field profiles, numerical computation examples are provided for a bianisotropic Bragg Array with central defect. It is shown that the location of the magnetoelectric tensor elements has a significant effect on the eigenstates of an equivalent isotropic (anisotropic) structure. One form of the magnetoelectric tensor (diagonal elements only) leads to the observation of merging states and the formation of exceptional points. The numerical approach presented can be implemented as an add-on to the familiar plane wave numerical technique.
基金sponsored by the National Natural Science Foundation of China(NSFC10902053)the Zijin Intelligent Program,Nanjing University of Science and Technology (2013_ZJ_0101)+1 种基金Qing Lan Project of Jiangsu province and the Fundamental Research Funds for the Central Universities (No.2012XQTR006)the State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology)(KFJJ15-07M) for their great support on the research work presented in this paper
文摘Plane wave generators(PWGs) are used to accelerate flyer plates to high velocities with their generated plane waves,which are widely used in the test of dynamic properties of materials.The traditional PWG is composed of two explosives with different detonation velocities.It is difficult to implement the related fabrication processes and control the generated waves due to its complicated structures.A simple plane wave generator is presented in this paper,which is composed of two identical cylindrical high explosive(HE) charges and an air-metal barrier.A theoretical model was established based on two different paths of the propagation of detonation waves,based on which the size of air-metal barrier was calculated for a given charge.The corresponding numerical simulations were also carried out by AUTODYN-2D based on the calculated results,which were used to compare with the theoretical calculations.A detonation wave with a flatness of 0.039 μs within the range of 70-percent diameter of the main charge was obtained through the simulations.
文摘The governing equations of a transversely isotropic dissipative medium are solved analytically to obtain the speeds of plane waves. The appropriate solutions satisfy the required boundary conditions at the stress-free surface to obtain the expressions of the reflection coefficients of reflected quasi-P (qP) and quasi-SV (qSV) waves in closed form for the incidence of qP and qSV waves. A particular model is chosen for numerical computation of these reflection coefficients for a certain range of the angle of incidence. The numerical values of these reflection coefficients are shown graphically against the angle of incidence for different values of initial stress parameter. The impact of initial stress parameter on the reflection coefficients is observed significantly.
