An Improved Ce/Se Scheme And Its Application To Detonation Propagation

An improved two-dimensional space-time conservation element and solution element ( CE/ SE) method with second-order accuracy is proposed, examined and extended to simulate the detonation propagations using detailed chemical reaction models. The numerical results of planar and cellular detonation are compared with corresponding results by the Chapman-Jouguet theory and experiments, and prove that the method is a new reliable way for numerical simulations of detonation propagation.

Case Study: Influence Of Morphological Changes On Flooding In Jingjiang River

The middle reach of the Yangtze River, customarily called the Jingjiang River, together with its diversion channels and Dongting Lake, form a large complicated drainage system. In the last five decades, significant geomorphological changes have occurred in the drainage system, including the shrinkage of diversion channels, contraction of Dongting Lake, changes in the rating curve at the Luoshan station, and cutoffs of the lower Jingjiang River. These changes are believed to be the cause of the occurrence of abnormal floods in the Jingjiang River. Qualitative analyses suggest that the first three factors aggravate the flood situation in the lower Jingjiang River, while the last factor seems beneficial for flood prevention. To quantitatively evaluate these conclusions, a finite-volume numerical model was constructed. A series of numerical simulations were carried out to test the individual and combined effects of the aforementioned four factors, and these simulations showed that high flood stages in the Jingjiang River clearly are related to the geomorphological changes.

Thermocapillary Flows In Liquid Bridges Of Molten Tin With Small Aspect Ratios

Numerical simulations were conducted to study thermocapillary flows in short half-zone liquid bridges of molten tin with Prandtl number Pr = 0.009, under ramped temperature difference. The spatio-temporal structures in the thermocapillary flows in short half-zone liquid bridges with aspect ratios As = 0.6, 0.8, and 1.0 were investigated. The first critical Marangoni numbers were compared with those predicted by linear stability analyses (LSA). The second critical Marangoni numbers for As = 0.6 and 0.8 were found to be larger than that for As = 1.0. The time evolutions of the thermocapillary flows exhibited unusual features such as a change in the azimuthal wave number during the three-dimensional stationary (non-oscillating) flow regime, a change in the oscillation mode during the three-dimensional oscillatory flow regime, and the decreasing and then increasing of amplitudes in a single oscillation mode. The effects of the ramping rate of the temperature difference on the flow modes and critical conditions were studied as well. In this paper, the experimental observability of the critical conditions was also discussed. (C) 2008 Elsevier Inc. All rights reserved.

A Pressure-Correction Method And Its Applications On An Unstructured Chimera Grid

In this paper, an unstructured Chimera mesh method is used to compute incompressible flow around a rotating body. To implement the pressure correction algorithm on unstructured overlapping sub-grids, a novel interpolation scheme for pressure correction is proposed. This indirect interpolation scheme can ensure a tight coupling of pressure between sub-domains. A moving-mesh finite volume approach is used to treat the rotating sub-domain and the governing equations are formulated in an inertial reference frame. Since the mesh that surrounds the rotating body undergoes only solid body rotation and the background mesh remains stationary, no mesh deformation is encountered in the computation. As a benefit from the utilization of an inertial frame, tensorial transformation for velocity is not needed. Three numerical simulations are successfully performed. They include flow over a fixed circular cylinder, flow over a rotating circular cylinder and flow over a rotating elliptic cylinder. These numerical examples demonstrate the capability of the current scheme in handling moving boundaries. The numerical results are in good agreement with experimental and computational data in literature. (C) 2007 Elsevier Ltd. All rights reserved.

Experimental Verification And Theoretical Analysis Of The Relationships Between Hardness, Elastic Modulus, And The Work Of Indentation

The relationship between hardness (H), reduced modulus (E-r), unloading work (W-u), and total work (W-t) of indentation is examined in detail experimentally and theoretically. Experimental study verifies the approximate linear relationship. Theoretical analysis confirms it. Furthermore, the solutions to the conical indentation in elastic-perfectly plastic solid, including elastic work (W-e), H, W-t, and W-u are obtained using Johnson's expanding cavity model and Lame solution. Consequently, it is found that the W-e should be distinguished from W-u, rather than their equivalence as suggested in ISO14577, and (H/E-r)/(W-u/W-t) depends mainly on the conical angle, which are also verified with numerical simulations. (C) 2008 American Institute of Physics.

Comparative Study Of The Influence Of Natural Convection On Directional Solidification Of Al-3.5 Wt% Ni And Al-7 Wt% Si Alloys

We present numerical simulations of thermosolutal convection for directional solidification of Al-3.5 wt% Ni and Al-7 wt% Si. Numerical results predict that fragmentation of dendrite arms resulting from dissolution could be favored in Al-7 wt% Si, but not in Al-3.5 wt% Ni. Corresponding experiments are in qualitative agreement with the numerical predictions. Distinguishing the two fragmentation mechanisms, namely dissolution and remelting, is critical during experiments on earth, when fluid flow is dominant. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Time Correlations Of Pressure In Isotropic Turbulence

