Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Cochran's statistical theorem for outer inverses of matrices and matrix quadratic forms

We extend the matrix version of Cochran's statistical theorem to outer inverses of a matrix. As applications, we investigate the Wishartness and independence of matrix quadratic forms for Kronecker product covariance structures.     

Matrix effects in the determination of butyltin compounds in environmental samples by GC AA after hydride generation

Experiments for the determination of mono-, di and tri-butyltin (MBT, DBT and TBT) by hydride generation/gas chromatography/atomic absorption spectrometry in various matrices (sediment, suspended matter, mussel, algae and water) have revealed that poor butyltin recoveries are obtained in sediments displaying high sulphur and hydrocarbon contents; very poor recoveries were also observed for TBT in sediments with high chlorophyll pigment contents as well as in algal samples. It was however not clear whether the hydride generatin was inhibited by these infering compounds, as was previously assumed in the case of hydrocarbons, or whether interferences affected the atomization rate. Further studies were performed to solve this problem in order to validate this method in the case of analyses of, for example, oil-contaminated sediment and algae. This paper presents the results obtained. It is concluded here that the poor recoveries were due to an inhibition of hydride generation rather than to interference at the atomization stage.     

Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

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.     

A note on an SOR-like method for augmented systems

Golub et al. (2001, BIT, 41, 71–85) gave a generalizedsuccessive over-relaxation method for the augmented systems.In this paper, the connection between the SOR-like method andthe preconditioned conjugate gradient (PCG) method for the augmentedsystems is investigated. It is shown that the PCG method isat least as accurate (fast) as the SOR-like method. Numericalexamples demonstrate that the PCG method is much faster thanthe SOR-like method.     

Morphology of Nylon-6 blends with styrenic polymers

The morphology of blends of styrenic polymers in a matrix of 75% Nylon-6 prepared in a Brabender Plasti-Corder was examined by scanning electron microscopy. Styrene/acrylonitrile copolymers (SAN) form smaller particles as the AN level increases owing to the corresponding decrease in the SAN–polyamide interfacial tension. Various styrenic polymers containing functional groups, maleic anhydride or oxazoline type, that can react with Nylon-6 during melt processing were added to the SAN phase which also led to a decrease in the particle size owing to the graft copolymer formed in situ. The effects of functional group type, amount of functional groups per chain, amount of functional polymer added, and the miscibility of the styrene/maleic anhydride (SMA) and SAN copolymers on the morphology of the styrenic phase in the Nylon-6 matrix are described. © 1992 John Wiley & Sons, Inc.     

The completely positive and doubly nonnegative completion problems

We prove the following. Let G be an undirected graph. Every partially specified symmetric matrix, the graph of whose specified entries is G and each of whose fully specified submatrices is completely positive (equal to BB^T for some entrywise nonnegative matrix B), may be completed to a completely positive matrix if and only if G is a block-clique graph (a chordal graph in which distinct maximal cliques overlap in at most one vertex). The same result holds for matrices that are doubly nonnegative (entrywise nonnegative and positive semidefinite).     

Synthesis and characterization of functional gradient copolymers of allyl methacrylate and butyl acrylate

Functional spontaneous gradient copolymers of allyl methacrylate (A) and butyl acrylate (B) were synthesized via atom transfer radical polymerization. The copolymerization reactions were carried out in toluene solutions at 100 °C with methyl 2‐bromopropionate as the initiator and copper bromide with N,N,N′,N″,N″‐pentamethyldiethylenetriamine as the catalyst system. Different aspects of the statistical reaction copolymerizations, such as the kinetic behavior, crosslinking density, and gel fraction, were studied. The gel data were compared with Flory's gelation theory, and the sol fractions of the synthesized copolymers were characterized by size exclusion chromatography and nuclear magnetic resonance spectroscopy. The copolymer composition, demonstrating the gradient character of the copolymers, and the microstructure were analyzed. The experimental data agreed well with data calculated with the Mayo–Lewis terminal model and Bernoullian statistics, with monomer reactivity ratios of 2.58 ± 0.37 and 0.51 ± 0.05 for A and B, respectively, an isotacticity parameter for A of 0.24, and a coisotacticity parameter of 0.33. © 2006 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 5304–5315, 2006     

Literal‐paraconsistent and literal‐paracomplete matrices

We introduce a family of matrices that define logics in which paraconsistency and/or paracompleteness occurs only at the level of literals, that is, formulas that are propositional letters or their iterated negations. We give a sound and complete axiomatization for the logic defined by the class of all these matrices, we give conditions for the maximality of these logics and we study in detail several relevant examples. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)