The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity mo...The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity model used for the conventional intersection method is based on homogeneous or laterally homogeneous media, which is too simple. In order to improve the accuracy, we have modified the existing intersection method. In the modified approach, the earthquake loci are not assumed to be circular or hyperbolic and calculation accuracy is improved using a minimum traveltime tree algorithm for tracing rays. The numerical model shows that the modified method can locate earthquakes in complex velocity models.展开更多
This paper introduces part of the content in the association standard,T/CAAM0002–2020 Nomenclature and Location of Acupuncture Points for Laboratory Animals Part 3:Mouse.This standard was released by the China Associ...This paper introduces part of the content in the association standard,T/CAAM0002–2020 Nomenclature and Location of Acupuncture Points for Laboratory Animals Part 3:Mouse.This standard was released by the China Association of Acupuncture and Moxibustion on May 15,2020,implemented on October 31,2020,and published by Standards Press of China.The standard was drafted by the Institute of Acupuncture and Moxibustion,China Academy of Chinese Medical Sciences,and the Nanjing University of Chinese Medicine.Principal draftsmen:Xiang-hong JING and Xing-bang HUA.Participating draftsmen:Wan-zhu BAI,Bin XU,Dong-sheng XU,Yi GUO,Tie-ming MA,Xin-jun WANG,and Sheng-feng LU.展开更多
The efficiency of an optimization method for acoustic emission/microseismic（AE/MS） source location is determined by the compatibility of its error definition with the errors contained in the input data.This compatib...The efficiency of an optimization method for acoustic emission/microseismic（AE/MS） source location is determined by the compatibility of its error definition with the errors contained in the input data.This compatibility can be examined in terms of the distribution of station residuals.For an ideal distribution,the input error is held at the station where it takes place as the station residual and the error is not permitted to spread to other stations.A comparison study of two optimization methods,namely the least squares method and the absolute value method,shows that the distribution with this character constrains the input errors and minimizes their impact,which explains the much more robust performance by the absolute value method in dealing with large and isolated input errors.When the errors in the input data are systematic and/or extreme in that the basic data structure is altered by these errors,none of the optimization methods are able to function.The only means to resolve this problem is the early detection and correction of these errors through a data screening process.An efficient data screening process is of primary importance for AE/MS source location.In addition to its critical role in dealing with those systematic and extreme errors,data screening creates a favorable environment for applying optimization methods.展开更多
This paper introduces part of the content in the association standard,T/CAAM0002–2020 Nomenclature and Location of Acupuncture Points for Laboratory Animals Part 2:Rat.This standard was released by the China Associat...This paper introduces part of the content in the association standard,T/CAAM0002–2020 Nomenclature and Location of Acupuncture Points for Laboratory Animals Part 2:Rat.This standard was released by the China Association of Acupuncture and Moxibustion on May 15,2020,implemented on October 31,2020,and published by Standards Press of China.The standard was drafted by the Institute of Acupuncture and Moxibustion,China Academy of Chinese Medical Sciences,and the Nanjing University of Chinese Medicine.Principal draftsmen:Xiang-hong JING and Xing-bang HUA.Participating draftsmen:Wan-Zhu BAI,Bin XU,Dong-sheng XU,Yi GUO,Tie-ming MA,Xin-jun WANG,and Sheng-feng LU.展开更多
In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of th...In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of the mappings which solves some variational inequality. The results presented extend the corresponding results of Shimizu and Takahashi IT. Shimizu, W. Takahashi, Strong convergence to common fixed point of families of nonexpansive mappings, J. Math. Anal. Appl. 211 （1997）, 71-83], and Yao and Chen [Y. Yao, R. Chert, Convergence to common fixed points of average mappings without commutativity assumption in Hilbert spaces, Nonlinear Analysis 67（2007）, 1758-1763].展开更多
Geo-interfaces refer to the contact surfaces between multiple media within geological strata,as well as the transition zones that regulate the migration of three-phase matter,changes in physical states,and the deforma...Geo-interfaces refer to the contact surfaces between multiple media within geological strata,as well as the transition zones that regulate the migration of three-phase matter,changes in physical states,and the deformation and stability of rock and soil masses.Owing to the combined effects of natural factors and human activities,geo-interfaces play crucial roles in the emergence,propagation,and triggering of geological disasters.Over the past three decades,the material point method(MPM)has emerged as a preferred approach for addressing large deformation problems and simulating soil-water-structure interactions,making it an ideal tool for analyzing geo-interface behaviors.In this review,we offer a systematic summary of the basic concepts,classifications,and main characteristics of the geo-interface,and provide a comprehensive overview of recent advances and developments in simulating geo-interface using the MPM.We further present a brief description of various MPMs for modeling different types of geo-interfaces in geotechnical engineering applications and highlight the existing limitations and future research directions.This study aims to facilitate innovative applications of the MPM in modeling complex geo-interface problems,providing a reference for geotechnical practitioners and researchers.展开更多
The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their ...The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report.展开更多
Fixture locating layout has a direct and influential impact on aeronautical thin-walled component(ATWC)manufacturing quality.The purpose is to develop a topological optimization method for ATWC fixture locating layout...Fixture locating layout has a direct and influential impact on aeronautical thin-walled component(ATWC)manufacturing quality.The purpose is to develop a topological optimization method for ATWC fixture locating layout to minimize the manufacturing deformation.Firstly,a topological optimization model that takes the stiffness of ATWC as the objective function and the volume of the locating structure as the constraint is established.Secondly,ATWC and the locating structure are regarded as an integrated entity,and the variable-density method based topological optimization approach is adopted for the optimization of the locating structure using ABAQUS topology optimization module(ATOM).Thirdly,through a subsequent model reconstruction referring to the obtained topological structure,the optimal fixture locating layout is achieved.Finally,a case study is conducted to verify the proposed method and the comparison results with firefly algorithm(FA)coupled with finite element analysis(FEA)indicate that the number and positions of the locators for ATWC can be optimized simultaneously and successfully by the proposed topological optimization model.展开更多
Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of...Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.展开更多
Tuina is a traditional Chinese treatment for sensory disturbances caused by peripheral nerve injury and related diseases. Our previous studies showed that tuina regulates relevant regions and indices of the spinal dor...Tuina is a traditional Chinese treatment for sensory disturbances caused by peripheral nerve injury and related diseases. Our previous studies showed that tuina regulates relevant regions and indices of the spinal dorsal horn using the Dian, Bo, and Rou method in Yinmen（BL37）, Yanglingquan（GB34）, and Weizhong（BL40）. Treatment prevents muscle atrophy, protects spinal cord neurons, and promotes sciatic nerve repair. The mechanisms of action of tuina for treating peripheral nerve injury remain poorly understood. This study established rat models of sciatic nerve injury using the crushing method. Rats received Chinese tuina in accordance with the principle of ＂Three Methods and Three Points,＂ once daily for 20 days. Tuina intervention reduced paw withdrawal latency and improved wet weight of the gastrocnemius muscle, as well as promoting morphological recovery of sciatic nerve fibers, Schwann cells, and axons. The protein expression levels of phospho-p38 mitogen-activated protein kinase, tumor necrosis factor-α, and interleukin-1β also decreased. These findings indicate that ＂Three Methods and Three Points＂ promoted morphological recovery and improved behavior of rats with peripheral nerve injury.展开更多
According to earthquake catalog records of Fujian Seismic Network, the Tnow method and the fourstation continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first ...According to earthquake catalog records of Fujian Seismic Network, the Tnow method and the fourstation continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations in each earthquake. It shows that the fourstation continuous location method can locate more seismic events than the Tnow method. By analyzing the results, it is concluded that the reason for this is that the Tnow method makes use of information from stations without being triggered, while some stations failed to be reflected in earthquake catalog because of discontinuous records or unclear records of seismic phases. For seismic events whose location results can be given, there is no obvious difference in location results of the two methods and positioning deviation of most seismic events is also not significant. For earthquakes outside the network, the positioning deviation may amplify as the epicentral distance enlarges, which may relate to the situation that the seismic stations are centered on one side of epicenter and the opening angle between seismic stations used for location and epicenter is small.展开更多
Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicit...Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicity and quick progress in the early iterations. In this work, to accelerate the convergence of the interior point method, few iterations of this generalized algorithm are applied to the Mehrotra’s heuristic, which determines the starting point for the interior point method in the PCx software. Computational experiments in a set of linear programming problems have shown that this approach reduces the total number of iterations and the running time for many of them, including large-scale ones.展开更多
The aim of this paper is to obtain the numerical solutions of generalized space-fractional Burgers' equations with initial-boundary conditions by the Jacobi spectral collocation method using the shifted Jacobi-Gau...The aim of this paper is to obtain the numerical solutions of generalized space-fractional Burgers' equations with initial-boundary conditions by the Jacobi spectral collocation method using the shifted Jacobi-Gauss-Lobatto collocation points. By means of the simplifed Jacobi operational matrix, we produce the diferentiation matrix and transfer the space-fractional Burgers' equation into a system of ordinary diferential equations that can be solved by the fourth-order Runge-Kutta method. The numerical simulations indicate that the Jacobi spectral collocation method is highly accurate and fast convergent for the generalized space-fractional Burgers' equation.展开更多
This paper presents a class of r-point (r+1)st-order A-stable one-block methods with damping at the infinite point (DIAOB r,r+1 ). Under the conditions of the same order, A-stability, operation count (at each iterativ...This paper presents a class of r-point (r+1)st-order A-stable one-block methods with damping at the infinite point (DIAOB r,r+1 ). Under the conditions of the same order, A-stability, operation count (at each iterative step) and storage space are the same as the methods in , the methods in the paper improve the stability in a neighborhood at the infinite point. And, by using the OOPI method , it possesses much faster rate of convergence for solving systems of nonlinear equations produced by the DIAOB r,r+1 .展开更多
In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this sche...In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this scheme seamlessly combines the material point method(MPM)for accurately tracking solid-phase deformation and heat transport,and the Eulerian finite element method(FEM)for effectively capturing fluid flow and heat advection-diffusion behavior.The proposed approach circumvents the substantial challenges posed by large nonlinear equation systems with the monolithic solution scheme.The staggered solution process strategically separates each physical field through explicit or implicit integration.The characteristic-based method is used to stabilize advection-dominated heat flows for efficient numerical implementation.Furthermore,a fractional step approach is employed to decompose fluid velocity and pressure,thereby suppressing pore pressure oscillation on the linear background grid.The fracturing initiation and propagation are simulated by a rate-dependent phase field model.Through a series of quasi-static and transient simulations,the exceptional performance and promising potential of the proposed model in addressing THM fracturing problems in poro-elastic media is demonstrated.展开更多
A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a...A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a background mesh. It is shown that, in a point collocation approach, the remesh problem because of the mesh distortion in FEM （finite element method） and the low efficiency in Galerkin-based meshfree method are avoided. The corrected kernel functions are introduced to the stabilization of free-surface boundary conditions. The solution of symmetric ring compression problem is compared with a conventional finite element solution, and reasonable results have been obtained.展开更多
This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-...This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-oped using the material point method.To reduce the computational cost of Monte Carlo simulations,response surface models are created as surrogate models for the material point system to approximate its dynamic behavior.An adaptive randomized greedy algorithm is employed to construct a sparse polynomial chaos expansion model with a fixed order,effectively balancing the accuracy and computational efficiency of the surrogate model.Based on the sparse polynomial chaos expansion,sensitivity analysis is conducted using the global finite difference and Sobol methods.Several examples of structural dynamics are provided to demonstrate the effectiveness of the proposed method in addressing structural dynamics problems.展开更多
The mathematical method of ZTD(zenith tropospheric delay)spatial prediction is important for precise ZTD derivation and real-time precise point positioning(PPP)augmentation.This paper analyses the performance of the p...The mathematical method of ZTD(zenith tropospheric delay)spatial prediction is important for precise ZTD derivation and real-time precise point positioning(PPP)augmentation.This paper analyses the performance of the popular optimal function coefficient(OFC),sphere cap harmonic analysis(SCHA),kriging and inverse distance weighting(IDW)interpolation in ZTD spatial prediction and Beidou satellite navigation system(BDS)-PPP augmentation over China.For ZTD spatial prediction,the average time consumption of the OFC,kriging,and IDW methods is less than 0.1 s,which is significantly better than that of the SCHA method(63.157 s).The overall ZTD precision of the OFC is 3.44 cm,which outperforms those of the SCHA(9.65 cm),Kriging(10.6 cm),and IDW(11.8 cm)methods.We confirmed that the low performance of kriging and IDW is caused by their weakness in modelling ZTD variation in the vertical direction.To mitigate such deficiencies,an elevation normalization factor(ENF)is introduced into the kriging and IDW models(kriging-ENF and IDW-ENF).The overall ZTD spatial prediction accuracies of IDW-ENF and kriging-ENF are 2.80 cm and 2.01 cm,respectively,which are both superior to those of the OFC and the widely used empirical model GPT3(4.92 cm).For BDS-PPP enhancement,the ZTD provided by the kriging-ENF,IDW-ENF and OFC as prior constraints can effectively reduce the convergence time.Compared with unconstrained BDS-PPP,our proposed kriging-ENF outperforms IDW-ENF and OFC by reducing the horizontal and vertical convergence times by approximately 13.2%and 5.8%in Ningxia and 30.4%and 7.84%in Guangdong,respectively.These results indicate that kriging-ENF is a promising method for ZTD spatial prediction and BDS-PPP enhancement over China.展开更多
基金This work is supported by the National Natural Science Foundation of China(40674044)the Special Foundation for Basic Professional Scientific Research (DQJB06A02)
文摘The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity model used for the conventional intersection method is based on homogeneous or laterally homogeneous media, which is too simple. In order to improve the accuracy, we have modified the existing intersection method. In the modified approach, the earthquake loci are not assumed to be circular or hyperbolic and calculation accuracy is improved using a minimum traveltime tree algorithm for tracing rays. The numerical model shows that the modified method can locate earthquakes in complex velocity models.
文摘This paper introduces part of the content in the association standard,T/CAAM0002–2020 Nomenclature and Location of Acupuncture Points for Laboratory Animals Part 3:Mouse.This standard was released by the China Association of Acupuncture and Moxibustion on May 15,2020,implemented on October 31,2020,and published by Standards Press of China.The standard was drafted by the Institute of Acupuncture and Moxibustion,China Academy of Chinese Medical Sciences,and the Nanjing University of Chinese Medicine.Principal draftsmen:Xiang-hong JING and Xing-bang HUA.Participating draftsmen:Wan-zhu BAI,Bin XU,Dong-sheng XU,Yi GUO,Tie-ming MA,Xin-jun WANG,and Sheng-feng LU.
文摘The efficiency of an optimization method for acoustic emission/microseismic（AE/MS） source location is determined by the compatibility of its error definition with the errors contained in the input data.This compatibility can be examined in terms of the distribution of station residuals.For an ideal distribution,the input error is held at the station where it takes place as the station residual and the error is not permitted to spread to other stations.A comparison study of two optimization methods,namely the least squares method and the absolute value method,shows that the distribution with this character constrains the input errors and minimizes their impact,which explains the much more robust performance by the absolute value method in dealing with large and isolated input errors.When the errors in the input data are systematic and/or extreme in that the basic data structure is altered by these errors,none of the optimization methods are able to function.The only means to resolve this problem is the early detection and correction of these errors through a data screening process.An efficient data screening process is of primary importance for AE/MS source location.In addition to its critical role in dealing with those systematic and extreme errors,data screening creates a favorable environment for applying optimization methods.
文摘This paper introduces part of the content in the association standard,T/CAAM0002–2020 Nomenclature and Location of Acupuncture Points for Laboratory Animals Part 2:Rat.This standard was released by the China Association of Acupuncture and Moxibustion on May 15,2020,implemented on October 31,2020,and published by Standards Press of China.The standard was drafted by the Institute of Acupuncture and Moxibustion,China Academy of Chinese Medical Sciences,and the Nanjing University of Chinese Medicine.Principal draftsmen:Xiang-hong JING and Xing-bang HUA.Participating draftsmen:Wan-Zhu BAI,Bin XU,Dong-sheng XU,Yi GUO,Tie-ming MA,Xin-jun WANG,and Sheng-feng LU.
基金the Thailand Research Fund for financial support under Grant BRG5280016
文摘In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of the mappings which solves some variational inequality. The results presented extend the corresponding results of Shimizu and Takahashi IT. Shimizu, W. Takahashi, Strong convergence to common fixed point of families of nonexpansive mappings, J. Math. Anal. Appl. 211 （1997）, 71-83], and Yao and Chen [Y. Yao, R. Chert, Convergence to common fixed points of average mappings without commutativity assumption in Hilbert spaces, Nonlinear Analysis 67（2007）, 1758-1763].
基金supported by the National Science Fund for Distinguished Young Scholars of China(Grant No.42225702)the National Natural Science Foundation of China(Grant Nos.42461160266 and 52379106).
文摘Geo-interfaces refer to the contact surfaces between multiple media within geological strata,as well as the transition zones that regulate the migration of three-phase matter,changes in physical states,and the deformation and stability of rock and soil masses.Owing to the combined effects of natural factors and human activities,geo-interfaces play crucial roles in the emergence,propagation,and triggering of geological disasters.Over the past three decades,the material point method(MPM)has emerged as a preferred approach for addressing large deformation problems and simulating soil-water-structure interactions,making it an ideal tool for analyzing geo-interface behaviors.In this review,we offer a systematic summary of the basic concepts,classifications,and main characteristics of the geo-interface,and provide a comprehensive overview of recent advances and developments in simulating geo-interface using the MPM.We further present a brief description of various MPMs for modeling different types of geo-interfaces in geotechnical engineering applications and highlight the existing limitations and future research directions.This study aims to facilitate innovative applications of the MPM in modeling complex geo-interface problems,providing a reference for geotechnical practitioners and researchers.
文摘The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report.
基金supported by the National Natural Science Foundation of China(No.51375396)the Shaanxi Science and Technology Innovation Project Plan,China(No.2016KTCQ01-50)
文摘Fixture locating layout has a direct and influential impact on aeronautical thin-walled component(ATWC)manufacturing quality.The purpose is to develop a topological optimization method for ATWC fixture locating layout to minimize the manufacturing deformation.Firstly,a topological optimization model that takes the stiffness of ATWC as the objective function and the volume of the locating structure as the constraint is established.Secondly,ATWC and the locating structure are regarded as an integrated entity,and the variable-density method based topological optimization approach is adopted for the optimization of the locating structure using ABAQUS topology optimization module(ATOM).Thirdly,through a subsequent model reconstruction referring to the obtained topological structure,the optimal fixture locating layout is achieved.Finally,a case study is conducted to verify the proposed method and the comparison results with firefly algorithm(FA)coupled with finite element analysis(FEA)indicate that the number and positions of the locators for ATWC can be optimized simultaneously and successfully by the proposed topological optimization model.
文摘Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.
基金supported by the National Natural Science Foundation of China,No.81373759the Natural Science Foundation of Beijing of China,No.7142097
文摘Tuina is a traditional Chinese treatment for sensory disturbances caused by peripheral nerve injury and related diseases. Our previous studies showed that tuina regulates relevant regions and indices of the spinal dorsal horn using the Dian, Bo, and Rou method in Yinmen（BL37）, Yanglingquan（GB34）, and Weizhong（BL40）. Treatment prevents muscle atrophy, protects spinal cord neurons, and promotes sciatic nerve repair. The mechanisms of action of tuina for treating peripheral nerve injury remain poorly understood. This study established rat models of sciatic nerve injury using the crushing method. Rats received Chinese tuina in accordance with the principle of ＂Three Methods and Three Points,＂ once daily for 20 days. Tuina intervention reduced paw withdrawal latency and improved wet weight of the gastrocnemius muscle, as well as promoting morphological recovery of sciatic nerve fibers, Schwann cells, and axons. The protein expression levels of phospho-p38 mitogen-activated protein kinase, tumor necrosis factor-α, and interleukin-1β also decreased. These findings indicate that ＂Three Methods and Three Points＂ promoted morphological recovery and improved behavior of rats with peripheral nerve injury.
文摘According to earthquake catalog records of Fujian Seismic Network, the Tnow method and the fourstation continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations in each earthquake. It shows that the fourstation continuous location method can locate more seismic events than the Tnow method. By analyzing the results, it is concluded that the reason for this is that the Tnow method makes use of information from stations without being triggered, while some stations failed to be reflected in earthquake catalog because of discontinuous records or unclear records of seismic phases. For seismic events whose location results can be given, there is no obvious difference in location results of the two methods and positioning deviation of most seismic events is also not significant. For earthquakes outside the network, the positioning deviation may amplify as the epicentral distance enlarges, which may relate to the situation that the seismic stations are centered on one side of epicenter and the opening angle between seismic stations used for location and epicenter is small.
文摘Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicity and quick progress in the early iterations. In this work, to accelerate the convergence of the interior point method, few iterations of this generalized algorithm are applied to the Mehrotra’s heuristic, which determines the starting point for the interior point method in the PCx software. Computational experiments in a set of linear programming problems have shown that this approach reduces the total number of iterations and the running time for many of them, including large-scale ones.
基金This work is supported by the National Natural Science Foundation of China(Grant Nos.11701358,11774218)。
文摘The aim of this paper is to obtain the numerical solutions of generalized space-fractional Burgers' equations with initial-boundary conditions by the Jacobi spectral collocation method using the shifted Jacobi-Gauss-Lobatto collocation points. By means of the simplifed Jacobi operational matrix, we produce the diferentiation matrix and transfer the space-fractional Burgers' equation into a system of ordinary diferential equations that can be solved by the fourth-order Runge-Kutta method. The numerical simulations indicate that the Jacobi spectral collocation method is highly accurate and fast convergent for the generalized space-fractional Burgers' equation.
文摘This paper presents a class of r-point (r+1)st-order A-stable one-block methods with damping at the infinite point (DIAOB r,r+1 ). Under the conditions of the same order, A-stability, operation count (at each iterative step) and storage space are the same as the methods in , the methods in the paper improve the stability in a neighborhood at the infinite point. And, by using the OOPI method , it possesses much faster rate of convergence for solving systems of nonlinear equations produced by the DIAOB r,r+1 .
基金supported by National Natural Science Foundation of China(Grant No.42377149)the Research Grants Council of Hong Kong(General Research Fund Project No.17202423).
文摘In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this scheme seamlessly combines the material point method(MPM)for accurately tracking solid-phase deformation and heat transport,and the Eulerian finite element method(FEM)for effectively capturing fluid flow and heat advection-diffusion behavior.The proposed approach circumvents the substantial challenges posed by large nonlinear equation systems with the monolithic solution scheme.The staggered solution process strategically separates each physical field through explicit or implicit integration.The characteristic-based method is used to stabilize advection-dominated heat flows for efficient numerical implementation.Furthermore,a fractional step approach is employed to decompose fluid velocity and pressure,thereby suppressing pore pressure oscillation on the linear background grid.The fracturing initiation and propagation are simulated by a rate-dependent phase field model.Through a series of quasi-static and transient simulations,the exceptional performance and promising potential of the proposed model in addressing THM fracturing problems in poro-elastic media is demonstrated.
基金the National Natural Science Foundation of China （No. 50275059）.
文摘A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a background mesh. It is shown that, in a point collocation approach, the remesh problem because of the mesh distortion in FEM （finite element method） and the low efficiency in Galerkin-based meshfree method are avoided. The corrected kernel functions are introduced to the stabilization of free-surface boundary conditions. The solution of symmetric ring compression problem is compared with a conventional finite element solution, and reasonable results have been obtained.
基金support from the National Natural Science Foundation of China(Grant Nos.52174123&52274222).
文摘This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-oped using the material point method.To reduce the computational cost of Monte Carlo simulations,response surface models are created as surrogate models for the material point system to approximate its dynamic behavior.An adaptive randomized greedy algorithm is employed to construct a sparse polynomial chaos expansion model with a fixed order,effectively balancing the accuracy and computational efficiency of the surrogate model.Based on the sparse polynomial chaos expansion,sensitivity analysis is conducted using the global finite difference and Sobol methods.Several examples of structural dynamics are provided to demonstrate the effectiveness of the proposed method in addressing structural dynamics problems.
基金co-supported by the National Nature Science Foundation of China(No.12303071)the Shanghai Science and Technology Plan Project,China(No.23YF1455500)+1 种基金the China Postdoctoral Science Foundation(No.2023M743653)Ministry of Industry and Information Technology of China through the High Precision Timing Service Project(No.TC220A04A-80)。
文摘The mathematical method of ZTD(zenith tropospheric delay)spatial prediction is important for precise ZTD derivation and real-time precise point positioning(PPP)augmentation.This paper analyses the performance of the popular optimal function coefficient(OFC),sphere cap harmonic analysis(SCHA),kriging and inverse distance weighting(IDW)interpolation in ZTD spatial prediction and Beidou satellite navigation system(BDS)-PPP augmentation over China.For ZTD spatial prediction,the average time consumption of the OFC,kriging,and IDW methods is less than 0.1 s,which is significantly better than that of the SCHA method(63.157 s).The overall ZTD precision of the OFC is 3.44 cm,which outperforms those of the SCHA(9.65 cm),Kriging(10.6 cm),and IDW(11.8 cm)methods.We confirmed that the low performance of kriging and IDW is caused by their weakness in modelling ZTD variation in the vertical direction.To mitigate such deficiencies,an elevation normalization factor(ENF)is introduced into the kriging and IDW models(kriging-ENF and IDW-ENF).The overall ZTD spatial prediction accuracies of IDW-ENF and kriging-ENF are 2.80 cm and 2.01 cm,respectively,which are both superior to those of the OFC and the widely used empirical model GPT3(4.92 cm).For BDS-PPP enhancement,the ZTD provided by the kriging-ENF,IDW-ENF and OFC as prior constraints can effectively reduce the convergence time.Compared with unconstrained BDS-PPP,our proposed kriging-ENF outperforms IDW-ENF and OFC by reducing the horizontal and vertical convergence times by approximately 13.2%and 5.8%in Ningxia and 30.4%and 7.84%in Guangdong,respectively.These results indicate that kriging-ENF is a promising method for ZTD spatial prediction and BDS-PPP enhancement over China.
