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.

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.

PREDICTION OF ACOUSTIC PERFORMANCE IN EXPANSION CHAMBER MUFFLERS WITH MEAN FLOW BY FINITE ELEMENT METHOD

The equation of wave propagation in a circular chamber with mean flow is obtained. Computational solution based on finite element method is employed to determine the transmission loss of expansive chamber. The effect of the mean flow and geometry (length of expansion chamber and expansion ratio)on acoustic attenuation performance is discussed, the predicted values of transmission loss of expansion chamber without and with mean flow are compared with those reported in the literature and they agree well. The accuracy of the prediction of transmission loss implies that finite element approximations are applicable to a lot of practical applications.

SPACE-TIME FINITE ELEMENT METHOD FOR SCHRDINGER EQUATION AND ITS CONSERVATION

Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.

Continuous finite element methods for Hamiltonian systems

By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems, and they both keep energy conservative. The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems. The numerical results are in agree-ment with theory.

Finite Element Methods for Coupled Stokes and Darcy Problems

We derived and analyzed a new numerical scheme for the coupled Stokes and Darcy problems by using H（div） conforming elements in the entire domain. The approach employs the mixed finite element method for the Darcy equations and a stabilized H（div） finite element method for the Stokes equations. Optimal error estimates for the fluid velocity and pressure are derived. The finite element solutions from the new scheme not only feature a full satisfaction of the continuity equation, which is highly demanded in scientific computing, but also satisfy the mass conservation.

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.

A Posteriori Error Estimates for Finite Element Methods for Systems of Nonlinear,Dispersive Equations

The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.

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.

Research on Analytical Method of Fatigue Characteristics of Soft Yoke Mooring System Based on Full-Scale Measurement

By focusing on the vulnerability of the structure of marine equipments, together with considering the randomness of meta-ocean load in statistics, a kind of analytical method of fatigue characteristics of marine structure based on full- scale and actual measurement of prototype is proposed. On the basis of short-term field measurement results of structural response, research is carried out on the fatigue analysis of hinge joints at the upper part of the Soft Yoke single point Mooring System （SYMS） by simultaneously monitoring the environmental load and considering the design criteria of offshore structure. Through analysis of finite element modeling, the time-histories of typical stress response are obtained, and then the assessment of fatigue damage at key hinge joints is conducted. The simulation results indicate that the proposed method can accurately analyze the fatigue damage of offshore engineering structure caused by the effect of wave load. The present analytical method of fatigue characteristics can be extended on other offshore engineering structures subjected to meta-ocean load.

Characteristics of permanent magnet linear synchronous motor fed by spwm inverter based on field-circuit coupled method

Presented field-circuit coupled adaptive time-stepping finite element method to study on permanent magnet linear synchronous motor (PMLSM) characteristics fed by SPWM voltage source inverter.In air-gap field where the direction or magnitude of the field is changing rapidly,the smallest elements are demanded due to high accuracy to use adaptive meshing technique.The co-simulation was used with the status space functions and time-step finite element functions,in which time-step of the status space functions was the smallest than finite element functions'.The magnitude relation of the normal elec- tromagnetic force and tangential electromagnetic force and the period were attained,and current curve was very abrupt at current zero area due to the bigger resistance and leak- age reactance,including main characteristics of motor voltage and velocity.The simulation results compare triumphantly with the experiments results.

A Conservative SAV-RRK Finite Element Method for the Nonlinear Schrodinger Equation

Abstract.In this paper,we propose,analyze and numerically validate a conservative finite element method for the nonlinear Schrodinger equation.A scalar auxiliary variable(SAV)is introduced to reformulate the nonlinear Schrodinger equation into an equivalent system and to transform the energy into a quadratic form.We use the standard continuous finite element method for the spatial discretization,and the relaxation Runge-Kutta method for the time discretization.Both mass and energy conservation laws are shown for the semi-discrete finite element scheme,and also preserved for the full-discrete scheme with suitable relaxation coefficient in the relaxation Runge-Kutta method.Numerical examples are presented to demonstrate the accuracy of the proposed method,and the conservation of mass and energy in long time simulations.

Calculation and Analysis of Acoustic Characteristics of Straight-Through Perforated Pipe Muffler Based on Multilayer Sound Absorbing Material