The time correlations of pressure modes in stationary isotropic turbulence are investigated under the Kraichnan and Tennekes "random sweeping" hypothesis. A simple model is obtained which predicts a universal form for the time correlations. It implies that the decorrelation process of pressure fluctuations in time is mainly dominated by the sweeping velocity, and the pressure correlations have the same decorrelation time scales as the velocity correlations. These results are verified using direct numerical simulations of isotropic turbulence at two moderate Reynolds numbers; the mode correlations collapse to the universal form when the time separations are scaled by wavenumber times the sweeping velocity, and the ratios of the correlation coefficients of pressure modes to those of velocity modes are approximately unity for the entire range of time separation. (c) 2008 American Institute of Physics.

Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed. to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number.

Space experiments of thermocapillary convection in two-liquid layers

In 1999, the space experiments on the Marangoni convection and thermocapillary convection in a system of two immiscible liquid layers in microgravity environment were conducted on board the Chinese scientific satellite SJ-5. A new system of two-layer liquids such as FC-70 liquid and paraffin was used successfully, with the paraffin melted in the space. Two different test-cells are subjected to a temperature gradient perpendicular or parallel to the interface to study the Marangoni convection and thermocapillary convection, respectively. The experimental data obtained in the first Chinese space experiment of fluid are presented. Two-dimensional numerical simulations of thermocapillary convections are carried out using SIMPLEC method A reasonable agreement between the experimental investigation and the numerical results is obtained.

Influence of threshold stress on the estimation of the Weibull statistics

The influence of threshold stress on the estimation of the Weibull statistics is discussed in terms of the Akaike information criterion. Numerical simulations show that, if sample data are limited in number and threshold stress is not too large, the two-parameter Weibull distribution is still a preferred choice. For example, the fit of strength data of glass and ceramics to the two- and three-parameter Weibull distributions is compared.

Coupling effects of heterogeneity and stress fluctuation on rupture

The coupling of mesoscopic strength distribution and stress fluctuation on damage evolution and rupture are examined. The numerical simulations show that there is only weak stress fluctuation at the initial damage stage when the mean field approximation is in effect. As the damage fraction becomes larger than the threshold value, the fluctuation is amplified significantly, and damage localization appears. The coupling between stress fluctuation, disordered heterogeneity and the damage localization may play an essential role in catastrophic rupture. (C) 2003 Elsevier Ltd. All rights reserved.

On dispersion-controlled principles for non-oscillatory shock-capturing schemes

The role of dispersions in the numerical solutions of hydrodynamic equation systems has been realized for long time. It is only during the last two decades that extensive studies on the dispersion-controlled dissipative (DCD) schemes were reported. The studies have demonstrated that this kind of the schemes is distinct from conventional dissipation-based schemes in which the dispersion term of the modified equation is not considered in scheme construction to avoid nonphysical oscillation occurring in shock wave simulations. The principle of the dispersion controlled aims at removing nonphysical oscillations by making use of dispersion characteristics instead of adding artificial viscosity to dissipate the oscillation as the conventional schemes do. Research progresses on the dispersion controlled principles are reviewed in this paper, including the exploration of the role of dispersions in numerical simulations, the development of the dispersion-controlled principles, efforts devoted to high-order dispersion-controlled dissipative schemes, the extension to both the finite volume and the finite element methods, scheme verification and solution validation, and comments on several aspects of the schemes from author's viewpoint.

Synchronization of modified Chua's circuit with x vertical bar x vertical bar function

This paper considers the chaos synchronization of the modified Chua's circuit with x vertical bar x vertical bar function. We firstly show that a couple of the modified Chua systems with different parameters and initial conditions can be synchronized using active control when the values of parameters both in drive system and response system are known aforehand. Furthermore, based on Lyapunov stability theory we propose an adaptive active control approach to make the states of two identical Chua systems with unknown constant parameters asymptotically synchronized. Moreover the designed controller is independent of those unknown parameters. Numerical simulations are given to validate the proposed synchronization approach.

The degree of predictability of earthquakes in several regions of China: statistical analysis of historical data

The stress release model, a stochastic version of the elastic-rebound theory, is applied to the historical earthquake data from three strong earthquake-prone regions of China, including North China, Southwest China, and the Taiwan seismic regions. The results show that the seismicity along a plate boundary (Taiwan) is more active than in intraplate regions (North and Southwest China). The degree of predictability or regularity of seismic events in these seismic regions, based on both the Akaike information criterion (AIC) and fitted sensitivity parameters, follows the order Taiwan, Southwest China, and North China, which is further identified by numerical simulations. (c) 2004 Elsevier Ltd. All rights reserved.