Fractional Steps Scheme to Approximate the Phase-Field Transition System with Non-Homogeneous Cauchy–Neumann Boundary Conditions

The paper presents a proof of the convergence for an iterative scheme of fractional steps type associated to the phase-field transition system (a nonlinear parabolic system) with non-homogeneous Cauchy–Neumann boundary conditions. The advantage of such method consists in simplifying the numerical computation necessary to be done in order to approximate the solution of a nonlinear parabolic system. On the basis of this approach, a numerical algorithm in the two dimensional case is introduced and an industrial implementation is made.     

On the multiple hyperbolic systems modelling phase transformation kinetics

We discuss Cahn’s time cone method modelling phase transformation kinetics. The model equation by the time cone method is an integral equation in the space-time region. First, we reduce it to a system of hyperbolic equations, and in the case of odd spatial dimensions, the reduced system is a multiple hyperbolic equation. Next, we propose a numerical method for such a hyperbolic system. By means of alternating direction implicit methods, numerical simulations for practical forward problems are implemented with satisfactory accuracy and efficiency. In particular, in the three dimensional case, our numerical method on the basis of reduced multiple hyperbolic equation is fast.     

Linear B-spline finite element method for the improved Boussinesq equation

In this paper, we develop and validate a numerical procedure for solving a class of initial boundary value problems for the improved Boussinesq equation. The finite element method with linear B-spline basis functions is used to discretize the nonlinear partial differential equation in space and derive a second order system involving only ordinary derivatives. It is shown that the coefficient matrix for the second order term in this system is invertible. Consequently, for the first time, the initial boundary value problem can be reduced to an explicit initial value problem to which many accurate numerical methods are readily applicable. Various examples are presented to validate this technique and demonstrate its capacity to simulate wave splitting, wave interaction and blow-up behavior.     

Well-posedness and finite element approximation of time dependent generalized bioconvective flow

We consider the finite element approximation of a time dependent generalized bioconvective flow. The underlying system of partial differential equations consists of incompressible Navier–Stokes type convection equations coupled with an equation describing the transport of micro-organisms. The viscosity of the fluid is assumed to be a function of the concentration of the micro-organisms. We show the existence and uniqueness of the weak solution of the system in two dimensions and construct numerical approximations based on the finite element method, for which we obtain error estimates. In addition, we conduct several numerical experiments to demonstrate the accuracy of the numerical method and perform simulations of the bioconvection pattern formations based on realistic model parameters to demonstrate the validity of the proposed numerical algorithm.     

三维数值流形方法的理论研究

在二维数值流形方法的基础上,对三维数值流形进行了理论研究．研究了三维覆盖位移函数,进行了三维数值流形的力学分析,给出了三维流形单元的刚度矩阵,详细推导了三维数值流形的Hammer积分及剖分规则,系统地研究三维数值流形的理论体系与数值实现方法．作为数值算例,给出了相应的悬臂梁的计算结果,计算结果表明算法的精度和计算效益较高．     

An optimization‐based domain decomposition method for numerical simulation of the incompressible Navier‐Stokes flows

This article is concerned about an optimization‐based domain decomposition method for numerical simulation of the incompressible Navier‐Stokes flows. Using the method, an classical domain decomposition problem is transformed into a constrained minimization problem for which the objective functional is chosen to measure the jump in the dependent variables across the common interfaces between subdomains. The Lagrange multiplier rule is used to transform the constrained optimization problem into an unconstrained one and that rule is applied to derive an optimality system from which optimal solutions may be obtained. The optimality system is also derived using “sensitivity” derivatives instead of the Lagrange multiplier rule. We consider a gradient‐type approach to the solution of domain decomposition problem. The results of some numerical experiments are presented to demonstrate the feasibility and applicability of the algorithm developed in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Numerical solution of nonlinear two-dimensional Fredholm integral equations of the second kind using Gauss product quadrature rules

The Gauss product quadrature rules and collocation method are applied to reduce the second-kind nonlinear two-dimensional Fredholm integral equations (FIE) to a nonlinear system of equations. The convergence of the proposed numerical method is proved under certain conditions on the kernel of the integral equation. An iterative method for approximating the solution of the obtained nonlinear system is provided and its convergence is proved. Also, some numerical examples are presented to show the efficiency and accuracy of the proposed method.     

Theoretical analysis of high speed spindle air bearings by a hybrid numerical method

To study the behavior of the high speed spindle air bearing (HSSAB) system, we conduct the research by means of a hybrid numerical method which combines the differential transformation method and the finite difference method in this paper. According to the results of the research, the flexible rotor center is found to include a complex dynamic behavior that comprises periodic, sub-harmonic and quasi-periodic responses. In addition, as the rotor mass and the bearing number are increased, there will be some changes taking place in the dynamic behavior of the bearing system. The results are proven to have no conflict with those of the other numerical methods, which enables an effective means in gaining insights into the nonlinear dynamics of HSSAB systems.     

Newton-like method with modification of the right-hand-side vector

This paper proposes a new Newton-like method which defines new iterates using a linear system with the same coefficient matrix in each iterate, while the correction is performed on the right-hand-side vector of the Newton system. In this way a method is obtained which is less costly than the Newton method and faster than the fixed Newton method. Local convergence is proved for nonsingular systems. The influence of the relaxation parameter is analyzed and explicit formulae for the selection of an optimal parameter are presented. Relevant numerical examples are used to demonstrate the advantages of the proposed method.

     

A new filled function method for an unconstrained nonlinear equation

In this paper, we transform an unconstrained system of nonlinear equations into a special optimization problem. A new filled function is constructed by employing the special properties of the transformed optimization problem. Theoretical and numerical properties of the proposed filled function are investigated and a solution of the algorithm is proposed. Under some conditions, we can find a solution or an approximate solution to the system of nonlinear equations in finite iterations. The implementation of the algorithm on six test problems is reported with satisfactory numerical results.     

