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.     

Validation of a thermal decomposition mechanism of formaldehyde by detection of CH2O and HCO behind shock waves

The thermal decomposition of formaldehyde was investigated behind shock waves at temperatures between 1675 and 2080 K. Quantitative concentration time profiles of formaldehyde and formyl radicals were measured by means of sensitive 174 nm VUV absorption (CH₂O) and 614 nm FM spectroscopy (HCO), respectively. The rate constant of the radical forming channel (1a), CH₂O + M → HCO + H + M, of the unimolecular decomposition of formaldehyde in argon was measured at temperatures from 1675 to 2080 K at an average total pressure of 1.2 bar, k_1a = 5.0 × 10¹⁵ exp(‐308 kJ mol^?1/RT) cm³ mol^?1 s^?1. The pressure dependence, the rate of the competing molecular channel (1b), CH₂O + M → H₂ + CO + M, and the branching fraction β = k_1a/(kA_1a + k_1b) was characterized by a two‐channel RRKM/master equation analysis. With channel (1b) being the main channel at low pressures, the branching fraction was found to switch from channel (1b) to channel (1a) at moderate pressures of 1–50 bar. Taking advantage of the results of two preceding publications, a decomposition mechanism with six reactions is recommended, which was validated by measured formyl radical profiles and numerous literature experimental observations. The mechanism is capable of a reliable prediction of almost all formaldehyde pyrolysis literature data, including CH₂O, CO, and H atom measurements at temperatures of 1200–3200 K, with mixtures of 7 ppm to 5% formaldehyde, and pressures up to 15 bar. Some evidence was found for a self‐reaction of two CH₂O molecules. At high initial CH₂O mole fractions the reverse of reaction (6), CH₂OH + HCO ? CH₂O + CH₂O becomes noticeable. The rate of the forward reaction was roughly measured to be k₆ = 1.5 × 10¹³ cm³ mol^?1 s^?1. © 2004 Wiley Periodicals, Inc. Int J Chem Kinet 36: 157–169 2004     

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     

Model based online diagnosis of unbalance and transverse fatigue crack in rotor systems

A model based technique for online identification of malfunctions in rotor systems is discussed. Presence of fault changes the dynamic behavior of the system. This change is taken into account by equivalent loads acting on the undamaged system model. Equivalent loads are fictitious forces and moments acting on the undamaged system model, which generate a dynamic behavior identical to that of the real damaged system. The mathematical representation of equivalent loads is referred to as Fault Model. The work focuses on developing a fault model for a transverse fatigue crack in shaft and testing it through simulated studies. The basic principle of the technique is validated for unbalance identification, through numerical simulations as well as by experiments on a real rotor system.     

Time resolved spectroscopic NMR imaging using hyperpolarized Xe

We have visualized the melting and dissolution processes of xenon (Xe) ice into different solvents using the methods of nuclear magnetic resonance (NMR) spectroscopy, imaging, and time resolved spectroscopic imaging by means of hyperpolarized ¹²⁹Xe. Starting from the initial condition of a hyperpolarized solid Xe layer frozen on top of an ethanol (ethanol/water) ice block we measured the Xe phase transitions as a function of time and temperature. In the pure ethanol sample, pieces of Xe ice first fall through the viscous ethanol to the bottom of the sample tube and then form a thin layer of liquid Xe/ethanol. The xenon atoms are trapped in this liquid layer up to room temperature and keep their magnetization over a time period of 11 min. In the ethanol/water mixture (80 vol%/20%), most of the polarized Xe liquid first stays on top of the ethanol/water ice block and then starts to penetrate into the pores and cracks of the ethanol/water ice block. In the final stage, nearly all the Xe polarization is in the gas phase above the liquid and trapped inside the pores. NMR spectra of homogeneous samples of pure ethanol containing thermally polarized Xe and the spectroscopic images of the melting process show that very high concentrations of hyperpolarized Xe (about half of the density of liquid Xe) can be stored or delivered in pure ethanol.     

Dynamics of capture and fusion in heavy ion collisions

The effect of entrance channel on decrease of the complete fusion cross sections and on the yield of reaction products are associated with the quasifission which depends on the mass asymmetry and shell structure of colliding nuclei. In reactions of massive projectile and target nuclei, the competition between complete fusion and quasifission appears at the stage of compound nucleus formation, in addition to the increase of the fission probability. It is shown that the yield of quasifission products may be symmetric or asymmetric in dependence on peculiarities of shell structure of reaction fragments. Marima of mass or charge distributions are connected with the peculiarities of shell structure of reaction fragments.