A hybrid finite-element–volume-of-fluid method for simulating free surface flows and interfaces

A numerical technique is developed for the simulation of free surface flows and interfaces. This technique combines the strength on the finite element method (FEM) in calculating the field variables for a deforming boundary and the versatility of the volume-of-fluid (VOF) technique in advection of the fluid interfaces. The advantage of the VOF technique is that it allows the simulation of interfaces with large deformations, including surface merging and breaking. However, its disadantage is that is solving the flow equations, it cannot resolve interfaces smaller than the cell size, since information on the subgrid scale is lost. Therefore the accuracy of the interface reconstruction and the treatment of the boundary conditions (i.e. viscous stresses and surface tension forces) become grid-size-dependent. On the other hand, the FEM with deforming interface mesh allows accurate implementation of the boundary conditions, but it cannot handle large surface deformations occurring in breaking and merging of liquid regions. Combining the two methods into a hybrid FEM-VOF method eliminates the major shortcomings of both. The outcome is a technique which can handle large surface deformations with accurate treatment of the boundary conditions. For illustration, two computational examples are presented, namely the instability and break-up of a capillary jet and the coalescence collision of two liquid drops.     

Algorithms for finding proper essential surfaces in 3-manifolds

In this paper, we present an algorithm which, for a given compact orientable irreducible boundary irreducible 3-manifold M, verifies whether M contains an essential orientable surface (possibly, with boundary), whose genus is at most N. The algorithm is based on Haken’s theory of normal surfaces, and on a trick suggested by Jaco and consisting in estimating the mean length of boundary curves in an unknown essential surface of a given genus in the given manifold.     

On error operators related to the arbitrary functions principle

The error on a real quantity Y due to the graduation of the measuring instrument may be asymptotically represented, when the graduation is regular and fines down, by a Dirichlet form on R whose square field operator does not depend on the probability law of Y as soon as this law possesses a continuous density. This feature is related to the “arbitrary functions principle” (Poincaré, Hopf). We give extensions of this property to Rd and to the Wiener space for some approximations of the Brownian motion. This gives new approximations of the Ornstein-Uhlenbeck gradient. These results apply to the discretization of some stochastic differential equations encountered in mechanics.     

On the invariant measures for the two-dimensional euler flow

It is proven that the canonical Gibbs measure associated with a gas of vortices of intensity ± converges, in the limitN, 0,Nconst, to a Gaussian measure, which is invariant for the two-dimensional Euler equation.On leave from Dipartimento di Matematica Università di Roma Tor Vergata Roma, Italy.On leave from Dipartimento di Matematica Università di Roma La Sapienza, Roma, Italy.     

Elliptic Reconstruction and a Posteriori Error Estimates for Fully Discrete Semilinear Parabolic Optimal Control Problems

This article studies a posteriori error analysis of fully discrete finite element approximations for semilinear parabolic optimal control problems. Based on elliptic reconstruction approach introduced earlier by Makridakis and Nochetto [25], a residual based a posteriori error estimators for the state, co-state and control variables are derived. The space discretization of the state and co-state variables is done by using the piecewise linear and continuous finite elements, whereas the piecewise constant functions are employed for the control variable. The temporal discretization is based on the backward Euler method. We derive a posteriori error estimates for the state, co-state and control variables in the $L^\infty(0,T;L^2(\Omega))$-norm. Finally, a numerical experiment is performed to illustrate the performance of the derived estimators.     

Lattice-BGK approach to simulating granular flows

Many continuum theories for granular flow produce an equation of motion for the fluctuating kinetic energy density (granular temperature) that accounts for the energy lost in inelastic collisions. Apart from the presence of an extra dissipative term, this equation is very similar in form to the usual temperature equation in hydrodynamics. It is shown how a lattice-kinetic model based on the Bhatnagar-Gross-Krook (BGK) equation that was previously derived for a miscible two-component fluid may be modified to model the continuum equations for granular flow. This is done by noting that the variable corresponding to the concentration of one species follows an equation that is essentially analogous to the granular temperature equation. A simulation of an unforced granular fluid using the modified model reproduces the phenomenon of clustering instability, namely the spontaneous agglomeration of particles into dense clusters, which occurs generically in all granular flows. The success of the continuum theory in capturing the gross features of this basic phenomenon is discussed. Some shear flow simulations are also presented.     

Characterization of the organization in shear layers via the Proper Orthogonal Decomposition

Experiments are performed in an incompressible plane turbulent mixing layer, using various hot wire rake configurations. From these experiments, the Proper Orthogonal Decomposition is applied for kernels where the space-time correlation tensor is evaluated over different spatial meshes and velocity components configurations. The resulting decompositions are then discussed in terms of characterization of the organization of the flow for various scalar or vectorial approaches of the POD. An incrtial range law is evidenced. The instantaneous contribution of the first modes of the POD to the organization of the flow is analyzed. A dynamical behavior for the organization of the flow is observed from the correlation between the first two modes contribution.     

Spreading of viscous droplets on a non viscous liquid

Some polymer melts (of high viscosity ) can wet completely the surface of a non miscible, simple liquid. We discuss here the laws of spreading for a macroscopic droplet of this type, when the internal friction of the droplet dominates. We predict a droplet radius increasing liket ^1/4 wheret is the spreading time, or equivalently a droplet curvature decreasing liket ^–1. The droplet should be surrounded by a precursor film, which is not discussed in the present note.     

飞行时间质谱仪中Ni^＋与酯分子束流的气相反应及机理

在一个石英反应室内，让激光解离的Ｎｉ＾＋与连续喷入的酯类化合物分子束流反应，生成的产物离子经飞行时间质谱仪检测，研究了Ｎｉ＾＋与乙酸乙酯、乙酸丙烯酯、乙酸乙烯酯、乙酸丁酯的气相反应。根据如下机理满意的解释了发生的反应：（１）Ｎｉ＾＋与酯形成激发态络合物；（２）β－Ｈ迁移至烃氧基的氧上形成Ｎｉ＾＋的双配体化合物离子；（３）失去一中性配体分子得Ｎｉ＾＋的单配体化合物离子。由酯的ΔＨ离和Ｄ（Ｎｉ＾＋－Ａ     

Cathode of a fuel cell with a solid polymer electrolyte: Calculating overall currents in the presence of a gas-diffusion layer

Mathematical apparatus, which makes it possible to perform calculations of the current-voltage characteristics of cathodes of fuel cells with a solid polymer electrolyte in conditions where there are present extraneous diffusion restrictions is proposed. In so doing, the partial pressure of oxygen and the absolute pressure of gas in the gas chamber may assume any values. First of all presented are the results of calculations of the current-voltage characteristics intrinsic to active layers of the air and oxygen cathodes, which are performed under the assumption that the extraneous diffusion restrictions are absent altogether. Thereafter, in the same conditions (at the same parameters that characterize the active layer of a cathode), obtained are results of a calculation of the current-voltage characteristics inherent in the air and oxygen cathodes in the presence of extraneous diffusion restrictions. Afterward there is performed an analysis of the way a gas-diffusion layer restricts the process of generation of current in a cathode and of what measures should be taken in order for the extraneous diffusion restrictions to become less significant.