A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element m...A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction-diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction-diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the high-gradient boundary layers.展开更多
In this paper, a semi-discrete defect-correction mixed finite element method (MFEM) for solving the non-stationary conduction-convection problems in two dimension is presented. In this method, we solve the nonlinear e...In this paper, a semi-discrete defect-correction mixed finite element method (MFEM) for solving the non-stationary conduction-convection problems in two dimension is presented. In this method, we solve the nonlinear equations with an added artificial viscosity term on a finite element grid and correct this solutions on the same grid using a linearized defect-correction technique. The stability and the error analysis are derived. The theory analysis shows that our method is stable and has a good convergence property.展开更多
Under the' assumption of linearization of the free-surface condition, making use of Green's function method and the convolution theorem, analytic solutions of perturbation velocity potentials which correspond ...Under the' assumption of linearization of the free-surface condition, making use of Green's function method and the convolution theorem, analytic solutions of perturbation velocity potentials which correspond to three dimensional unsteady thickness problem and lifting problem caused respectively by arbitrary motions of a body and a hydrofoil beneath the water surface can be achieved in the closed form, In general, the whole perturbation velocity potential consists of three terms, namely φ=φ1+φ2+φ3 , where φ1 denotes the induced velocity potential of the surface singularity distribution in an unbounded fluid, φ2 denotes its mirror image and φ3 denotes that of wave formation which includes the memory effect of the action of the singularity distribution. Utilizing the polynomial expansion of sin[(t-τ)] , the similarity between φ2 and φ3 is discovered and thus a simpler differential relation between them is obtained. Applying this relation, the amount of work in calculation of φ3 which is the most time-consuming one will be reduced significantly. It is favorable not only for dealing with unsteady wave- making problems but also for solving the steady ones in virtue of evading a major difficulty which has to be encountered during the evaluation of an improper inte- gral containing a singularity in the Green's function. The limitation of this new technique turns out to be its slower convergence as the Froude number is lower.展开更多
文摘A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction-diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction-diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the high-gradient boundary layers.
基金supported by National Natural Science Foundation of China (Grant No.10971166)the National Basic Research Program of China (Grant No. 2005CB321703)
文摘In this paper, a semi-discrete defect-correction mixed finite element method (MFEM) for solving the non-stationary conduction-convection problems in two dimension is presented. In this method, we solve the nonlinear equations with an added artificial viscosity term on a finite element grid and correct this solutions on the same grid using a linearized defect-correction technique. The stability and the error analysis are derived. The theory analysis shows that our method is stable and has a good convergence property.
文摘Under the' assumption of linearization of the free-surface condition, making use of Green's function method and the convolution theorem, analytic solutions of perturbation velocity potentials which correspond to three dimensional unsteady thickness problem and lifting problem caused respectively by arbitrary motions of a body and a hydrofoil beneath the water surface can be achieved in the closed form, In general, the whole perturbation velocity potential consists of three terms, namely φ=φ1+φ2+φ3 , where φ1 denotes the induced velocity potential of the surface singularity distribution in an unbounded fluid, φ2 denotes its mirror image and φ3 denotes that of wave formation which includes the memory effect of the action of the singularity distribution. Utilizing the polynomial expansion of sin[(t-τ)] , the similarity between φ2 and φ3 is discovered and thus a simpler differential relation between them is obtained. Applying this relation, the amount of work in calculation of φ3 which is the most time-consuming one will be reduced significantly. It is favorable not only for dealing with unsteady wave- making problems but also for solving the steady ones in virtue of evading a major difficulty which has to be encountered during the evaluation of an improper inte- gral containing a singularity in the Green's function. The limitation of this new technique turns out to be its slower convergence as the Froude number is lower.
