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.展开更多
This study focuses on determining the second-order irregular wave loads in the time domain without using the Inverse Fast Fourier Transform(IFFT).Considering the substantial displacement effects that Floating Offshore...This study focuses on determining the second-order irregular wave loads in the time domain without using the Inverse Fast Fourier Transform(IFFT).Considering the substantial displacement effects that Floating Offshore Wind Turbine(FOWT)support structures undergo when subjected to wave loads,the time-domain wave method is more suitable,while the frequency-domain method requiring IFFT cannot be used for moving bodies.Nonetheless,the computational challenges posed by the considerable computer time requirements of the time-domain wave method remain a significant obstacle.Thus,the paper incorporates various numerical schemes,including parallel computing and extrapolation of wave forces during specific time steps to improve overall efficiency.Despite the effectiveness of these schemes,the computational difficulties associated with the time-domain wave method persist.This study then proposes an innovative approach utilizing different randomnumbers in distinct segments,significantly reducing the computation of second-order wave loads.This random number interpolation ensures a smooth curve transition between two segments,emphasizingminimizing errors near the end of the first segment.Numerical analyses demonstrate substantial decreases in total computer time for FOWT structural analyses while maintaining consistent steel design results.The proposed method is uncomplicated,requiring only a simple subprogram modification in a conventional wave load computer program.展开更多
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.展开更多
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.展开更多
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.展开更多
Microcracks are common in compact bone,but their continued propagation can lead to macroscopic fractures.These microcracks cannot be visualized radiographically,necessitating alternative noninvasive methods to identif...Microcracks are common in compact bone,but their continued propagation can lead to macroscopic fractures.These microcracks cannot be visualized radiographically,necessitating alternative noninvasive methods to identify excessive microcracking and prevent fractures.In this study,terahertz time-domain spectroscopy(THz-TDS)was used to examine bone interiors near cracks resulting from loading in bovine tibia samples.Various loading configurations,such as impact,quasi-static loading,and fatigue loading,known to induce different types of micro-scale damage,were applied.The values of refractive index and absorption coefficient of the bone samples were then determined from the THz-TDS spectra acquired before loading and after fracture.The study revealed that different loading configurations led to varying terahertz optical coefficients associated with various types of bone fractures.Specifically,the refractive index notably increased under fatigue loading but remained relatively stable during quasi-static bending.The absorption coefficient of bone decreased only under fatigue loading.Furthermore,samples were subjected to axial and radial impacts without sustaining damage.Results indicated that in the undamaged state,the change in refractive index was smaller compared to after impact failure,while the change in absorption coefficient remained consistent after failure.Under radial impact loading,changes in refractive index and absorption coefficient were significantly more pronounced than under axial loading.Prior to loading,the measured value of refractive index was 2.72±0.11,and the absorption coefficient was 6.33±0.09 mm^(−1)at 0.5 THz.展开更多
In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible plat...In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible platform was performed.Relevant hydrodynamic parameters were obtained using the retardation function method of three-dimensional frequency-domain potential flow theory.The results of the hydrodynamic analysis were highly consistent with the test findings,verifying the accuracy of the multifloating hydrodynamic coupling analysis,and key hydrodynamic parameters were solved for different water depths and the coupling effect.According to the obtained results,the hydrodynamic influence was the largest in shallow waters when the coupling effect was considered.Furthermore,the coupled motion equation combined with viscous damping,fender system,and mooring system was established,and the hydrodynamics,floating body motion,and dynamic response of the fender system were analyzed.Motion analysis revealed good agreement among the surge,sway,and yaw motions of the two floating bodies.However,when the wave period reached 10 s,the motion of the two floating bodies showed severe shock,and a relative motion was also observed.Therefore,excessive constraints should be added between the two floating bodies during construction to ensure construction safety.The numerical analysis and model test results of the semisubmersible platform and HYSY 229 barge at a water depth of 42 m and sea conditions of 0°,45°,and 90° were in good agreement,and the error was less than 5%.The maximum movement of the HYSY 229 barge reached 2.61 m in the sway direction,whereas that of the semisubmersible platform was 2.11 m.During construction,excessive constraints should be added between the two floating bodies to limit their relative movement and ensure construction safety.展开更多
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.展开更多
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.展开更多
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.展开更多
The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is pr...The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is prone to body freedomflutter(BFF),which is a result of coupling of the rigid body short-periodmodewith 1st wing bendingmode.Accurate prediction of the BFF characteristics is helpful to reflect the attitude changes of the vehicle intuitively and design the active flutter suppression control law.Instead of using the rigid body mode,this work simulates the rigid bodymotion of the model by using the six-degree-of-freedom(6DOF)equation.A dynamicmesh generation strategy particularly suitable for BFF simulation of free flying aircraft is developed.An accurate Computational Fluid Dynamics/Computational Structural Dynamics/six-degree-of-freedom equation(CFD/CSD/6DOF)-based BFF prediction method is proposed.Firstly,the time-domain CFD/CSD method is used to calculate the static equilibrium state of the model.Based on this state,the CFD/CSD/6DOF equation is solved in time domain to evaluate the structural response of themodel.Then combinedwith the variable stiffnessmethod,the critical flutter point of the model is obtained.This method is applied to the BFF calculation of a flyingwing model.The calculation results of the BFF characteristics of the model agree well with those fromthe modalmethod andNastran software.Finally,the method is used to analyze the influence factors of BFF.The analysis results show that the flutter speed can be improved by either releasing plunge constraint or moving the center ofmass forward or increasing the pitch inertia.展开更多
We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular...We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.展开更多
To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of t...To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.展开更多
基金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.
基金funded by National Science and Technology Council,grant number NSTC 113-2223-E-006-014.
文摘This study focuses on determining the second-order irregular wave loads in the time domain without using the Inverse Fast Fourier Transform(IFFT).Considering the substantial displacement effects that Floating Offshore Wind Turbine(FOWT)support structures undergo when subjected to wave loads,the time-domain wave method is more suitable,while the frequency-domain method requiring IFFT cannot be used for moving bodies.Nonetheless,the computational challenges posed by the considerable computer time requirements of the time-domain wave method remain a significant obstacle.Thus,the paper incorporates various numerical schemes,including parallel computing and extrapolation of wave forces during specific time steps to improve overall efficiency.Despite the effectiveness of these schemes,the computational difficulties associated with the time-domain wave method persist.This study then proposes an innovative approach utilizing different randomnumbers in distinct segments,significantly reducing the computation of second-order wave loads.This random number interpolation ensures a smooth curve transition between two segments,emphasizingminimizing errors near the end of the first segment.Numerical analyses demonstrate substantial decreases in total computer time for FOWT structural analyses while maintaining consistent steel design results.The proposed method is uncomplicated,requiring only a simple subprogram modification in a conventional wave load computer program.
基金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.
基金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.
基金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 Nos.11972247 and 12372080).
文摘Microcracks are common in compact bone,but their continued propagation can lead to macroscopic fractures.These microcracks cannot be visualized radiographically,necessitating alternative noninvasive methods to identify excessive microcracking and prevent fractures.In this study,terahertz time-domain spectroscopy(THz-TDS)was used to examine bone interiors near cracks resulting from loading in bovine tibia samples.Various loading configurations,such as impact,quasi-static loading,and fatigue loading,known to induce different types of micro-scale damage,were applied.The values of refractive index and absorption coefficient of the bone samples were then determined from the THz-TDS spectra acquired before loading and after fracture.The study revealed that different loading configurations led to varying terahertz optical coefficients associated with various types of bone fractures.Specifically,the refractive index notably increased under fatigue loading but remained relatively stable during quasi-static bending.The absorption coefficient of bone decreased only under fatigue loading.Furthermore,samples were subjected to axial and radial impacts without sustaining damage.Results indicated that in the undamaged state,the change in refractive index was smaller compared to after impact failure,while the change in absorption coefficient remained consistent after failure.Under radial impact loading,changes in refractive index and absorption coefficient were significantly more pronounced than under axial loading.Prior to loading,the measured value of refractive index was 2.72±0.11,and the absorption coefficient was 6.33±0.09 mm^(−1)at 0.5 THz.
基金the National Natural Science Foundation of China(No.U20A20328).
文摘In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible platform was performed.Relevant hydrodynamic parameters were obtained using the retardation function method of three-dimensional frequency-domain potential flow theory.The results of the hydrodynamic analysis were highly consistent with the test findings,verifying the accuracy of the multifloating hydrodynamic coupling analysis,and key hydrodynamic parameters were solved for different water depths and the coupling effect.According to the obtained results,the hydrodynamic influence was the largest in shallow waters when the coupling effect was considered.Furthermore,the coupled motion equation combined with viscous damping,fender system,and mooring system was established,and the hydrodynamics,floating body motion,and dynamic response of the fender system were analyzed.Motion analysis revealed good agreement among the surge,sway,and yaw motions of the two floating bodies.However,when the wave period reached 10 s,the motion of the two floating bodies showed severe shock,and a relative motion was also observed.Therefore,excessive constraints should be added between the two floating bodies during construction to ensure construction safety.The numerical analysis and model test results of the semisubmersible platform and HYSY 229 barge at a water depth of 42 m and sea conditions of 0°,45°,and 90° were in good agreement,and the error was less than 5%.The maximum movement of the HYSY 229 barge reached 2.61 m in the sway direction,whereas that of the semisubmersible platform was 2.11 m.During construction,excessive constraints should be added between the two floating bodies to limit their relative movement and ensure construction safety.
基金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.
基金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 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.
基金This work was supported by the National Natural Science Foundation of China(No.11872212)and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions.
文摘The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is prone to body freedomflutter(BFF),which is a result of coupling of the rigid body short-periodmodewith 1st wing bendingmode.Accurate prediction of the BFF characteristics is helpful to reflect the attitude changes of the vehicle intuitively and design the active flutter suppression control law.Instead of using the rigid body mode,this work simulates the rigid bodymotion of the model by using the six-degree-of-freedom(6DOF)equation.A dynamicmesh generation strategy particularly suitable for BFF simulation of free flying aircraft is developed.An accurate Computational Fluid Dynamics/Computational Structural Dynamics/six-degree-of-freedom equation(CFD/CSD/6DOF)-based BFF prediction method is proposed.Firstly,the time-domain CFD/CSD method is used to calculate the static equilibrium state of the model.Based on this state,the CFD/CSD/6DOF equation is solved in time domain to evaluate the structural response of themodel.Then combinedwith the variable stiffnessmethod,the critical flutter point of the model is obtained.This method is applied to the BFF calculation of a flyingwing model.The calculation results of the BFF characteristics of the model agree well with those fromthe modalmethod andNastran software.Finally,the method is used to analyze the influence factors of BFF.The analysis results show that the flutter speed can be improved by either releasing plunge constraint or moving the center ofmass forward or increasing the pitch inertia.
基金supported by the Young Scientists Fund of the National Natural Science Foundation of China(No.41304082)the China Postdoctoral Science Foundation(No.2016M590731)+2 种基金the Young Scientists Fund of the Natural Science Foundation of Hebei Province(No.D2014403011)the Program for Young Excellent Talents of Higher Education Institutions of Hebei Province(No.BJ2016046)the Geological survey project of China Geological Survey(No.1212011121197)
文摘We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.
基金supported by the National High Technology Research and Development Program (863 Program)(2009AA06Z108)
文摘To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.
