CHARACTERISTIC GALERKIN METHOD FOR CONVECTION-DIFFUSION EQUATIONS AND IMPLICIT ALGORITHM USING PRECISE INTEGRATION

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.

The Characteristic Function Method and Its Application to (1 + 1)-Dimensional Dispersive Long Wave Equation

In this paper, the characteristic function method is applied to seek traveling wave solutions of nonlinear partial differential equations in a unified way. We consider the Wu-Zhang equation (which describes (1 + 1)-dimensional disper-sive long wave). The equations governing the wave propagation consist of a pair of non linear partial differential equations. The characteristic function method reduces the system of nonlinear partial differential equations to a system of nonlinear ordinary differential equations which is solved via the shooting method, coupled with Rungekutta scheme. The results include kink-profile solitary wave solutions, periodic wave solutions and rational solutions. As an illustrative example, the properties of some soliton solutions for Wu-Zhang equation are shown by some figures.

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION (Ⅱ)——THE DISCRETE-TIME CASE ALONG CHARACTERISTICS

The mixed finite element(MFE) methods for a shallow water equation system consisting of water dynamics equations,silt transport equation,and the equation of bottom topography change were derived.A fully discrete MFE scheme for the discrete_time along characteristics is presented and error estimates are established.The existence and convergence of MFE solution of the discrete current velocity,elevation of the bottom topography,thickness of fluid column,and mass rate of sediment is demonstrated.

A Characteristics-Mix Stabilized Finite Element Method for Variable Density Incompressible Navier-Stokes Equations

This paper describes a characteristics-mix finite element method for the computation of incompressible Navi-er-Stokes equations with variable density. We have introduced a mixed scheme which combines a characteristics finite element scheme for treating the mass conservation equation and a finite element method to deal with the momentum equation and the divergence free constraint. The proposed method has a lot of attractive computational properties: parameter-free, very flexible, and averting the difficulties caused by the original equations. The stability of the method is proved. Finally, several numerical experiments are given to show that this method is efficient for variable density incompressible flows problem.

A Second Order Characteristic Mixed Finite Element Method for Convection Diffusion Reaction Equations

A combined approximate scheme is defined for convection-diffusion-reaction equations. This scheme is constructed by two methods. Standard mixed finite element method is used for diffusion term. A second order characteristic finite element method is presented to handle the material derivative term, that is, the time derivative term plus the convection term. The stability is proved and the L2-norm error estimates are derived for both the scalar unknown variable and its flux. The scheme is of second order accuracy in time increment, symmetric, and unconditionally stable.

Two-grid method for characteristic mixed finite-element solutions of nonlinear convection-diffusion equations

A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlinear systems by two steps. The error analysis shows that the two-grid scheme combined with the characteristic mixed finite-element method can decrease numerical oscillation caused by dominated convections and solve nonlinear advection-dominated diffusion problems efficiently.

Matlab Methods for Roots of Characteristic Equation of Cylindrical Bessel Equation Eigenvalue Problems on General Interval

Root of characteristic equation for cylindrical Bessel equation eigenvalue problems on general interval is of great real physical importance at engineering and physical. First, the characteristic equation of cylindrical Bessel equation eigenvalue problem on general interval is given, second, by mean of compared method, we obtaining roots of characteristic equation with Matlab program is discussed.

Theoretical framework for geoacoustic inversion by adjoint method

Traditional geoacoustic inversions are generally solved by matched-field processing in combination with metaheuristic global searching algorithms which usually need massive computations. This paper proposes a new physical framework for geoacoustic retrievals. A parabolic approximation of wave equation with non-local boundary condition is used as the forward propagation model. The expressions of the corresponding tangent linear model and the adjoint operator are derived, respectively, by variational method. The analytical expressions for the gradient of the cost function with respect to the control variables can be formulated by the adjoint operator, which in turn can be used for optimization by the gradient-based method.

