A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

Robust handling of non‐linear constraints for GA optimization of aerodynamic shapes

A new approach to the robust handling of non‐linear constraints for GAs (genetic algorithms) optimization is proposed. A specific feature of the approach consists of the change in the conventional search strategy by employing search paths which pass through both feasible and infeasible points (contrary to the traditional approach where only feasible points may be included in a path). The method (driven by full Navier–Stokes computations) was applied to the problem of multiobjective optimization of aerodynamic shapes subject to various geometrical and aerodynamic constraints. The results demonstrated that the method retains high robustness of conventional GAs while keeping CFD computational volume to an acceptable level, which allowed the algorithm to be used in a demanding engineering environment. Copyright © 2004 John Wiley & Sons, Ltd.     

Reversible nonlinear conduction in high‐density polyethylene/acetylene carbon black composites at various ambient temperatures

The reversible nonlinear conduction (RNC) in of high‐density polyethylene/acetylene carbon black composites with different degrees of crosslinking was studied above room temperature and below the melting point of high‐density polyethylene (HDPE). The experimental current density‐electric field strength curves can be overlapped onto a master curve, suggesting that the microscopic mechanisms for the appearance of RNC exist regardless of the ambient temperature and the crosslinking degree of the HDPE matrix. The relationship between the crossover current density and the linear conductivity can be explained in the framework of the dynamic random‐resistor‐network model. According to these results, two electron‐tunneling models are suggested to interpret the microscopic conduction behavior. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1212–1217, 2004     

Dynamic surface tension measurements of surfactant solutions using the maximum bubble pressure method – limits of applicability

One of the essential differences in the design of bubble pressure tensiometers consists in the geometry of the measuring capillaries. To reach extremely short adsorption times of milliseconds and below, the so-called deadtime of the capillaries must be of the order of some 10 ms. In particular, for concentrated surfactant solutions, such as micellar solutions, short deadtimes are needed to minimize the initial surfactant load of the generated bubbles. A theoretical model is derived and confirmed by experiments performed for a wide range of experimental conditions, mainly in respect to variations in deadtime and bubble volume.     

A note on the necessary conditions for the algebraic Riccati equation

We derive some useful and easily computable necessary conditionsfor the existence of a positive semi-definite solution to thealgebraic Ricatti equation (ARE). A motivating example is givento highlight the usefulness of the conditions for controllerand observer designs for nonlinear systems. Further, an upperbound on the trace of the solution to the ARE is also derived.     

Nonlinear optical properties of regioregular main‐chain polyesters

We have prepared new polyesters containing quadratic, nonlinear optical (NLO) active chromophores covalently incorporated into the main chain. In these polymers, the sequence of the chromophore units along the main chain is rigorously head to tail. All the polyesters are processable, both in the melt and in solution. For one polyester, a full second‐order NLO characterization has been performed. An out‐of‐resonance d₃₃ coefficient of 21 pm/V at 1368 nm has been measured. © 2007 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 2719–2725, 2007     

Studies on calcium‐containing poly(urethane ether)s

The calcium salt of mono(hydroxyethoxyethyl)phthalate [Ca(HEEP)₂] was synthesized by the reaction of diethylene glycol, phthalic anhydride, and calcium acetate. Calcium‐containing poly(urethane ether)s (PUEs) were synthesized by the reaction of hexamethylene diisocyanate (HMDI) or tolylene 2,4‐diisocyanate (TDI) with a mixture of Ca(HEEP)₂ and poly(ethylene glycol) (PEG₃₀₀ or PEG₄₀₀) with di‐n‐butyltin dilaurate as a catalyst. A series of calcium‐containing PUEs of different compositions were synthesized with Ca(HEEP)₂/PEG₃₀₀ (or PEG₄₀₀)/diisocyanate (HMDI or TDI) molar ratios of 2:2:4, 3:1:4, and 1:3:4 so that the coating properties of the PUEs could be studied. Blank PUEs without calcium‐containing ionic diols were also prepared by the reaction of PEG₃₀₀ or PEG₄₀₀ with HMDI or TDI. The PUEs were well characterized by Fourier transform infrared, ¹H and ¹³C NMR, solid‐state cross‐polarity/magic‐angle‐spinning ¹³C NMR, viscosity, solubility, and X‐ray diffraction studies. The thermal properties of the polymers were also studied with thermogravimetric analysis and differential scanning calorimetry. The PUEs were applied as top coats on acrylic‐coated leather, and their physicomechanical properties were also studied. The coating properties of PUEs, such as the tensile strength, elongation at break, tear strength, water vapor permeability, flexing endurance, cold crack resistance, abrasion resistance, color fastness, and adhesive strength, were better than the standard values. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 41: 2865–2878, 2003     

Extension of a combined analytical/numerical initial value problem solver for unsteady periodic flow

Here we describe analytical and numerical modifications that extend the Differential Reduced Ejector/ mixer Analysis (DREA), a combined analytical/numerical, multiple species ejector/mixing code developed for preliminary design applications, to apply to periodic unsteady flow. An unsteady periodic flow modelling capability opens a range of pertinent simulation problems including pulse detonation engines (PDE), internal combustion engine ICE applications, mixing enhancement and more fundamental fluid dynamic unsteadiness, e.g. fan instability/vortex shedding problems. Although mapping between steady and periodic forms for a scalar equation is a classical problem in applied mathematics, we will show that extension to systems of equations and, moreover, problems with complex initial conditions are more challenging. Additionally, the inherent large gradient initial condition singularities that are characteristic of mixing flows and that have greatly influenced the DREA code formulation, place considerable limitations on the use of numerical solution methods. Fortunately, using the combined analytical–numerical form of the DREA formulation, a successful formulation is developed and described. Comparison of this method with experimental measurements for jet flows with excitation shows reasonable agreement with the simulation. Other flow fields are presented to demonstrate the capabilities of the model. As such, we demonstrate that unsteady periodic effects can be included within the simple, efficient, coarse grid DREA implementation that has been the original intent of the DREA development effort, namely, to provide a viable tool where more complex and expensive models are inappropriate. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Numerical analysis for the periodic response of nonsmooth dynamic systems of high-dimension

This paper deals with the numerical solution of high dimensional dynamic systems with nonsmooth characteristics, such as structural elastoplasticity, mechanical clearance and dry friction. As a stable periodic response of the system attracts its transient response, it may be extracted from the transient response somehow. The paper suggests a novel iteration scheme to extrapolate the periodic response from a short time history of its transient response by curve fitting. Compared with the current schemes such as shooting and incremental harmonic balancing, the present scheme makes full use of the information within the transient response and the system characters, then its computation efficiency has increased by an order and the numerical convergence depending on the initial iteration has greatly improved.The subject supported by the Chinese Foundation of Aeronautical Science