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.

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.

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.

Analysis of Frequency Characteristics and Sensitivity of Compliant Mechanisms

Based on a modified pseudo-rigid-body model,the frequency characteristics and sensitivity of the large-deformation compliant mechanism are studied.Firstly,the pseudo-rigid-body model under the static and kinetic conditions is modified to enable the modified pseudo-rigid-body model to be more suitable for the dynamic analysis of the compliant mechanism.Subsequently,based on the modified pseudo-rigid-body model,the dynamic equations of the ordinary compliant four-bar mechanism are established using the analytical mechanics.Finally,in combination with the finite element analysis software ANSYS,the frequency characteristics and sensitivity of the compliant mechanism are analyzed by taking the compliant parallel-guiding mechanism and the compliant bistable mechanism as examples.From the simulation results,the dynamic characteristics of compliant mechanism are relatively sensitive to the structure size,section parameter,and characteristic parameter of material on mechanisms.The results could provide great theoretical significance and application values for the structural optimization of compliant mechanisms,the improvement of their dynamic properties and the expansion of their application range.

THE MODIFIED METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OPERATOR-SPLITTING PROCEDURES FOR COMPRESSIBLE MULTICOMPONENT DISPLACEMENT PROBLEM

For the three-dimensional compressible multicomponent displacement problem we put forward the modified method of characteristics with finite element operator-splitting procedures and make use of operator-splitting,characteristic method,calculus of variations,energy method,negative norm estimate,two kinds of test functions and the theory of prior estimates and techniques.Optimal order estimates in L^2 norm are derived for the error in the approximate solution.These methods have been successfully used in oil-gas resources estimation,enhanced oil recovery simulation and seawater intrusion numerical simulation.

A SPLIT-CHARACTERISTIC FINITE ELEMENT MODEL FOR 1-D UNSTEADY FLOWS

An efficient and accurate solution algorithm was proposed for 1-D unsteady flow problems widely existing in hydraulic engineering. Based on the split-characteristic finite element method, the numerical model with the Saint-Venant equations of 1-D unsteady flows was established. The assembled f＇mite element equations were solved with the tri-diagonal matrix algorithm. In the semi-implicit and explicit scheme, the critical time step of the method was dependent on the space step and flow velocity, not on the wave celerity. The method was used to eliminate the restriction due to the wave celerity for the computational analysis of unsteady open-channel flows. The model was verified by the experimental data and theoretical solution and also applied to the simulation of the flow in practical river networks. It shows that the numerical method has high efficiency and accuracy and can be used to simulate 1-D steady flows, and unsteady flows with shock waves or flood waves. Compared with other numerical methods, the algorithm of this method is simpler with higher accuracy, less dissipation, higher computation efficiency and less computer storage.

A CHARACTERISTIC FINITE ELEMENT METHOD FOR CONSTRAINED CONVECTION-DIFFUSION-REACTION OPTIMAL CONTROL PROBLEMS

In this paper, we develop a priori error estimates for the solution of constrained convection-diffusion-reaction optimal control problems using a characteristic finite element method. The cost functional of the optimal control problems consists of three parts： The first part is about integration of the state over the whole time interval, the second part refers to final-time state, and the third part is a regularization term about the control. We discretize the state and co-state by piecewise linear continuous functions, while the control is approximated by piecewise constant functions. Pointwise inequality function constraints on the control are considered, and optimal a L2-norm priori error estimates are obtained. Finally, we give two numerical examples to validate the theoretical analysis.

Theory and Application of Characteristic Finite Element Domain Decomposition Procedures for Coupled System of Dynamics of Fluids in Porous Media

For a coupled system of multiplayer dynamics of fluids in porous media, the characteristic finite element domain decomposition procedures applicable to parallel arithmetic are put forward. Techniques such as calculus of variations, domain decomposition, characteristic method, negative norm estimate, energy method and the theory of prior estimates are adopted. Optimal order estimates in L^2 norm are derived for the error in the approximate solution.