EXPLICITLY SOLVABLE CHARACTERISTIC METHOD FOR CONVECTION DIFFUSION PROBLEMS

In this paper, a new method for numerically solving nonlinear convection-dominated diffusion problems is devised and analysed. The discrete time approximations with time stepping along charactcristics are cstablished and solved in spaces posscssing reproducing kernel functions. At each time step, the exact solution of the approximate problem is given by explicit expression. The computational advantage of this method is that the schemes are absolutely stable, and are explicitly solvable as well. The stability and error estimates are derived. Some numerical results are given.

Highly-Efficient Aerodynamic Optimal Design of Rotor Airfoil Using Viscous Adjoint Method

In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.

Characteristic-based finite volume scheme for 1D Euler equations

In this paper, a high-order finite-volume scheme is presented for the one- dimensional scalar and inviscid Euler conservation laws. The Simpson＇s quadrature rule is used to achieve high-order accuracy in time. To get the point value of the Simpson＇s quadrature, the characteristic theory is used to obtain the positions of the grid points at each sub-time stage along the characteristic curves, and the third-order and fifth-order central weighted essentially non-oscillatory （CWENO） reconstruction is adopted to estimate the cell point values. Several standard one-dimensional examples are used to verify the high-order accuracy, convergence and capability of capturing shock.

TWO-GRID METHOD FOR CHARACTERISTICS FINITE-ELEMENT SOLUTION OF 2D NONLINEAR CONVECTION-DOMINATED DIFFUSION PROBLEM

For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid and the fine-grid solution can be obtained in a single linear step. For the nonlinear convection-dominated diffusion equation, this method can not only stabilize the numerical oscillation but also accelerate the convergence and improve the computational efficiency. The error analysis demonstrates if the mesh sizes between coarse-grid and fine-grid satisfy the certain relationship, the two-grid solution and the characteristics finite-element solution have the same order of accuracy. The numerical is more efficient than that of characteristics example confirms that the two-grid method finite-element method.

ON DIRECT METHOD OF SOLUTION FOR A CLASS OF SINGULAR INTEGRAL EQUATIONS

In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.

A Computational Study with Finite Difference Methods for Second Order Quasilinear Hyperbolic Partial Differential Equations in Two Independent Variables

In this paper we consider the numerical method of characteristics for the numerical solution of initial value problems (IVPs) for quasilinear hyperbolic Partial Differential Equations, as well as the difference scheme Central Time Central Space (CTCS), Crank-Nicolson scheme, ω scheme and the method of characteristics for the numerical solution of initial and boundary value prob-lems for the one-dimension homogeneous wave equation. The initial deriva-tive condition is approximated by different second order difference quotients in order to examine which gives more accurate numerical results. The local truncation error, consistency and stability of the difference schemes CTCS, Crank-Nicolson and ω are also considered.

Optimal Trajectory of Underwater Manipulator Using Adjoint Variable Method for Reducing Drag

In order to decrease the fluid drag on an underwater robot manipulator, an optimal trajectory method based on the variational method is presented. By introducing the adjoint variables, which are Lagrange multipliers, we formulate a Lagrange function under certain constraints related to the target angle, target angular velocity, and dynamic equation of the robot manipulator. The state equation (the partial differentiation of the Lagrange function with respect to the state variables), adjoint equation (the partial differentiation of the Lagrange function with respect to the adjoint variables), and sensitivity equation (the partial differentiation of the Lagrange function with respect to torques) can be derived from the stationary conditions of the Lagrange function. Using the state equation, we can calculate the state variables (angles, angular velocities, and angular acceleration) at every time step in the forward time direction. These state variables are stored as data at every time step. Next, by using the adjoint equation, we can calculate the adjoint variables by using these state variables at every time step in the backward time direction. These adjoint variables are stored as data at every time step. Third, the sensitivity equation is calculated by using both the state variables and the adjoint variables. Finally, the optimal trajectory of the manipulator is obtained using the sensitivities. The proposed method is applied to the problem of two-link manipulators. It can obtain the optimal drag reduction trajectory of the manipulator under the constraints mentioned above.

