3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic m...3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.展开更多
For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study prop...For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.展开更多
In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibi...In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.展开更多
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.展开更多
The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy...The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy absorption and its distribution in homogeneousand layered homogenous lossy dielectric spheres are presented,and a comparison of these resultswith analytical solution is given.The calculation is carried out for dielectric cylinder on conduct-ing ground as well,and the results are compared with the image theory.All the computationsshew that the finite-difference time-domain method can give satisfactory results.展开更多
The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces o...The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
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.展开更多
This article presents four techniques for assessing verticality:the plumb line approach,the total station distance technique,the three-point centering method,and the centroid method.Given the significant error associa...This article presents four techniques for assessing verticality:the plumb line approach,the total station distance technique,the three-point centering method,and the centroid method.Given the significant error associated with the total station horizontal distance technique when measuring circular piers,this paper proposes the centroid method.This method calculates verticality by determining the coordinates of the center points at both ends of the pier.Experimental findings indicate that the centroid method achieves accuracy in measuring the verticality of circular piers comparable to the three-point centering method,while offering a faster inspection process.Furthermore,the paper explores the concept of composite verticality and validates the effectiveness of the centroid method in measuring composite verticality and its practical applications through comparative experiments.展开更多
This article aims to enhance seismic hazard assessment methods for Kazakhstan’s seismotectonic conditions.It combines probabilistic seismic hazard analysis(PSHA),ground motion simulation,sitespecific geological and g...This article aims to enhance seismic hazard assessment methods for Kazakhstan’s seismotectonic conditions.It combines probabilistic seismic hazard analysis(PSHA),ground motion simulation,sitespecific geological and geotechnical data analysis,and seismic scenario analysis to develop Probabilistic General Seismic Zoning(GSZ)maps for Kazakhstan and Probabilistic Seismic Microzoning maps for Almaty.These maps align with Eurocode 8 principles,incorporating seismic intensity and engineering parameters like peak ground acceleration(PGA).The new procedure,applied in national projects,has resulted in GSZ maps for the country,seismic microzoning maps for Almaty,and detailed seismic zoning maps for East Kazakhstan.These maps,part of a regulatory document,guide earthquake-resistant design and construction.They offer a comprehensive assessment of seismic hazards,integrating traditional Medvedev-Sponheuer-Karnik(MSK-64)intensity scale points with quantitative parameters like peak ground acceleration.This innovative approach promises to advance methods for quantifying seismic hazards in specific regions.展开更多
The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the M...The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.展开更多
With deep mining of coal mines, prospecting multilayer water-filled goaf has become a new content that results from geophysical exploration in coalfields. The central loop transient electromagnetic （TEM） method is f...With deep mining of coal mines, prospecting multilayer water-filled goaf has become a new content that results from geophysical exploration in coalfields. The central loop transient electromagnetic （TEM） method is favorable for prospecting conductive layers because of the coupling relationship between its field structure and formation. However, the shielding effect of conductive overburden would not only require a longer observation time when prospecting the same depth but also weaken the anomalous response of underlying layers. Through direct time domain numerical simulation and horizontal layered earth forward modeling, this paper estimates the length of observation time required to prospect the target, and the distinguishable criterion of multilayer water-filled goal is presented with observation error according to the effect of noise on observation data. The observed emf curves from Dazigou Coal Mine, Shanxi Province can distinguish multilayer water-filled goaf. In quantitative inversion interpretation of observed curves, using electric logging data as initial parameters restrains the equivalence caused by coal formation thin layers. The deduced three-layer and two-layer water-filled goals are confirmed by the drilling hole. The result suggests that when observation time is long enough and with the anomalous situation of underlying layers being greater than the observation error, the use of the central loop TEM method to orosoect a multilaver water-filled goaf is feasible.展开更多
BACKGROUND Although the number of patients who need central venous ports for permanent vascular access is increasing,there is still no“gold standard”for the implantation technique.AIM To identify the implantation te...BACKGROUND Although the number of patients who need central venous ports for permanent vascular access is increasing,there is still no“gold standard”for the implantation technique.AIM To identify the implantation technique that should be favored.METHODS Two hundred central venous port-implanted patients in a tertiary hospital were retrospectively evaluated.Patients were assigned into two groups according to the access method.The first group comprised patients whose jugular veins were used,and the second group comprised patients whose subclavian veins were used.Groups were evaluated regarding age,sex,application side,primary diagnosis,active follow-up period in the hospital,chemotherapy agents administered,number of complications,and the Clavien-Dindo severity score.The distribution of the variables was tested with the Kolmogorov-Smirnov test and the Mann-Whitney U test.Theχ^(2) test was used to analyze the variables.RESULTS There was no statistically significant difference between the groups regarding age,sex,side,number of chemotherapy drugs,and duration of port usage(P>0.05).Only 2 patients in group 1 had complications,whereas in group 2 we observed 19 patients with complications(P<0.05).No port occlusion was found in group 1,but the catheters of 4 patients were occluded in group 2.One port was infected in group 1 compared to three infected ports in group 2.Two port ruptures,two pneumothorax,one revision due to a mechanical problem,one tachyarrhythmia during implantation,and four suture line problems were also recorded in group 2 patients.We also showed that it would be sufficient to evaluate and wash ports once every 2 mo.CONCLUSION Our results robustly confirm that the jugular vein route is safer than the subclavian vein approach for central venous port implantation.展开更多
The prediction of central bursting defects in the rod extrusion process through conical dies using the upper bound analysisis investigated. A kinematically admissible velocity field, including the radial and angular v...The prediction of central bursting defects in the rod extrusion process through conical dies using the upper bound analysisis investigated. A kinematically admissible velocity field, including the radial and angular velocity components, is proposed. A newcriterion is presented to predict the occurrence of the central bursting defects. Parameter bobt, which represents the risk probability ofcracking, is proposed. It is calculated using the shape of the boundary at the entrance by minimizing the total power dissipationduring the extrusion process. When bobt is equal to or greater than bcr, central bursting occurs. Furthermore, the quantitativerelationships between central bursting defects and process parameters （semi die angle, reduction in area and frictional factor） arestudied. The results show that the central bursting defects are affected primarily by the reduction in area and the friction factor. Thepresented criterion is verified by comparing with the FEM simulation data and the results of the published paper.展开更多
The relativistic Hartree-Bogoliubov(RHB)theory is a powerful tool for describing exotic nuclei near drip lines.The key technique is to solve the RHB equation in the coordinate space to obtain the quasi-particle states...The relativistic Hartree-Bogoliubov(RHB)theory is a powerful tool for describing exotic nuclei near drip lines.The key technique is to solve the RHB equation in the coordinate space to obtain the quasi-particle states.In this paper,we solve the RHB equation with the Woods-Saxon-type mean-field and Delta-type pairing-field potentials by using the finite-difference method(FDM).We inevitably obtain spurious states when using the common symmetric central difference formula(CDF)to construct the Hamiltonian matrix,which is similar to the problem resulting from solving the Dirac equation with the same method.This problem is solved by using the asymmetric difference formula(ADF).In addition,we show that a large enough box is necessary to describe the continuum quasi-particle states.The canonical states obtained by diagonalizing the density matrix constructed by the quasi-particle states are not particularly sensitive to the box size.Part of the asymptotic wave functions can be improved by applying the ADF in the FDM compared to the shooting method with the same box boundary condition.展开更多
Despite the CaCO<sub>3</sub> estimation using titration method was not reliable, but up to the present time, some soil laboratories in Sudan still used this method. The objective of this study was to compa...Despite the CaCO<sub>3</sub> estimation using titration method was not reliable, but up to the present time, some soil laboratories in Sudan still used this method. The objective of this study was to compare and assess the results of calcimetric and titrimetric methods of quantitative estimation for soil calcium carbonate of different soils in Sudan. 26 soil samples from five soil profiles were collected from different climatological and ecological regions in central Sudan. CaCO<sub>3</sub> equivalent was estimated using calcimeter and titration methods in order to find accurate, rapid and suitable method for soils of Sudan. The results revealed that there are no significant differences between calcimeter and titration methods for calcium carbonate estimation in all studied samples except in samples from Gedaref area. We concluded that when the Calcimeter method used for CaCO<sub>3</sub> estimation, the differences between one person and another in detecting titration end point would be avoided, rapid and accurate results would be obtained compared to titration method. Additionally, time would be saved;fewer amounts of chemicals would be used. From this study, we highly recommend using calcimeter method for CaCO<sub>3</sub> estimation for soils of Sudan.展开更多
[Objectives] To optimize the extraction process of total flavonoids in stems of Mallotus apelta. [Methods]On the basis of singlefactor test,with volume fraction of ethanol,extraction time and ratio of solvent as indep...[Objectives] To optimize the extraction process of total flavonoids in stems of Mallotus apelta. [Methods]On the basis of singlefactor test,with volume fraction of ethanol,extraction time and ratio of solvent as independent variables,the content of total flavonoids as dependent variables,the completely secondary response surface regression fitting was conducted on the independent and dependent variables,and the Response Surface Method was used to optimize the optimum extraction process of total flavonoids in Mallotus apelta stems and predict the optimum process. [Results] The optimum extraction process of total flavonoids in Mallotus apelta was determined as follows: ethanol concentration of 71. 5%; extraction time of 154. 6 min; solid-liquid ratio of 1∶19. 2; total flavonoids content of 7. 060 mg/g; fitted binomial squared correlation coefficient R^2= 0. 8751.[Conclusions]Composite Design/Response Surface Method could be used in the extraction process optimization of total flavonoids in Mallotus apelta stems,the mathematical model established had high prediction accuracy,the method was simple and operability was good.展开更多
So far, Schrodinger equation with central potential has been solved in different methods but solving this equation with non-central potentials is less dealt with. Solving such equations are way more difficult and comp...So far, Schrodinger equation with central potential has been solved in different methods but solving this equation with non-central potentials is less dealt with. Solving such equations are way more difficult and complicated and a certain and limited number of non-central potentials can be solved. In this paper, we introduce one of the solvable kinds of such potentials and we will use NU method for solving Schrodinger equation and then by using this method we have calculated particular figures of its energy and function.展开更多
基金The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04National Natural Science Foundation of China(42074132,42274154).
文摘3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.
基金National Natural Science Foundation of China under Grant Nos.51978213 and 51778190the National Key Research and Development Program of China under Grant Nos.2017YFC0703605 and 2016YFC0701106。
文摘For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.
基金the National Natural Science Foundation of China(No.40774056)Program of Excellent Team in Harbin Institute of Technology
文摘In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.
基金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.
文摘The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy absorption and its distribution in homogeneousand layered homogenous lossy dielectric spheres are presented,and a comparison of these resultswith analytical solution is given.The calculation is carried out for dielectric cylinder on conduct-ing ground as well,and the results are compared with the image theory.All the computationsshew that the finite-difference time-domain method can give satisfactory results.
基金Taif University Researchers Supporting Project Number(TURSP-2020/164),Taif University,Taif,Saudi Arabia.
文摘The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
基金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 article presents four techniques for assessing verticality:the plumb line approach,the total station distance technique,the three-point centering method,and the centroid method.Given the significant error associated with the total station horizontal distance technique when measuring circular piers,this paper proposes the centroid method.This method calculates verticality by determining the coordinates of the center points at both ends of the pier.Experimental findings indicate that the centroid method achieves accuracy in measuring the verticality of circular piers comparable to the three-point centering method,while offering a faster inspection process.Furthermore,the paper explores the concept of composite verticality and validates the effectiveness of the centroid method in measuring composite verticality and its practical applications through comparative experiments.
基金The work was carried out in the framework of earmarked funding“Assessment of seismic hazard of territories of Kazakhstan on modern scientific and methodological basis”,programme code number F.0980.Source of funding-Ministry of Science and Higher Education of the Republic of Kazakhstan.
文摘This article aims to enhance seismic hazard assessment methods for Kazakhstan’s seismotectonic conditions.It combines probabilistic seismic hazard analysis(PSHA),ground motion simulation,sitespecific geological and geotechnical data analysis,and seismic scenario analysis to develop Probabilistic General Seismic Zoning(GSZ)maps for Kazakhstan and Probabilistic Seismic Microzoning maps for Almaty.These maps align with Eurocode 8 principles,incorporating seismic intensity and engineering parameters like peak ground acceleration(PGA).The new procedure,applied in national projects,has resulted in GSZ maps for the country,seismic microzoning maps for Almaty,and detailed seismic zoning maps for East Kazakhstan.These maps,part of a regulatory document,guide earthquake-resistant design and construction.They offer a comprehensive assessment of seismic hazards,integrating traditional Medvedev-Sponheuer-Karnik(MSK-64)intensity scale points with quantitative parameters like peak ground acceleration.This innovative approach promises to advance methods for quantifying seismic hazards in specific regions.
基金supported by the Natural Science Foundation of China(No.41374078)Geological Survey Projects of Ministry of Land and Resources of China(No.12120113086100 and 12120113101300)
文摘The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.
基金supported by the National Science Foundation of China(No.41374129)Science and Technology Project of Shanxi Province(No.20100321066)Research and Development Project of National Major Scientifi c Research Equipment(No.ZDYZ2012-1-05-04)
文摘With deep mining of coal mines, prospecting multilayer water-filled goaf has become a new content that results from geophysical exploration in coalfields. The central loop transient electromagnetic （TEM） method is favorable for prospecting conductive layers because of the coupling relationship between its field structure and formation. However, the shielding effect of conductive overburden would not only require a longer observation time when prospecting the same depth but also weaken the anomalous response of underlying layers. Through direct time domain numerical simulation and horizontal layered earth forward modeling, this paper estimates the length of observation time required to prospect the target, and the distinguishable criterion of multilayer water-filled goal is presented with observation error according to the effect of noise on observation data. The observed emf curves from Dazigou Coal Mine, Shanxi Province can distinguish multilayer water-filled goaf. In quantitative inversion interpretation of observed curves, using electric logging data as initial parameters restrains the equivalence caused by coal formation thin layers. The deduced three-layer and two-layer water-filled goals are confirmed by the drilling hole. The result suggests that when observation time is long enough and with the anomalous situation of underlying layers being greater than the observation error, the use of the central loop TEM method to orosoect a multilaver water-filled goaf is feasible.
文摘BACKGROUND Although the number of patients who need central venous ports for permanent vascular access is increasing,there is still no“gold standard”for the implantation technique.AIM To identify the implantation technique that should be favored.METHODS Two hundred central venous port-implanted patients in a tertiary hospital were retrospectively evaluated.Patients were assigned into two groups according to the access method.The first group comprised patients whose jugular veins were used,and the second group comprised patients whose subclavian veins were used.Groups were evaluated regarding age,sex,application side,primary diagnosis,active follow-up period in the hospital,chemotherapy agents administered,number of complications,and the Clavien-Dindo severity score.The distribution of the variables was tested with the Kolmogorov-Smirnov test and the Mann-Whitney U test.Theχ^(2) test was used to analyze the variables.RESULTS There was no statistically significant difference between the groups regarding age,sex,side,number of chemotherapy drugs,and duration of port usage(P>0.05).Only 2 patients in group 1 had complications,whereas in group 2 we observed 19 patients with complications(P<0.05).No port occlusion was found in group 1,but the catheters of 4 patients were occluded in group 2.One port was infected in group 1 compared to three infected ports in group 2.Two port ruptures,two pneumothorax,one revision due to a mechanical problem,one tachyarrhythmia during implantation,and four suture line problems were also recorded in group 2 patients.We also showed that it would be sufficient to evaluate and wash ports once every 2 mo.CONCLUSION Our results robustly confirm that the jugular vein route is safer than the subclavian vein approach for central venous port implantation.
文摘The prediction of central bursting defects in the rod extrusion process through conical dies using the upper bound analysisis investigated. A kinematically admissible velocity field, including the radial and angular velocity components, is proposed. A newcriterion is presented to predict the occurrence of the central bursting defects. Parameter bobt, which represents the risk probability ofcracking, is proposed. It is calculated using the shape of the boundary at the entrance by minimizing the total power dissipationduring the extrusion process. When bobt is equal to or greater than bcr, central bursting occurs. Furthermore, the quantitativerelationships between central bursting defects and process parameters （semi die angle, reduction in area and frictional factor） arestudied. The results show that the central bursting defects are affected primarily by the reduction in area and the friction factor. Thepresented criterion is verified by comparing with the FEM simulation data and the results of the published paper.
基金Supported by the National Natural Science Foundation of China(11775119,2175109)the Natural Science Foundation of Tianjin,China(19JCYBJC30800)。
文摘The relativistic Hartree-Bogoliubov(RHB)theory is a powerful tool for describing exotic nuclei near drip lines.The key technique is to solve the RHB equation in the coordinate space to obtain the quasi-particle states.In this paper,we solve the RHB equation with the Woods-Saxon-type mean-field and Delta-type pairing-field potentials by using the finite-difference method(FDM).We inevitably obtain spurious states when using the common symmetric central difference formula(CDF)to construct the Hamiltonian matrix,which is similar to the problem resulting from solving the Dirac equation with the same method.This problem is solved by using the asymmetric difference formula(ADF).In addition,we show that a large enough box is necessary to describe the continuum quasi-particle states.The canonical states obtained by diagonalizing the density matrix constructed by the quasi-particle states are not particularly sensitive to the box size.Part of the asymptotic wave functions can be improved by applying the ADF in the FDM compared to the shooting method with the same box boundary condition.
文摘Despite the CaCO<sub>3</sub> estimation using titration method was not reliable, but up to the present time, some soil laboratories in Sudan still used this method. The objective of this study was to compare and assess the results of calcimetric and titrimetric methods of quantitative estimation for soil calcium carbonate of different soils in Sudan. 26 soil samples from five soil profiles were collected from different climatological and ecological regions in central Sudan. CaCO<sub>3</sub> equivalent was estimated using calcimeter and titration methods in order to find accurate, rapid and suitable method for soils of Sudan. The results revealed that there are no significant differences between calcimeter and titration methods for calcium carbonate estimation in all studied samples except in samples from Gedaref area. We concluded that when the Calcimeter method used for CaCO<sub>3</sub> estimation, the differences between one person and another in detecting titration end point would be avoided, rapid and accurate results would be obtained compared to titration method. Additionally, time would be saved;fewer amounts of chemicals would be used. From this study, we highly recommend using calcimeter method for CaCO<sub>3</sub> estimation for soils of Sudan.
基金Supported by Chinese Medicine Science and Technology Project of Guangxi Administration of Traditional Chinese Medicine(GZLC14-31)Science and Technology Research and Development Program of Guilin Bureau of Technology(20130403-4)+1 种基金Guangxi"2011 Collaborative Innovation Center"-Zhuang and Yao Medicine Collaborative Innovation Center(Gui201320)the Autonomous Region-Level College Students’ Innovation and Entrepreneurship Training Program(201710601082)
文摘[Objectives] To optimize the extraction process of total flavonoids in stems of Mallotus apelta. [Methods]On the basis of singlefactor test,with volume fraction of ethanol,extraction time and ratio of solvent as independent variables,the content of total flavonoids as dependent variables,the completely secondary response surface regression fitting was conducted on the independent and dependent variables,and the Response Surface Method was used to optimize the optimum extraction process of total flavonoids in Mallotus apelta stems and predict the optimum process. [Results] The optimum extraction process of total flavonoids in Mallotus apelta was determined as follows: ethanol concentration of 71. 5%; extraction time of 154. 6 min; solid-liquid ratio of 1∶19. 2; total flavonoids content of 7. 060 mg/g; fitted binomial squared correlation coefficient R^2= 0. 8751.[Conclusions]Composite Design/Response Surface Method could be used in the extraction process optimization of total flavonoids in Mallotus apelta stems,the mathematical model established had high prediction accuracy,the method was simple and operability was good.
文摘So far, Schrodinger equation with central potential has been solved in different methods but solving this equation with non-central potentials is less dealt with. Solving such equations are way more difficult and complicated and a certain and limited number of non-central potentials can be solved. In this paper, we introduce one of the solvable kinds of such potentials and we will use NU method for solving Schrodinger equation and then by using this method we have calculated particular figures of its energy and function.
