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.     

Spectral Representation Theory for Dielectric Behavior of Nonspherical Cell Suspensions

Recent experiments revealed that the dielectric dispersion spectrum of fission yeast cells in a suspension was mainly composed of two sub-dispersions. The low-frequency sub-dispersion depended on the cell length, while the high-frequency one was independent of it. The cell shape effect was simulated by an ellipsoidal cell model but the comparison between theory and experiment was far from being satisfactory. Prompted by the discrepancy, we proposed the use of spectral representation to analyze more realistic cell models. We adopted a shell-spheroidal model to analyze the effects of the cell membrane. It is found that the dielectric property of the cell membrane has only a minor effect on the dispersion magnitude ratio and the characteristic frequency ratio. We further included the effect of rotation of dipole induced by an external electric field, and solved the dipole-rotation spheroidal model in the spectral representation.Good agreement between theory and experiment has been obtained.     

Physical aging of a polyetherimide: Creep and DSC measurements

Creep and differential scanning calorimetry (DSC) measurements have been used to study the physical aging behavior of a polyetherimide. Isothermal aging temperatures ranged from 160°C to T_g with aging times ranging from 10 min to 8 days. The only measurable effect of physical aging on the short-time creep curves is a shift of the creep compliance to longer times. Andrade plots of the compliance versus the cube root of time are linear at short times with the slope β decreasing with increasing aging time to a constant value once equilibrium is reached. Log β³ is related directly to the degree to which the creep curves shift to longer times with physical aging, and is used in this work as a measure of physical aging. A reduced curve of log β³ versus log aging time is obtained for the aging temperatures investigated by appropriate vertical and horizontal shifts. The enthalpy change during aging increases linearly with the logarithm of the aging time, t_a, leveling off at equilibrium at values which increase with decreasing aging temperature. Hence, both nonequilibrium and equilibrium temperature shift factors can be calculated from the DSC data. Good agreement is observed between the equilibrium temperature shift factors obtained from the creep and DSC data. The temperature dependence of the nonequilibrium temperature shift factors is found to be an order of magnitude smaller than that of the equilibrium shift factors. The time scales to reach equilibrium for enthalpy and for mechanical measurements are found to be the same within experimental error. © 1995 John Wiley & Sons, Inc.     

Numerical simulation on an interaction of a vortex street with an elliptical leading edge using an unstructured grid

An algorithm for a time accurate incompressible Navier–Stokes solver on an unstructured grid is presented. The algorithm uses a second order, three‐point, backward difference formula for the physical time marching. For each time step, a divergence free flow field is obtained based on an artificial compressibility method. An implicit method with a local time step is used to accelerate the convergence for the pseudotime iteration. To validate the code, an unsteady laminar flow over a circular cylinder at a Reynolds number of 200 is calculated. The results are compared with available experimental and numerical data and good agreements are achieved. Using the developed unsteady code, an interaction of a Karman vortex street with an elliptical leading edge is simulated. The incident Karman vortex street is generated by a circular cylinder located upstream. A clustering to the path of the vortices is achieved easily due to flexibility of an unstructured grid. Details of the interaction mechanism are analysed by investigating evolutions of vortices. Characteristics of the interactions are compared for large‐ and small‐scale vortex streets. Different patterns of the interaction are observed for those two vortex streets and the observation is in agreement with experiment. Copyright © 2004 John Wiley & Sons, Ltd.     

Space–time integrated least squares: a time‐marching approach

A time‐marching formulation is derived from the space–time integrated least squares (STILS) method for solving a pure hyperbolic convection equation and is numerically compared to various known methods. Copyright © 2004 John Wiley & Sons, Ltd.     

AN ENCOMPASSING APPROACH TO MODELING FISHERY DYNAMICS: MODELING DYNAMIC NONLINEAR SYSTEMS

Uncertainty about the nature and significance of nonlinearities and the manner in which dynamics affect future realizations makes model specification the most difficult aspect of modeling dynamic systems. By interpreting several popular fishery models as subcases of a nesting dynamic Taylor series approximation, we isolate the specification differences between these models in a way that accounts for commonalities. On the argument that the differences due to alternative nonlinear forms are likely to be small compared to more mundane considerations such as delay difference and general dynamic lag specification, we propose an alternative model that uses the terms from the first order approximation common to all models combined with a data-based determination of the appropriate lags using the methods of state space time series analysis. Finally, the success of the alternative models is judged in an application to Pacific halibut data.     

Orientation and relaxation of polymer–clay solutions studied by rheology and small-angle neutron scattering

The influence of shear on viscoelastic solutions of poly(ethylene oxide) (PEO) and clay [montmorillonite, i.e., Cloisite NA+ (CNA)] was investigated with rheology and small-angle neutron scattering (SANS). The steady-state viscosity and SANS were used to measure the shear-induced orientation and relaxation of the polymer and clay platelets. Anisotropic scattering patterns developed at much lower shear rates than in pure clay solutions. The scattering anisotropy saturated at low shear rates, and the CNA clay platelets aligned with the flow, with the surface normal parallel to the gradient direction. The cessation of shear led to partial and slow randomization of the CNA platelets, whereas extremely fast relaxation was observed for laponite (LRD) platelets. These PEO–CNA networklike solutions were compared with previously reported PEO–LRD networks, and the differences and similarities, with respect to the shear orientation, relaxation, and polymer–clay interactions, were examined. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 3102–3112, 2004     

狭义相对论效应──运动时钟变慢的计算机模拟演示实验

利用精炼的C语言进行程序设计，定量模拟狭义相结论的运动时钟变慢，并对处于两种不同时空观的时钟运行进行对比，使得时间膨胀效应变得直观。     

Electrical properties of polypyrrole/p‐InP structure

The polypyrrole/p‐InP structure has been fabricated by the electrochemical polymerization of the organic polypyrrole onto the p‐InP substrate. The current–voltage (I–V), capacitance–voltage (C–V), and capacitance–frequency (C–f) characteristics of the PPy/p‐InP structure have been determined at room temperature. The structure showed nonideal I–V behavior with the ideality factor and the barrier height 1.48 and 0.69 eV respectively. C–f measurements of the structure have been carried out using the Schottky capacitance spectroscopy technique and it has been seen that there is a good agreement between the experimental and theoretical values. Also, it has been seen that capacitance almost show a plateau up to a certain value of frequency, after which, the capacitance decreases. The higher values of capacitance at low frequencies were attributed to the excess capacitance resulting from the interface states in equilibrium with the p‐InP that can follow the a.c. signal. The interface state density N_ss and relaxation time τ of the structure were determined from C–f characteristics. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 1572–1579, 2006     

Studies on solution cast perfluorocarbonsulfonic acid ionomers

The thermal properties of perfluorocarbonsulfonic acid ionomer films cast from solution and their temperature-dependent far infrared spectra have been investigated. In addition to the endotherms and spectral changes associated with the loss of solvent as the films were heated, a significant exotherm has been observed and assigned to the relaxation associated with rearrangements resulting in partially crystallized phases in insoluble thermally treated membranes. © 1993 John Wiley & Sons, Inc.