Hydrodynamic Lubrication of Elastic Foil Gas Bearing Using Over Relaxation Iteration Method and Non-Dimensional Equation

The purpose is to accurately predict the performance of foil bearing and achieve accurate results in the design of foil bearing structure.A new type of foil bearing with surface microstructure is used as experimental material.First,the lubrication mechanism of elastic foil gas bearing is analyzed.Then,the numerical solution process of the static bearing capacity and friction torque is analyzed,including the discretization of the governing equation of rarefied gas pressure based on the non-dimensional modified Reynolds equation and the over relaxation iteration method,the grid planning within the calculation range,the static solution of boundary parameters and static solution of the numerical process.Finally,the solution program is analyzed.The experimental data in National Aeronautics and Space Administration(NASA)public literature are compared with the simulation results of this exploration,so as to judge the accuracy of the calculation process.The results show that under the same static load,the difference between the minimum film thickness calculated and the test results is not obvious;when the rotor speed of the bearing is 60000 r/min,the influence of the boundary slip effect increases with the increase of the micro groove depth on the flat foil surface;when the eccentricity or the micro groove depth of the bearing increases,the bearing capacity will be strengthened.When the eccentricity is 6µm and 14µm,the viscous friction torque of the new foil bearing increases significantly with the increase of the depth of the foil micro groove,but when the eccentricity is 22µm,the viscous friction torque does not change with the change of the depth of the foil micro groove.It shows that the bearing capacity and performance of foil bearing are improved.

Goal-Oriented Anisotropic hp-Adaptive Discontinuous Galerkin Method for the Euler Equations

We deal with the numerical solution of the compressible Euler equations with the aid of the discontinuous Galerkin(DG)method with focus on the goal-oriented error estimates and adaptivity.We analyse the adjoint consistency of the DG scheme where the adjoint problem is not formulated by the differentiation of the DG form and the target functional but using a suitable linearization of the nonlinear forms.Furthermore,we present the goal-oriented anisotropic hp-mesh adaptation method for the Euler equations.The theoretical results are supported by numerical experiments.

一个解特定二阶驻定方程的新方法——特征伴随方程法

本文对文献[1]中提出的一类二阶驻定方程G″/G=∑^(m)_(j=0)P_(j)(G′/G)^(j)的求法进行探讨,提出了全新的求解方法——特征伴随方程法,通过该求法得到这类方程当m取不同非负整数(本文仅讨论m=2)时的通解,并相应给出该方程作为求解非线性偏微分方程的辅助方程时所需要的G′/G的解析表达式;同时给出一个作为常微分方程中驻定方程的应用实例,以及该方程作为非线性偏微分方程的辅助方程时利用扩展G′/G展开法的求解例子,通过该实例给出方程的精确行波解。

不同前倾角度倾转旋翼噪声数值计算分析

基于自由尾迹结合FW-H(Ffowcs Williams-Hawkings)方程的方法建立了倾转旋翼气动噪声计算模型,并采用倾转旋翼模型噪声试验数据验证了计算分析方法的有效性。选取典型过渡路径,进行考虑配平的倾转双旋翼气动噪声特性计算,获得了旋翼桨叶剖面非定常气动载荷以及不同测点气动噪声等计算结果,分析了倾转旋翼在不同前倾角下噪声指向性和噪声声压级的变化。结果表明:由于双旋翼噪声在传播中的叠加和抵消,倾转双旋翼和孤立单旋翼的噪声指向特性存在较大的不同;倾转旋翼噪声随前倾角增加总体上呈现先增加后减小的变化趋势,在前倾角30°附近噪声最强;不同前倾角下噪声声压级和指向性的变化与旋翼桨尖马赫数、气动载荷和桨盘角度等多种因素相关。