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.展开更多
Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modu...Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modulus and combining the dynamic skeleton curve and the damping degradation coefficient,the constitutive equation of the logarithmic dynamic skeleton can be obtained,which considers the damping effect in a soil dynamics problem.Based on the finite difference method and the multi-transmitting boundary condition,a 1D site seismic response analysis program called Soilresp1D has been developed herein and used to analyze the time-domain seismic response in three types of sites.At the same time,this study also provides numerical simulation results based on the hyperbolic constitutive model and the equivalent linear method.The results verify the rationality of the new soil dynamic constitutive model.It can analyze the mucky soil site nonlinear seismic response,reflecting the deformation characteristics and damping effect of the silty soil.The hysteresis loop area is more extensive,and the residual strain is evident.展开更多
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.展开更多
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.展开更多
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.展开更多
Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and ...Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.展开更多
An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge...An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the divole source.展开更多
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.展开更多
Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the...Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the broken spatial translational symmetry.However,the intensity of first-order Raman scattering at crystal edges has been rarely explored,although the mechanical stress and edge characteristics have been thoroughly studied by the Raman peak shift and the spectral features of the edge-related Raman modes.Here,by taking Ga As crystal with a well-defined edge as an example,we reveal the intensity enhancement of Raman-active modes and the emergence of Raman-forbidden modes under specific polarization configurations at the edge.This is attributed to the presence of a hot spot at the edge due to the redistributed electromagnetic fields and electromagnetic wave propagations of incident laser and Raman signal near the edge,which are confirmed by the finite-difference time-domain simulations.Spatially-resolved Raman intensities of both Raman-active and Raman-forbidden modes near the edge are calculated based on the redistributed electromagnetic fields,which quantitatively reproduce the corresponding experimental results.These findings offer new insights into the intensity enhancement of Raman scattering at crystal edges and present a new avenue to manipulate light–matter interactions of crystal by manufacturing various types of edges and to characterize the edge structures in photonic and optoelectronic devices.展开更多
During normal de-ballasting operations for floating docks,each ballast pump independently manages a specific group of ballast tanks.However,when a pump malfunctions,a connection valve between the two groups of ballast...During normal de-ballasting operations for floating docks,each ballast pump independently manages a specific group of ballast tanks.However,when a pump malfunctions,a connection valve between the two groups of ballast water systems is opened.This allows the adjacent pump to serve as a helper pump,simultaneously controlling two groups of ballast water systems.This study explores a full-scale floating dock’s dynamic behaviours during the de-ballasting operations under this situation through a numerical model.In the developed numerical model,the dock is described as a six-degree-of-freedom rigid body which is subjected to hydrostatic,hydrodynamic,and mooring loads.A hydraulic model of the piping network of the malfunctioning pump and the helper pump is proposed.A modified P-controller regulates opening angles of all tank valves for minimal pitch and roll.Two configurations of the floating dock,i.e.,a single floating dock and a floating dock with an onboard vessel,are considered.The numerical results show that the optimal helper pumps can be identified regarding the pumps’total de-ballasting capacity and the dock’s stability.The most severe scenarios can be determined in term of the dock’s maximum draught differences caused by its roll and pitch.The observed maximum draught differences remain small relative to the dock’s width,indicating the effectiveness of employing helper pumps and the proposed automatic ballast control strategy for one-pump malfunction scenarios.展开更多
Ultrasonic Lamb waves undergo complex mode conversion and diffraction at non-penetrating defects, such as plate corrosion and cracks. Lamb wave imaging has a resolution limit due to the guided wave dispersion characte...Ultrasonic Lamb waves undergo complex mode conversion and diffraction at non-penetrating defects, such as plate corrosion and cracks. Lamb wave imaging has a resolution limit due to the guided wave dispersion characteristics and Rayleigh criterion limitations. In this paper, a full convolutional network is designed to segment and reconstruct the received signals, enabling the automatic identification of target modalities. This approach eliminates clutter and mode conversion interference when calculating direct and accompanying acoustic fields in time-domain topological energy(TDTE) imaging.Subsequently, the measured accompanying acoustic field is reversed for adaptive focusing on defects and enhance the imaging quality. To circumvent the limitations of the Rayleigh criterion, the direct acoustic field and the accompanying acoustic field were fused to characterize the pixel distribution in the imaging region, achieving Lamb wave super-resolution imaging. Experimental results indicate that compared to the sign coherence factor-total focusing method(SCF-TFM),the proposed method achieves a 31.41% improvement in lateral resolution and a 29.53% increase in signal-to-noise ratio for single-blind-hole defects. In the case of multiple-blind-hole defects with spacings greater than the Rayleigh criterion resolution limit, it exhibits a 27.23% enhancement in signal-to-noise ratio. On the contrary, when the defect spacings are relatively smaller than the limit, this method has a higher resolution limit than SCF-TFM in super-resolution imaging.展开更多
Investigating thermal transport mechanisms at the interface between phase change materials(PCMs)and high thermally conductive fillers has become increasingly significant in developing phase change energy storage techn...Investigating thermal transport mechanisms at the interface between phase change materials(PCMs)and high thermally conductive fillers has become increasingly significant in developing phase change energy storage technologies.This study explores the interfacial thermal transport between a representative PCM,erythritol,and various fillers,including crystalline(Si C,Si_(3)N_(4))and amorphous(Si O_(2))nanoparticles,using molecular dynamics(MD)simulations.Additionally,time-domain thermoreflectance(TDTR)experiments were performed to quantify the interfacial thermal conductance between erythritol and the three types of fillers,yielding values of 50.1,40.0,and25.6 MW m^(–2)K^(-1).These results align well with the trends observed in the simulations.Furthermore,the underlying mechanisms of interfacial heat transfer were analyzed by examining the phonon density of states,overlap energy,and interaction energy.This research provides innovative insights into nanoscale interfacial thermal transport in composite PCMs.This could lead to significant advancements in thermal management technologies,particularly in developing more efficient thermal energy storage systems.展开更多
Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotrop...Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic （TDAEM） method is still confined to one dimensional （1D） cases, and the corresponding three-dimensional （3D） numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell＇s equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.展开更多
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.展开更多
Copper sulfate pentahydrate is investigated by terahertz time-domain spectroscopy. It is shown that the terahertz absorption coefficients are correlated with the particle size of the samples, as well as the heating ra...Copper sulfate pentahydrate is investigated by terahertz time-domain spectroscopy. It is shown that the terahertz absorption coefficients are correlated with the particle size of the samples, as well as the heating rates of the ambient temperature. Furthermore, the water molecules of copper sulfate pentahydrate can be quantitatively characterized due to the high sensitivity of the terahertz wave to water molecules. Based on such results, the status of water incorporated in mineral opal is also characterized using terahertz time-domain spectroscopy. It indicates that terahertz technology can be considered as an efficient method to detect the dehydration of minerals.展开更多
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.展开更多
The path loss analysis model based on 5 rays in vacancy indoor environment is proposed. The relationship between multipath overlapping and the path loss is analyzed mathematically. Time-domain technique is introduced ...The path loss analysis model based on 5 rays in vacancy indoor environment is proposed. The relationship between multipath overlapping and the path loss is analyzed mathematically. Time-domain technique is introduced to compute reflection coefficient in a very short time interval. A 5 rays path loss calculation method, which is satisfactory accurate, is developed. 5 typical environments are involved to analyze and generalize the common path loss characteristics in vacancy indoor environment. The simulation result shows that the path loss can be characterized as 3 zones with different path loss exponent as distance between transmitter and receiver increasing.展开更多
We report a broadband terahertz time-domain spectroscopy(THz-TDS)which enables twenty vibrational modes of adenosine nucleoside to be resolved in a wide frequency range of 1-20 THz.The observed spectroscopic features ...We report a broadband terahertz time-domain spectroscopy(THz-TDS)which enables twenty vibrational modes of adenosine nucleoside to be resolved in a wide frequency range of 1-20 THz.The observed spectroscopic features of adenosine are in good agreement with the published spectra obtained using Fourier transform infrared spectroscopy(FTIR)and Raman spectroscopy.This much extended bandwidth leads to enhanced material characterization capability as it provides spectroscopic information on both intra-and inter-molecular vibrations.In addition,we also report a low-cost frequency modulation continuous wave(FMCW)imaging system which has a fast measurement speed of 40000 waveforms per second.Cross-sectional imaging capability through cardboard has also been demonstrated using its excellent penetration capability at a frequency range of 76-81 GHz.We anticipate that the integration of these two complementary imaging technologies would be highly desirable for many real-world applications because it provides both spectroscopic discrimination and penetration capabilities in a single instrument.展开更多
Simulating the coupled motions of multiple bodies in the time domain is a complex problem because of the strong hydrodynamic interactions and coupled effect of various mechanical connectors. In this study, we investig...Simulating the coupled motions of multiple bodies in the time domain is a complex problem because of the strong hydrodynamic interactions and coupled effect of various mechanical connectors. In this study, we investigate the hydrodynamic responses of three barges moored side-by-side in a floatover operation in the frequency and time domains. In the frequency domain, the damping lid method is adopted to improve the overestimated hydrodynamic coefficients calculated from conventional potential flow theory. A time-domain computing program based on potential flow theory and impulse theory is compiled for analyses that consider multibody hydrodynamic interactions and mechanical effects from lines and fenders. Correspondingly, an experiment is carried out for comparison with the numerical results. All statistics, time series, and power density spectra from decay and irregular wave tests are in a fairly good agreement.展开更多
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.展开更多
基金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.
基金Major Program of the National Natural Science Foundation of China under Grant No.52192675 and the 111 Project of China under Grant No.D21001。
文摘Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modulus and combining the dynamic skeleton curve and the damping degradation coefficient,the constitutive equation of the logarithmic dynamic skeleton can be obtained,which considers the damping effect in a soil dynamics problem.Based on the finite difference method and the multi-transmitting boundary condition,a 1D site seismic response analysis program called Soilresp1D has been developed herein and used to analyze the time-domain seismic response in three types of sites.At the same time,this study also provides numerical simulation results based on the hyperbolic constitutive model and the equivalent linear method.The results verify the rationality of the new soil dynamic constitutive model.It can analyze the mucky soil site nonlinear seismic response,reflecting the deformation characteristics and damping effect of the silty soil.The hysteresis loop area is more extensive,and the residual strain is evident.
基金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.
基金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.
基金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 National Natural Science Foundation of China (Grant No. 40974039)High-Tech Research and Development Program of China (Grant No.2006AA06205)Leading Strategic Project of Science and Technology, Chinese Academy of Sciences (XDA08020500)
文摘Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.
基金supported by Chinese National Programs for Fundamental Research and Development(No.2012CB416605)the National Natural Science Foundation of China(No.41174090)Development Project of National Key Scientific Equipment(No.ZDYZ2012-1-05-04)
文摘An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the divole source.
基金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.
基金Project supported by the National Key Research and Development Program of China(Grant No.2023YFA1407000)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDB0460000)+4 种基金the National Natural Science Foundation of China(Grant Nos.12322401,12127807,and 12393832)CAS Key Research Program of Frontier Sciences(Grant No.ZDBS-LY-SLH004)Beijing Nova Program(Grant No.20230484301)Youth Innovation Promotion Association,Chinese Academy of Sciences(Grant No.2023125)CAS Project for Young Scientists in Basic Research(Grant No.YSBR-026)。
文摘Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the broken spatial translational symmetry.However,the intensity of first-order Raman scattering at crystal edges has been rarely explored,although the mechanical stress and edge characteristics have been thoroughly studied by the Raman peak shift and the spectral features of the edge-related Raman modes.Here,by taking Ga As crystal with a well-defined edge as an example,we reveal the intensity enhancement of Raman-active modes and the emergence of Raman-forbidden modes under specific polarization configurations at the edge.This is attributed to the presence of a hot spot at the edge due to the redistributed electromagnetic fields and electromagnetic wave propagations of incident laser and Raman signal near the edge,which are confirmed by the finite-difference time-domain simulations.Spatially-resolved Raman intensities of both Raman-active and Raman-forbidden modes near the edge are calculated based on the redistributed electromagnetic fields,which quantitatively reproduce the corresponding experimental results.These findings offer new insights into the intensity enhancement of Raman scattering at crystal edges and present a new avenue to manipulate light–matter interactions of crystal by manufacturing various types of edges and to characterize the edge structures in photonic and optoelectronic devices.
基金The article is a result of joined research performed during the project:“A Floating Dock Digital Twin towards Efficient,Safer and Autonomous Docking Operations”-NOR/POLNOR/DigiFloDock/0009/2019-00 which is cofinanced by the programme“Applied research”under the Norwegian Financial Mechanisms 2014-2021 POLNOR 2019-Digital and Industry.
文摘During normal de-ballasting operations for floating docks,each ballast pump independently manages a specific group of ballast tanks.However,when a pump malfunctions,a connection valve between the two groups of ballast water systems is opened.This allows the adjacent pump to serve as a helper pump,simultaneously controlling two groups of ballast water systems.This study explores a full-scale floating dock’s dynamic behaviours during the de-ballasting operations under this situation through a numerical model.In the developed numerical model,the dock is described as a six-degree-of-freedom rigid body which is subjected to hydrostatic,hydrodynamic,and mooring loads.A hydraulic model of the piping network of the malfunctioning pump and the helper pump is proposed.A modified P-controller regulates opening angles of all tank valves for minimal pitch and roll.Two configurations of the floating dock,i.e.,a single floating dock and a floating dock with an onboard vessel,are considered.The numerical results show that the optimal helper pumps can be identified regarding the pumps’total de-ballasting capacity and the dock’s stability.The most severe scenarios can be determined in term of the dock’s maximum draught differences caused by its roll and pitch.The observed maximum draught differences remain small relative to the dock’s width,indicating the effectiveness of employing helper pumps and the proposed automatic ballast control strategy for one-pump malfunction scenarios.
基金Project supported by the National Natural Science Foundation of China (Grant No. 12174085)the Key Research and Development Project of Changzhou, Jiangsu Province, China (Grant No. CE20235054)the Postgraduate Research and Practice Innovation Program of Jiangsu Province, China (Grant No. KYCX24 0833)。
文摘Ultrasonic Lamb waves undergo complex mode conversion and diffraction at non-penetrating defects, such as plate corrosion and cracks. Lamb wave imaging has a resolution limit due to the guided wave dispersion characteristics and Rayleigh criterion limitations. In this paper, a full convolutional network is designed to segment and reconstruct the received signals, enabling the automatic identification of target modalities. This approach eliminates clutter and mode conversion interference when calculating direct and accompanying acoustic fields in time-domain topological energy(TDTE) imaging.Subsequently, the measured accompanying acoustic field is reversed for adaptive focusing on defects and enhance the imaging quality. To circumvent the limitations of the Rayleigh criterion, the direct acoustic field and the accompanying acoustic field were fused to characterize the pixel distribution in the imaging region, achieving Lamb wave super-resolution imaging. Experimental results indicate that compared to the sign coherence factor-total focusing method(SCF-TFM),the proposed method achieves a 31.41% improvement in lateral resolution and a 29.53% increase in signal-to-noise ratio for single-blind-hole defects. In the case of multiple-blind-hole defects with spacings greater than the Rayleigh criterion resolution limit, it exhibits a 27.23% enhancement in signal-to-noise ratio. On the contrary, when the defect spacings are relatively smaller than the limit, this method has a higher resolution limit than SCF-TFM in super-resolution imaging.
基金supported by the National Natural Science Foundation of China(Nos.52222602,and52236006)the Fundamental Research Funds for the Central Universities(Nos.FRF-EYIT-23-05,and FRF-TP-22-001C1)+1 种基金Noncommunicable Chronic Diseases-National Science and Technology Major Project(No.2023ZD0500902)the member of the Youth Innovation Promotion Association Foundation of CAS,China(No.2023310)。
文摘Investigating thermal transport mechanisms at the interface between phase change materials(PCMs)and high thermally conductive fillers has become increasingly significant in developing phase change energy storage technologies.This study explores the interfacial thermal transport between a representative PCM,erythritol,and various fillers,including crystalline(Si C,Si_(3)N_(4))and amorphous(Si O_(2))nanoparticles,using molecular dynamics(MD)simulations.Additionally,time-domain thermoreflectance(TDTR)experiments were performed to quantify the interfacial thermal conductance between erythritol and the three types of fillers,yielding values of 50.1,40.0,and25.6 MW m^(–2)K^(-1).These results align well with the trends observed in the simulations.Furthermore,the underlying mechanisms of interfacial heat transfer were analyzed by examining the phonon density of states,overlap energy,and interaction energy.This research provides innovative insights into nanoscale interfacial thermal transport in composite PCMs.This could lead to significant advancements in thermal management technologies,particularly in developing more efficient thermal energy storage systems.
基金financially supported by National Nonprofit institute Research Grant of IGGE(Nos.AS2017J06,AS2017Y04,and AS2016J10)Survey on coastal area for airborne magnetic method of UNV in Jiangsu(No.DD20160151-03)+3 种基金Key National Research Project of China(No.2017YFC0601900)Key Program of National Natural Science Foundation of China(No.41530320)Natural Science Foundation(No.41274121)China Natural Science Foundation for Young Scientists(No.41404093)
文摘Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic （TDAEM） method is still confined to one dimensional （1D） cases, and the corresponding three-dimensional （3D） numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell＇s equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant No.61805214)the Fundamental Research Funds for the Central Universities,China(Grant No.2652017142)
文摘Copper sulfate pentahydrate is investigated by terahertz time-domain spectroscopy. It is shown that the terahertz absorption coefficients are correlated with the particle size of the samples, as well as the heating rates of the ambient temperature. Furthermore, the water molecules of copper sulfate pentahydrate can be quantitatively characterized due to the high sensitivity of the terahertz wave to water molecules. Based on such results, the status of water incorporated in mineral opal is also characterized using terahertz time-domain spectroscopy. It indicates that terahertz technology can be considered as an efficient method to detect the dehydration of minerals.
基金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.
基金This project was supported by the key programof the Nationed Natural Science Foundation of China (60432040)
文摘The path loss analysis model based on 5 rays in vacancy indoor environment is proposed. The relationship between multipath overlapping and the path loss is analyzed mathematically. Time-domain technique is introduced to compute reflection coefficient in a very short time interval. A 5 rays path loss calculation method, which is satisfactory accurate, is developed. 5 typical environments are involved to analyze and generalize the common path loss characteristics in vacancy indoor environment. The simulation result shows that the path loss can be characterized as 3 zones with different path loss exponent as distance between transmitter and receiver increasing.
基金Project supported by the Royal Society and Natural Science Foundation of China(NSFC)International Exchanges Cost Share(IEC\NSFC\181415).
文摘We report a broadband terahertz time-domain spectroscopy(THz-TDS)which enables twenty vibrational modes of adenosine nucleoside to be resolved in a wide frequency range of 1-20 THz.The observed spectroscopic features of adenosine are in good agreement with the published spectra obtained using Fourier transform infrared spectroscopy(FTIR)and Raman spectroscopy.This much extended bandwidth leads to enhanced material characterization capability as it provides spectroscopic information on both intra-and inter-molecular vibrations.In addition,we also report a low-cost frequency modulation continuous wave(FMCW)imaging system which has a fast measurement speed of 40000 waveforms per second.Cross-sectional imaging capability through cardboard has also been demonstrated using its excellent penetration capability at a frequency range of 76-81 GHz.We anticipate that the integration of these two complementary imaging technologies would be highly desirable for many real-world applications because it provides both spectroscopic discrimination and penetration capabilities in a single instrument.
基金financially supported by Lloyd’s Register Foundation(LRF),a UK-registered charity and sole shareholder of Lloyd’s Register Group Ltd.the Youth Innovation Fund of State Key Laboratory of Ocean Engineering(Grant No.GKZD010059-21)
文摘Simulating the coupled motions of multiple bodies in the time domain is a complex problem because of the strong hydrodynamic interactions and coupled effect of various mechanical connectors. In this study, we investigate the hydrodynamic responses of three barges moored side-by-side in a floatover operation in the frequency and time domains. In the frequency domain, the damping lid method is adopted to improve the overestimated hydrodynamic coefficients calculated from conventional potential flow theory. A time-domain computing program based on potential flow theory and impulse theory is compiled for analyses that consider multibody hydrodynamic interactions and mechanical effects from lines and fenders. Correspondingly, an experiment is carried out for comparison with the numerical results. All statistics, time series, and power density spectra from decay and irregular wave tests are in a fairly good agreement.
基金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.
