Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<discrete cells was proposed and named the discrete fluid dynamics (DFD) algorithm. Analysis shows that the basic conservative equations for discrete cells are the Euler equations, NS- and diffusion parabolized (DP) NS equations. In this paper, a new multiscale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

Numerical Analyses of Discrete Gust Response for an Aircraft

Based on Navier-Stokes equations and structural and flight dynamic equations of motion, dynamic responses in vertical discrete gust flow perturbation are investigated for a supersonic transport model. A tightly coupled method was developed by subiterations between aerodynamic equations and dynamic equations of motion. First, under the assumption of rigid-body and single freedom of motion in the vertical plunging, the results of a direct-coupling method are compared with the results of quasi-steady model method. Then, gust responses for the one-minus-cosine gust profile arc analyzed with two freedoms of motion in plunging and pitching for the airplane configurations with and without the consideration of structural deformation.

Magnetic properties and magnetic domain structure of bulk glass forming Nd60Al10Fe20Co10 alloy

The transition from hard to soft magnetic behaviour with increasing quenching rate is shown for Nd60WAl10Fe20Co10 melt-spun ribbons with different thickness. Microstructure and magnetic domain structure of ribbons were studied by magnetic force microscopy (MFM). Particle sizes < 5 nm decreasing gradually with increasing quenching rate were deduced from topographic images which differ from large-scale magnetic domains with a periodicity of about 350 nm in all ribbons irrespective the coercivity. This indicates that the magnetic properties of the alloy are governed by interaction of small magnetic particles. It is concluded that the presence of short-range-ordered structures with a local ordering similar to the Al metastable Nd-Fe binary phase is responsible for the hard magnetic properties in samples subjected to relatively low quenching rate.

Domain structure of hard magnetic NdAlFeCo bulk metallic glass

Magnetic domain structure of hard magnetic Nd60Al10Fe20Co10 bulk metallic glass (BMG) has been studied by using magnetic force microscopy. In the magnetic force images it is shown that the exchange interaction type magnetic domains with a period of about 360 nm do exist in the BMG, which is believed to be associated with the appearance of hard-magnetic properties in this system. As the scale of the magnetic domain is much larger than the size of the short-range ordered atomic clusters existing in the BMG, it is believed that the large areas of magnetic contrast are actually a collection of a group of clusters aligned in parallel by strong exchange coupling interaction. After fully crystallization, the BMG exhibits paramagnetism. No obvious magnetic contrast is observed in the magnetic force images of fully crystallized samples, except for a small quantity of ferromagnetic crystalline phase with low coercivity and an average size of 900 nm.

Characterization of a Laser-Discrete Quenched Steel Substrate/Chromium System by Dissolving Coatings

A laser-discrete quenched steel (LDQS) substrate/as-deposited chromium (top high-contraction (HC) and underlying low-contraction (LC) chromium) system was investigated by dissolving coatings in order to reveal the mechanism that the service life of the coated parts is largely improved using the hybrid technique of laser pre-quenching plus chromium post-depositing. It was found that the surface characteristics of the substrate, LC and HC chromium layer can be simultaneously revealed owing to the dissolution edge effect of chromium coatings. Moreover, the periodical gradient morphologies of the LDQS substrate are clearly shown: the surfaces of laser transformation-hardened regions are rather smooth; a lot of fine micro-holes exist in the transition zones; there are many micro-dimples in the original substrate. Furthermore, the novel method of dissolving coatings with sharp interfaces may be used to reveal the structural features of a substrate/coating system, explore the effect of the substrate on the initial microstructure and morphologies of coatings, and check the quality of the coated-parts.

Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.

Domain structure of a Nd60Al10Fe20Co10 bulk metallic glass

Magnetic domain structure of Nd60Al10Fe20Co10 bulk metallic glass (BMG) has been studied by using magnetic-force microscopy. In the magnetic-force images it is shown that the exchange-interaction-type magnetic domains with a period of about 360 nm do exist in the BMG, which is believed to be associated with the appearance of hard-magnetic properties in this system. The existence of the large-scale domains demonstrates that the magnetic moments of a great deal of short-scale ordered atomic clusters in the BMG have been aligned by exchange coupling. Annealing at 715 K leads to partial crystallization of the BMG. However, the exchange coupling is stronger in the annealed sample, which is considered to arise from the increase of transition-metal concentration in the amorphous phase due to the precipitation of Nd crystalline phase.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Size Effects in Micro-bending Deformation of MEMS Devices Based on the Discrete Dislocation Theory

In the present research, the discrete dislocation theory is used to analyze the size effect phenomena for the MEMS devices undergoing micro-bending load. A consistent result with the experimental one in literature is obtained. In order to check the effectiveness to use the discrete dislocation theory in predicting the size effect, both the basic version theory and the updated one are adopted simultaneously. The normalized stress-strain relations of the material are obtained for different plate thickness or for different obstacle density. The prediction results are compared with experimental results.

A time domain computation method for dynamic behavior of mooring system

A quasi-steady time domain method is developed for the prediction of dynamic behavior of a mooring system under the environmental disturbances, such as regular or irregular waves, winds and currents. The mooring forces are obtained in a static sense at each instant. The dynamic feature of the mooring cables can be obtained by incorporating the extended 3-D lumped-mass method with the known ship motion history. Some nonlinear effects, such as the influence of the instantaneous change of the wetted hull surface on the hydrostatic restoring forces and Froude-Krylov forces, are included. The computational results show a satisfactory agreement with the experimental ones.

Degradation failure features of chromium-plated gun barrels with a laser-discrete-quenched substrate

The effect of substrate laser-discrete quenching on the degradation failure of chromium-plated gun barrels was metallurgically investigated. The results show that substrate laser-discrete quenching changes the failure patterns of chromium coatings during firing, and some periodic through-thickness cracks in the fired chromium coatings are justly located at original substrate zones between two adjacent laser-quenched tracks. Moreover, chromium coatings and the laser-quenched zones on the substrate are simultaneously degraded in microstructure and property during firing. Furthermore, the periodic structure of the laser-discrete-quenched steel (LDQS) substrate near the breech remains after firing, and the hardness of the fired laser-quenched zones is still higher than that of original substrates. The specific failure features were utilized to illustrate the mechanism of the extended service life of chromium-plated gun barrels with the LDQS substrate. (c) 2007 Elsevier B.V All rights reserved.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.