首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
We investigate the numerical behavior of a fourth-order accurate method. The results are compared with those of a second-order method. We have verified the rate of convergence for the numerical solution. As test cases the simple hyperbolic model equationu 1+u x =0 and the two-dimensional Euler equations over backward-facing step have been used. The fourth-order method has been implemented on a dataparallel computer, and the difference operators have been designed to minimize the bandwidth. We also derive boundary modified, semidefinite artificial viscosity operators of arbitrary order of accuracy. The viscosity operators are presented in a form that is particularly well-suited for the implementation on dataparallel computers.  相似文献   

2.
In this paper a dissipative exponentially-fitted method for the numerical integration of the Schrödinger equation and related problems is developed. The method is called dissipative since is a nonsymmetric multistep method. An application to the the resonance problem of the radial Schrödinger equation and to other well known related problems indicates that the new method is more efficient than the corresponding classical dissipative method and other well known methods. Based on the new method and the method of Raptis and Cash a new variable-step method is obtained. The application of the new variable-step method to the coupled differential equations arising from the Schrödinger equation indicates the power of the new approach.  相似文献   

3.
Several interesting problems in neuroscience are of multiscale type, i.e. possess different temporal and spatial scales that cannot be disregarded. Such characteristics impose severe burden to numerical simulations since the need to resolve small scale features pushes the computational costs to unreasonable levels. Classical numerical methods that do not resolve the small scales suffer from spurious oscillations and lack of precision.This paper presents an innovative numerical method of multiscale type that ameliorates these maladies. As an example we consider the case of a cable equation modeling heterogeneous dendrites. Our method is not only easy to parallelize, but it is also nodally exact, i.e., it matches the values of the exact solution at every node of the discretization mesh, for a class of problems.To show the validity of our scheme under different physiological regimes, we describe how the model behaves whenever the dendrites are thin or long, or the longitudinal conductance is small. We also consider the case of a large number of synapses and of large or low membrane conductance.  相似文献   

4.
In this paper a singularly perturbed fourth-order ordinary differential equation is considered. The differential equation is transformed into a coupled system of singularly perturbed equations. A hybrid finite difference scheme on a Vulanovi?–Shishkin mesh is used to discretize the system. This hybrid difference scheme is a combination of a non-equidistant generalization of the Numerov scheme and the central difference scheme based on the relation between the local mesh widths and the perturbation parameter. We will show that the scheme is maximum-norm stable, although the difference scheme may not satisfy the maximum principle. The scheme is proved to be almost fourth-order uniformly convergent in the discrete maximum norm. Numerical results are presented for supporting the theoretical results.  相似文献   

5.
6.
In this paper we define a new accurate fast implicit method for the finite difference solution of the two dimensional parabolic partial differential equations with first level condition, which may be obtained by any other method. The stability region is discussed. The suggested method is considered as an accelerating technique for the implicit finite difference scheme, which is used to find the first level condition. The obtained results are compared with some famous finite difference schemes and it is in satisfactory agreement with the exact solution.  相似文献   

7.
For a Bose-Einstein Condensate placed in a rotating trap and confined in the z-axis, a multisymplectic difference scheme was constructed to investigate the evolution of vortices in this paper. First, we look for a steady state solution of the imaginary time G-P equation. Then, we numerically study the vortices's development in real time, starting with the solution in imaginary time as initial value.  相似文献   

8.
We study a finite difference continuation (FDC) method for computing energy levels and wave functions of Bose-Einstein condensates (BEC), which is governed by the Gross-Pitaevskii equation (GPE). We choose the chemical potential λ as the continuation parameter so that the proposed algorithm can compute all energy levels of the discrete GPE. The GPE is discretized using the second-order finite difference method (FDM), which is treated as a special case of finite element methods (FEM) using the piecewise bilinear and linear interpolatory functions. Thus the mathematical theory of FEM for elliptic eigenvalue problems (EEP) also holds for the Schrödinger eigenvalue problem (SEP) associated with the GPE. This guarantees the existence of discrete numerical solutions for the ground-state as well as excited-states of the SEP in the variational form. We also study superconvergence of FDM for solution derivatives of parameter-dependent problems (PDP). It is proved that the superconvergence O(ht) in the discrete H1 norm is achieved, where t=2 and t=1.5 for rectangular and polygonal domains, respectively, and h is the maximal boundary length of difference grids. Moreover, the FDC algorithm can be implemented very efficiently using a simplified two-grid scheme for computing energy levels of the BEC. Numerical results are reported for the ground-state of two-coupled NLS defined in a large square domain, and in particular, for the second-excited state solutions of the 2D BEC in a periodic potential.  相似文献   

