A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

Dispersion analysis of the least‐squares finite‐element shallow‐water system

The frequency or dispersion relation for the least‐squares mixed formulation of the shallow‐water equations is analysed. We consider the use of different approximation spaces corresponding to co‐located and staggered meshes, respectively. The study includes the effect of Coriolis, and the dispersion properties are compared analytically and graphically with those of the mixed Galerkin formulation. Numerical solutions of a test problem to simulate slow Rossby modes illustrate the theoretical results. Copyright © 2003 John Wiley & Sons, Ltd.     

Diffusion via splitting and remeshing via merging in vortex methods

The technique of splitting a fat vortex element (with a core width larger than some threshold) into some thin ones in order to fix the convergence problem of the core‐spreading vortex methods is convenient and efficient. In particular, it keeps the method purely Lagrangian. In the present investigation, the splitting process is further viewed as part of the physical diffusion process. A new splitting method in which several weaker child vortices surround a thinned but still strong parent vortex is proposed. It is found that because of the survival of the parent vortex, the error arising from the splitting events can be largely reduced. The computational amount on the other hand is kept reasonably large by merging similar and close‐by vortices. The merging scheme designed herein not only involves fewer restrictions but also allows merging vortices of opposite rotations through the viewpoint of remeshing. The validity and accuracy of these techniques, proposed particularly for simulations undergoing lots of splitting and merging events, are verified by successfully simulating the interactions between two Burgers vortices under an external straining field. Copyright © 2005 John Wiley & Sons, Ltd.     

On the use of the Galerkin method for 3D numerical modelling of the general circulation: the South Atlantic eperiment

A linear three‐dimensional hydrodynamical numerical model, with the application of the Galerkin Method for the vertical dependence, is here presented. The spherical coordinate system is used, in order to allow large‐scale simulations. The equations and mathematical development of the model are shown in detail, together with the boundary and initial conditions, and the sequence of equations' solution. The model is applied to the South Atlantic Ocean, for estimating typical seasonal circulations, and the results are summarized in maps of currents at surface and 1000 m depth, and in transport values of the Brazil Current between 30°S and 40°S. Copyright © 2002 John Wiley & Sons, Ltd.     

Meshless simulation of equilibrium swelling/deswelling of pH‐sensitive hydrogels

Hydrogels have been widely used in microelectromechanical systems (MEMS) and Bio‐MEMS devices. In this article, the equilibrium swelling/deswelling of the pH‐stimulus cylindrical hydrogel in the microchannel is studied and simulated by the meshless method. The multi‐field coupling model, called multi‐effect‐coupling pH‐stimulus (MECpH) model, is presented and used to describe the chemical field, electric field, and the mechanical field involved in the problem. The partial differential equations (PDEs) describing these three fields are either nonlinear or coupled together. This multi‐field coupling and high nonlinear characteristics produce difficulties for the conventional numerical methods (e.g., the finite element method or the finite difference method), so an alternative—meshless method is developed to discretize the PDEs, and the efficient iteration technique is adopted to solve the nonlinear problem. The computational results for the swelling/deswelling diameter of the hydrogel under the different pH values are firstly compared with experimental results, and they have a good agreement. The influences of other parameters on the mechanical properties of the hydrogel are also investigated in detail. It is shown that the multi‐field coupling model and the developed meshless method are efficient, stable, and accurate for simulation of the properties of the stimuli‐sensitive hydrogel. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 326–337, 2006     

Interval parameter perturbation method for evaluating the bounds on natural frequencies of structures with interval parameters

For structural parameters with uncertainties, interval mathematics can, in the case where the probabilistic distribution density of uncertain variables is unavailable, deal with the influence of uncertainties in structural parameters on the response of structures. In order to evaluate the region containing natural frequencies of structures with interval parameters, the interval parameter perturbation method is presented in this paper. The advantage of the present method is its computational efficiency in evaluating the region containing natural frequencies. A numerical example is used to illustrate the efficiency of the method proposed. The project is supported by National Youth Natural Science Foundation of China and National Outstanding Youth Science Foundation of China.     

苯胺气相N－烷基化反应研究Ⅱ．苯胺和乙醇常压催化合成N－乙基苯胺的工业应用

苯胺气相Ｎ－烷基化反应研究Ⅱ．苯胺和乙醇常压催化合成Ｎ－乙基苯胺的工业应用陈骏如，李瑞祥，韩银仙，李贤均（四川大学化学系成都６１００６４）Ｎ－乙基苯胺是制造染料、塑料、炸药和医药等许多化工产品的重要中间体，其合成方法过去多采用无机酸催化剂，如Ｈ＿２Ｓ...     

AN INVESTIGATION ON THE ab initio CALCULATION OF INTERMOLECULAR INTERACTION——NON-BONDED INTERACTION IN THE HYDROGEN BOND COMPLEX (HF)_2

An ab initio calculation of the hydrogen bond complex (HF)_2 is given with the 6-311 G~(**) basis set, according to which the potential surface around the balance point of the distancès and the orientations between two HF molecules is obtained. The atomic charges in the system are calculated with the PD/LSF method (potential-derived/least-square-fitting method) and then an analysis of the hydrogen bond interaction between two HF molecules is given with the (exp-6-1) potential function, by means of which it is shown that the main interaction between them is not an electro-static but a charge transfer one. The potential curve between two HF molecules is like a Morse function.     

A Robust and Efficient Method for the Computation of Equilibrium Composition in Gaseous Mixtures

This paper presents in detail a robust, efficient and accurate methodology for the computation of equilibrium composition in gaseous mixtures. The methodology is founded on the concept of the chemical basis, which is defined and formalized using a powerful matricial approach. The method is specially designed to be general, thus providing basic thermodynamic data in several areas, such as combustion, plasma chemistry and, more generally speaking, computational fluid dynamics. The performance of the method is given in terms of CPU usage and the computed results are compared with those in the published literature. The method is shown to yield results of very high quality in terms of accuracy and smoothness.     

Kinetic model ofCO oxidation on nonuniform surface analyzed regardingCOads spillover, 2. homotopic method. Effect ofm1 surface portion

The homotopic method has been used to analyze the kinetic model of three-stageCO oxidation on two nonuniform surface patches conjugated byCO ₂ spillover. Diagrams of steady states depending on the portion of surface patchm ₁ at various temperatures and pressure ratiosP(O ₂)/P(CO) have been constructed. The ratios of different type patches corresponding to the maximum overall reaction rate have been found.