A linear,symmetric and energy-conservative scheme for the space-fractional Klein–Gordon–Schrödinger equations

In this paper, we propose an efficient numerical scheme for the space-fractional Klein–Gordon–Schrödinger (SFKGS) equations. Motivated by the “Invariant Energy Quadratization” (IEQ) approach, we introduce two auxiliary variables to transform the SFKGS system into a new equivalent system in which the time derivative is discretized by the Crank–Nicolson method, and the space discretization is based on the Fourier spectral method. Consequently, the numerical scheme shares two good features. The first feature is that the nonlinear terms are treated semi-explicitly and a linear symmetric system is solved at each time step. The second feature is the energy conservation at the discrete level. These two advantages are proved by the theoretical analysis and illustrated by a given numerical example.     

Modularization modeling and simulation of turbine test rig main test system

Comprehensive applications of modularization modeling method have proven its effectiveness and versatility in system simulation field. This paper establishes the modularization numerical model of a turbine test rig main test system by using a finite volume numerical system developed. The simulation study based on an experiment is conducted. The comparison with available experimental data indicates that the general trends of simulation curves are in agreement with test curves and that there is obvious thermal stratification phenomenon at different positions along combustion gas flow direction. Accordingly, it can be concluded that the analysis of experimental data is reasonable and the established numerical system is effective. It is also found that the modeling of valve spool throttling and the modeling of components-wall heat transfer are two key factors of affecting simulation accuracy.     

A Dynamical Method for Solving the Obstacle Problem

In this paper, we consider the unilateral obstacle problem, trying to find the numerical solution and coincidence set. We construct an equivalent format of the original problem and propose a method with a second-order in time dissipative system for solving the equivalent format. Several numerical examples are given to illustrate the effectiveness and stability of the proposed algorithm. Convergence speed comparisons with existent numerical algorithm are also provided and our algorithm is fast.     

Auxiliary Equations Approach for the Stochastic Unsteady Navier-Stokes Equations with Additive Random Noise

This paper presents a Martingale regularization method for the stochastic Navier-Stokes equations with additive noise. The original system is split into two equivalent parts, the linear stochastic Stokes equations with Martingale solution and the stochastic modified Navier-Stokes equations with relatively-higher regularities. Meanwhile, a fractional Laplace operator is introduced to regularize the noise term. The stability and convergence of numerical scheme for the pathwise modified Navier-Stokes equations are proved.The comparisons of non-regularized and regularized noises for the Navier-Stokes system are numerically presented to further demonstrate the efficiency of our numerical scheme.     

带损伤弹性反问题的数值分析

考虑一类由椭圆性方程和热传导方程共同来刻画的准静态弹性模型,通过给定观测值来反演边界的牵引力.首先构造一个凸目标泛函,并引入Tikhonov正则化方法,使之极小化得到一个稳定的近似解.再用有限元离散求解,导出误差估计.最后,用数值例子说明算法的可行性和有效性.     

Laguerre polynomial approach for solving linear delay difference equations

This paper presents a numerical method for the approximate solution of mth-order linear delay difference equations with variable coefficients under the mixed conditions in terms of Laguerre polynomials. The aim of this article is to present an efficient numerical procedure for solving mth-order linear delay difference equations with variable coefficients. Our method depends mainly on a Laguerre series expansion approach. This method transforms linear delay difference equations and the given conditions into matrix equation which corresponds to a system of linear algebraic equation. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments and performed on the computer algebraic system Maple.     

Some iterative methods for the solution of a symmetric indefinite KKT system

This paper is concerned with the numerical solution of a Karush–Kuhn–Tucker system. Such symmetric indefinite system arises when we solve a nonlinear programming problem by an Interior-Point (IP) approach. In this framework, we discuss the effectiveness of two inner iterative solvers: the method of multipliers and the preconditioned conjugate gradient method. We discuss the implementation details of these algorithms in an IP scheme and we report the results of a numerical comparison on a set of large scale test-problems arising from the discretization of elliptic control problems. This research was supported by the Italian Ministry for Education, University and Research (MIUR), FIRB Project RBAU01JYPN.     

The Runge‐Kutta DG finite element method and the KFVS scheme for compressible flow simulations

This article presents a new type of second‐order scheme for solving the system of Euler equations, which combines the Runge‐Kutta discontinuous Galerkin (DG) finite element method and the kinetic flux vector splitting (KFVS) scheme. We first discretize the Euler equations in space with the DG method and then the resulting system from the method‐of‐lines approach will be discretized using a Runge‐Kutta method. Finally, a second‐order KFVS method is used to construct the numerical flux. The proposed scheme preserves the main advantages of the DG finite element method including its flexibility in handling irregular solution domains and in parallelization. The efficiency and effectiveness of the proposed method are illustrated by several numerical examples in one‐ and two‐dimensional spaces. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A new approach to grid generation

An entirely new approach to numerical grid generation, the deformation method, is presented. Each point of an existing grid is moved, in an inter-related manner, to a new position according to a system of n ordinary differential equations (n=spacial dimension). The resulting grid has prescribed mesh sizes.     

Numerical study of homotopy‐perturbation method applied to Burgers equation in fluid

In this article, the problem of Burgers equation is presented and the homotopy perturbation method (HPM) is employed to compute an approximation to the solution of the system of nonlinear differential equations governing on the problem. Comparison is made between the HPM and Exact solutions. The obtained solutions, in comparison with the exact solutions, admit a remarkable accuracy. A clear conclusion can be drawn from the numerical results that the HPM provides highly accurate numerical solutions for nonlinear differential equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010