A CHARACTERISTIC MIXED FINITE ELEMENT TWO-GRID METHOD FOR COMPRESSIBLE MISCIBLE DISPLACEMENT PROBLEM

A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium.The concentration equation is treated by a mixed finite element method with characteristics(CMFEM)and the pressure equation is treated by a parabolic mixed finite element method(PMFEM).Two-grid algorithm is considered to linearize nonlinear coupled system of two parabolic partial differential equations.Moreover,the L q error estimates are conducted for the pressure,Darcy velocity and concentration variables in the two-grid solutions.Both theoretical analysis and numerical experiments are presented to show that the two-grid algorithm is very effective.

Characteristics of motorized spindle supported by active magnetic bearings

A motorized spindle supported by active magnetic bearings（AMBs） is generally used for ultra-high-speed machining. Iron loss of radial AMB is very great owing to high rotation speed, and it will cause severe thermal deformation. The problem is particularly serious on the occasion of large power application, such as all electric aero-engine. In this study, a prototype motorized spindle supported by five degree-of-freedom AMBs is developed. Homopolar and heteropolar AMBs are independently adopted as radial bearings. The influences of the two types of radial AMBs on the dynamic characteristics of the motorized spindle are comparatively investigated by theoretical analysis, test modal analysis and actual operation of the system. The iron loss of the two types of radial AMBs is analyzed by finite element software and verified through run-down experiments of the system. The results show that the structures of AMB have less influence on the dynamic characteristics of the motorized spindle. However, the homopolar structure can effectively reduce the iron loss of the radial AMB and it is useful for improving the overall performance of the motorized spindle.

The MFE-CFE-GFE Method for the Fully Coupled Quasi-Static Thermo-Poroelastic Problem

In this work,we consider a combined finite element method for fully coupled nonlinear thermo-poroelastic model problems.The mixed finite element(MFE)method is used for the pressure,the characteristics finite element(CFE)method is used for the temperature,and the Galerkin finite element(GFE)method is used for the elastic displacement.The semi-discrete and fully discrete finite element schemes are established and the stability of this method is presented.We derive error estimates for the pressure,temperature and displacement.Several numerical examples are presented to confirm the accuracy of the method.

GALERKIN ALTERNATING-DIRECTION METHODS FOR NONRECTANGULAR REGIONS FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

For the transient behavior of a semiconductor device, the modified method of characteristics with alternating-direction finite element procedures for nonrectangular region is put forward. Some techniques, such as calculus of variations, isoparametric transformation，patch approximation, operator-splitting, characteristic method, symmetrical reflection,energy method, negative norm estimate and a prior estimates and techniques, are employed. In the nonrectangular region case, optimal order estimates in L^2 norm are derived for the error in the approximation solution. Thus the well-known theoretical problem has been thoroughly and completely solved.

A posteriori error estimates for optimal control problems constrained by convection-diffusion equations

We propose a characteristic finite element evolutionary type convection-diffusion optimal control discretization of problems. Non- divergence-free velocity fields and bilateral inequality control constraints are handled. Then some residual type a posteriori error estimates are analyzed for the approximations of the control, the state, and the adjoint state. Based on the derived error estimators, we use them as error indicators in developing efficient multi-set adaptive meshes characteristic finite element algorithm for such optimal control problems. Finally, one numerical example is given to check the feasibility and validity of multi-set adaptive meshes refinements.

The Alternating-Direction Schemes and Numerical Analysis for the Three-dimensional Seawater Intrusion Simulation

Both numerical simulation and theoretical analysis of seawater intrusion in coastal regions are of great theoretical importance in environmental sciences. The mathematical model can be described as a coupled system of three dimensional nonlinear partial differential equations with initial-boundary value problems. In this paper, according to the actual conditions of molecular and three-dimensional characteristic of the problem, we construct the characteristic finite element alternating-direction schemes which can be divided into three continuous one-dimensional problems. By making use of tensor product algorithm, and priori estimation theory and techniques, the optimal order estimates in H1 norm are derived for the error in the approximate solution.