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.     

Closed-form solutions for a rectangular layer with central crack under anti-plane shear stress

We consider the problem of determining the stress distributionin a finite rectangular elastic layer containing a Griffithcrack which is opened by internal shear stress acting alongthe length of the crack. The mode III crack is assumed to belocated in the middle plane of the rectangular layer. The followingtwo problems are considered: (A) the central crack is perpendicularto the two fixed lateral surfaces and parallel to the othertwo stress-free surfaces; (B) all the lateral surfaces of therectangular layer are clamped and the central crack is parallelto the two lateral surfaces. By using Fourier transformations,we reduce the solution of each problem to the solution of dualintegral equations with sine kernels and a weight function whichare solved exactly. Finally, we derive closed-form expressionsfor the stress intensity factor at the tip of the crack andthe numerical values for the stress intensity factor at theedges of the cracks are presented in the form of tables.     

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.     

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.     

An improved numerical simulator for multiphase flow in porous media

In basin modelling the thermodynamics of a multicomponent multiphase fluid flux are computationally too expensive when derived from an equation of state and the Gibbs equality constraints. In this article we present a novel implicit molar mass formulation technique using binary mixture thermodynamics. The two proposed solution methods, with and without cross derivative terms between components, are based on a preconditioned Newton‐GMRES scheme for each time‐step with analytical computation of the derivatives. These new algorithms reduce significantly the numerical effort for the computation of the molar masses, and we illustrate the behavior of these methods with numerical computations. 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.     

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