9.
In the present work, a numerical study has been carried out for the singularly perturbed generalized Burgers-Huxley equation using a three-step Taylor-Galerkin finite element method. A Burgers-Huxley equation represents the traveling wave phenomena. In singular perturbed problems, a very small positive parameter, ?, called the singular perturbation parameter is multiplied with the highest order derivative term. As this parameter tends towards zero, the problem exhibits boundary layers. The traditional methods fail to capture the boundary layers when ? becomes very small. In this paper a three-step Taylor-Galerkin finite element method is used to capture the boundary layers. The method is third-order accurate and has inbuilt upwinding. Stability analysis has been carried out and the numerical results show that the method is efficient in capturing the boundary layers.  相似文献   

10.
The present paper is devoted to the development of a new scheme to solve the one-dimensional time-dependent Burgers' equation locally on sub-domains, using similarity reductions for partial differential equations. Each sub-domain is divided into three grid points. The ordinary differential equation deduced from the similarity reduction can be integrated and is then used to approximate the flux vector in the Burgers' equation. The arbitrary constants in the analytical solution of the similarity equation can be determined in terms of the dependent variables at the grid points in each sub-domain. This approach eliminates the difficulties associated with boundary conditions for the similarity reductions over the whole solution domain. Numerical results are obtained for two different test cases and are compared with other numerical results.  相似文献   

11.
In topology optimization, elements without any contribution to the improvement of the objective function vanish by decrease of density of the design parameter. This easily causes a singular stiffness matrix. To avoid the numerical breakdown caused by this singularity, conventional optimization techniques employ additional procedures. These additional procedures, however, raise some problems. On the other hand, convergence of Krylov subspace methods for singular systems have been studied recently. Through subsequent studies, it has been revealed that the conjugate gradient method (CGM) does not converge to the local optimal solution in some singular systems but in those satisfying certain condition, while the conjugate residual method (CRM) yields converged solutions in any singular systems. In this article, we show that a local optimal solution for topology optimization is obtained by using the CRM and the CGM as a solver of the equilibrium equation in the structural analysis, even if the stiffness matrix becomes singular. Moreover, we prove that the CGM, without any additional procedures, realizes convergence to a local optimal solution in that case. Computer simulation shows that the CGM gives almost the same solutions obtained by the CRM in the case of the two-bar truss problem.  相似文献   

12.
An innovative application focused on the segmentation of decay zones from images of stone materials is presented. The adopted numerical approach to extract decay regions from the color images of monuments provides a tool that helps experts analyze degraded regions by contouring them. In this way even if the results of the proposed procedure depend on the evaluation of experts, the approach can be a contribution to improving the efficiency of the boundary detection process. The segmentation is a process that allows an image to be divided into disjoint zones so that partitioned zones contain homogeneous characteristics. The numerical method, used to segment color images, is based on the theory of interface evolution, which is described by the eikonal equation. We adopted the fast marching technique to solve the upwind finite difference approximation of the eikonal equation. The fast marching starts from a seed point in the region of interest and generates a front which evolves according to a specific speed function until the boundary of the region is identified. We describe the segmentation results obtained with two speed functions, attained by the image gradient computation and color information about the object of interest. Moreover, we present the extension of the working modality of the method by introducing the possibility to extract the regions not only in a local way but also in a global way on the entire image. In this case, in order to improve the segmentation efficiency the application of the fast marching technique starts with more seed points defined as seed regions. The study case concerns the impressive remains of the Roman Theatre in the city of Aosta (Italy). In the image segmentation process the color space LabLab is utilized.  相似文献   

13.
This paper presents selected approximation techniques, typical for the meshless finite difference method (MFDM), although applied to the finite element method (FEM). Finite elements with standard or hierarchical shape functions are coupled with higher order meshless schemes, based upon the correction terms of a simple difference operator. Those terms consist of higher order derivatives, which are evaluated by means of the appropriate formulas composition as well as a numerical solution, which corresponds to the primary interpolation order, assigned to element shape functions. Correction terms modify the right-hand sides of algebraic FE equations only, yielding an iterative procedure. Therefore, neither re-generation of the stiffness matrix nor introduction of any additional nodes and/or degrees of freedom is required. Such improved FE-MFD solution approach allows for the optimal application of advantages of both methods, for instance, a high accuracy of the nodal FE solution and a derivatives’ super-convergence phenomenon at arbitrary domain points, typical for the meshless FDM. Existing and proposed higher order techniques, applied in the FEM, are compared with each other in terms of the solution accuracy, algorithm efficiency and computational complexity.In order to examine the considered algorithms, numerical results of several two-dimensional benchmark elliptic problems are presented. Both the accuracy of a solution and the solution’s derivatives as well as their convergence rates, evaluated on irregular and structured meshes as well as arbitrarily irregular adaptive clouds of nodes, are taken into account.  相似文献   

