Uniform convergence of Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity

In this article we prove uniform convergence estimates for the recently developed Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity. These multigrid methods are defined in terms of the “Galerkin approach,” where quadratic forms over coarse grids are constructed using the quadratic form on the finest grid and iterated coarse‐to‐fine intergrid transfer operators. Previously, uniform estimates were obtained for problems with full elliptic regularity, whereas these estimates are derived with less than full elliptic regularity here. Applications to the nonconforming P₁, rotated Q₁, and Wilson finite elements are analyzed. The result applies to the mixed method based on finite elements that are equivalent to these nonconforming elements. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 203–217, 2002; DOI 10.1002/num.10004     

OPTIMAL INSTAR PARASITIZATION IN A STAGE STRUCTURED HOST-PARASITOID MODEL

A stage structured host-parasitoid model is derived and the equilibria studied. It is shown under what conditions the parasitoid controls an exponentially growing host in the sense that a coexistence equilibrium exists. Furthermore, for host populations whose inherent growth rate is not too large it is proved that in order to minimize the adult host equilibrium level it is necessary that the parasitoids attack only one of the larval stages. It is also proved in this case that the minimum adult host equilibrium level is attained when the parasitoids attack that larval stage which also maximizes the expected number of emerging adult parasitoid per larva at equilibrium. Numerical simulations tentatively indicate that the first conclusion remains in general valid for the model. However, numerical studies also show that it is not true in general that the optimal strategy will maximize the number of emerging adult parasitoid per larva at equilibrium.     

On disjoint configurations in infinite graphs

For a graph A and a positive integer n, let nA denote the union of n disjoint copies of A; similarly, the union of ?₀ disjoint copies of A is referred to as ?₀A. It is shown that there exist (connected) graphs A and G such that nA is a minor of G for all n??, but ?₀A is not a minor of G. This supplements previous examples showing that analogous statements are true if, instead of minors, isomorphic embeddings or topological minors are considered. The construction of A and G is based on the fact that there exist (infinite) graphs G₁, G₂,… such that G_i is not a minor of G_j for all i ≠ j. In contrast to previous examples concerning isomorphic embeddings and topological minors, the graphs A and G presented here are not locally finite. The following conjecture is suggested: for each locally finite connected graph A and each graph G, if nA is a minor of G for all n ? ?, then ?₀A is a minor of G, too. If true, this would be a far‐reaching generalization of a classical result of R. Halin on families of disjoint one‐way infinite paths in graphs. © 2002 Wiley Periodicals, Inc. J Graph Theory 39: 222–229, 2002; DOI 10.1002/jgt.10016     

Dynamics of liquid membranes. II: Adaptive finite difference methods

Two domain-adaptive finite difference methods are presented and applied to study the dynamic response of incompressible, inviscid, axisymmetric liquid membranes subject to imposed sinusoidal pressure oscillations. Both finite difference methods map the time-dependent physical domain whose downstream boundary is unknown onto a fixed computational domain. The location of the unknown time-dependent downstream boundary of the physical domain is determined from the continuity equation and results in an integrodifferential equation which is non-linearly coupled with the partial differential equations which govern the conservation of mass and linear momentum and the radius of the liquid membrane. One of the finite difference methods solves the non-conservative form of the governing equations by means of a block implicit iterative method. This method possesses the property that the Jacobian matrix of the convection fluxes has an eigenvalue of algebraic multiplicity equal to four and of geometric multiplicity equal to one. The second finite difference procedure also uses a block implicit iterative method, but the governing equations are written in conservation law form and contain an axial velocity which is the difference between the physical axial velocity and the grid speed. It is shown that these methods yield almost identical results and are more accurate than the non-adaptive techniques presented in Part I. It is also shown that the actual value of the pressure coefficient determined from linear analyses can be exceeded without affecting the stability and convergence of liquid membranes if the liquid membranes are subjected to sinusoidal pressure variations of sufficiently high frequencies.     

Equivalence of u–p and ζ–ψ formulations of the time-dependent Navier–Stokes equations

This paper deals with the non-stationary incompressible Navier--Stokes equations for two-dimensional flows expressed in terms of the velocity and pressure and of the vorticity and streamfunction. The equivalence of the two formulations is demonstrated, both formally and rigorously, by virtue of a condition of compatibility between the boundary and initial values of the normal component of velocity. This condition is shown to be the only compatibility condition necessary to allow for solutions of a minimal regularity, namely H¹ for the velocity, as in most current numerical schemes relying on spatial discretizations of local type.     

Different acute effects of oral and intratracheal administration of disodium arsenate and gallium arsenide on heme synthesis in rats

