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 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.展开更多
In this study,we analyze the convergence of the finite difference method on non-uniform grids and provide examples to demonstrate its effectiveness in approximating fractional differential equations involving the frac...In this study,we analyze the convergence of the finite difference method on non-uniform grids and provide examples to demonstrate its effectiveness in approximating fractional differential equations involving the fractional Laplacian.By utilizing non-uniform grids,it becomes possible to achieve higher accuracy and improved resolution in specific regions of interest.Overall,our findings indicate that finite difference approximation on non-uniform grids can serve as a dependable and efficient tool for approximating fractional Laplacians across a diverse array of applications.展开更多
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.展开更多
Compared with other migration methods, reverse-time migration is based on a precise wave equation, not an approximation, and performs extrapolation in the depth domain rather than the time domain. It is highly accurat...Compared with other migration methods, reverse-time migration is based on a precise wave equation, not an approximation, and performs extrapolation in the depth domain rather than the time domain. It is highly accurate and not affected by strong subsurface structure complexity and horizontal velocity variations. The difference method based on triangular grids maintains the simplicity of the difference method and the precision of the finite element method. It can be used directly for forward modeling on models with complex top surfaces and migration without statics preprocessing. We apply a finite difference method based on triangular grids for post-stack reverse-time migration for the first time. Tests on model data verify that the combination of the two methods can achieve near-perfect results in application.展开更多
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.展开更多
D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration a...D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.展开更多
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 new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite differen...A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.展开更多
The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux ...The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.展开更多
In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite differ...In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.展开更多
Characteristic finite difference fractional step schemes are put forward. The electric potential equation is described by a seven-point finite difference scheme, and the electron and hole concentration equations are t...Characteristic finite difference fractional step schemes are put forward. The electric potential equation is described by a seven-point finite difference scheme, and the electron and hole concentration equations are treated by a kind of characteristic finite difference fractional step methods. The temperature equation is described by a fractional step method. Thick and thin grids are made use of to form a complete set. Piecewise threefold quadratic interpolation, symmetrical extension, calculus of variations, commutativity of operator product, decomposition of high order difference operators and prior estimates are also made use of. Optimal order estimates in l2 norm are derived to determine the error of the approximate solution. The well-known problem is thorongley and completely solred.展开更多
Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly use...Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly used mesh-free FD methods and can accurately simulate seismic wave propagation in the non-rectangular computational domain.In this paper,we propose a perfectly matched layer(PML)boundary condition for a meshfree FD solution of the elastic wave equation,which can be applied to the boundaries of the non-rectangular velocity model.The performance of the PML is,however,severely reduced for near-grazing incident waves and low-frequency waves.We thus also propose the complexfrequency-shifted PML(CFS-PML)boundary condition for a mesh-free FD solution of the elastic wave equation.For two PML boundary conditions,we derive unsplit time-domain expressions by constructing auxiliary differential equations,both of which require less memory and are easy for programming.Numerical experiments demonstrate that these two PML boundary conditions effectively eliminate artificial boundary reflections in mesh-free FD simulations.When compared with the PML boundary condition,the CFS-PML boundary condition results in better absorption for near-grazing incident waves and evanescent waves.展开更多
The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of c...The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity 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.展开更多
The dynamic calculations of slender marine risers, such as Finite Element Method （FEM） or Modal Expansion Solution Method （MESM）, are mainly for the slender structures with their both ends hinged to the surface an...The dynamic calculations of slender marine risers, such as Finite Element Method （FEM） or Modal Expansion Solution Method （MESM）, are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation （FDA） method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bemoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser＇s free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.展开更多
基金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 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 Spanish MINECO through Juan de la Cierva fellow-ship FJC2021-046953-I.
文摘In this study,we analyze the convergence of the finite difference method on non-uniform grids and provide examples to demonstrate its effectiveness in approximating fractional differential equations involving the fractional Laplacian.By utilizing non-uniform grids,it becomes possible to achieve higher accuracy and improved resolution in specific regions of interest.Overall,our findings indicate that finite difference approximation on non-uniform grids can serve as a dependable and efficient tool for approximating fractional Laplacians across a diverse array of applications.
基金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.
基金sponsored by National Natural Science Foundation(40474041)National Symposium of 863(2006AA06Z206)+1 种基金National Symposium of 973(2007CB209605)CNPC Geophysical Key Laboratory of the China University of Petroleum (East China) Research Department
文摘Compared with other migration methods, reverse-time migration is based on a precise wave equation, not an approximation, and performs extrapolation in the depth domain rather than the time domain. It is highly accurate and not affected by strong subsurface structure complexity and horizontal velocity variations. The difference method based on triangular grids maintains the simplicity of the difference method and the precision of the finite element method. It can be used directly for forward modeling on models with complex top surfaces and migration without statics preprocessing. We apply a finite difference method based on triangular grids for post-stack reverse-time migration for the first time. Tests on model data verify that the combination of the two methods can achieve near-perfect results in application.
文摘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.
基金This project is sponsored by the Specialized Prophasic Basic Research of the"973"Programme,contract No:2001cca02300
文摘D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.
基金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.
基金the National Natural Science Foundation of China
文摘A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.
基金the National Key R&D Program of China(No.2018YFE0180900).
文摘The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.
文摘In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.
基金This work is supported by the Major State Basic Research Program of China (19990328), the National Tackling Key Problem Program, the National Science Foundation of China (10271066 and 0372052), and the Doctorate Foundation of the Ministry of Education of China (20030422047).
文摘Characteristic finite difference fractional step schemes are put forward. The electric potential equation is described by a seven-point finite difference scheme, and the electron and hole concentration equations are treated by a kind of characteristic finite difference fractional step methods. The temperature equation is described by a fractional step method. Thick and thin grids are made use of to form a complete set. Piecewise threefold quadratic interpolation, symmetrical extension, calculus of variations, commutativity of operator product, decomposition of high order difference operators and prior estimates are also made use of. Optimal order estimates in l2 norm are derived to determine the error of the approximate solution. The well-known problem is thorongley and completely solred.
基金supported by the National Science and Technology Major Project(2016ZX05006-002)the National Natural Science Foundation of China(Nos.41874153,41504097)
文摘Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly used mesh-free FD methods and can accurately simulate seismic wave propagation in the non-rectangular computational domain.In this paper,we propose a perfectly matched layer(PML)boundary condition for a meshfree FD solution of the elastic wave equation,which can be applied to the boundaries of the non-rectangular velocity model.The performance of the PML is,however,severely reduced for near-grazing incident waves and low-frequency waves.We thus also propose the complexfrequency-shifted PML(CFS-PML)boundary condition for a mesh-free FD solution of the elastic wave equation.For two PML boundary conditions,we derive unsplit time-domain expressions by constructing auxiliary differential equations,both of which require less memory and are easy for programming.Numerical experiments demonstrate that these two PML boundary conditions effectively eliminate artificial boundary reflections in mesh-free FD simulations.When compared with the PML boundary condition,the CFS-PML boundary condition results in better absorption for near-grazing incident waves and evanescent waves.
文摘The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity 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.
基金supported and sponsored jointly by the National Natural Science Foundation of China(Grand Nos.51009092 and 50909061)Doctoral Foundation of the Ministry of Education of China(Grand No.20090073120013)the National High Technology Research and Development Program of China(863Program,Grand No.2008AA092301-1)
文摘The dynamic calculations of slender marine risers, such as Finite Element Method （FEM） or Modal Expansion Solution Method （MESM）, are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation （FDA） method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bemoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser＇s free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.
