Atomistic phenomena in dense fluid shock waves

The shock structure problem is one of the classical problems of fluid mechanics and at least for non-reacting dilute gases it has been considered essentially solved. Here we present a few recent findings, to show that this is not the case. There are still new physical effects to be discovered provided that the numerical technique is general enough to not rule them out a priori. While the results have been obtained for dense fluids, some of the effects might also be observable for shocks in dilute gases.     

Approximation of transient temperatures in complex geometries using fractional derivatives

It is shown that time-dependent temperatures in a transient, conductive system can be approximately modeled by a fractional-order differential equation, the order of which depends on the Biot number. This approximation is particularly suitable for complex shapes for which a first-principles approach is too difficult or computationally time-consuming. Analytical solutions of these equations can be written in terms of the Mittag-Leffler function. The approximation is especially useful if a suitable fractional-order controller is to be designed for the system.     

Analytical formulas for the velocity field induced by an infinitely thin vortex ring

Three different analytical solutions are presented for a potential vortex ring using three different streamfunctions. Verification studies confirm that all three approaches are valid. It is found that the solution obtained using the Biot–Savart law is the most efficient method due to its simplicity. It is shown that all analytical results are accurate to within machine accuracy and sample calculations are included. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

Grid‐free surface vorticity method applied to flow induced vibration of flexible cylinders

In order to study cross flow induced vibration of heat exchanger tube bundles, a new fluid–structure interaction model based on surface vorticity method is proposed. With this model, the vibration of a flexible cylinder is simulated at Re=2.67 × 10⁴, the computational results of the cylinder response, the fluid force, the vibration frequency, and the vorticity map are presented. The numerical results reproduce the amplitude‐limiting and non‐linear (lock‐in) characteristics of flow‐induced vibration. The maximum vibration amplitude as well as its corresponding lock‐in frequency is in good agreement with experimental results. The amplitude of vibration can be as high as 0.88D for the case investigated. As vibration amplitude increases, the amplitude of the lift force also increases. With enhancement of vibration amplitude, the vortex pattern in the near wake changes significantly. This fluid–structure interaction model is further applied to simulate flow‐induced vibration of two tandem cylinders and two side‐by‐side cylinders at similar Reynolds number. Promising and reasonable results and predictions are obtained. It is hopeful that with this relatively simple and computer time saving method, flow induced vibration of a large number of flexible tube bundles can be successfully simulated. Copyright © 2004 John Wiley & Sons, Ltd.     

Damage and fracture evaluation of granular composite materials by digital image correlation method

This paper presents the applications of digital image correlation technique to the mesoscopic damage and fracture study of some granular based composite materials including steelfiber reinforced concrete, sandstone and crystal-polymer composite. The deformation fields of the composite materials resulted from stress localization were obtained by the correlation computation of the surface images with loading steps and thus the related damage prediction and fracture parameters were evaluated. The correlation searching could be performed either directly based on the gray levels of the digital images or from the wavelet transform (WT) coefficients of the transform spectrum. The latter was developed by the authors and showed higher resolution and sensitivity to the singularity detection. Because the displacement components came from the rough surfaces of the composite materials without any coats of gratings or fringes of optical interferometry, both surface profiles and the deformation fields of the composites were visualized which was helpful to compare each other to analyze the damage of those heterogeneous materials. The project supported by the National Natural Science Foundation of China (10125211 and 10072002), the Scientific Committee of Yunnan Province for the Program of Steel Fiber Reinforced Concrete, and the Institute of Chemical Materials, CAEP at Mianyang     

Comments on ‘A New Model for Viscous Dissipation in Porous Media Across a Range of Permeability Values’, Transport in Porous Media 53, 117–122, 2003

Transport in Porous Media -     

Properties of backscattered particles produced by protons of energy above 1 TeV in an iron absorber

Backscattered-particle production is studied by means of a detailed simulation of cascade processes in a dense medium. The energy dependence of the albedo and the spatial and angular distributions of various components of this flux are analyzed.     

Resonant destabilization of phonon microwave stimulated emission in an acoustic quantum oscillator (phaser) with periodically modulated pumping

Two qualitatively different types of resonant destabilization of phonon stimulated emission (SE) are discovered in experiments where a 9-GHz multimode ruby laser is periodically modulated (the electromagnetic pump frequency is 23 GHz). In the case of deep pump modulation at low modulation frequencies (ω=70–200 Hz, where ω_m is the modulation frequency), a fast random alternation of microwave phonon SE modes is observed. This destabilization range corresponds to relaxation resonance in optical lasers. Outside the relaxation resonance range (at ω_m≈10 Hz), the other type of resonant destabilization of stationary phonon SE is observed. This destabilization shows up as very slow regular self-detunings of the microwave SE spectra. The period of such self-organized motions depends significantly on ω_m and changes by several orders of magnitude when ω_m varies within several percent. The second type of SE resonant destabilization is explained in terms of antiphase energy exchange between modes in a modulated phaser.