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 test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

Fractional Derivatives Applied to Phase-Space Reconstructions

The concept and application of phase-space reconstructions are reviewed. Fractional derivatives are then proposed for the purpose of reconstructing dynamics from a single observed time history. A procedure is presented in which the fractional derivatives of time series data are obtained in the frequency domain. The method is applied to the Lorenz system. The ability of the method to unfold the data is assessed by the method of global false nearest neighbors. The reconstructed data is used to compute recurrences and correlation dimensions. The reconstruction is compared to the commonly used method of delays in order to assess the choice of reconstruction parameters, and also the quality of results.     

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.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. 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.     

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.     

Synthesis and characterization of oxidized cellulose

Samples of oxidized cellulose (OC) with various carboxyl contents and degrees of crystallinity were obtained by the oxidation of native and mercerized cellulose with a solution of nitrogen(IV) oxide in CCl₄. A detailed characterization of these OC samples was performed. The effect of oxidation conditions (concentration of N₂O₄ in the solution and oxidation time) and starting cellulose material on OC characteristics (carboxyl, carbonyl and nitrogen content, degree of crystallinity and polymerization, surface area and swelling, and acidic properties) was investigated. Reactivity in the oxidation process was higher in mercerized cellulose than in native cellulose. The action of dilute solutions (10–15%) of N₂O₄ did not affect the degree of crystallinity of cellulose samples. Under these conditions, the oxidation took place mainly in amorphous regions and on the surface of crystallites. Oxidation in a concentrated (40%) N₂O₄ solution led to the destruction of crystallites, which increased the surface area and swelling of cellulose in water. The surface area and the swelling of OC samples increased with a decrease in the index of crystallinity. The acidic properties of OC were shown to increase with an increase of swelling in water. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4785–4791, 2004