PERIODIC ENVIRONMENTS AND PERIODIC HARVESTING

ABSTRACT. We consider autonomous population models under periodic harvesting and population models in periodic environments and seek conditions under which there is an asymptotically stable periodic solution.     

Study of the dynamic structure factor of strong glasses

The dynamic structure factors (DSFs) of several strong glasses (SF6, SF10, BK7, SUPRASIL) measured by Brillouin light scattering spectroscopy are reported. Spectra have been collected, at and above room temperature, at two scattering angles, θ=90° and θ=180° corresponding to exchanged wavevector q values ranging from 0.0256 to 0.0448 nm⁻¹. In particular we find that the isotropic spectral lineshapes are in all cases well described by the simple hydrodynamic theory of an amorphous solid. The width of the Brillouin peaks are found to be consistent with the predicted q² dependence at both investigated temperatures. This damping is however found to account only partially for the strong asymmetry of the Brillouin line clearly visible on a logarithmic intensity scale. As a matter of fact there is an excess intensity in the very low frequency plateau underlying the central component. The height of this plateau and hence the entire lineshape is well reproduced if a relaxation process is taken into account in the hydrodynamic equations. Owing to the intense elastic scattering we are able to determine unambiguously only the ratio between amplitude and characteristic time of this process which quantifies the sound dispersion to be of the order of a few percent in all samples. The temperature dependence of the parameters indicates that this relaxation cannot be attributed to thermally activated relaxation phenomena. These general findings favorably compare with molecular dynamics simulation results on similar systems.     

Combined supercritical-fluid chromatography/mass spectrometry in the analysis of diuron in plasma using on-line phase-system switching

On-line sample pretreatment by means of the phase-system switching approach is an interesting technique for the analysis of aqueous samples, e.g., plasma, by means of supercritical-fluid chromatography. In order to analyse plasma samples the following analytical procedure is used. The plasma sample is injected on to a short precolumn, which is washed with water and subsequently dried with nitrogen. Next, the solutes are desorbed with the supercritical mobile phase, analysed with packed-column supercritical-fluid chromatography and detected with either a UV detector or a mass spectrometer, equipped with a moving-belt interface. The herbicide diuron is selected as a test compound to study the feasibility of this approach. Using a selective detector the procedure is sufficiently sensitive to detect diuron in plasma, but not appropriate to detect the diuron metabolites in a post-mortem plasma sample. These have been identified with liquid chromatography/mass spectrometry. The detection limit of diuron in plasma using the procedure described is about 30 ng/mL.     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

A transform method for linear evolution PDEs on a finite interval

We study initial boundary value problems for linear scalar evolutionpartial differential equations, with spatial derivatives ofarbitrary order, posed on the domain {t > 0, 0 < x <L}. We show that the solution can be expressed as an integralin the complex k-plane. This integral is defined in terms ofan x-transform of the initial condition and a t-transform ofthe boundary conditions. The derivation of this integral representationrelies on the analysis of the global relation, which is an algebraicrelation defined in the complex k-plane coupling all boundaryvalues of the solution. For particular cases, such as the case of periodic boundaryconditions, or the case of boundary value problems for even-orderPDEs, it is possible to obtain directly from the global relationan alternative representation for the solution, in the formof an infinite series. We stress, however, that there existinitial boundary value problems for which the only representationis an integral which cannot be written as an infinite series.An example of such a problem is provided by the linearized versionof the KdV equation. Similarly, in general the solution of odd-orderlinear initial boundary value problems on a finite intervalcannot be expressed in terms of an infinite series.     

A numerical and experimental investigation of the interactions between a non‐uniform planar array of incompressible free jets

The interaction between multiple incompressible air jets has been studied numerically and experimentally. The numerical predictions have been first validated using experimental data for a single jet configuration. The spreading features of five unequal jets in the configuration of one larger central jet surrounded by four smaller equi‐distant jets, have been studied, for different lateral spacing ratios of 1.5, 2.0 and 2.5 and a central jet Reynolds number of 1.24×105 (corresponding to a Mach number of 0.16). Flow of five equal jets has also been simulated, for the sake of comparison. The jet interactions commence at an axial distance of about 3–4 diameters and complete by an axial distance of about 10 diameters for the lowest spacing ratio of 1.5. For larger spacing ratios, the length required for the start and completion of jet interaction increase. Peripheral jets bend more towards the central jet and merge at a smaller distance, when their sizes are smaller than that of the central jet. The entrainment ratio for multiple jets is higher than that for a single jet. Excellent agreement is observed between the experimental data and theoretical predictions for both mean flow field and turbulent quantities, at regions away from the jet inlet. The potential core length and initial jet development, however, are not predicted very accurately due to differences in the assumed and actual velocity profiles at the jet inlet. Copyright © 2004 John Wiley & Sons, Ltd.     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.