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 method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order sy...The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.展开更多
The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the ...The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.展开更多
Nonlinear response of the driven Duffing oscillator to periodic or quasi-periodic signals has been well studied. In this paper, we investigate the nonlinear response of the driven Duffing oscillator to non-periodic, m...Nonlinear response of the driven Duffing oscillator to periodic or quasi-periodic signals has been well studied. In this paper, we investigate the nonlinear response of the driven Duffing oscillator to non-periodic, more specifically, chaotic time series. Through numerical simulations, we find that the driven Duffing oscillator can also show regular nonlinear response to the chaotic time series with different degree of chaos as generated by the same chaotic series generating model, and there exists a relationship between the state of the driven Duffing oscillator and the chaoticity of the input signal of the driven Duffing oscillator. One real-world and two artificial chaotic time series are used to verify the new feature of Duffing oscillator. A potential application of the new feature of Duffing oscillator is also indicated.展开更多
A fast Time Domain Integral Equation(TDIE) solver is presented for analysis of transient scattering from electrically large conducting complex objects.The numerical process of Marching-On-in-Time(MOT) method based TDI...A fast Time Domain Integral Equation(TDIE) solver is presented for analysis of transient scattering from electrically large conducting complex objects.The numerical process of Marching-On-in-Time(MOT) method based TDIE encounters high computational cost and exorbitant memory requirements.A group-style accelerated method-Plane Wave Time Domain(PWTD) algorithm,which permits rapid evaluation of transient wave field generated by temporally bandlimited sources,is employed to reduce the computational cost of MOT-based TDIE solvers.An efficient compressed storage technique for sparse matrix is adopted to decrease the enormous memory requirements of MOT.The scheme of the Multi-Level PWTD(MLPWTD)-enhanced MOT with compressed storage for sparse matrix is presented for analysis of transient scattering from electrically large complex objects in this paper.The numerical simulation results demonstrate the validity and efficiency of the presented scheme.展开更多
The time-history responses of the surface were obtained for a linear elastic half-plane including regularly distributed enormous embedded circular cavities subjected to propagating obliquely incident plane SH-waves. A...The time-history responses of the surface were obtained for a linear elastic half-plane including regularly distributed enormous embedded circular cavities subjected to propagating obliquely incident plane SH-waves. An advanced numerical approach named half-plane time-domain boundary element method（BEM）, which only located the meshes around the cavities, was used to create the model. By establishing the modified boundary integral equation（BIE）independently for each cavity and forming the matrices, the final coupled equation was solved step-by-step in the timedomain to obtain the boundary values. The responses were developed for a half-plane with 512 cavities. The amplification patterns were also obtained to illustrate the frequencydomain responses for some cases. According to the results,the presence of enormous cavities affects the scattering and diffraction of the waves arrived to the surface. The introduced method can be recommended for geotechnical/mechanical engineers to model structures in the fields of earthquake engineering and composite materials.展开更多
Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method （IBEM）, and the effec...Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method （IBEM）, and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA （peak ground acceleration） and its PRS （peak response spectrum） can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and E1 Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.展开更多
Generated by an ideal sinusoidal motion of the vertical plate, the simplest linear solution in time domain for two-dimensional regular waves is derived. The solution describes the propagation process of the plane prog...Generated by an ideal sinusoidal motion of the vertical plate, the simplest linear solution in time domain for two-dimensional regular waves is derived. The solution describes the propagation process of the plane progressive wave with a front, and will approach the linear steady- state solution as the oscillation time of the plate approaches infinity. The solution presented in this paper can be used to provide an incident wave model with analytical expression for solving the problems of diffraction and response of floating bodies in time domain.展开更多
The optimal plane for measurement of the right ventricular (RV) volumes by real-time three-dimensional echocardiography (RT3DE) was determined and the feasibility and accuracy of RT3DE in studying RV systolic function...The optimal plane for measurement of the right ventricular (RV) volumes by real-time three-dimensional echocardiography (RT3DE) was determined and the feasibility and accuracy of RT3DE in studying RV systolic function was assessed. RV “Full volume” images were acquired by RT3DE in 22 healthy subjects. RV end-diastolic volumes (RVEDV) and end-systolic volumes (RVESV) were outlined using apical biplane, 4-plane, 8-plane, 16-plane offline separately. RVSV and RVEF were calculated. Meanwhile tricuspid annual systolic excursion (TASE) was measured by M-mode echo. LVSV was outlined by 2-D echo according to the biplane Simpson's rule. The results showed: (1) There was a good correlation between RVSV measured from series planes and LVSV from 2-D echo (r=0.73; r=0.69; r=0.63; r=0.66, P<0.25—0.0025); (2) There were significant differences between RVEDV in biplane and those in 4-, 8-, 16-plane (P<0.001). There was also difference between RV volume in 4-plane and that in 8-plane (P<0.05), but there was no significant difference between RV volume in 8-plane and that in 16-plane (P>0.05); (3) Inter-observers and intro-observers variability analysis showed that there were close agreements and relations for RV volumes (r=0.986, P<0.001; r=0.93, P<0.001); (4) There was a significantly positive correlation of TASE to RVSV and RVEF from RT3DE (r=0.83; r=0.90). So RV volume measures with RT3DE are rapid, accurate and reproducible. In view of RV's complex shape, apical 8-plane method is better in clinical use. It may allow early detection of RV systolic function.展开更多
Room-temperature terahertz （THz） detectors indicate a great potential in the imaging application because of their real-time, compact bulk, and wide spectral band responding characteristics. THz detectors with differ...Room-temperature terahertz （THz） detectors indicate a great potential in the imaging application because of their real-time, compact bulk, and wide spectral band responding characteristics. THz detectors with different dimensions based on a micro-bridge structure have been designed and fabricated to get optimized micro-bolometer parameters from the test results of membrane deformation. A nanostructured titanium （Ti） thin film absorber is integrated in the micro-bridge structure of the VOx micro-bolometer by a combined process of magnetron sputtering and reactive ion etching （RIE）, and its improvement of THz absorption is verified by an optical characteristics mesurement. Continuous-wave THz detection and imaging are demonstrated by using a 2.52 THz far infrared CO2 laser and a 320x240 vanadium oxide micro-bolometer focal plane array with an optimized cell structure. With this detecting system, THz imaging of metal concealed in a wiping cloth and an envelope is demonstrated, respectively.展开更多
基金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.
基金supported by the National Natural Science Foundation of China(Grant Nos.60931002 and 61101064)the Universities Natural Science Foundation of Anhui Province,China(Grant Nos.KJ2011A002 and 1108085J01)
文摘The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11304074,61475042,and 11274088)the Natural Science Foundation of Hebei Province,China(Grant Nos.A2015202320 and GCC2014048)the Key Subject Construction Project of Hebei Province University,China
文摘The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.
基金supported by the National Natural Science Foundation of China (Grant Nos 40574051 and 40774054)
文摘Nonlinear response of the driven Duffing oscillator to periodic or quasi-periodic signals has been well studied. In this paper, we investigate the nonlinear response of the driven Duffing oscillator to non-periodic, more specifically, chaotic time series. Through numerical simulations, we find that the driven Duffing oscillator can also show regular nonlinear response to the chaotic time series with different degree of chaos as generated by the same chaotic series generating model, and there exists a relationship between the state of the driven Duffing oscillator and the chaoticity of the input signal of the driven Duffing oscillator. One real-world and two artificial chaotic time series are used to verify the new feature of Duffing oscillator. A potential application of the new feature of Duffing oscillator is also indicated.
文摘A fast Time Domain Integral Equation(TDIE) solver is presented for analysis of transient scattering from electrically large conducting complex objects.The numerical process of Marching-On-in-Time(MOT) method based TDIE encounters high computational cost and exorbitant memory requirements.A group-style accelerated method-Plane Wave Time Domain(PWTD) algorithm,which permits rapid evaluation of transient wave field generated by temporally bandlimited sources,is employed to reduce the computational cost of MOT-based TDIE solvers.An efficient compressed storage technique for sparse matrix is adopted to decrease the enormous memory requirements of MOT.The scheme of the Multi-Level PWTD(MLPWTD)-enhanced MOT with compressed storage for sparse matrix is presented for analysis of transient scattering from electrically large complex objects in this paper.The numerical simulation results demonstrate the validity and efficiency of the presented scheme.
文摘The time-history responses of the surface were obtained for a linear elastic half-plane including regularly distributed enormous embedded circular cavities subjected to propagating obliquely incident plane SH-waves. An advanced numerical approach named half-plane time-domain boundary element method（BEM）, which only located the meshes around the cavities, was used to create the model. By establishing the modified boundary integral equation（BIE）independently for each cavity and forming the matrices, the final coupled equation was solved step-by-step in the timedomain to obtain the boundary values. The responses were developed for a half-plane with 512 cavities. The amplification patterns were also obtained to illustrate the frequencydomain responses for some cases. According to the results,the presence of enormous cavities affects the scattering and diffraction of the waves arrived to the surface. The introduced method can be recommended for geotechnical/mechanical engineers to model structures in the fields of earthquake engineering and composite materials.
基金supported by National Natural Science Foundation of China under grant No. 50978183Tianjin Key Project for Applied Basic Research under grant No. 12JCZDJC29000
文摘Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method （IBEM）, and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA （peak ground acceleration） and its PRS （peak response spectrum） can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and E1 Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.
基金This study is financially supported by the National Natural Science Foundation of China
文摘Generated by an ideal sinusoidal motion of the vertical plate, the simplest linear solution in time domain for two-dimensional regular waves is derived. The solution describes the propagation process of the plane progressive wave with a front, and will approach the linear steady- state solution as the oscillation time of the plate approaches infinity. The solution presented in this paper can be used to provide an incident wave model with analytical expression for solving the problems of diffraction and response of floating bodies in time domain.
文摘The optimal plane for measurement of the right ventricular (RV) volumes by real-time three-dimensional echocardiography (RT3DE) was determined and the feasibility and accuracy of RT3DE in studying RV systolic function was assessed. RV “Full volume” images were acquired by RT3DE in 22 healthy subjects. RV end-diastolic volumes (RVEDV) and end-systolic volumes (RVESV) were outlined using apical biplane, 4-plane, 8-plane, 16-plane offline separately. RVSV and RVEF were calculated. Meanwhile tricuspid annual systolic excursion (TASE) was measured by M-mode echo. LVSV was outlined by 2-D echo according to the biplane Simpson's rule. The results showed: (1) There was a good correlation between RVSV measured from series planes and LVSV from 2-D echo (r=0.73; r=0.69; r=0.63; r=0.66, P<0.25—0.0025); (2) There were significant differences between RVEDV in biplane and those in 4-, 8-, 16-plane (P<0.001). There was also difference between RV volume in 4-plane and that in 8-plane (P<0.05), but there was no significant difference between RV volume in 8-plane and that in 16-plane (P>0.05); (3) Inter-observers and intro-observers variability analysis showed that there were close agreements and relations for RV volumes (r=0.986, P<0.001; r=0.93, P<0.001); (4) There was a significantly positive correlation of TASE to RVSV and RVEF from RT3DE (r=0.83; r=0.90). So RV volume measures with RT3DE are rapid, accurate and reproducible. In view of RV's complex shape, apical 8-plane method is better in clinical use. It may allow early detection of RV systolic function.
基金supported by the National Science Funds for Creative Research Groups of China under Grant No.61421002National High Technology Research and Development Program under Grant No.2015AA8123014
文摘Room-temperature terahertz （THz） detectors indicate a great potential in the imaging application because of their real-time, compact bulk, and wide spectral band responding characteristics. THz detectors with different dimensions based on a micro-bridge structure have been designed and fabricated to get optimized micro-bolometer parameters from the test results of membrane deformation. A nanostructured titanium （Ti） thin film absorber is integrated in the micro-bridge structure of the VOx micro-bolometer by a combined process of magnetron sputtering and reactive ion etching （RIE）, and its improvement of THz absorption is verified by an optical characteristics mesurement. Continuous-wave THz detection and imaging are demonstrated by using a 2.52 THz far infrared CO2 laser and a 320x240 vanadium oxide micro-bolometer focal plane array with an optimized cell structure. With this detecting system, THz imaging of metal concealed in a wiping cloth and an envelope is demonstrated, respectively.