The acute influences of arsenic compounds on the metabolism of porphyrins and heme in various organs of rats after oral or intratracheal administration of disodium arsenate (Na₂HAsO₄) and gallium arsenide (GaAs) were examined and compared. For the oral administration experiments, 21 or 84 mg of Na₂HAsO₄, or 2 or 4 g of GaAs, per cm³ saline per kg body weight of each animal was administered to Jcl: Wistar male rats and the organs were removed after exsanguination from the vein of the right axilla under anesthesia with ether, 16 h after administration. In the case of intratracheal administration, rats given 8.2 or 16.4 mg of Na₂HAsO₄, or 0.2 or 0.4 g GaAs per cm³ saline per kg body weight were examined under the same experimental conditions as for the administration route. Increase in the body weight of rats was suppressed after intratracheal administration of the two arsenic compounds. In these rats the hematocrit value increased significantly. These changes were not shown by the orally administered rats. Elevation in δ-aminolevulinate synthase (ALA-S, EC 2.3.1.37) activity in erythroblasts by Na₂HAsO₄ was much higher after intratracheal administration than after oral administration. Suppression in the activities of δ-aminolevulinate dehydratase (ALA-D, EC 4.2.1.24) and porphobilinogen deaminase (PBG-D, EC 4.3.1.8) in peripheral erythrocytes by Na₂HAsO₄ and GaAs were stronger by intratracheal administration than by the oral route. Influences of GaAs on the activity of PBG-D in rat liver were shown to be more effective by oral administration than by the intratracheal route. Oral administration of Na₂HAsO₄ and GaAs had a stronger suppression effect on the activities of ALA-D and PBG-D in rat kidney. It seems from these results that the different extents of the influence of arsenic compounds might depend on the routes of intake.     

Optimal control of unsteady compressible viscous flows

The control of complex, unsteady flows is a pacing technology for advances in fluid mechanics. Recently, optimal control theory has become popular as a means of predicting best case controls that can guide the design of practical flow control systems. However, most of the prior work in this area has focused on incompressible flow which precludes many of the important physical flow phenomena that must be controlled in practice including the coupling of fluid dynamics, acoustics, and heat transfer. This paper presents the formulation and numerical solution of a class of optimal boundary control problems governed by the unsteady two‐dimensional compressible Navier–Stokes equations. Fundamental issues including the choice of the control space and the associated regularization term in the objective function, as well as issues in the gradient computation via the adjoint equation method are discussed. Numerical results are presented for a model problem consisting of two counter‐rotating viscous vortices above an infinite wall which, due to the self‐induced velocity field, propagate downward and interact with the wall. The wall boundary control is the temporal and spatial distribution of wall‐normal velocity. Optimal controls for objective functions that target kinetic energy, heat transfer, and wall shear stress are presented along with the influence of control regularization for each case. Copyright © 2002 John Wiley & Sons, Ltd.     

Optimality conditions and duality theory for minimizing sums of the largest eigenvalues of symmetric matrices

The sum of the largest eigenvalues of a symmetric matrix is a nonsmooth convex function of the matrix elements. Max characterizations for this sum are established, giving a concise characterization of the subdifferential in terms of a dual matrix. This leads to a very useful characterization of the generalized gradient of the following convex composite function: the sum of the largest eigenvalues of a smooth symmetric matrix-valued function of a set of real parameters. The dual matrix provides the information required to either verify first-order optimality conditions at a point or to generate a descent direction for the eigenvalue sum from that point, splitting a multiple eigenvalue if necessary. Connections with the classical literature on sums of eigenvalues and eigenvalue perturbation theory are discussed. Sums of the largest eigenvalues in the absolute value sense are also addressed.This paper is dedicated to Phil Wolfe on the occasion of his 65th birthday.The work of this author was supported by the National Science Foundation under grants CCR-8802408 and CCR-9101640.The work of this author was supported in part during a visit to Argonne National Laboratory by the Applied Mathematical Sciences subprogram of the Office of Energy Research of the U.S. Department of Energy under contract W-31-109-Eng-38, and in part during a visit to the Courant Institute by the U.S. Department of Energy under Contract DEFG0288ER25053.     

Differential–algebraic equations in multibody system modeling

Numerical methods for the efficient integration of both stiff and nonstiff equations of motion of multibody systems having the form of differential-algebraic equations (DAE) of index 3 are discussed. Linear multi-step ABM and BDF methods are considered for the non-iterational integration of nonstiff DAE. The Park method is proposed for integration of stiff equations. This revised version was published online in June 2006 with corrections to the Cover Date.     

Synthesis of functional polymers—Vinylidene fluoride based fluorinated copolymers and terpolymers bearing bromoaromatic side‐group

The radical co‐ and terpolymerization of 4‐[(α,β,β‐trifluorovinyl)oxy]bromo benzene (TFVOBB) with 1,1‐difluoroethylene (or vinylidene fluoride, VDF, or VF₂), hexafluoropropene (HFP), perfluoromethyl vinyl ether (PMVE), and chlorotrifluroroethylene (CTFE) is presented. Although TFVOBB could be thermocyclodimerized, it could not homopolymerize under radical initiation. TFVOBB could be copolymerized in solution under a radical initiator with VDF or CTFE comonomers, while its copolymerization with HFP or PMVE were unsuccessful. The terpolymerization of TFVOBB with VDF and HFP, or VDF and PMVE, or VDF and CTFE also led to original fluorinated terpolymers bearing bromoaromatic side‐groups. The conditions of co‐ and terpolymerization were optimized in terms of the nature of the radical initiators, and of the nature of solvents (fluorinated or nonhalogenated). Various monomer concentrations in the co‐ and terpolymers were assessed by ¹⁹F and ¹H‐NMR spectroscopy. The thermal and physico chemical properties were also studied. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5077–5097, 2004