14.
通过将原方程变换为对流扩散方程,将所得方程的对流项采用四阶组合紧致迎风格式离散,扩散项采用四阶对称紧致格式离散之后,对得到的空间半离散格式采用四阶龙格库塔方法进行时间推进,得到了一种求解非定常对流扩散反应问题的高精度方法,其收敛阶为O(h4+τ4).经数值实验并与文献结果进行对比,表明该格式适用于对流占优问题的数值模拟,验证了格式的良好性能.  相似文献   

15.
In this paper we analyze and implement a second-order-in-time numerical scheme for the three-dimensional phase field crystal (PFC) equation. The numerical scheme was proposed in Hu et al. (2009), with the unique solvability and unconditional energy stability established. However, its convergence analysis remains open. We present a detailed convergence analysis in this article, in which the maximum norm estimate of the numerical solution over grid points plays an essential role. Moreover, we outline the detailed multigrid method to solve the highly nonlinear numerical scheme over a cubic domain, and various three-dimensional numerical results are presented, including the numerical convergence test, complexity test of the multigrid solver and the polycrystal growth simulation.  相似文献   

16.
主要研究了如何使用迭代法求解阶差分方程,利用Java小程序内嵌网页,实现了多种激励下的输出序列(包括零输入响应、零状态响应、全响应)波形的动态绘制。界面交互性强,动态性好,丰富了差分方程的教学,对信号与系统课程的网上教学作了很好的尝试。  相似文献   

17.
In this paper we have discussed a general method of deriving high order multilevel schemes for the heat conduction equation in higher dimensions. The application of the three level parabolic schemes for the solution of the steady state Dirichlet problem in two and three dimensions is discussed by introducing two parameters. The error bounds given by Hadjidimos are further refined which drastically reduces the theoretical estimate of the number of iteration cycles required for a given accuracy η. We have also defined a new set of iteration parameters. A numerical example is computed and it is seen that the present methods produce accurate results.  相似文献   

18.
Spectral preconditioners are based on the fact that the convergence rate of the Krylov subspace methods is improved if the eigenvalues of the smallest magnitude of the system matrix are ‘removed’. In this paper, two preconditioning strategies are studied to solve a set of linear systems associated with the numerical integration of the time-dependent neutron diffusion equation. Both strategies can be implemented using the matrix–vector product as the main operation and succeed at reducing the total number of iterations needed to solve the set of systems.  相似文献   

19.
Fractional partial differential equations (PDEs) provide a powerful and flexible tool for modeling challenging phenomena including anomalous diffusion processes and long-range spatial interactions, which cannot be modeled accurately by classical second-order diffusion equations. However, numerical methods for space-fractional PDEs usually generate dense or full stiffness matrices, for which a direct solver requires O(N3) computations per time step and O(N2) memory, where N is the number of unknowns. The significant computational work and memory requirement of the numerical methods makes a realistic numerical modeling of three-dimensional space-fractional diffusion equations computationally intractable.Fast numerical methods were previously developed for space-fractional PDEs on multidimensional rectangular domains, without resorting to lossy compression, but rather, via the exploration of the tensor-product form of the Toeplitz-like decompositions of the stiffness matrices. In this paper we develop a fast finite difference method for distributed-order space-fractional PDEs on a general convex domain in multiple space dimensions. The fast method has an optimal order storage requirement and almost linear computational complexity, without any lossy compression. Numerical experiments show the utility of the method.  相似文献   

20.
The fourth-order compact approximation for the spatial second-derivative and several linearized approaches, including the time-lagging method of Zhang et al. (1995), the local-extrapolation technique of Chang et al. (1999) and the recent scheme of Dahlby et al. (2009), are considered in constructing fourth-order linearized compact difference (FLCD) schemes for generalized NLS equations. By applying a new time-lagging linearized approach, we propose a symmetric fourth-order linearized compact difference (SFLCD) scheme, which is shown to be more robust in long-time simulations of plane wave, breather, periodic traveling-wave and solitary wave solutions. Numerical experiments suggest that the SFLCD scheme is a little more accurate than some other FLCD schemes and the split-step compact difference scheme of Dehghan and Taleei (2010). Compared with the time-splitting pseudospectral method of Bao et al. (2003), our SFLCD method is more suitable for oscillating solutions or the problems with a rapidly varying potential.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号