A Nonconforming Arbitrary Quadrilateral Finite Element Method for Approximating Maxwell's Equations

The main aim of this paper is to provide convergence analysis of Quasi-Wilson nonconforming finite element to Maxwell's equations under arbitrary quadrilateral meshes.The error estimates are derived,which are the same as those for conforming elements under conventional regular meshes.     

Critical percolation exploration path and SLE 6: a proof of convergence

It was argued by Schramm and Smirnov that the critical site percolation exploration path on the triangular lattice converges in distribution to the trace of chordal SLE ₆. We provide here a detailed proof, which relies on Smirnov’s theorem that crossing probabilities have a conformally invariant scaling limit (given by Cardy’s formula). The version of convergence to SLE ₆ that we prove suffices for the Smirnov–Werner derivation of certain critical percolation crossing exponents and for our analysis of the critical percolation full scaling limit as a process of continuum nonsimple loops. Research of Charles M.Newman was partially supported by the US NSF under grants DMS-01-04278 and DMS-06-06696.     

An algebraic multigrid solver for Navier-Stokes problems

An efficient numerical method is presented for solving the equations of motion for viscous fluids. The equations are discretized on the basis of unstructured finite element meshes and then solved by direct iteration. Advective fluxes are temporarily fixed at each iteration to provide a linearized set of coupled equations which are then also solved by iteration using a fully implicit algebraic multigrid (AMG) scheme. A rapid convergence to machine accuracy is achieved that is almost mesh-independent. The scaling of computing time with mesh size is therefore close to the optimum.     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

ContributionsTriangular embeddings of Kn−Km with unboundedly large m

The author has proposed methods of constructing index 2 and 3 current graphs generating triangular embeddings of graphs K_n−K_m with unboundedly large m (as n increases). As a result, triangular embeddings of graphs of many families of graphs K_n−K_m with unboundedly large m were constructed. The paper gives a survey of these results and a short explanation of the methods.     

Maximum-norm error analysis of a non-monotone FEM for a singularly perturbed reaction-diffusion problem

A non-monotone FEM discretization of a singularly perturbed one-dimensional reaction-diffusion problem whose solution exhibits strong layers is considered. The method is shown to be maximum-norm stable although it is not inverse monotone. Both a priori and a posteriori error bounds in the maximum norm are derived. The a priori result can be used to deduce uniform convergence of various layer-adapted meshes proposed in the literature. Numerical experiments complement the theoretical results. AMS subject classification (2000)  65L10, 65L50, 65L60     

An $L^{\infty}$ Second Order Cartesian Method for 3D Anisotropic Interface Problems

A second order accurate method in the infinity norm is proposed for general three dimensional anisotropic elliptic interface problems in which the solution and its derivatives, the coefficients, and source terms all can have finite jumps across one or several arbitrary smooth interfaces. The method is based on the 2D finite element-finite difference (FE-FD) method but with substantial differences in method derivation, implementation, and convergence analysis. One of challenges is to derive 3D interface relations since there is no invariance anymore under coordinate system transforms for the partial differential equations and the jump conditions. A finite element discretization whose coefficient matrix is a symmetric semi-positive definite is used away from the interface; and the maximum preserving finite difference discretization whose coefficient matrix part is an M-matrix is constructed at irregular elements where the interface cuts through. We aim to get a sharp interface method that can have second order accuracy in the point-wise norm. We show the convergence analysis by splitting errors into several parts. Nontrivial numerical examples are presented to confirm the convergence analysis.     

High throughput continuous titration based on a flow ratiometry controlled with feedback-based variable triangular waves and subsequent fixed triangular waves

We propose a new approach for high-throughput continuous titration based on a flow ratiometry. The method was applied to potentiometric titrations of acids and bases. A base solution, the flow rate F_B of which was varied in response to controller output voltage V_c, was merged with an acid solution under constant total (acid + base) flow rate. Downstream, the pH of the mixed solution was measured with a flow-through glass electrode. Initially, V_c, and thus F_B was increased linearly. At the instant the equivalence point (EP) was sensed, the ramp direction of V_c was reversed from upward to downward. The direction was reversed to upward when EP was sensed again. Such the feedback-based operation gives a triangular waveform of V_c, because there is a delay corresponding mainly to the transit time of merged solutions to reach the sensor. The value of V_c that gives EP composition, V_E, was estimated by averaging the most recent maximum and minimum values of V_c. Next, fixed triangular waves of V_c was used to control F_B. The amplitude and the scan rate of the waves were fixed narrower and faster, respectively, than those in the feedback-based operation in order to improve the throughput rate. The EP can be located as long as the scan range covers V_E. These automated processes limited the titration to just the narrow range around EP, and thus realized extremely high throughput rate of maximally 17.6 titrations per minute (=3.4 s per titration) at R.S.D. = 0.35%.     

On ω-limit Sets of Triangular Maps on the Unit Cube

We consider Sturm–Liouville eigenvalue problems of second order with arbitary seperated boundary conditions and perform a suitabl discretization of them. The obtained discrete Sturm–Liouville eigenvalue problems are examined and the asymptotic behavior of their eigenvalue as the norm of the partition tends to zero is investigated.