This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the...This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.展开更多
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.展开更多
Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open...Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain（ECT-CFDTD） method is used to simulate the wave propagation inside a perfect electric conductor（PEC） waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.展开更多
An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D tra...An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.展开更多
Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional fin...Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.展开更多
An Improved Locally Conformal Finite-Difference Time-Domain (ILC-FDTD) method is presented in this paper, which is used to analyze the edge inclined slots penetrating adjacent broadwalls of a finite wall thickness wav...An Improved Locally Conformal Finite-Difference Time-Domain (ILC-FDTD) method is presented in this paper, which is used to analyze the edge inclined slots penetrating adjacent broadwalls of a finite wall thickness waveguide. ILC-FDTD not only removes tile instability of the original locally conformal FDTD algorithm, but also improves the computational accuracy by locally modifying magnetic field update equations and the virtual iterative electric fields accordlng to the complexity of tile slot fringe fields. The mutual coupling between two edge inclined slots can also be analyzed by ILC-FDTD effectively.展开更多
Metallic ring-shaped nanotube arrays are proposed and its optical transmission properties are studied by using finite-difference time-domain （FDTD） method. Compared with the transmission spectra of conventional circ...Metallic ring-shaped nanotube arrays are proposed and its optical transmission properties are studied by using finite-difference time-domain （FDTD） method. Compared with the transmission spectra of conventional circular nanotube arrays, two photonic band gaps are emerged in the transmission spectra offing-shaped nanotube arrays, the two band gaps and transmission spectra are adjusted by the length, inner radius, intertube spacing and the dielectric constants of the core and embedding medium, and magnitude modification, redshift and blueshift of the resonance modes are observed. A metallic ring-shaped nanotube arrays for subwavelength band-stop filter in the range of visible light can be achieved. To understand its physical origin, field-interference mechanism was suggested by the field distributions. The proposed nanostructures and results may have great potential applications in subwavelength near-field optics.展开更多
The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy...The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy absorption and its distribution in homogeneousand layered homogenous lossy dielectric spheres are presented,and a comparison of these resultswith analytical solution is given.The calculation is carried out for dielectric cylinder on conduct-ing ground as well,and the results are compared with the image theory.All the computationsshew that the finite-difference time-domain method can give satisfactory results.展开更多
To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE met...To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.展开更多
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.展开更多
A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedur...A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.展开更多
Guided waves are generally considered as a powerful approach for crack detection in structures,which are commonly investigated using the finite element method(FEM).However,the traditional FEM has many disadvantages in...Guided waves are generally considered as a powerful approach for crack detection in structures,which are commonly investigated using the finite element method(FEM).However,the traditional FEM has many disadvantages in solving wave propagation due to the strict requirement of mesh density.To tackle this issue,this paper proposes an efficient time-domain spectral finite element method(SFEM)to analyze wave propagation in cracked structures,in which the breathing crack is modeled by definiiig the spectral gap element.Moreover,novel orthogonal polynomials and Gauss-Lobatto-Legendre quadrature rules are adopted to construct the spectral element.Meanwhile,a separable hard contact is utilized to simulate the breathing behavior.Finally,a comparison of the numerical results between the FEM and the SFEM is conducted to demonstrate the high efficiency and accuracy of the proposed method.Based on the developed SFEM,the nonlinear features of waves and influence of the incident mode are also studied in detail,which provides a helpful guide for a physical understanding of the wave propagation behavior in structures with breathing cracks.展开更多
A method for extracting optical parameters of plastics materials based on terahertz time domain spectroscopy is presented. The transmission-type Terahertz Time-Domain Spectroscopy(THz TDS) system is adopted to detect ...A method for extracting optical parameters of plastics materials based on terahertz time domain spectroscopy is presented. The transmission-type Terahertz Time-Domain Spectroscopy(THz TDS) system is adopted to detect the refractive index and extinction coefficient on different plastic materials. Then the corresponding spectral information is obtained by Fourier transform of the terahertz time domain waveform of the sampling points, including the corresponding amplitude and phase information of the waveform. The optical parameter extraction model is built. By using the simplex optimization method, the curves of the refractive index and extinction coefficient for the plastic material are obtained. The experimental samples are made of different plastic parallel plate materials. The experimental results show that the optimization of optical parameters can improve their extraction accuracy, and the error of refractive index is ±0.005. Extraction technology with the simplex optimization method of optical parameter based on THz TDS can help to extract the optical parameters of engineering plastics. It is of great significance for the research of terahertz nondestructive testing.展开更多
In this work, we treat scattering objects, water, surface and bottom in a truly unified manner in a parallel finitedifference time-domain （FDTD） scheme, which is suitable for distributed parallel computing in a mess...In this work, we treat scattering objects, water, surface and bottom in a truly unified manner in a parallel finitedifference time-domain （FDTD） scheme, which is suitable for distributed parallel computing in a message passing interface （MPI） programming environment. The algorithm is implemented on a cluster-based high performance computer system. Parallel computation is performed with different division methods in 2D and 3D situations. Based on analysis of main factors affecting the speedup rate and parallel efficiency, data communication is reduced by selecting a suitable scheme of task division. A desirable scheme is recommended, giving a higher speedup rate and better efficiency. The results indicate that the unified parallel FDTD algorithm provides a solution to the numerical computation of acoustic scattering.展开更多
To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed...To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.展开更多
Stochastic heat conduction and thermal stress analysis of structures has received considerable attention in recent years.The propagation of uncertain thermal environments will lead to stochastic variations in temperat...Stochastic heat conduction and thermal stress analysis of structures has received considerable attention in recent years.The propagation of uncertain thermal environments will lead to stochastic variations in temperature fields and thermal stresses.Therefore,it is reasonable to consider the variability of thermal environments while conducting thermal analysis.However,for ambient thermal excitations,only stationary random processes have been investigated thus far.In this study,the highly efficient explicit time-domain method(ETDM)is proposed for the analysis of non-stationary stochastic transient heat conduction and thermal stress problems.The explicit time-domain expressions of thermal responses are first constructed for a thermoelastic body.Then the statistical moments of thermal displacements and stresses can be directly obtained based on the explicit expressions of thermal responses.A numerical example involving non-stationary stochastic internal heat generation rate is investigated.The accuracy and efficiency of the proposed method are validated by comparison with the Monte-Carlo simulation.展开更多
In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibi...In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.展开更多
The hydroelastic behavior of a moored oil storage vessel subjected to arbitrary time-dependent external loads,which include wind,waves,and currents with different incident directions,is investigated with the time-doma...The hydroelastic behavior of a moored oil storage vessel subjected to arbitrary time-dependent external loads,which include wind,waves,and currents with different incident directions,is investigated with the time-domain modal expansion method.First,the water boundary integral equations on the body surface of a quarter model,which can be obtained via the free-surface Green’s function method,are established.Then,the time-dependent elastic deflection of the moored oil storage vessel is expressed by a superposition of modal functions and corresponding modal amplitudes,and a Galerkin scheme is applied to derive the linear system of equations for the modal amplitudes.The second-order linear differential equations for modal amplitudes are solved via the fourth-order Runge−Kutta method.The present model is validated against existing frequency domain results for a truncated cylinder and a VLFS.Numerical calculations for the moored oil storage vessel are then conducted to obtain the time series of various modal amplitudes and elastic displacements of the measurement points and the corresponding spectra with different incident directions.展开更多
A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spheric...A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.展开更多
Based on chromatographic theory, the moment method and the time-domain fitting analysis were applied to measure and evaluate the adsorption equilibrium constant and mass transfer properties (axial dispersion coefficie...Based on chromatographic theory, the moment method and the time-domain fitting analysis were applied to measure and evaluate the adsorption equilibrium constant and mass transfer properties (axial dispersion coefficient and effective intra-particle diffusivity) for toluene and p-dichlorobenzene on silica gel adsorbent in the subcritical and supercritical CO2. An apparatus based on supercritical fluid chromatography was established and the experiments were performed at temperatures of 298.15-318.15 K and pressures of 7.5-17.8 MPa. The two methods have been compared. The results show that for the systems studied here the moment method can give reasonable values for both adsorption equilibrium constant and mass transfer properties, but the time-domain analysis only can obtain the adsorption equilibrium constant. The dependence of adsorption equilibrium constant and mass transfer properties on temperature and pressure was investigated.展开更多
基金Project supported by Tianjin Research Program Application Foundation and Advanced Technology,China(Grant No.15JCQNJC01100)
文摘This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.
基金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 No.61231003)
文摘Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain（ECT-CFDTD） method is used to simulate the wave propagation inside a perfect electric conductor（PEC） waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.
基金supported by National Key Research and Development Program of China(No.2022YFF0800602)Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(No.2022B1212010002)Shenzhen Science and Technology Program(No.KQTD20170810111725321).
文摘Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.
基金Supported by the National Natural Science Fbundation of China(No.69931030)
文摘An Improved Locally Conformal Finite-Difference Time-Domain (ILC-FDTD) method is presented in this paper, which is used to analyze the edge inclined slots penetrating adjacent broadwalls of a finite wall thickness waveguide. ILC-FDTD not only removes tile instability of the original locally conformal FDTD algorithm, but also improves the computational accuracy by locally modifying magnetic field update equations and the virtual iterative electric fields accordlng to the complexity of tile slot fringe fields. The mutual coupling between two edge inclined slots can also be analyzed by ILC-FDTD effectively.
基金Projects(11164007,61275174)supported by the National Natural Science Foundation of ChinaProject(20100162110068)supported by the Doctoral Program of Higher Education of China+1 种基金Project(20132BAB212007)supported by the Natural Science Foundation of Jiangxi Province,ChinaProject(GJJ11107)supported by Scientific Foundation of Jiangxi Education Department,China
文摘Metallic ring-shaped nanotube arrays are proposed and its optical transmission properties are studied by using finite-difference time-domain （FDTD） method. Compared with the transmission spectra of conventional circular nanotube arrays, two photonic band gaps are emerged in the transmission spectra offing-shaped nanotube arrays, the two band gaps and transmission spectra are adjusted by the length, inner radius, intertube spacing and the dielectric constants of the core and embedding medium, and magnitude modification, redshift and blueshift of the resonance modes are observed. A metallic ring-shaped nanotube arrays for subwavelength band-stop filter in the range of visible light can be achieved. To understand its physical origin, field-interference mechanism was suggested by the field distributions. The proposed nanostructures and results may have great potential applications in subwavelength near-field optics.
文摘The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy absorption and its distribution in homogeneousand layered homogenous lossy dielectric spheres are presented,and a comparison of these resultswith analytical solution is given.The calculation is carried out for dielectric cylinder on conduct-ing ground as well,and the results are compared with the image theory.All the computationsshew that the finite-difference time-domain method can give satisfactory results.
基金Project supported by the National Natural Science Foundation of China(Grant No.51025622)
文摘To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.
基金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 project was financially supported by the National Natural Science Foundation of China under the Grant No. 19732004 the National Science Fund for Distinguished Young Scholars under the Grant No. 50029002
文摘A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.
基金the National Natural Sclenee Foundation of China(Grant No.51704222)China Pastdoctoral Science Foundation(Grant No.2018M633570)Fundamental Research Funds for the Cemtal Unveritiee(Grant No.3102017090004).
文摘Guided waves are generally considered as a powerful approach for crack detection in structures,which are commonly investigated using the finite element method(FEM).However,the traditional FEM has many disadvantages in solving wave propagation due to the strict requirement of mesh density.To tackle this issue,this paper proposes an efficient time-domain spectral finite element method(SFEM)to analyze wave propagation in cracked structures,in which the breathing crack is modeled by definiiig the spectral gap element.Moreover,novel orthogonal polynomials and Gauss-Lobatto-Legendre quadrature rules are adopted to construct the spectral element.Meanwhile,a separable hard contact is utilized to simulate the breathing behavior.Finally,a comparison of the numerical results between the FEM and the SFEM is conducted to demonstrate the high efficiency and accuracy of the proposed method.Based on the developed SFEM,the nonlinear features of waves and influence of the incident mode are also studied in detail,which provides a helpful guide for a physical understanding of the wave propagation behavior in structures with breathing cracks.
基金National defense technical basic research project,Terahertz detection technology and application research on ceramic matrix composites(JSZL2015411C002)
文摘A method for extracting optical parameters of plastics materials based on terahertz time domain spectroscopy is presented. The transmission-type Terahertz Time-Domain Spectroscopy(THz TDS) system is adopted to detect the refractive index and extinction coefficient on different plastic materials. Then the corresponding spectral information is obtained by Fourier transform of the terahertz time domain waveform of the sampling points, including the corresponding amplitude and phase information of the waveform. The optical parameter extraction model is built. By using the simplex optimization method, the curves of the refractive index and extinction coefficient for the plastic material are obtained. The experimental samples are made of different plastic parallel plate materials. The experimental results show that the optimization of optical parameters can improve their extraction accuracy, and the error of refractive index is ±0.005. Extraction technology with the simplex optimization method of optical parameter based on THz TDS can help to extract the optical parameters of engineering plastics. It is of great significance for the research of terahertz nondestructive testing.
基金Project supported by the National Defense Laboratory Foundation (Grant No.51444020103QT0601)the Shanghai Leading Academic Discipline Project (Grant No.T0102)
文摘In this work, we treat scattering objects, water, surface and bottom in a truly unified manner in a parallel finitedifference time-domain （FDTD） scheme, which is suitable for distributed parallel computing in a message passing interface （MPI） programming environment. The algorithm is implemented on a cluster-based high performance computer system. Parallel computation is performed with different division methods in 2D and 3D situations. Based on analysis of main factors affecting the speedup rate and parallel efficiency, data communication is reduced by selecting a suitable scheme of task division. A desirable scheme is recommended, giving a higher speedup rate and better efficiency. The results indicate that the unified parallel FDTD algorithm provides a solution to the numerical computation of acoustic scattering.
基金Supported by the National Natural Science Foundation of China under (Grant No.107 72040,50709005 and 50921001)the Major National Science and Technology Projects of China under (Grant No.2008ZX05026-02)the Open Fund of State Key Laboratory of Ocean Engineering
文摘To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.
基金funded by the National Natural Science Foundation of China (51678252)the Guangzhou Science and Technology Project (201804020069)
文摘Stochastic heat conduction and thermal stress analysis of structures has received considerable attention in recent years.The propagation of uncertain thermal environments will lead to stochastic variations in temperature fields and thermal stresses.Therefore,it is reasonable to consider the variability of thermal environments while conducting thermal analysis.However,for ambient thermal excitations,only stationary random processes have been investigated thus far.In this study,the highly efficient explicit time-domain method(ETDM)is proposed for the analysis of non-stationary stochastic transient heat conduction and thermal stress problems.The explicit time-domain expressions of thermal responses are first constructed for a thermoelastic body.Then the statistical moments of thermal displacements and stresses can be directly obtained based on the explicit expressions of thermal responses.A numerical example involving non-stationary stochastic internal heat generation rate is investigated.The accuracy and efficiency of the proposed method are validated by comparison with the Monte-Carlo simulation.
基金the National Natural Science Foundation of China(No.40774056)Program of Excellent Team in Harbin Institute of Technology
文摘In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.
基金financially supported by the Department of Natural Resources of Guangdong Province(Grant No.[2024]31)the National Natural Science Foundation of China(Grant No.52071145)+1 种基金the Natural Science Foundation of Guangdong Province,China(Grant No.2022B1515020071)the Fundamental Research Funds for the Central Universities(Grant No.2023ZYGXZR029).
文摘The hydroelastic behavior of a moored oil storage vessel subjected to arbitrary time-dependent external loads,which include wind,waves,and currents with different incident directions,is investigated with the time-domain modal expansion method.First,the water boundary integral equations on the body surface of a quarter model,which can be obtained via the free-surface Green’s function method,are established.Then,the time-dependent elastic deflection of the moored oil storage vessel is expressed by a superposition of modal functions and corresponding modal amplitudes,and a Galerkin scheme is applied to derive the linear system of equations for the modal amplitudes.The second-order linear differential equations for modal amplitudes are solved via the fourth-order Runge−Kutta method.The present model is validated against existing frequency domain results for a truncated cylinder and a VLFS.Numerical calculations for the moored oil storage vessel are then conducted to obtain the time series of various modal amplitudes and elastic displacements of the measurement points and the corresponding spectra with different incident directions.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61301063 and 41305017)
文摘A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.
基金Supported by research grant 02KJB530002 from Jiangsu Provincial Committee of Education.
文摘Based on chromatographic theory, the moment method and the time-domain fitting analysis were applied to measure and evaluate the adsorption equilibrium constant and mass transfer properties (axial dispersion coefficient and effective intra-particle diffusivity) for toluene and p-dichlorobenzene on silica gel adsorbent in the subcritical and supercritical CO2. An apparatus based on supercritical fluid chromatography was established and the experiments were performed at temperatures of 298.15-318.15 K and pressures of 7.5-17.8 MPa. The two methods have been compared. The results show that for the systems studied here the moment method can give reasonable values for both adsorption equilibrium constant and mass transfer properties, but the time-domain analysis only can obtain the adsorption equilibrium constant. The dependence of adsorption equilibrium constant and mass transfer properties on temperature and pressure was investigated.