Using the multi-physical field simulation software COMSOL,the acoustic characteristics of the multilayer sound absorbing material straight-through perforated pipe muffler are studied by the finite element method.The results show that the finite element calculation of the multilayer sound absorbing material straight-through the perforated pipe muffler agrees well with the experimental measurement results.The reliability of the finite element method for studying the acoustic performance of the straight-through perforated pipe muffler with multilayer sound absorbing materials is shown.Furthermore,the influence of some structural parameters of porous sound absorbing material and micro-perforated plate on the acoustic performance of the multilayer sound absorbing material straight-through perforated pipe muffler is analyzed.The muffler based on multilayer sound absorbing material is different from the traditional muffler.After applying the multilayer sound absorbing material to the straight-through perforated pipe muffler,the transmission loss value greatly increases,which provides new ideas and directions for future research on the muffler.

Cell Area and Strut Distribution Changes of Bent Coronary Stents:A Finite Element Analysis

Coronary stents are metal coils or mesh tubes delivered to blocked vessels through catheters, which are expanded by balloons to reopen and scaffold target vessels. Recently,special drugs are carried by stents (drug-eluting stents) to further reduce in-stent restenosis rate after stenting procedure. However,continual study on biomechanical characteristics of stents is necessary for better interactions between stents and tissue, or to provide a more suitable drug loading platform for drug-eluting stents. The purpose of this paper is to show how finite element methods can be used to study cell area and strut distribution changes of bent coronary stents. A same bending deformation was applied to two commercial coronary stent models by a rigid curved vessel. Results show that the stent design influenced the changes of cell area and strut distribution under bending situation. The stent with links had more cell area changes at outer curvature, and the stent with peak-peak (><) strut design could have strut contact and overlapping at inner curvature. In conclusion, this finite element method can be used to study and compare cell area and strut distribution changes of bent stents,and to provide a convenient tool for designers in testing and improving biomechanical characteristics of new stents.

考虑激励幅值慢变的卷绕头转子系统动力学建模

卷绕头转子系统结构之间的耦合振动和时频变激励使振动问题十分复杂,为了更准确地研究其动力学特性,根据转子动力学理论,使用有限元方法建立了考虑多盘质量慢变的卷绕头非线性转子系统整机动力学模型,运用Newmark-β数值积分方法对其动力学特性求解,得到仅存在初始不平衡量以及卷装不良带来的激励频率、幅值慢变下的非线性动力学行为。通过分析轴心轨迹、时域波形和幅频特性曲线等结果,发现卷绕头转子系统对因卷装不良引起的激励幅值慢变较为敏感,出现新的频率成分。进一步对比分析了不同节点处的振动频谱,其结果可为卷绕头转子系统动力学设计和故障识别提供一定理论支持。

脉冲风洞天平基础动力特征的简化计算方法

为了研究脉冲风洞天平基础在脉冲气动荷载作用下的动力响应特征,以某脉冲风洞为例,选择某典型气动荷载作用形式,建立了脉冲风洞天平基础在气动荷载作用下竖直方向、水平方向位移以及回转角度等动力特征的简化计算方法,并采用数值模拟的方法验证其可靠性.研究结果表明:天平基础在典型气动荷载作用下,竖直向最大振幅0.00175 mm,频率7.94 Hz,水平向最大振幅0.00283 mm,频率7.94 Hz,回转角度最大振幅0.00034°,频率7.94 Hz,典型气动荷载对天平基础振动影响较小,也未产生共振现象;基础振动最大振幅随气动荷载增大而增大,基础振动频率随气动荷载频率增大而增大;在气动荷载不变的条件下,基础振动的最大振幅与频率随尺寸的变大而变小,基础振动的最大振幅也随地基土性质的增强而减小,但地基土性质的变化对基础振动频率无影响.

基于功率流有限元法的涡旋压缩机振动特性研究

针对涡旋压缩机内部振动传递路径复杂,采用传统的传递函数法往往不能真实地描述振动产生与传递的本质问题,引入功率流有限元法从能量角度分析振动传递特性,进行涡旋压缩机内部振动产生、传递与空间分布的可视化研究。在分析涡旋压缩机的动力拓扑关系基础上建立激励源,提取与振动响应模型并进行仿真分析;利用实体试验和响应分析结果对比进行模型有效性分析;在验证后的模型上定义8个能量输出面,分析谐波与激励作用下上支撑和下支撑的振动能量传递特性,并且计算节点的功率流矢量确定振源的空间分布。结果表明,涡旋压缩机曲轴在运转到280°时,振动频率在1 000 Hz与1 890 Hz时振动能量传递率较高,上支撑能量分布均匀并集于自身,下支撑能量集中于支脚并向外传递。研究方法和研究成果对于涡旋压缩机振动分布及传递、振动控制和优化设计等方面具有一定借鉴意义。