THE COUPLING OF NATURAL BOUNDARY ELEMENT AND FINITE ELEMENT METHOD FOR 2D HYPERBOLIC EQUATIONS

In this paper, we investigate the coupling of natural boundary element and finite element methods of exterior initial boundary value problems for hyperbolic equations. The governing equation is first discretized in time, leading to a time-step scheme, where an exterior elliptic problem has to be solved in each time step. Second, a circular artificial boundary FR consisting of a circle of radius R is introduced, the original problem in an unbounded domain is transformed into the nonlocal boundary value problem in abounded subdomain. And the natural integral equation and the Poisson integral formula are obtained in the infinite domainΩ2 outside circle of radius R. The coupled variational formulation is given. Only the function itself, not its normal derivative at artificial boundary ΓR, appears in the variational equation, so that the unknown numbers are reducedand the boundary element stiffness matrix has a few different elements. Such a coupled method is superior to the one based on direct boundary element method. This paper discusses finite element discretization for variational problem and its corresponding numerical technique, and the convergence for the numerical solutions. Finally, the numerical example is presented to illustrate feasibility and efficiency of this method.     

藻类植物中钙、镁、铁、锰、铜和锌含量分析

本文采用火焰原子吸收光谱法测定了七种藻类植物中钙 ,镁 ,铁 ,锰 ,铜 ,锌六种元素的含量。结果表明 ,藻类植物中含有较丰富的钙 ,镁 ,铁 ,锰 ,锌元素 ,其含量由大到小顺序为Ca >Mg >Fe >Zn >Mn >Cu ,尤以钙的含量最为丰富 ,铜含量较少。本实验结果为探讨藻类植物作为资源开发利用提供了有用数据。     

Duality Theory in Fuzzy Linear Programming Problems with Fuzzy Coefficients

The concept of fuzzy scalar (inner) product that will be used in the fuzzy objective and inequality constraints of the fuzzy primal and dual linear programming problems with fuzzy coefficients is proposed in this paper. We also introduce a solution concept that is essentially similar to the notion of Pareto optimal solution in the multiobjective programming problems by imposing a partial ordering on the set of all fuzzy numbers. We then prove the weak and strong duality theorems for fuzzy linear programming problems with fuzzy coefficients.     

UPWIND FINITE VOLUME SCHEMES FOR SEMICONDUCTOR DEVICE

The mathematical model of semiconductor devices is described by the initial boundary value problem of a system of three nonlinear partial differential equations. One equation in elliptic form is for the electrostatic potential; two equations of convection-dominated diffusion type are for the electron and hole concentrations. Finite volume element procedure are put forward for the electrostatic potential, while upwind     

Structural optimization design for large mirror

The system design must be started from the concept with low cost and high performance. In this point, the topological shape of the structure is very important in the view of the structural rigidities and light-weight design.In this paper, the optimization methodology is presented in the design stage of the large optical mirror. We obtain the optimal layout through the topology optimization and then design the details through the size or shape optimization for structural rigidity.     

Mathematical and numerical analysis of a simplified time-dependent viscoelastic flow

A time-dependent model corresponding to an Oldroyd-B viscoelastic fluid is considered, the convective terms being disregarded. Global existence in time is proved in Banach spaces provided the data are small enough, using the implicit function theorem and a maximum regularity property for a three fields Stokes problem. A finite element discretization in space is then proposed. Existence of the numerical solution is proved for small data, so as a priori error estimates, using again an implicit function theorem. Supported by the Swiss National Science Foundation. Fellowship PBEL2–114311.     

Modeling hyperelastic behavior of rubber: A novel invariant‐based and a review of constitutive models

A constitutive phenomenological model completing the Gent‐Thomas concept is carried out to formulate laws governing the hyperelastic behavior of incompressible rubber materials. It is shown that the phenomenological Gent‐Thomas model (1958) and the constrained chain model (1992) give similar precise results at small to moderate deformation. On the other hand, comparisons of the outcome of the proposed model with that of the molecular model from the combined concepts of Flory‐Erman and Boyce‐Arruda (2000), and with those of the phenomenological models of Ogden (1982), Yeoh‐Fleming (1997), Pucci‐Saccomandi (2002) and Beda (2005) are made. Residual inconveniences raised by attractive continuum models in rubber elasticity literature have been successfully overcome. Results from both the statistical and phenomenological mechanics concepts are compared with the data of some useful classical materials (rubbers of Treloar, Rivlin‐Saunders, Pak‐Flory and Yeoh‐Fleming). The results permit one to see salient equivalence of the two theories for a more reliable prediction of stress‐stretch response for all states of any mode of deformation. A complete and exhaustive analysis of the Mooney plot that combines small and very large extension‐compression has been quite essential in assessing the validity of models. A method of identification of material parameters is presented and data of the simple tension suffice for the determination of the parameter values. It is shown that the ordinary identification procedures, such as the usual least squares, a very much used numerical method in materials investigation, can be unsuitable in some cases of hyperelastic modeling. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1713–1732, 2007     

Finite element computations of two-dimensional arterial flow in the presence of a transverse magnetic field

A finite element solution of the Navier-Stokes equations for steady flow under the magnetic effect through a double-branched two-dimensional section of a three-dimensional model of the canine aorta is discussed. The numerical scheme involves transforming the physical co-ordinates to a curvilinear boundary-fitted co-ordinate system. The shear stress at the wall is calculated for a Reynolds number of 1000 with the branch-to-main aortic flow rate ratio as a parameter. The results are compared with earlier works involving experimental data and found to be in reasonable qualitative agreement. The steady flow, shear stress and branch flow under the effect of a magnetic field have been discussed in detail.     

Speed of adaptation in structured populations

Adaptation of populations takes place with the occurrence and subsequent fixation of mutations that confer some selective advantage to the individuals which acquire it. For this reason, the study of the process of fixation of advantageous mutations has a long history in the population genetics literature. Particularly, the previous investigations aimed to find out the main evolutionary forces affecting the strength of natural selection in the populations. In the current work, we investigate the dynamics of fixation of beneficial mutations in a subdivided population. The subpopulations (demes) can exchange migrants among their neighbors, in a migration network which is assumed to have either a random graph or a scale-free topology. We have observed that the migration rate drastically affects the dynamics of mutation fixation, despite of the fact that the probability of fixation is invariant on the migration rate, accordingly to Maruyama's conjecture. In addition, we have noticed a topological dependence of the adaptive evolution of the population when clonal interference becomes effective.     

热传导型半导体瞬态问题特征变网格有限元法及其分析

热传导型半导体器件的瞬时状态由四个方程的非线性偏微分方程组的初边值问题所决定，其中电子位势方程是椭圆型的，电子和空穴浓度方程是对流扩散型的，温度方程为热传导型的。本文提出解这类问题的特征变网格有限元法，并进行了理论分析，在一定条件下，得到了某种意义下的最佳L^2误差估计结果。