A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

A complete boundary integral formulation for compressible Navier–Stokes equations

A complete boundary integral formulation for compressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for wall pressure and wall skin friction of two‐dimensional compressible laminar viscous flow around airfoils are in good agreement with field numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.     

Finite element computations of two-dimensional arterial flow in the presence of a transverse magnetic field

A finite element solution of the Navier-Stokes equations for steady flow under the magnetic effect through a double-branched two-dimensional section of a three-dimensional model of the canine aorta is discussed. The numerical scheme involves transforming the physical co-ordinates to a curvilinear boundary-fitted co-ordinate system. The shear stress at the wall is calculated for a Reynolds number of 1000 with the branch-to-main aortic flow rate ratio as a parameter. The results are compared with earlier works involving experimental data and found to be in reasonable qualitative agreement. The steady flow, shear stress and branch flow under the effect of a magnetic field have been discussed in detail.     

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     

Numerical analysis of finite amplitude motion of waves and a moored floating body under severe storm conditions

The motion of a moored floating body under the action of wave forces, which is influenced by fluid forces, shape of the floating body and mooring forces, should be analysed as a complex coupled motion system. Especially under severe storm conditions or resonant motion of the floating body it is necessary to consider finite amplitude motions of the waves, the floating body and the mooring lines as well as non-linear interactions of these finite amplitude motions. The problem of a floating body has been studied on the basis of linear wave theory by many researchers. However, the finite amplitude motion under a correlated motion system has rarely been taken into account. This paper presents a numerical method for calculating the finite amplitude motion when a floating body is moored by non-linear mooring lines such as chains and cables under severe storm conditions.     

THE MODELING AND ESTIMATION OF THE SHORT-RUN SUPPLY OF ROUNDWOOD UNDER NONLINEAR INCOME TAXES

This paper deals with the theory and estimation of the short-run supply of roundwood under nonlinear (progressive) income taxation. The theoretical part of the paper uses the two-period consumption harvesting model to introduce nonlinear income taxes into the forest management problem and presents a simple method to determine the optimal supply of roundwood in this case. In the empirical part of the paper, the model is specified in a way suitable for econometric analysis. The estimation method is the maximum likelihood method, and the estimation results seem to support the theoretical model. Finally, some results from model simulation imply that income taxation may have a nonnegligible impact on the supply of roundwood.     

非线性光学材料──水草酸铵

一水草酸铵（简称ＡＯＭ）是一种优良的电光晶体￣［１］。文献￣［１，２］报道了该晶体的晶体结构，化学式为（ＮＨ＿４）＿２Ｃ＿２Ｏ＿２·Ｈ＿２Ｏ，属正交晶系，Ｐ２２２空间群，晶胞参数ａ＝０．８０３５ｎｍ，ｂ＝１。０３１ｎｍ，ｃ＝０．３８０１ｎｍ，ｚ＝２。本文拟报道ＡＯＭ单晶生长，热化学分析，透过波段以及非线性光学性质的初步研究。     

FRICTIONAL DISSIPATION AND NONLINEAR BAROTROPIC INSTABILITY

Based on the nonlinear quasi-geostrophic barotropic vorticity equation with Ekman friction, the criteria for nonlinear barotropic stability of the zonal basic flows are derived using Serrin-Joseph energy approach and through total energy, total enstrophy and their linear combination, separately, in terms of variational principle. Since the new transformation for Euler equation is utilized, the estimation of eigenvalue is more accurate, and the previous results of the author are improved very well.