An operator-splitting Galerkin/SUPG finite element method for population balance equations : stability and convergence

We present a heterogeneous finite element method for the solution of a high-dimensional population balance equation, which depends both the physical and the internal property coordinates. The proposed scheme tackles the two main difficulties in the finite element solution of population balance equation: (i) spatial discretization with the standard finite elements, when the dimension of the equation is more than three, (ii) spurious oscillations in the solution induced by standard Galerkin approximation due to pure advection in the internal property coordinates. The key idea is to split the high-dimensional population balance equation into two low-dimensional equations, and discretize the low-dimensional equations separately. In the proposed splitting scheme, the shape of the physical domain can be arbitrary, and different discretizations can be applied to the low-dimensional equations. In particular, we discretize the physical and internal spaces with the standard Galerkin and Streamline Upwind Petrov Galerkin (SUPG) finite elements, respectively. The stability and error estimates of the Galerkin/SUPG finite element discretization of the population balance equation are derived. It is shown that a slightly more regularity, i.e. the mixed partial derivatives of the solution has to be bounded, is necessary for the optimal order of convergence. Numerical results are presented to support the analysis.

Modified Newton, rank-1 Broyden update and rank-2 BFGS update methods in helicopter trim: A comparative study

Nonlinear equations in mathematical physics and engineering are solved by linearizing the equations and forming various iterative procedures, then executing the numerical simulation. For strongly nonlinear problems, the solution obtained in the iterative process can diverge due to numerical instability. As a result, the application of numerical simulation for strongly nonlinear problems is limited. Helicopter aeroelasticity involves the solution of systems of nonlinear equations in a computationally expensive environment. Reliable solution methods which do not need Jacobian calculation at each iteration are needed for this problem. In this paper, a comparative study is done by incorporating different methods for solving the nonlinear equations in helicopter trim. Three different methods based on calculating the Jacobian at the initial guess are investigated. (C) 2011 Elsevier Masson SAS. All rights reserved.

Effect of leachate recirculation and extent of degradation on the stability of bioreactor landfill slopes

The stability of a bioreactor landfill slope is influenced by the quantity and method of leachate recirculation as well as on the degree of decomposition. Other factors include properties variation of waste material and geometrical configurations, i.e., height and slope of landfills. Conventionally, the stability of slopes is evaluated using factor of safety approach, in which the variability in the engineering properties of MSW is not considered directly and stability issues are resolved from past experiences and good engineering judgments. On the other hand, probabilistic approach considers variability in mathematical framework and provides stability in a rational manner that helps in decision making. The objective of the present study is to perform a parametric study on the stability of a bioreactor landfill slope in probabilistic framework considering important influencing factors, such as, variation in MSW properties, amount of leachate recirculation, and age of degradation, in a systematic manner. The results are discussed in the light of existing relevant regulations, design and operation issues.

Role of different factors affecting interdiffusion in Cu(Ga) and Cu(Si) solid solutions

A detailed diffusion study was carried out on Cu(Ga) and Cu(Si) solid solutions in order to assess the role of different factors in the behaviour of the diffusing components. The faster diffusing species in the two systems, interdiffusion, intrinsic and impurity diffusion coefficients, are determined to facilitate the discussion. It was found that Cu was more mobile in the Cu-Si system, whereas Ga was the faster diffusing species in the Cu-Ga system. In both systems, the interdiffusion coefficients increased with increasing amount of solute (e.g. Si or Ga) in the matrix (Cu). Impurity diffusion coefficients for Si and Ga in Cu, found out by extrapolating interdiffusion coefficient data to zero composition of the solute, were both higher than the Cu tracer diffusion coefficient. These observed trends in diffusion behaviour could be rationalized by considering: (i) formation energies and concentration of vacancies, (ii) elastic moduli (indicating bond strengths) of the elements and (iii) the interaction parameters and the related thermodynamic factors. In summary, we have shown here that all the factors introduced in this paper should be considered simultaneously to understand interdiffusion in solid solutions. Otherwise, some of the aspects may look unusual or even impossible to explain.