To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical c...To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.展开更多
This study investigates the complex heat transfer dynamics inmultilayer bifacial photovoltaic(bPV)solar modules under spectrally resolved solar irradiation.A novel numericalmodel is developed to incorporate internal h...This study investigates the complex heat transfer dynamics inmultilayer bifacial photovoltaic(bPV)solar modules under spectrally resolved solar irradiation.A novel numericalmodel is developed to incorporate internal heat generation resulting from optical absorption,grounded in the physical equations governing light-matter interactions within the module’smultilayer structure.The model accounts for reflection and transmission at each interface between adjacent layers,as well as absorption within individual layers,using the wavelength-dependent dielectric properties of constituent materials.These properties are used to calculate the spectral reflectance,transmittance,and absorption coefficients,enabling precise quantification of internal heat sources from irradiance incidents on both the front and rear surfaces of the module.The study further examines the influence of irradiance reflection on thermal behavior,evaluates the thermal impact of various supporting materials placed beneath the module,and analyzes the role of albedo in modifying heat distribution.By incorporating spectrally resolved heat generation across each layer often simplified or omitted in conventional models,the proposed approach enhances physical accuracy.The transient heat equation is solved using a one-dimensional finite difference(FD)method to produce detailed temperature profiles under multiple operating scenarios,including Standard Test Conditions(STC),Bifacial Standard Test Conditions(BSTC),Normal Operating Cell Temperature(NOCT),and Bifacial NOCT(BNOCT).The results offer valuable insights into the interplay between optical and thermal phenomena in bifacial systems,informing the design and optimization of more efficient photovoltaic technologies.展开更多
In response to the issue of fuzzy matching and association when optical observation data are matched with the orbital elements in a catalog database,this paper proposes a matching and association strategy based on the...In response to the issue of fuzzy matching and association when optical observation data are matched with the orbital elements in a catalog database,this paper proposes a matching and association strategy based on the arcsegment difference method.First,a matching error threshold is set to match the observation data with the known catalog database.Second,the matching results for the same day are sorted on the basis of target identity and observation residuals.Different matching error thresholds and arc-segment dynamic association thresholds are then applied to categorize the observation residuals of the same target across different arc-segments,yielding matching results under various thresholds.Finally,the orbital residual is computed through orbit determination(OD),and the positional error is derived by comparing the OD results with the orbit track from the catalog database.The appropriate matching error threshold is then selected on the basis of these results,leading to the final matching and association of the fuzzy correlation data.Experimental results showed that the correct matching rate for data arc-segments is 92.34% when the matching error threshold is set to 720″,with the arc-segment difference method processing the results of an average matching rate of 97.62% within 8 days.The remaining 5.28% of the fuzzy correlation data are correctly matched and associated,enabling identification of orbital maneuver targets through further processing and analysis.This method substantially enhances the efficiency and accuracy of space target cataloging,offering robust technical support for dynamic maintenance of the space target database.展开更多
In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving th...In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.展开更多
A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficien...A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficients(k)are constructed based on FDM.The rock bridge area was divided through k-means algorithm and the optimal number of clusters was determined by sum of squared errors(SSE)and elbow method.The influence of maximum principal stress and stress change rate as clustering indexes on the clustering results of rock bridges was compared by using Euclidean distance.The results show that using stress change rate as clustering index is more effective.When the joint coalescence coefficient is less than 0.6,there is no significant stress concentration in the middle area of adjacent joints,that is,no generation of rock bridge.In addition,the range of rock bridge is affected by the coalescence coefficient(k),the relative position of joints and the parameters of weak interlayer.展开更多
For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study prop...For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.展开更多
In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolso...In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.展开更多
In this paper,we give improved error estimates for linearized and nonlinear CrankNicolson type finite difference schemes of Ginzburg-Landau equation in two dimensions.For linearized Crank-Nicolson scheme,we use mathem...In this paper,we give improved error estimates for linearized and nonlinear CrankNicolson type finite difference schemes of Ginzburg-Landau equation in two dimensions.For linearized Crank-Nicolson scheme,we use mathematical induction to get unconditional error estimates in discrete L^(2)and H^(1)norm.However,it is not applicable for the nonlinear scheme.Thus,based on a‘cut-off’function and energy analysis method,we get unconditional L^(2)and H^(1)error estimates for the nonlinear scheme,as well as boundedness of numerical solutions.In addition,if the assumption for exact solutions is improved compared to before,unconditional and optimal pointwise error estimates can be obtained by energy analysis method and several Sobolev inequalities.Finally,some numerical examples are given to verify our theoretical analysis.展开更多
Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is a...Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is analyzed. The dynamic model uses Euler-Bernoulli beam theory for cantilever beam, in which the driving force (or driving moment) of the electromagnetic system is taken into account, and the contact force between moving contact and stationary contact is simulated by the Kelvin-Voigt vis-coelastic...展开更多
An improved finite difference method （FDM）is described to solve existing problems such as low efficiency and poor convergence performance in the traditional method adopted to derive the pressure distribution of aero...An improved finite difference method （FDM）is described to solve existing problems such as low efficiency and poor convergence performance in the traditional method adopted to derive the pressure distribution of aerostatic bearings. A detailed theoretical analysis of the pressure distribution of the orifice-compensated aerostatic journal bearing is presented. The nonlinear dimensionless Reynolds equation of the aerostatic journal bearing is solved by the finite difference method. Based on the principle of flow equilibrium, a new iterative algorithm named the variable step size successive approximation method is presented to adjust the pressure at the orifice in the iterative process and enhance the efficiency and convergence performance of the algorithm. A general program is developed to analyze the pressure distribution of the aerostatic journal bearing by Matlab tool. The results show that the improved finite difference method is highly effective, reliable, stable, and convergent. Even when very thin gas film thicknesses （less than 2 Win）are considered, the improved calculation method still yields a result and converges fast.展开更多
Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with ...Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with theoretical method, the finite difference method has been verified to be feasible by a case study. It is found that under seismic loading, the dynamic response of anchorage system is synchronously fluctuated with the seismic vibration. The change of displacement amplitude of material points is slight, and comparatively speaking, the displacement amplitude of the outside point is a little larger than that of the inside point, which shows amplification effect of surface. While the axial force amplitude transforms considerably from the inside to the outside. It increases first and reaches the peak value in the intersection between the anchoring section and free section, then decreases slowly in the free section. When considering damping effect of anchorage system, the finite difference method can reflect the time attenuation characteristic better, and the calculating result would be safer and more reasonable than the dynamic steady-state theoretical method. What is more, the finite difference method can be applied to the dynamic response analysis of harmonic and seismic random vibration for all kinds of anchor, and hence has a broad application prospect.展开更多
An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for ...An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.展开更多
Thermal stress simulation can provide a scientific reference to eliminate defects such as crack,residual stress centralization and deformation etc.,caused by thermal stress during casting solidification.To study the t...Thermal stress simulation can provide a scientific reference to eliminate defects such as crack,residual stress centralization and deformation etc.,caused by thermal stress during casting solidification.To study the thermal stress distribution during casting process,a unilateral thermal-stress coupling model was employed to simulate 3D casting stress using Finite Difference Method(FDM),namely all the traditional thermal-elastic-plastic equations are numerically and differentially discrete.A FDM/FDM numerical simulation system was developed to analyze temperature and stress fields during casting solidification process.Two practical verifications were carried out,and the results from simulation basically coincided with practical cases.The results indicated that the FDM/FDM stress simulation system can be used to simulate the formation of residual stress,and to predict the occurrence of hot tearing.Because heat transfer and stress analysis are all based on FDM,they can use the same FD model,which can avoid the matching process between different models,and hence reduce temperature-load transferring errors.This approach makes the simulation of fluid flow,heat transfer and stress analysis unify into one single model.展开更多
A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is ...A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.展开更多
With the finite difference method to calculate the temperature distribution in mass concrete structures, the solution precision will increase with a smaller step size, at the cost of computational time. In view of the...With the finite difference method to calculate the temperature distribution in mass concrete structures, the solution precision will increase with a smaller step size, at the cost of computational time. In view of the basic characteristics of the finite difference method, a simple yet powerful improvement is introduced. By multiplying the adiabatic temperature function with a correction factor, the precision of the solution can be assured without an increase in the computation time. In addition, the correction rules for three types of commonly used concrete hydration formulas are investigated.展开更多
In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical ...In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical reaction,say zi,is between 0 and 1,withΣz_(i)=1.Due to the lack of maximum-principle,most of the previous bound-preserving technique cannot be applied directly.To preserve those bounds,we will use the positivity-preserving technique to all the zi'is and enforceΣz_(i)=1 by constructing conservative schemes,thanks to conservative time integrations and consistent numerical fluxes in the system.Moreover,detonation is an extreme singular mode of flame propagation in premixed gas,and the model contains a significant stiff source.It is well known that for hyperbolic equations with stiff source,the transition points in the numerical approximations near the shocks may trigger spurious shock speed,leading to wrong shock position.Intuitively,the high-order weighted essentially non-oscillatory(WENO)scheme,which can suppress oscillations near the discontinuities,would be a good choice for spatial discretization.However,with the nonlinear weights,the numerical fluxes are no longer“consistent”,leading to nonconservative numerical schemes and the bound-preserving technique does not work.Numerical experiments demonstrate that,without further numerical techniques such as subcell resolutions,the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.展开更多
A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditiona...A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.展开更多
In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald...In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.展开更多
In this study, we examined the varia- tions between the difference method and the enzyme- hydrolyzed casein method for determining endogenous amino acid loss and the true amino acid digestibility in ducks fed normal p...In this study, we examined the varia- tions between the difference method and the enzyme- hydrolyzed casein method for determining endogenous amino acid loss and the true amino acid digestibility in ducks fed normal protein-containing diets. These methods were compared to the nitrogen-free （N-free） diet method. The difference method was based on soy- bean meal as the only protein source, with the experi- mental diets containing crude protein levels at 15% and 20%. The enzyme-hydrolyzed casein method was based on enzyme-hydrolyzed casein meal as the pro- tein source, with the experimental diet containing a crude protein level of 17.5%. The N-free diet was prepared with starches and paper fibers. In each meth- od,64 Tianfu meat drakes （7-weeks-old） with an av- erage body weight of 2.77±0.16 kg were used and divided into four groups, and fed four different diets. Each group contained four replicates of four drakes and they were force fed trial diets according to the Sirbald method for detecting their apparent amino aciddigestibility, endogenous amino acid loss and true a- mino acid digestibility. The results demonstrated that using the difference, enzyme-hydrolyzed casein and N-free diet methods, endogenous amino acid losses were 0. 9946,1. 2243 and 0. 9297 mg/g dry matter in- take （ DMI）, respectively. The true amino acid digest- ibility measured by the difference method was 88.93 %±4.43 %. Using the enzyme-hydrolyzed ca- sein method with two dietary crude protein levels of 15% and 20%, the digestibility was 91.15%±4.33% and 91.97%±4. 16%, respectively, and by the N-free diet methods with two dietary crude protein levels of 15% and 20% ,it was 88.55%±4.29% and 88.82 %±4.61%, respectively. The results suggested that when the dietary protein level was 15% to 20 %, the true amino acid digestibility and endogenous ami- no acid loss as determined by the difference method was more accurate than the values determined by the enzyme-hydrolyzed casein method.展开更多
The nonlinear Schr?dinger equation with a Dirac delta potential is considered in this paper.It is noted that the equation can be transformed into an equation with a drift-admitting jump.Then following the procedure pr...The nonlinear Schr?dinger equation with a Dirac delta potential is considered in this paper.It is noted that the equation can be transformed into an equation with a drift-admitting jump.Then following the procedure proposed in Chen and Deng(2018 Phys.Rev.E 98033302),a new second-order finite difference scheme is developed,which is justified by numerical examples.展开更多
基金Supported by Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)+3 种基金National Natural Science Foundation of China(12301556)Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)Basic Research Plan of Shanxi Province(202203021211129)。
文摘To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.
文摘This study investigates the complex heat transfer dynamics inmultilayer bifacial photovoltaic(bPV)solar modules under spectrally resolved solar irradiation.A novel numericalmodel is developed to incorporate internal heat generation resulting from optical absorption,grounded in the physical equations governing light-matter interactions within the module’smultilayer structure.The model accounts for reflection and transmission at each interface between adjacent layers,as well as absorption within individual layers,using the wavelength-dependent dielectric properties of constituent materials.These properties are used to calculate the spectral reflectance,transmittance,and absorption coefficients,enabling precise quantification of internal heat sources from irradiance incidents on both the front and rear surfaces of the module.The study further examines the influence of irradiance reflection on thermal behavior,evaluates the thermal impact of various supporting materials placed beneath the module,and analyzes the role of albedo in modifying heat distribution.By incorporating spectrally resolved heat generation across each layer often simplified or omitted in conventional models,the proposed approach enhances physical accuracy.The transient heat equation is solved using a one-dimensional finite difference(FD)method to produce detailed temperature profiles under multiple operating scenarios,including Standard Test Conditions(STC),Bifacial Standard Test Conditions(BSTC),Normal Operating Cell Temperature(NOCT),and Bifacial NOCT(BNOCT).The results offer valuable insights into the interplay between optical and thermal phenomena in bifacial systems,informing the design and optimization of more efficient photovoltaic technologies.
基金supported by National Natural Science Foundation of China(12273080).
文摘In response to the issue of fuzzy matching and association when optical observation data are matched with the orbital elements in a catalog database,this paper proposes a matching and association strategy based on the arcsegment difference method.First,a matching error threshold is set to match the observation data with the known catalog database.Second,the matching results for the same day are sorted on the basis of target identity and observation residuals.Different matching error thresholds and arc-segment dynamic association thresholds are then applied to categorize the observation residuals of the same target across different arc-segments,yielding matching results under various thresholds.Finally,the orbital residual is computed through orbit determination(OD),and the positional error is derived by comparing the OD results with the orbit track from the catalog database.The appropriate matching error threshold is then selected on the basis of these results,leading to the final matching and association of the fuzzy correlation data.Experimental results showed that the correct matching rate for data arc-segments is 92.34% when the matching error threshold is set to 720″,with the arc-segment difference method processing the results of an average matching rate of 97.62% within 8 days.The remaining 5.28% of the fuzzy correlation data are correctly matched and associated,enabling identification of orbital maneuver targets through further processing and analysis.This method substantially enhances the efficiency and accuracy of space target cataloging,offering robust technical support for dynamic maintenance of the space target database.
基金supported by the National Natural Science Foundation of China(12201228,12171047)the Fundamental Research Funds for the Central Universities(3034011102)supported by National Key R&D Program of China(2020YFA0713701).
文摘In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.
基金supported by the National Natural Science Foundation of China(No.42277175)Guangxi Emergency Management Department 2024 Innovation and Technology Research Project,China(No.2024GXYJ006)+2 种基金Hunan Provincial Department of Natural Resources Geological Exploration Project,China(No.2023ZRBSHZ056)The First National Natural Disaster Comprehensive Risk Survey in Hunan Province,China(No.2022-70)Guizhou Provincial Major Scientific and Technological Program,China(No.2023-425).
文摘A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficients(k)are constructed based on FDM.The rock bridge area was divided through k-means algorithm and the optimal number of clusters was determined by sum of squared errors(SSE)and elbow method.The influence of maximum principal stress and stress change rate as clustering indexes on the clustering results of rock bridges was compared by using Euclidean distance.The results show that using stress change rate as clustering index is more effective.When the joint coalescence coefficient is less than 0.6,there is no significant stress concentration in the middle area of adjacent joints,that is,no generation of rock bridge.In addition,the range of rock bridge is affected by the coalescence coefficient(k),the relative position of joints and the parameters of weak interlayer.
基金National Natural Science Foundation of China under Grant Nos.51978213 and 51778190the National Key Research and Development Program of China under Grant Nos.2017YFC0703605 and 2016YFC0701106。
文摘For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.
基金supported by the Key Laboratory of Road Construction Technology and Equipment(Chang’an University,No.300102253502)the Natural Science Foundation of Shandong Province of China(GrantNo.ZR2022YQ06)the Development Plan of Youth Innovation Team in Colleges and Universities of Shandong Province(Grant No.2022KJ140).
文摘In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.
基金Supported by the National Natural Science Foundation of China(Grant No.11571181)the Research Start-Up Foundation of Nantong University(Grant No.135423602051).
文摘In this paper,we give improved error estimates for linearized and nonlinear CrankNicolson type finite difference schemes of Ginzburg-Landau equation in two dimensions.For linearized Crank-Nicolson scheme,we use mathematical induction to get unconditional error estimates in discrete L^(2)and H^(1)norm.However,it is not applicable for the nonlinear scheme.Thus,based on a‘cut-off’function and energy analysis method,we get unconditional L^(2)and H^(1)error estimates for the nonlinear scheme,as well as boundedness of numerical solutions.In addition,if the assumption for exact solutions is improved compared to before,unconditional and optimal pointwise error estimates can be obtained by energy analysis method and several Sobolev inequalities.Finally,some numerical examples are given to verify our theoretical analysis.
文摘Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is analyzed. The dynamic model uses Euler-Bernoulli beam theory for cantilever beam, in which the driving force (or driving moment) of the electromagnetic system is taken into account, and the contact force between moving contact and stationary contact is simulated by the Kelvin-Voigt vis-coelastic...
基金The National Natural Science Foundation of China(No50475073,50775036)the High Technology Research Program of Jiangsu Province(NoBG2006035)
文摘An improved finite difference method （FDM）is described to solve existing problems such as low efficiency and poor convergence performance in the traditional method adopted to derive the pressure distribution of aerostatic bearings. A detailed theoretical analysis of the pressure distribution of the orifice-compensated aerostatic journal bearing is presented. The nonlinear dimensionless Reynolds equation of the aerostatic journal bearing is solved by the finite difference method. Based on the principle of flow equilibrium, a new iterative algorithm named the variable step size successive approximation method is presented to adjust the pressure at the orifice in the iterative process and enhance the efficiency and convergence performance of the algorithm. A general program is developed to analyze the pressure distribution of the aerostatic journal bearing by Matlab tool. The results show that the improved finite difference method is highly effective, reliable, stable, and convergent. Even when very thin gas film thicknesses （less than 2 Win）are considered, the improved calculation method still yields a result and converges fast.
基金Projects(51308273,41372307,41272326) supported by the National Natural Science Foundation of ChinaProjects(2010(A)06-b) supported by Science and Technology Fund of Yunan Provincial Communication Department,China
文摘Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with theoretical method, the finite difference method has been verified to be feasible by a case study. It is found that under seismic loading, the dynamic response of anchorage system is synchronously fluctuated with the seismic vibration. The change of displacement amplitude of material points is slight, and comparatively speaking, the displacement amplitude of the outside point is a little larger than that of the inside point, which shows amplification effect of surface. While the axial force amplitude transforms considerably from the inside to the outside. It increases first and reaches the peak value in the intersection between the anchoring section and free section, then decreases slowly in the free section. When considering damping effect of anchorage system, the finite difference method can reflect the time attenuation characteristic better, and the calculating result would be safer and more reasonable than the dynamic steady-state theoretical method. What is more, the finite difference method can be applied to the dynamic response analysis of harmonic and seismic random vibration for all kinds of anchor, and hence has a broad application prospect.
基金supported by the National Natural Science Foundation of China(Nos.11171193 and11371229)the Natural Science Foundation of Shandong Province(No.ZR2014AM033)the Science and Technology Development Project of Shandong Province(No.2012GGB01198)
文摘An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.
基金supported by the National Natural Science Foundation of China (No.50805056)New Century Excellent Talents in University (No.NCET-09-0396)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education (2009)
文摘Thermal stress simulation can provide a scientific reference to eliminate defects such as crack,residual stress centralization and deformation etc.,caused by thermal stress during casting solidification.To study the thermal stress distribution during casting process,a unilateral thermal-stress coupling model was employed to simulate 3D casting stress using Finite Difference Method(FDM),namely all the traditional thermal-elastic-plastic equations are numerically and differentially discrete.A FDM/FDM numerical simulation system was developed to analyze temperature and stress fields during casting solidification process.Two practical verifications were carried out,and the results from simulation basically coincided with practical cases.The results indicated that the FDM/FDM stress simulation system can be used to simulate the formation of residual stress,and to predict the occurrence of hot tearing.Because heat transfer and stress analysis are all based on FDM,they can use the same FD model,which can avoid the matching process between different models,and hence reduce temperature-load transferring errors.This approach makes the simulation of fluid flow,heat transfer and stress analysis unify into one single model.
基金supported by the Yunnan Provincial Applied Basic Research Program of China(No. KKSY201207019)
文摘A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.
基金the Shanghai Excellent Young Professor Candidate Program (Grant No.04YQHB139)
文摘With the finite difference method to calculate the temperature distribution in mass concrete structures, the solution precision will increase with a smaller step size, at the cost of computational time. In view of the basic characteristics of the finite difference method, a simple yet powerful improvement is introduced. By multiplying the adiabatic temperature function with a correction factor, the precision of the solution can be assured without an increase in the computation time. In addition, the correction rules for three types of commonly used concrete hydration formulas are investigated.
基金the National Natural Science Foundation of China under Grant Number NSFC 11801302Tsinghua University Initiative Scientific Research Program.Yang Yang is supported by the NSF Grant DMS-1818467.
文摘In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical reaction,say zi,is between 0 and 1,withΣz_(i)=1.Due to the lack of maximum-principle,most of the previous bound-preserving technique cannot be applied directly.To preserve those bounds,we will use the positivity-preserving technique to all the zi'is and enforceΣz_(i)=1 by constructing conservative schemes,thanks to conservative time integrations and consistent numerical fluxes in the system.Moreover,detonation is an extreme singular mode of flame propagation in premixed gas,and the model contains a significant stiff source.It is well known that for hyperbolic equations with stiff source,the transition points in the numerical approximations near the shocks may trigger spurious shock speed,leading to wrong shock position.Intuitively,the high-order weighted essentially non-oscillatory(WENO)scheme,which can suppress oscillations near the discontinuities,would be a good choice for spatial discretization.However,with the nonlinear weights,the numerical fluxes are no longer“consistent”,leading to nonconservative numerical schemes and the bound-preserving technique does not work.Numerical experiments demonstrate that,without further numerical techniques such as subcell resolutions,the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.
基金Project(Z110803)supported by the State Key Laboratory of Geomechanics and Geotechnical Engineering,ChinaProject(2008AA062303)supported by the National High Technology Research and Development Program of China
文摘A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.
文摘In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.
文摘In this study, we examined the varia- tions between the difference method and the enzyme- hydrolyzed casein method for determining endogenous amino acid loss and the true amino acid digestibility in ducks fed normal protein-containing diets. These methods were compared to the nitrogen-free （N-free） diet method. The difference method was based on soy- bean meal as the only protein source, with the experi- mental diets containing crude protein levels at 15% and 20%. The enzyme-hydrolyzed casein method was based on enzyme-hydrolyzed casein meal as the pro- tein source, with the experimental diet containing a crude protein level of 17.5%. The N-free diet was prepared with starches and paper fibers. In each meth- od,64 Tianfu meat drakes （7-weeks-old） with an av- erage body weight of 2.77±0.16 kg were used and divided into four groups, and fed four different diets. Each group contained four replicates of four drakes and they were force fed trial diets according to the Sirbald method for detecting their apparent amino aciddigestibility, endogenous amino acid loss and true a- mino acid digestibility. The results demonstrated that using the difference, enzyme-hydrolyzed casein and N-free diet methods, endogenous amino acid losses were 0. 9946,1. 2243 and 0. 9297 mg/g dry matter in- take （ DMI）, respectively. The true amino acid digest- ibility measured by the difference method was 88.93 %±4.43 %. Using the enzyme-hydrolyzed ca- sein method with two dietary crude protein levels of 15% and 20%, the digestibility was 91.15%±4.33% and 91.97%±4. 16%, respectively, and by the N-free diet methods with two dietary crude protein levels of 15% and 20% ,it was 88.55%±4.29% and 88.82 %±4.61%, respectively. The results suggested that when the dietary protein level was 15% to 20 %, the true amino acid digestibility and endogenous ami- no acid loss as determined by the difference method was more accurate than the values determined by the enzyme-hydrolyzed casein method.
基金supported by the National Natural Science Foundation of China(Grant Nos.11601517,61772542).
文摘The nonlinear Schr?dinger equation with a Dirac delta potential is considered in this paper.It is noted that the equation can be transformed into an equation with a drift-admitting jump.Then following the procedure proposed in Chen and Deng(2018 Phys.Rev.E 98033302),a new second-order finite difference scheme is developed,which is justified by numerical examples.
