The dynamic parameters for internal solitary waves(ISWs)derived from the extended Korteweg-de Vries(eKdV)equation play an important role in the understanding and prediction of ISWs.The spatiotemporal variations of the...The dynamic parameters for internal solitary waves(ISWs)derived from the extended Korteweg-de Vries(eKdV)equation play an important role in the understanding and prediction of ISWs.The spatiotemporal variations of the dynamic parameters of the ISWs in the northern South China Sea(SCS)were studied based on the reanalysis of long-term temperature and salinity datasets.The results for spectrum analysis show that there are definite geographical differences for the periodic variation of the parameters:in shallow water,all parameters vary with a wave period of one year,while in deep water wave components of the parameters at other frequencies exist.Using wavelet analysis,the wavelet power spectral densities in deep water exhibited an inter-annual variation pattern.For example,the wave component of the dispersion coefficient with a wave period of about half a year reached its power peak once every two years.Based on previous work,this inter-annual variation pattern was deduced to be caused by dynamic processes.In further work on the regulatory mechanisms,empirical orthogonal function(EOF)decomposition was performed.It was found that the modes of the dispersion coefficient have different geographical distributions,explaining the reason why the wave components in different frequencies appeared in different locations.The numerical simulation results confirm that the variations in the parameters of the ISWs derived from the eKdV equation could affect the waveforms significantly because of changes in the polarity of the ISWs.Therefore,the periodic variations of the dynamic parameters are related to the geographical location because of dynamic processes operating.展开更多
This paper focuses on studying the integration method of a generalized Birkhoffian system.The method of variation on parameters for the dynamical equations of a generalized Birkhoffian system is presented.The procedur...This paper focuses on studying the integration method of a generalized Birkhoffian system.The method of variation on parameters for the dynamical equations of a generalized Birkhoffian system is presented.The procedure for solving the problem can be divided into two steps:the first step,a system of auxiliary equations is constructed and its general solution is given;the second step,the parameters are varied,and the solution of the problem is obtained by using the properties of generalized canonical transformation.The method of variation on parameters for the generalized Birkhoffian system is of universal significance,and we take a nonholonomic system and a nonconservative system as examples to illustrate the application of the results of this paper.展开更多
The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(R...The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.展开更多
As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accura...As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.展开更多
The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-...The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-fidelity models has been challenging due to the prohibitive computational costs.This paper presents an efficient parallel algorithm tailored for HVNM based on the Message Passing Interface standard.The algorithm evenly distributes the response matrix sets among processors during the matrix formation process,thus enabling independent construction without communication.Once the formation tasks are completed,a collective operation merges and shares the matrix sets among the processors.For the solution process,the problem domain is decomposed into subdomains assigned to specific processors,and the red-black Gauss-Seidel iteration is employed within each subdomain to solve the response matrix equation.Point-to-point communication is conducted between adjacent subdomains to exchange data along the boundaries.The accuracy and efficiency of the parallel algorithm are verified using the KAIST and JRR-3 test cases.Numerical results obtained with multiple processors agree well with those obtained from Monte Carlo calculations.The parallelization of HVNM results in eigenvalue errors of 31 pcm/-90 pcm and fission rate RMS errors of 1.22%/0.66%,respectively,for the 3D KAIST problem and the 3D JRR-3 problem.In addition,the parallel algorithm significantly reduces computation time,with an efficiency of 68.51% using 36 processors in the KAIST problem and 77.14% using 144 processors in the JRR-3 problem.展开更多
For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for s...For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for solving the dynamical equations of the constrained Birkhoffian system is provided.First the differential equations of motion for the constrained Birkhoffian system as well as for the corresponding free Birkhoffian system are established.Secondly a system of auxiliary equations is constructed and the general solution of the equations is found.Finally by varying the parameters and utilizing the properties of the generalized canonical transformation of the Birkhoffian system the solution of the problem can be obtained.The proposed method reveals the inherent relationship between the solution of a free Birkhoffian system and that of a constrained Birkhoffian system. The research results are of universal significance which can be further used in a variety of constrained mechanical systems such as non-conservative systems and nonholonomic systems etc.展开更多
Synthetic aperture radar(SAR)and wave spectrometers,crucial in microwave remote sensing,play an essential role in monitoring sea surface wind and wave conditions.However,they face inherent limitations in observing sea...Synthetic aperture radar(SAR)and wave spectrometers,crucial in microwave remote sensing,play an essential role in monitoring sea surface wind and wave conditions.However,they face inherent limitations in observing sea surface phenomena.SAR systems,for instance,are hindered by an azimuth cut-off phenomenon in sea surface wind field observation.Wave spectrometers,while unaffected by the azimuth cutoff phenomenon,struggle with low azimuth resolution,impacting the capture of detailed wave and wind field data.This study utilizes SAR and surface wave investigation and monitoring(SWIM)data to initially extract key feature parameters,which are then prioritized using the extreme gradient boosting(XGBoost)algorithm.The research further addresses feature collinearity through a combined analysis of feature importance and correlation,leading to the development of an inversion model for wave and wind parameters based on XGBoost.A comparative analysis of this model with ERA5 reanalysis and buoy data for of significant wave height,mean wave period,wind direction,and wind speed reveals root mean square errors of 0.212 m,0.525 s,27.446°,and 1.092 m/s,compared to 0.314 m,0.888 s,27.698°,and 1.315 m/s from buoy data,respectively.These results demonstrate the model’s effective retrieval of wave and wind parameters.Finally,the model,incorporating altimeter and scatterometer data,is evaluated against SAR/SWIM single and dual payload inversion methods across different wind speeds.This comparison highlights the model’s superior inversion accuracy over other methods.展开更多
A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO mode...A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.展开更多
The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inv...The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace's equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace's equation. Then, the 'L-Curve' principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.展开更多
The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized funct...The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized functional was established, and the functional was solved by the sensitivity coefficient and Newtonaphson iteration method. Moreover, the orthogonal experimental design was used to estimate the appropriate initial value and variation domain of each variable to decrease the number of iteration and improve the identification accuracy and efficiency. It illustrated a detailed case of AlSiTMg sand mold casting and the temperature measurement experiment was done. The physical properties of sand mold and the interracial heat transfer coefficient were identified at the meantime. The results indicated that the new regularization method was efficient in overcoming the ill-posedness of the inverse heat conduction problem and improving the stability and accuracy of the solutions.展开更多
The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows ...The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows through an aquifer,its composition and temperature may variation dependent on the aquifer condition through which it flows.Thus,hydrologic investigations can also provide useful information about the subsurface geology of a region.But because such studies investigate processes that follow under the Earth's shallow,obtaining the information necessary to answer these questions is not continuously easy.Springs,which discharge groundwater table directly,afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity(T)and storativity(S)are vital for the evaluation of groundwater resources.There are several methods to estimate the accurate aquifer parameters(i.e.hydrograph analysis,pumping test,etc.).In initial days,these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory.The simultaneous information on the hydraulic behavior of the well(borehole)that provides on this method,the reservoir and the reservoir boundaries,are important for efficient aquifer and well data management and analysis.The most common in-situ test is pumping test performed on wells,which involves the measurement of the fall and increase of groundwater level with respect to time.The alteration in groundwater level(drawdown/recovery)is caused due to pumping of water from the well.Theis(1935)was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer.It is essential to know the transmissivity(T=Kb,where b is the aquifer thickness;pumping flow rate,Q=TW(dh/dl)flow through an aquifer)and storativity(confined aquifer:S=bS_s,unconfined:S=S_y),for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer.The determination of aquifer's parameters is an important basis for groundwater resources evaluation,numerical simulation,development and protection as well as scientific management.For determining aquifer's parameters,pumping test is a main method.A case study shows that these techniques have been fast speed and high correctness.The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.展开更多
Based on the nonlinear characiers of the discrete problems of some ellipticalvariational inequalities, this paper presents a numerical iterative method, the schemesof which are pithy and converge rapidly The new metho...Based on the nonlinear characiers of the discrete problems of some ellipticalvariational inequalities, this paper presents a numerical iterative method, the schemesof which are pithy and converge rapidly The new method possesses a high efficiency. insolving such applied engineering problems as obstacle problems and .free boundary.problems arising in fluid lubrications.展开更多
In this paper,we investigate pseudomonotone and Lipschitz continuous variational inequalities in real Hilbert spaces.For solving this problem,we propose a new method that combines the advantages of the subgradient ext...In this paper,we investigate pseudomonotone and Lipschitz continuous variational inequalities in real Hilbert spaces.For solving this problem,we propose a new method that combines the advantages of the subgradient extragradient method and the projection contraction method.Some very recent papers have considered different inertial algorithms which allowed the inertial factor is chosen in[0;1].The purpose of this work is to continue working in this direction,we propose another inertial subgradient extragradient method that the inertial factor can be chosen in a special case to be 1.Under suitable mild conditions,we establish the weak convergence of the proposed algorithm.Moreover,linear convergence is obtained under strong pseudomonotonicity and Lipschitz continuity assumptions.Finally,some numerical illustrations are given to confirm the theoretical analysis.展开更多
There is substantial individual variation in the growth rates of sea cucumber Apostiehopus japonicus individuals. This necessitates additional work to grade the seed stock and lengthens the production period. We evalu...There is substantial individual variation in the growth rates of sea cucumber Apostiehopus japonicus individuals. This necessitates additional work to grade the seed stock and lengthens the production period. We evaluated the influence of three culture methods (free-mixed, isolated-mixed, isolated-alone) on individual variation in growth and assessed the relationship between feeding, energy conversion efficiency, and individual growth variation in individually cultured sea cucumbers. Of the different culture methods, animals grew best when reared in the isolated-mixed treatment (i.e., size classes were held separately), though there was no difference in individual variation in growth between rearing treatment groups. The individual variation in growth was primarily attributed to genetic factors. The difference in food conversion efficiency caused by genetic differences among individuals was thought to be the origin of the variance. The level of individual growth variation may be altered by interactions among individuals and environmental heterogeneity. Our results suggest that, in addition to traditional seed grading, design of a new kind of substrate that changes the spatial distribution of sea cucumbers would effectively enhance growth and reduce individual variation in growth of sea cucumbers in culture.展开更多
Lane-Emden type equation is a nonlinear differential equation appears in many fields such as stellar structure, radioactive cooling and modeling of clusters of galaxies. In this work, this equation is investigated usi...Lane-Emden type equation is a nonlinear differential equation appears in many fields such as stellar structure, radioactive cooling and modeling of clusters of galaxies. In this work, this equation is investigated using a semi-analytical method called the Variation of parameters method with an auxiliary parameter. In the applied technique, an unknown auxiliary parameter is inserted in Variation of Parameters Method to solve some special cases of these equations. The used algorithm is easy to implement and very effective. The obtained solutions are also fairly accurate.展开更多
This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODE...This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODEs).A unification of the concepts underlying these three methods is attempted by considering a very general iterative algorithm for VIM.It is found that all the three methods can be regarded as special cases of using a very general matrix of Lagrange multipliers in the iterative algorithm of VIM.The global variational iteration method is briefly reviewed,and further recast into a Local VIM,which is much more convenient and capable of predicting long term complex dynamic responses of nonlinear systems even if they are chaotic.展开更多
In remote sensing sea surface temperature (SST), the traditional fusion method is used to compute the dot product of a subjective weight vector with a satellite measurement vector, while the result requires validati...In remote sensing sea surface temperature (SST), the traditional fusion method is used to compute the dot product of a subjective weight vector with a satellite measurement vector, while the result requires validation by field measurement. However, field measurement that relative to the satellite measurement is very sparse, many information may not be verified. A relative objective weight vector is constructed by using the limited field measurement, which is based on coefficient of variation method. And then it make an application of the data fusion by the weighted average method in the SST data. fuse SST data with the weighted average method. In this way, some posteriori information can be added to the fusion process. The model reduces the dependence on verification, and some of the satellite measurement can be handled without corresponding to the field measurement, and the fusion result matches transfer errors theory.展开更多
The rotation of the Earth's liquid core creates the Nearly Diurnal Free Wobble (NDFW). It is one of the problems of researching the Earth's liquid core for us to retrieve the parameters of the Earth's Free...The rotation of the Earth's liquid core creates the Nearly Diurnal Free Wobble (NDFW). It is one of the problems of researching the Earth's liquid core for us to retrieve the parameters of the Earth's Free Core Nutation (FCN),from the Earth's liquid core resonance of gravity tide waves on the diurnal frequency band. Since 1987, some scientists in many different countries have come to calculate the parameters of FCN by using the observational data of gravity tide waves on the diumal frequency band. They basically followed the Stacking method, which needed five diurnal waves for the resolution. In this paper, authors introduced the Aky-Btk value method with clearly geometrical meaning as a new method, which only requires three very high signal-noise-ratio waves o1, K1 and P1 for the resolution. Authors chose the observational data of the three superconducting gravimeter stations respectively located in Cantley of Canada, Wuhan of China and Brussels of Belgium, to compute the parameters of FCN. It was the first time that the observational results of the parameters of FCN obtained from gravitational tide were in accord with the parameters of FCN gained from VLBI since 1987.展开更多
The physical-cover-oriented variational principle of nonlinear numerical manifold method (NNMM) for the analysis of plastical problems is put forward according to the displacement model and the characters of numerical...The physical-cover-oriented variational principle of nonlinear numerical manifold method (NNMM) for the analysis of plastical problems is put forward according to the displacement model and the characters of numerical manifold method (NMM). The theoretical calculating formulations and the controlling equation of NNMM are derived. As an example, the plate with a hole in the center is calculated and the results show that the solution precision and efficiency of NNMM are agreeable.展开更多
Among the methods used for evaluating the potential hydraulic erodibility of rock,the most common methods are those based on the correlation between the force of flowing water and the capacity of a rock to resist eros...Among the methods used for evaluating the potential hydraulic erodibility of rock,the most common methods are those based on the correlation between the force of flowing water and the capacity of a rock to resist erosion,such as Annandale’s and Pells’methods.The capacity of a rock to resist erosion is evaluated based on erodibility indices that are determined from specific geomechanical parameters of a rock mass.These indices include unconfined compressive strength(UCS)of rock,rock block size,joint shear strength,a block’s shape and orientation relative to the direction of flow,joint openings,and the nature of the surface to be potentially eroded.However,it is difficult to determine the relevant geomechanical parameters for evaluating the hydraulic erodibility of rock.The assessment of eroded unlined spillways of dams has shown that the capacity of a rock to resist erosion is not accurately evaluated.Using more than 100 case studies,we develop a method to determine the relevant geomechanical parameters for evaluating the hydraulic erodibility of rock in unlined spillways.The UCS of rock is found not to be a relevant parameter for evaluating the hydraulic erodibility of rock.On the other hand,we find that the use of three-dimensional(3D)block volume measurements,instead of the block size factor used in Annandale’s method,improves the rock block size estimation.Furthermore,the parameter representing the effect of a rock block’s shape and orientation relative to the direction of flow,as considered in Pells’method,is more accurate than the parameter adopted by Annandale’s method.展开更多
基金Supported by the Hunan Provincial Science Fund for Distinguished Young Scholars(No.2023JJ10053)the National Natural Science Foundation of China(No.42276205)。
文摘The dynamic parameters for internal solitary waves(ISWs)derived from the extended Korteweg-de Vries(eKdV)equation play an important role in the understanding and prediction of ISWs.The spatiotemporal variations of the dynamic parameters of the ISWs in the northern South China Sea(SCS)were studied based on the reanalysis of long-term temperature and salinity datasets.The results for spectrum analysis show that there are definite geographical differences for the periodic variation of the parameters:in shallow water,all parameters vary with a wave period of one year,while in deep water wave components of the parameters at other frequencies exist.Using wavelet analysis,the wavelet power spectral densities in deep water exhibited an inter-annual variation pattern.For example,the wave component of the dispersion coefficient with a wave period of about half a year reached its power peak once every two years.Based on previous work,this inter-annual variation pattern was deduced to be caused by dynamic processes.In further work on the regulatory mechanisms,empirical orthogonal function(EOF)decomposition was performed.It was found that the modes of the dispersion coefficient have different geographical distributions,explaining the reason why the wave components in different frequencies appeared in different locations.The numerical simulation results confirm that the variations in the parameters of the ISWs derived from the eKdV equation could affect the waveforms significantly because of changes in the polarity of the ISWs.Therefore,the periodic variations of the dynamic parameters are related to the geographical location because of dynamic processes operating.
基金supported by the National Natural Science Foundation of China (10972151)
文摘This paper focuses on studying the integration method of a generalized Birkhoffian system.The method of variation on parameters for the dynamical equations of a generalized Birkhoffian system is presented.The procedure for solving the problem can be divided into two steps:the first step,a system of auxiliary equations is constructed and its general solution is given;the second step,the parameters are varied,and the solution of the problem is obtained by using the properties of generalized canonical transformation.The method of variation on parameters for the generalized Birkhoffian system is of universal significance,and we take a nonholonomic system and a nonconservative system as examples to illustrate the application of the results of this paper.
基金supported by the National Key Research and Development Program of China(No.2023YFC2907600)the National Natural Science Foundation of China(Nos.42077267,42277174 and 52074164)+2 种基金the Natural Science Foundation of Shandong Province,China(No.ZR2020JQ23)the Opening Project of State Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology(No.KFJJ21-02Z)the Fundamental Research Funds for the Central Universities,China(No.2022JCCXSB03).
文摘The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.
基金supported by the National Natural Science Foundation of China(Grant No.42277165)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant No.CUGCJ1821)the National Overseas Study Fund(Grant No.202106410040).
文摘As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.
基金supported by the National Key Research and Development Program of China(No.2020YFB1901900)the National Natural Science Foundation of China(Nos.U20B2011,12175138)the Shanghai Rising-Star Program。
文摘The heterogeneous variational nodal method(HVNM)has emerged as a potential approach for solving high-fidelity neutron transport problems.However,achieving accurate results with HVNM in large-scale problems using high-fidelity models has been challenging due to the prohibitive computational costs.This paper presents an efficient parallel algorithm tailored for HVNM based on the Message Passing Interface standard.The algorithm evenly distributes the response matrix sets among processors during the matrix formation process,thus enabling independent construction without communication.Once the formation tasks are completed,a collective operation merges and shares the matrix sets among the processors.For the solution process,the problem domain is decomposed into subdomains assigned to specific processors,and the red-black Gauss-Seidel iteration is employed within each subdomain to solve the response matrix equation.Point-to-point communication is conducted between adjacent subdomains to exchange data along the boundaries.The accuracy and efficiency of the parallel algorithm are verified using the KAIST and JRR-3 test cases.Numerical results obtained with multiple processors agree well with those obtained from Monte Carlo calculations.The parallelization of HVNM results in eigenvalue errors of 31 pcm/-90 pcm and fission rate RMS errors of 1.22%/0.66%,respectively,for the 3D KAIST problem and the 3D JRR-3 problem.In addition,the parallel algorithm significantly reduces computation time,with an efficiency of 68.51% using 36 processors in the KAIST problem and 77.14% using 144 processors in the JRR-3 problem.
基金The National Natural Science Foundation of China(No.10972151,11272227)
文摘For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for solving the dynamical equations of the constrained Birkhoffian system is provided.First the differential equations of motion for the constrained Birkhoffian system as well as for the corresponding free Birkhoffian system are established.Secondly a system of auxiliary equations is constructed and the general solution of the equations is found.Finally by varying the parameters and utilizing the properties of the generalized canonical transformation of the Birkhoffian system the solution of the problem can be obtained.The proposed method reveals the inherent relationship between the solution of a free Birkhoffian system and that of a constrained Birkhoffian system. The research results are of universal significance which can be further used in a variety of constrained mechanical systems such as non-conservative systems and nonholonomic systems etc.
基金The project supported by Key Laboratory of Space Ocean Remote Sensing and Application,Ministry of Natural Resources under contract No.2023CFO016the National Natural Science Foundation of China under contract No.61931025+1 种基金the Innovation Fund Project for Graduate Student of China University of Petroleum(East China)the Fundamental Research Funds for the Central Universities under contract No.23CX04042A.
文摘Synthetic aperture radar(SAR)and wave spectrometers,crucial in microwave remote sensing,play an essential role in monitoring sea surface wind and wave conditions.However,they face inherent limitations in observing sea surface phenomena.SAR systems,for instance,are hindered by an azimuth cut-off phenomenon in sea surface wind field observation.Wave spectrometers,while unaffected by the azimuth cutoff phenomenon,struggle with low azimuth resolution,impacting the capture of detailed wave and wind field data.This study utilizes SAR and surface wave investigation and monitoring(SWIM)data to initially extract key feature parameters,which are then prioritized using the extreme gradient boosting(XGBoost)algorithm.The research further addresses feature collinearity through a combined analysis of feature importance and correlation,leading to the development of an inversion model for wave and wind parameters based on XGBoost.A comparative analysis of this model with ERA5 reanalysis and buoy data for of significant wave height,mean wave period,wind direction,and wind speed reveals root mean square errors of 0.212 m,0.525 s,27.446°,and 1.092 m/s,compared to 0.314 m,0.888 s,27.698°,and 1.315 m/s from buoy data,respectively.These results demonstrate the model’s effective retrieval of wave and wind parameters.Finally,the model,incorporating altimeter and scatterometer data,is evaluated against SAR/SWIM single and dual payload inversion methods across different wind speeds.This comparison highlights the model’s superior inversion accuracy over other methods.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 90111011 and 10471039), the National Key Basic Research Special Foundation of China (Grant Nos 2003CB415101-03 and 2004CB418304), the Key Basic Research Foundation of the Chinese Academy of Sciences (Grant No KZCX3-SW-221) and in part by E-Institutes of Shanghai Municipal Education Commission (Grant No N.E03004).
文摘A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.
基金Project supported by the National Natural Science Foundation of China(Grant No.41175025)
文摘The simplified linear model of Grad-Shafranov (GS) reconstruction can be reformulated into an inverse boundary value problem of Laplace's equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace's equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace's equation. Then, the 'L-Curve' principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.
文摘The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized functional was established, and the functional was solved by the sensitivity coefficient and Newtonaphson iteration method. Moreover, the orthogonal experimental design was used to estimate the appropriate initial value and variation domain of each variable to decrease the number of iteration and improve the identification accuracy and efficiency. It illustrated a detailed case of AlSiTMg sand mold casting and the temperature measurement experiment was done. The physical properties of sand mold and the interracial heat transfer coefficient were identified at the meantime. The results indicated that the new regularization method was efficient in overcoming the ill-posedness of the inverse heat conduction problem and improving the stability and accuracy of the solutions.
文摘The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows through an aquifer,its composition and temperature may variation dependent on the aquifer condition through which it flows.Thus,hydrologic investigations can also provide useful information about the subsurface geology of a region.But because such studies investigate processes that follow under the Earth's shallow,obtaining the information necessary to answer these questions is not continuously easy.Springs,which discharge groundwater table directly,afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity(T)and storativity(S)are vital for the evaluation of groundwater resources.There are several methods to estimate the accurate aquifer parameters(i.e.hydrograph analysis,pumping test,etc.).In initial days,these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory.The simultaneous information on the hydraulic behavior of the well(borehole)that provides on this method,the reservoir and the reservoir boundaries,are important for efficient aquifer and well data management and analysis.The most common in-situ test is pumping test performed on wells,which involves the measurement of the fall and increase of groundwater level with respect to time.The alteration in groundwater level(drawdown/recovery)is caused due to pumping of water from the well.Theis(1935)was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer.It is essential to know the transmissivity(T=Kb,where b is the aquifer thickness;pumping flow rate,Q=TW(dh/dl)flow through an aquifer)and storativity(confined aquifer:S=bS_s,unconfined:S=S_y),for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer.The determination of aquifer's parameters is an important basis for groundwater resources evaluation,numerical simulation,development and protection as well as scientific management.For determining aquifer's parameters,pumping test is a main method.A case study shows that these techniques have been fast speed and high correctness.The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.
文摘Based on the nonlinear characiers of the discrete problems of some ellipticalvariational inequalities, this paper presents a numerical iterative method, the schemesof which are pithy and converge rapidly The new method possesses a high efficiency. insolving such applied engineering problems as obstacle problems and .free boundary.problems arising in fluid lubrications.
基金funded by the University of Science,Vietnam National University,Hanoi under project number TN.21.01。
文摘In this paper,we investigate pseudomonotone and Lipschitz continuous variational inequalities in real Hilbert spaces.For solving this problem,we propose a new method that combines the advantages of the subgradient extragradient method and the projection contraction method.Some very recent papers have considered different inertial algorithms which allowed the inertial factor is chosen in[0;1].The purpose of this work is to continue working in this direction,we propose another inertial subgradient extragradient method that the inertial factor can be chosen in a special case to be 1.Under suitable mild conditions,we establish the weak convergence of the proposed algorithm.Moreover,linear convergence is obtained under strong pseudomonotonicity and Lipschitz continuity assumptions.Finally,some numerical illustrations are given to confirm the theoretical analysis.
基金Supported by the National Natural Science Foundation of China(No.41106134)the National Marine Public Welfare Research Project of China(No.201305043)+1 种基金the National High Technology Research and Development Program of China(863 Program)(No.2012AA10A412)the Agriculture Science Technology Achievement Transformation Fund(No.2012GB24910656)
文摘There is substantial individual variation in the growth rates of sea cucumber Apostiehopus japonicus individuals. This necessitates additional work to grade the seed stock and lengthens the production period. We evaluated the influence of three culture methods (free-mixed, isolated-mixed, isolated-alone) on individual variation in growth and assessed the relationship between feeding, energy conversion efficiency, and individual growth variation in individually cultured sea cucumbers. Of the different culture methods, animals grew best when reared in the isolated-mixed treatment (i.e., size classes were held separately), though there was no difference in individual variation in growth between rearing treatment groups. The individual variation in growth was primarily attributed to genetic factors. The difference in food conversion efficiency caused by genetic differences among individuals was thought to be the origin of the variance. The level of individual growth variation may be altered by interactions among individuals and environmental heterogeneity. Our results suggest that, in addition to traditional seed grading, design of a new kind of substrate that changes the spatial distribution of sea cucumbers would effectively enhance growth and reduce individual variation in growth of sea cucumbers in culture.
文摘Lane-Emden type equation is a nonlinear differential equation appears in many fields such as stellar structure, radioactive cooling and modeling of clusters of galaxies. In this work, this equation is investigated using a semi-analytical method called the Variation of parameters method with an auxiliary parameter. In the applied technique, an unknown auxiliary parameter is inserted in Variation of Parameters Method to solve some special cases of these equations. The used algorithm is easy to implement and very effective. The obtained solutions are also fairly accurate.
文摘This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODEs).A unification of the concepts underlying these three methods is attempted by considering a very general iterative algorithm for VIM.It is found that all the three methods can be regarded as special cases of using a very general matrix of Lagrange multipliers in the iterative algorithm of VIM.The global variational iteration method is briefly reviewed,and further recast into a Local VIM,which is much more convenient and capable of predicting long term complex dynamic responses of nonlinear systems even if they are chaotic.
基金Project supported by the National Natural Science Foundation of China(Grant No.40976108)the Shanghai Leading Academic Discipline Project(Grant No.J50103)
文摘In remote sensing sea surface temperature (SST), the traditional fusion method is used to compute the dot product of a subjective weight vector with a satellite measurement vector, while the result requires validation by field measurement. However, field measurement that relative to the satellite measurement is very sparse, many information may not be verified. A relative objective weight vector is constructed by using the limited field measurement, which is based on coefficient of variation method. And then it make an application of the data fusion by the weighted average method in the SST data. fuse SST data with the weighted average method. In this way, some posteriori information can be added to the fusion process. The model reduces the dependence on verification, and some of the satellite measurement can be handled without corresponding to the field measurement, and the fusion result matches transfer errors theory.
文摘The rotation of the Earth's liquid core creates the Nearly Diurnal Free Wobble (NDFW). It is one of the problems of researching the Earth's liquid core for us to retrieve the parameters of the Earth's Free Core Nutation (FCN),from the Earth's liquid core resonance of gravity tide waves on the diurnal frequency band. Since 1987, some scientists in many different countries have come to calculate the parameters of FCN by using the observational data of gravity tide waves on the diumal frequency band. They basically followed the Stacking method, which needed five diurnal waves for the resolution. In this paper, authors introduced the Aky-Btk value method with clearly geometrical meaning as a new method, which only requires three very high signal-noise-ratio waves o1, K1 and P1 for the resolution. Authors chose the observational data of the three superconducting gravimeter stations respectively located in Cantley of Canada, Wuhan of China and Brussels of Belgium, to compute the parameters of FCN. It was the first time that the observational results of the parameters of FCN obtained from gravitational tide were in accord with the parameters of FCN gained from VLBI since 1987.
文摘The physical-cover-oriented variational principle of nonlinear numerical manifold method (NNMM) for the analysis of plastical problems is put forward according to the displacement model and the characters of numerical manifold method (NMM). The theoretical calculating formulations and the controlling equation of NNMM are derived. As an example, the plate with a hole in the center is calculated and the results show that the solution precision and efficiency of NNMM are agreeable.
基金the Natural Sciences and Engineering Research Council of Canada(Grant No.498020-16)Hydro-Quebec(NC525700)the Mitacs Accelerate Program(Grant Ref.IT10008)
文摘Among the methods used for evaluating the potential hydraulic erodibility of rock,the most common methods are those based on the correlation between the force of flowing water and the capacity of a rock to resist erosion,such as Annandale’s and Pells’methods.The capacity of a rock to resist erosion is evaluated based on erodibility indices that are determined from specific geomechanical parameters of a rock mass.These indices include unconfined compressive strength(UCS)of rock,rock block size,joint shear strength,a block’s shape and orientation relative to the direction of flow,joint openings,and the nature of the surface to be potentially eroded.However,it is difficult to determine the relevant geomechanical parameters for evaluating the hydraulic erodibility of rock.The assessment of eroded unlined spillways of dams has shown that the capacity of a rock to resist erosion is not accurately evaluated.Using more than 100 case studies,we develop a method to determine the relevant geomechanical parameters for evaluating the hydraulic erodibility of rock in unlined spillways.The UCS of rock is found not to be a relevant parameter for evaluating the hydraulic erodibility of rock.On the other hand,we find that the use of three-dimensional(3D)block volume measurements,instead of the block size factor used in Annandale’s method,improves the rock block size estimation.Furthermore,the parameter representing the effect of a rock block’s shape and orientation relative to the direction of flow,as considered in Pells’method,is more accurate than the parameter adopted by Annandale’s method.