In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs...In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.展开更多
One-dimensional heat equation was solved for different higher-order finite difference schemes, namely, forward time and fourth-order centered space explicit method, backward time and fourth-order centered space implic...One-dimensional heat equation was solved for different higher-order finite difference schemes, namely, forward time and fourth-order centered space explicit method, backward time and fourth-order centered space implicit method, and fourth-order implicit Crank-Nicolson finite difference method. Higher-order schemes have complexity in computing values at the neighboring points to the boundaries. It is required there a specification of the values of field variables at some points exterior to the domain. The complexity was incorporated using Hicks approximation. The convergence and stability analysis was also computed for those higher-order finite difference explicit and implicit methods in case of solving a one dimensional heat equation. The obtained numerical results were compared with exact solutions. It is found that backward time and fourth-order centered space implicit scheme along with Hicks approximation performed well over the other mentioned higher-order approaches.展开更多
A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spheric...A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
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.展开更多
In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved dire...In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly. The Poisson equation is approximated by fourth-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results.展开更多
An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solutio...An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.展开更多
Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been...Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.展开更多
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.展开更多
In order to address the complex uncertainties caused by interfacing between the fuzziness and randomness of the safety problem for embankment engineering projects, and to evaluate the safety of embankment engineering ...In order to address the complex uncertainties caused by interfacing between the fuzziness and randomness of the safety problem for embankment engineering projects, and to evaluate the safety of embankment engineering projects more scientifically and reasonably, this study presents the fuzzy logic modeling of the stochastic finite element method （SFEM） based on the harmonious finite element （HFE） technique using a first-order approximation theorem. Fuzzy mathematical models of safety repertories were introduced into the SFEM to analyze the stability of embankments and foundations in order to describe the fuzzy failure procedure for the random safety performance function. The fuzzy models were developed with membership functions with half depressed gamma distribution, half depressed normal distribution, and half depressed echelon distribution. The fuzzy stochastic mathematical algorithm was used to comprehensively study the local failure mechanism of the main embankment section near Jingnan in the Yangtze River in terms of numerical analysis for the probability integration of reliability on the random field affected by three fuzzy factors. The result shows that the middle region of the embankment is the principal zone of concentrated failure due to local fractures. There is also some local shear failure on the embankment crust. This study provides a referential method for solving complex multi-uncertainty problems in engineering safety analysis.展开更多
Based on the Reddy's theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadr...Based on the Reddy's theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadrature (DQ) method of nonlinear analysis to the problem was presented. New DQ approach, presented by Wang and Bert (DQWB), is extended to handle the multiple boundary conditions of plates. The techniques were also further extended to simplify nonlinear computations. The numerical convergence and comparison of solutions were studied. The results show that the DQ method presented is very reliable and valid. Moreover, the influences of geometric and material parameters as well as the transverse shear deformations on nonlinear bending were investigated. Numerical results show the influence of the shear deformation on the static bending of orthotropic moderately thick plate is significant.展开更多
Models of the coupled nonlinear Schr<span style="white-space:nowrap;">ödinger equations submit various critical physical phenomena with a typical equation for optical fibres with linear refrac...Models of the coupled nonlinear Schr<span style="white-space:nowrap;">ödinger equations submit various critical physical phenomena with a typical equation for optical fibres with linear refraction. In this article, we will presuppose the Compact Finite Difference method with Runge-Kutta of order 4 (explicit) method, which is sixth-order and fourth-order in space and time respectively, to solve coupled nonlinear Schr<span style="white-space:nowrap;">ödinger equations. Many methods used to solve coupled nonlinear Schr<span style="white-space:nowrap;">ödinger equations are second order in time and need to use extra-technique to rise up to fourth-order as Richardson Extrapolation technique. The scheme obtained is immediately fourth-order in one step. This approach is a conditionally stable method. The conserved quantities and the exact single soliton solution indicate the competence and accuracy of the article’s suggestion schemes. Furthermore, the article discusses the two solitons interaction dynamics.展开更多
This paper looks at the development of a class of Exponential Compact Higher Order (ECHO) schemes and attempts to comprehend their behaviour by introducing different combinations of discrete source function and its de...This paper looks at the development of a class of Exponential Compact Higher Order (ECHO) schemes and attempts to comprehend their behaviour by introducing different combinations of discrete source function and its derivatives. The characteristic analysis is performed for one-dimensional schemes to understand the efficiency of the scheme and a similar analysis has been introduced for higher dimensional schemes. Finally, the developed schemes are used to solve several example problems and compared the error norms and rates of convergence.展开更多
This paper uses the Taylor expansion to seek an approximate Korteweg- de Vries equation （KdV） solution to a higher-order traffic flow model with sufficiently large diffusion. It demonstrates the validity of the appr...This paper uses the Taylor expansion to seek an approximate Korteweg- de Vries equation （KdV） solution to a higher-order traffic flow model with sufficiently large diffusion. It demonstrates the validity of the approximate KdV solution considering all the related parameters to ensure the physical boundedness and the stability of the solution. Moreover, when the viscosity coefficient depends on the density and velocity of the flow, the wave speed of the KdV solution is naturally related to either the first or the second characteristic field. The finite element method is extended to solve the model and examine the stability and accuracy of the approximate KdV solution.展开更多
A mathematical model comprising of nonlinear reaction, diffusion, and convection mechanisms seen in natural and anthropogenic processes is numerically investigated here. It is proposed that a higher order numerical sc...A mathematical model comprising of nonlinear reaction, diffusion, and convection mechanisms seen in natural and anthropogenic processes is numerically investigated here. It is proposed that a higher order numerical scheme of finite difference method be used in conjunction with an iterative approach in order to solve the nonlinear one dimensional convection-diffusion-reaction equation. To account for the wide variety of physical characteristics and boundary conditions, an iterative approach is presented that yields a reliable and precise solution every time. We examined the accuracy and operational efficiency of two distinct finite difference approaches. The efficiency of the system is determined by comparing the estimated results to the appropriate analytical solution by adhering to established norms. Coherence and convergence were analyzed for each approach. The simulation results demonstrate the efficacy and accuracy of these methods in solving nonlinear convection- diffusion-reaction equations. Convection-diffusion-reaction equation modeling is critical for employing the offered results in heat and mass transport processes.展开更多
In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order a...In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order accuracy, while the exponential-sum-approximation (ESA) is used to approximate the variable-order Caputo fractional derivative in the temporal direction, and a novel spatial sixth-order hybrid ESA-CCD method is implemented successfully. Finally, the accuracy of the proposed method is verified by numerical experiments.展开更多
Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, et...Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, etc. Therefore, it is of great significance to study fractional-order time-delay differential equations. In this paper, we discuss a finite volume element method for a class of fractional-order neutral time-delay differential equations. By introducing an intermediate variable, the fourth-order problem is transformed into a system of equations consisting of two second-order partial differential equations. The L1 formula is used to approximate the time fractional order derivative terms, and the finite volume element method is used in space. A fully discrete format of the equations is established, and we prove the existence, uniqueness, convergence and stability of the solution. Finally, the validity of the format is verified by numerical examples.展开更多
This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-d...This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.展开更多
In this paper,we develop twofinite difference weighted essentially non-oscillatory(WENO)schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi(DP)andµ-Degasperis-Procesi(µDP)equations,whi...In this paper,we develop twofinite difference weighted essentially non-oscillatory(WENO)schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi(DP)andµ-Degasperis-Procesi(µDP)equations,which contain nonlinear high order derivatives,and possibly peakon solutions or shock waves.By introducing auxiliary variable(s),we rewrite the DP equation as a hyperbolic-elliptic system,and theµDP equation as afirst order system.Then we choose a linearfinite difference scheme with suitable order of accuracy for the auxiliary variable(s),and twofinite difference WENO schemes with unequal-sized sub-stencils for the primal variable.One WENO scheme uses one large stencil and several smaller stencils,and the other WENO scheme is based on the multi-resolution framework which uses a se-ries of unequal-sized hierarchical central stencils.Comparing with the classical WENO scheme which uses several small stencils of the same size to make up a big stencil,both WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil and enjoy the freedom of arbitrary positive linear weights.Another advantage is that thefinal reconstructed polynomial on the target cell is a polynomial of the same de-gree as the polynomial over the big stencil,while the classicalfinite difference WENO reconstruction can only be obtained for specific points inside the target interval.Nu-merical tests are provided to demonstrate the high order accuracy and non-oscillatory properties of the proposed schemes.展开更多
A high fidelity flow simulation for complex geometries for high Reynolds number(Re)flow is still very challenging,requiring a more powerful HPC system.However,the development of HPC with traditional CPU architecture s...A high fidelity flow simulation for complex geometries for high Reynolds number(Re)flow is still very challenging,requiring a more powerful HPC system.However,the development of HPC with traditional CPU architecture suffers bottlenecks due to its high power consumption and technical difficulties.Heterogeneous architecture computation is raised to be a promising solution to the challenges of HPC development.GPU accelerating technology has been utilized in low order scheme CFD solvers on the structured grid and high order scheme solvers on unstructured meshes.The high-order finite difference methods on structured grids possess many advantages,e.g.,high efficiency,robustness,and low storage.However,the strong dependence among points for a high-order finite difference scheme still limits its application on the GPU platform.In the present work,we propose a set of hardware-aware technology to optimize data transfer efficiency between CPU and GPU,as well as communication efficiency among GPUs.An in-house multi-block structured CFD solver with high order finite difference methods on curvilinear coordinates is ported onto the GPU platform and obtains satisfying performance with a speedup maximum of around 2000x over a single CPU core.This work provides an efficient solution to apply GPU computing in CFD simulation with specific high order finite difference methods on current GPU heterogeneous computers.The test shows that significant accelerating effects can be achieved for different GPUs.展开更多
文摘In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.
文摘One-dimensional heat equation was solved for different higher-order finite difference schemes, namely, forward time and fourth-order centered space explicit method, backward time and fourth-order centered space implicit method, and fourth-order implicit Crank-Nicolson finite difference method. Higher-order schemes have complexity in computing values at the neighboring points to the boundaries. It is required there a specification of the values of field variables at some points exterior to the domain. The complexity was incorporated using Hicks approximation. The convergence and stability analysis was also computed for those higher-order finite difference explicit and implicit methods in case of solving a one dimensional heat equation. The obtained numerical results were compared with exact solutions. It is found that backward time and fourth-order centered space implicit scheme along with Hicks approximation performed well over the other mentioned higher-order approaches.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61301063 and 41305017)
文摘A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
基金supported 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.
文摘In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly. The Poisson equation is approximated by fourth-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results.
文摘An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.
文摘Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.50379046)the Doctoral Fund of the Ministry of Education of China(Grant No.A50221)
文摘In order to address the complex uncertainties caused by interfacing between the fuzziness and randomness of the safety problem for embankment engineering projects, and to evaluate the safety of embankment engineering projects more scientifically and reasonably, this study presents the fuzzy logic modeling of the stochastic finite element method （SFEM） based on the harmonious finite element （HFE） technique using a first-order approximation theorem. Fuzzy mathematical models of safety repertories were introduced into the SFEM to analyze the stability of embankments and foundations in order to describe the fuzzy failure procedure for the random safety performance function. The fuzzy models were developed with membership functions with half depressed gamma distribution, half depressed normal distribution, and half depressed echelon distribution. The fuzzy stochastic mathematical algorithm was used to comprehensively study the local failure mechanism of the main embankment section near Jingnan in the Yangtze River in terms of numerical analysis for the probability integration of reliability on the random field affected by three fuzzy factors. The result shows that the middle region of the embankment is the principal zone of concentrated failure due to local fractures. There is also some local shear failure on the embankment crust. This study provides a referential method for solving complex multi-uncertainty problems in engineering safety analysis.
文摘Based on the Reddy's theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadrature (DQ) method of nonlinear analysis to the problem was presented. New DQ approach, presented by Wang and Bert (DQWB), is extended to handle the multiple boundary conditions of plates. The techniques were also further extended to simplify nonlinear computations. The numerical convergence and comparison of solutions were studied. The results show that the DQ method presented is very reliable and valid. Moreover, the influences of geometric and material parameters as well as the transverse shear deformations on nonlinear bending were investigated. Numerical results show the influence of the shear deformation on the static bending of orthotropic moderately thick plate is significant.
文摘Models of the coupled nonlinear Schr<span style="white-space:nowrap;">ödinger equations submit various critical physical phenomena with a typical equation for optical fibres with linear refraction. In this article, we will presuppose the Compact Finite Difference method with Runge-Kutta of order 4 (explicit) method, which is sixth-order and fourth-order in space and time respectively, to solve coupled nonlinear Schr<span style="white-space:nowrap;">ödinger equations. Many methods used to solve coupled nonlinear Schr<span style="white-space:nowrap;">ödinger equations are second order in time and need to use extra-technique to rise up to fourth-order as Richardson Extrapolation technique. The scheme obtained is immediately fourth-order in one step. This approach is a conditionally stable method. The conserved quantities and the exact single soliton solution indicate the competence and accuracy of the article’s suggestion schemes. Furthermore, the article discusses the two solitons interaction dynamics.
文摘This paper looks at the development of a class of Exponential Compact Higher Order (ECHO) schemes and attempts to comprehend their behaviour by introducing different combinations of discrete source function and its derivatives. The characteristic analysis is performed for one-dimensional schemes to understand the efficiency of the scheme and a similar analysis has been introduced for higher dimensional schemes. Finally, the developed schemes are used to solve several example problems and compared the error norms and rates of convergence.
基金supported by the National Natural Science Foundation of China(Nos.11072141 and11272199)the National Basic Research Program of China(No.2012CB725404)+2 种基金the Shanghai Program for Innovative Research Team in Universitiesthe Research Grants Council of the Hong KongSpecial Administrative Region,China(No.HKU7184/10E)the National Research Foundationof Korea(MEST)(No.NRF-2010-0029446)
文摘This paper uses the Taylor expansion to seek an approximate Korteweg- de Vries equation （KdV） solution to a higher-order traffic flow model with sufficiently large diffusion. It demonstrates the validity of the approximate KdV solution considering all the related parameters to ensure the physical boundedness and the stability of the solution. Moreover, when the viscosity coefficient depends on the density and velocity of the flow, the wave speed of the KdV solution is naturally related to either the first or the second characteristic field. The finite element method is extended to solve the model and examine the stability and accuracy of the approximate KdV solution.
文摘A mathematical model comprising of nonlinear reaction, diffusion, and convection mechanisms seen in natural and anthropogenic processes is numerically investigated here. It is proposed that a higher order numerical scheme of finite difference method be used in conjunction with an iterative approach in order to solve the nonlinear one dimensional convection-diffusion-reaction equation. To account for the wide variety of physical characteristics and boundary conditions, an iterative approach is presented that yields a reliable and precise solution every time. We examined the accuracy and operational efficiency of two distinct finite difference approaches. The efficiency of the system is determined by comparing the estimated results to the appropriate analytical solution by adhering to established norms. Coherence and convergence were analyzed for each approach. The simulation results demonstrate the efficacy and accuracy of these methods in solving nonlinear convection- diffusion-reaction equations. Convection-diffusion-reaction equation modeling is critical for employing the offered results in heat and mass transport processes.
文摘In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order accuracy, while the exponential-sum-approximation (ESA) is used to approximate the variable-order Caputo fractional derivative in the temporal direction, and a novel spatial sixth-order hybrid ESA-CCD method is implemented successfully. Finally, the accuracy of the proposed method is verified by numerical experiments.
文摘Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, etc. Therefore, it is of great significance to study fractional-order time-delay differential equations. In this paper, we discuss a finite volume element method for a class of fractional-order neutral time-delay differential equations. By introducing an intermediate variable, the fourth-order problem is transformed into a system of equations consisting of two second-order partial differential equations. The L1 formula is used to approximate the time fractional order derivative terms, and the finite volume element method is used in space. A fully discrete format of the equations is established, and we prove the existence, uniqueness, convergence and stability of the solution. Finally, the validity of the format is verified by numerical examples.
基金supported by the National Natural Science Foundation of China (Grants 11571223, 51404160)Shanxi Province Science Foundation for Youths (Grant 2014021025-1)
文摘This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.
基金supported by National Natural Science Foundation of China(Grant No.12071455)supported by National Natural Science Foundation of China(Grant No.11871428)。
文摘In this paper,we develop twofinite difference weighted essentially non-oscillatory(WENO)schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi(DP)andµ-Degasperis-Procesi(µDP)equations,which contain nonlinear high order derivatives,and possibly peakon solutions or shock waves.By introducing auxiliary variable(s),we rewrite the DP equation as a hyperbolic-elliptic system,and theµDP equation as afirst order system.Then we choose a linearfinite difference scheme with suitable order of accuracy for the auxiliary variable(s),and twofinite difference WENO schemes with unequal-sized sub-stencils for the primal variable.One WENO scheme uses one large stencil and several smaller stencils,and the other WENO scheme is based on the multi-resolution framework which uses a se-ries of unequal-sized hierarchical central stencils.Comparing with the classical WENO scheme which uses several small stencils of the same size to make up a big stencil,both WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil and enjoy the freedom of arbitrary positive linear weights.Another advantage is that thefinal reconstructed polynomial on the target cell is a polynomial of the same de-gree as the polynomial over the big stencil,while the classicalfinite difference WENO reconstruction can only be obtained for specific points inside the target interval.Nu-merical tests are provided to demonstrate the high order accuracy and non-oscillatory properties of the proposed schemes.
基金National Numerical Windtunnel Project,the National Natural Science Foundation projects(91952103,11772323,11621202)Fundamental Research Funds for the Central Universities.
文摘A high fidelity flow simulation for complex geometries for high Reynolds number(Re)flow is still very challenging,requiring a more powerful HPC system.However,the development of HPC with traditional CPU architecture suffers bottlenecks due to its high power consumption and technical difficulties.Heterogeneous architecture computation is raised to be a promising solution to the challenges of HPC development.GPU accelerating technology has been utilized in low order scheme CFD solvers on the structured grid and high order scheme solvers on unstructured meshes.The high-order finite difference methods on structured grids possess many advantages,e.g.,high efficiency,robustness,and low storage.However,the strong dependence among points for a high-order finite difference scheme still limits its application on the GPU platform.In the present work,we propose a set of hardware-aware technology to optimize data transfer efficiency between CPU and GPU,as well as communication efficiency among GPUs.An in-house multi-block structured CFD solver with high order finite difference methods on curvilinear coordinates is ported onto the GPU platform and obtains satisfying performance with a speedup maximum of around 2000x over a single CPU core.This work provides an efficient solution to apply GPU computing in CFD simulation with specific high order finite difference methods on current GPU heterogeneous computers.The test shows that significant accelerating effects can be achieved for different GPUs.
