Fractional Brownian Heavy Traffic Approximations of Multiclass Feedforward Queueing Networks

We consider multiclass feedforward queueing networks under first in first out and priority service disciplines driven by long-range dependent arrival and service time processes. We show that in critical loading the normalized workload, queue length and sojourn time processes can converge to a multi-dimensional reflected fractional Brownian motion. This weak heavy traffic approximation is deduced from a deterministic pathwise approximation of the network behavior close to constant critical load in terms of the solution of a Skorokhod problem. Since we model the doubly infinite time interval, our results directly cover the stationary case.AMS subject classification: primary 90B15, secondary 60K25, 68M20     

A thermodynamic approach with internal variables using Lagrange formalism. Part I: General framework

We present some reflections on the application of the Lagrangian formalism for continuous media locally uniform subjected to internal irreversible evolutions. The Lagrangian density, defined as the time derivative of a non-equilibrium thermodynamic potential, [Thermodynamics of Relaxation Processes using Internal variables within a Lagrange-formalism. P. Germain’s Anniversary Volume 2000. Contiuum Thermomechanics: the Art and Science of Modeling Matter’s Behaviour, 2000], contains all the symmetry properties of the system. The generalised Lagrange co-ordinates correspond to the state and internal variables of the time derivative of the generalised Gibbs potential. The latter being used within the framework of the De Donder’s method, must also account for the memory effect of the physical medium.This first part is devoted to the thermodynamic framework called the distribution of non-linear relaxations approach (DNLR) developed by C. Cunat on the basis of the generalised Gibbs’ relation.     

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.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. 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.     

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.     

Modeling Media with Oriented Structures

Heterogeneity is typically the result of space variability of soil parameters at different scales. Soil anisotropy may be defined as the spatial persistence in some direction only, across coarse-grid elements, of heterogeneous structures with different characteristic lengths in different directions. One can account for the effect of these structures by upscaling soil properties. Analyzing flow in a strongly anisotropic structured soil at different scales evidences how transverse dispersion reduces to a subscale process, leading to mixing within the conductive structures.     

On the hysteresis of vortex shedding in a periodically-varying flow

Efforts are made to explore the hysteresis characteristics of vortex shedding in a pipe flow, whose velocity varies periodically in time. Results obtained show that during acceleration of the flow, the vortex strength tends to be stronger, whereas during deceleration of the flow, the situation is reversed. As reconstructed from the velocity signals measured at a point in the flow field, the shed vortex arrays appear to possess uneven vortex strengths in response to periodically-varying incoming flows. Furthermore, in the hysteresis range, the streamwise spacings between the vortices appear to be unequal.     

Negative luminescence in semiconductors: a retrospective view

Twenty-five years ago, we introduced the phenomenon of negative luminescence (NL) into semiconductor physics. This paper provides an overview of work conducted to develop this fundamental concept. Initially, we consider the first-principle approach to radiation interaction with basic matter and the major properties of NL. Then we describe the problems of NL direct measurements in homogeneous materials and structures. Finally, we emphasize the use of NL approach in applications involving devices for infrared (IR) wavelength (3–12 μm) high-temperature (300–400 K) optoelectronics. Our subjects will include NL IR emitting diodes, radiative coolers, IR dynamic scene simulators, light up-conversion devices, and the Stealth effect in IR.     

The exponential power law: partial wetting kinetics and dynamic contact angles

An equation for the kinetics of partial drop spreading is proposed. This equation was empirically derived from experimental data for the spreading kinetics of partially wetting liquids in terms of the wet area versus time. The equation has the form of an exponential power law (EPL), and transforms into the well-known power law for complete wetting, when the equilibrium contact angle approaches zero. The EPL fits very well available experimental data. To lend additional support to the validity of this generalized equation, it will be demonstrated that when it is transformed to present the dynamic contact angle (DCA), it fits very well DCA experimental data for other wetting processes, such as capillary flow and tape coating.