Effect of non-uniform magnetic field on crystal growth by floating-zone method in microgravity

The magnetic damping effect of the non-uniform magnetic field on the floating-zone crystal growth process in microgravity is studied by numerical simulation. The results show that the non-uniform magnetic field with designed configuration can effectively reduce the flow near the free surface and then in the melt zone. At the same time, the designed magnetic field can improve the impurity concentration non-uniformity along the solidification interface. The primary principles of the magnetic field configuration design are also discussed.

Comparison of the Dynamic Response of a Saturated Soil Between Two Models

By comparing the dynamic responses of saturated soil to Biot's and Yamamoto's models, the properties of the two models have be pointed out. First of all, an analysis has been made for energy loss of each model from the basic equations. Then the damping of elastic waves in coarse sand and fine sand with loading frequency and soil's parameters have been calculated and the representation of viscous friction and Coulomb friction in the two models has been concluded. Finally, the variations of loading wave damping and stress phase angles with water depth and soil's parameters have been obtained as loading waves range in ocean waves.

Source Model for Blasting Vibration

By analyzing and comparing the experimental data, the point source moment theory and the cavity theory, it is concluded that the vibrating signals away from the blasting explosive come mainly from the natural vibrations of the geological structures near the broken blasting area. The source impulses are not spread mainly by the inelastic properties (such as through media damping, as believed to be the case by many researchers) of the medium in the propagation pass, but by this structure. Then an equivalent source model for the blasting vibrations of a fragmenting blasting is proposed, which shows the important role of the impulse of the source's time function under certain conditions. For the purpose of numerical simulation, the model is realized in FEM, The finite element results are in good agreement with the experimental data.

Dynamic behavior of a cylindrical crack in a functionally graded interlayer under torsional loading

In this paper the problem of a cylindrical crack located in a functionally graded material (FGM) interlayer between two coaxial elastic dissimilar homogeneous cylinders and subjected to a torsional impact loading is considered. The shear modulus and the mass density of the FGM interlayer are assumed to vary continuously between those of the two coaxial cylinders. This mixed boundary value problem is first reduced to a singular integral equation with a Cauchy type kernel in the Laplace domain by applying Laplace and Fourier integral transforms. The singular integral equation is then solved numerically and the dynamic stress intensity factor (DSIF) is also obtained by a numerical Laplace inversion technique. The DSIF is found to rise rapidly to a peak and then reduce and tend to the static value almost without oscillation. The influences of the crack location, the FGM interlayer thickness and the relative magnitudes of the adjoining material properties are examined. It is found among others that, by increasing the FGM gradient, the DSIF can be greatly reduced.

Oxyhydrogen combustion and detonation driven shock tube

The performance of combustion driver ignited by multi-spark plugs distributed along axial direction has been analysed and tested. An improved ignition method with three circumferential equidistributed ignitors at main diaphragm has been presented, by which the produced incident shock waves have higher repeatability, and better steadiness in the pressure, temperature and velocity fields of flow behind the incident shock, and thus meets the requirements of aerodynamic experiment. The attachment of a damping section at the end of the driver can eliminate the high reflection pressure produced by detonation wave, and the backward detonation driver can be employed to generate high enthalpy and high density test flow. The incident shock wave produced by this method is well repeated and with weak attenuation. The reflection wave caused by the contracted section at the main diaphragm will weaken the unfavorable effect of rarefaction wave behind the detonation wave, which indicates that the forward detonation driver can be applied in the practice. For incident shock wave of identical strength, the initial pressure of the forward detonation driver is about 1 order of magnitude lower than that of backward detonation.

Numerical simulation of standing solitons and their interaction

Standing soliton was studied by numerical simulation of ifs governing equation, a cubic Schrodiger equation with a complex conjugate term, which was derived by Miles and was accepted. The value of linear damping in Miles equation was studied. Calculations showed that linear damping effects strongly on the formation of a standing soliton and Laedke and Spatschek stable condition is only a necessary condition, but not a sufficient one. The interaction of two standing solitons was simulated. Simulations showed that the interaction pattern depends on system parameters. Calculations for the different initial condition and its development indicated that a stable standing soliton can be fanned only for proper initial disturbance, otherwise the disturbance will disappear or develop into several solitons.

Experimental Study on Thermocapillary Convection

In the present paper, the experimental studies on thermocapillary convection are reviewed. The author's interest is mainly focused on the onset of oscillatory thermocapillary convection, the features of oscillatory flow pattern, and the critical Marangoni number related with temperature and free surface oscillation. The coordinated measurement in a microgravity environment of a drop shaft is also addressed.

Experimental Study on the Transition Process to the Oscillatory Thermocapillary Convections in a Floating Half Zone

The transition process of the thermocapillary convection from a steady and axisymmetric mode to the oscillatory mode in a liquid bridge with a fixed aspect ratio and varied volume ratio was studied experimentally. To ensure the surface tension to play an important role in the ground-based experiment, the geometrical configuration of the liquid bridge was so designed that the associated dynamic Bond number Bd ≈ 1. The velocity fields were measured by Particle Image Velocimetry (PIV) technique to effectively distinguish the different flow modes during the transition period in the experiments. Our experiments showed that as the temperature difference increased the slender and fat bridges presented quite different features on the evolution in their flow feature: for the former the thermocapillary convection transformed from a steady and axisymmetric pattern directly into an oscillatory one; but for the latter a transition flow status, characterized by an axial asymmetric steady convection, appeared before reaching the oscillatory mode. Experimental observations agree with the results of numerical simulations and it is obvious that the volume of liquid bridge is a sensitive geometric parameter. In addition, at the initial stage of the oscillation, for the former a rotating oscillatory convection with azimuthal wave number m = 1 was observed while for the latter a pulsating oscillatory pattern with azimuthal wave number m = 2 emerged, and then with further increase of the temperature difference, the pulsating oscillatory convection with azimuthal wave number m = 2 evolved into a rotating oscillatory pattern with azimuthal wave number m = 2.

Fourth order accurate compact scheme with group velocity control (GVC)

For solving complex flow field with multi-scale structure higher order accurate schemes are preferred. Among high order schemes the compact schemes have higher resolving efficiency. When the compact and upwind compact schemes are used to solve aerodynamic problems there are numerical oscillations near the shocks. The reason of oscillation production is because of non-uniform group velocity of wave packets in numerical solutions. For improvement of resolution of the shock a parameter function is introduced in compact scheme to control the group velocity. The newly developed method is simple. It has higher accuracy and less stencil of grid points.

A simulation model of a floating half zone

A simulation model of a floating half zone was suggested by steady numerical simulation and experiment respectively, in the previous papers [Q.S. Chen, W.R. Hu, Int. J. Mass Heat Transfer 40 (1997) 757; J.H. Han, Y. Ar, R. Zhou, W.R. Hue, Int. J. Mass Heat Transfer 40 (1997) 2671]. In the present paper, the simulation model is studied by using the method of unsteady and three-dimensional numerical simulation, and the transient process from steady convection to oscillatory convection is especially analyzed. Comparison of onsets of oscillation for both simulation model and the usual model were obtained, and the results show that the critical Marangoni number of the simulation model is obviously smaller than that of the usual model for the same slender liquid bridge. This implies that the usual model of a floating half zone gives a lower estimation on the onset of oscillation for floating zone convection.

Simulations of oil spread on a water surface

Numerical simulation of an oil slick spreading on still and wavy surfaces is described in this paper. The so-called sigma transformation is used to transform the time-varying physical domain into a fixed calculation domain for the water wave motions and, at the same time, the continuity equation is changed into an advection equation of wave elevation. This evolution equation is discretized by the forward time and central space scheme, and the momentum equations by the projection method. A damping zone is set up in front of the outlet boundary coupled with a Sommerfeld-Orlanski condition at that boundary to minimize the wave reflection. The equations for the oil slick are depth-averaged and coupled with the water motions when solving numerically. As examples, sinusoidal and solitary water waves, the oil spread on a smooth plane and on still and wavy water surfaces are calculated to examine the accuracy of simulating water waves by Navier-Stokes equations, the effect of damping zone on wave reflection and the precise structures of oil spread on waves.

Two bifurcation transitions of the floating half zone convection in a fat liquid bridge of larger Pr

The transient process of the thermocapillary convection was obtained for the large Pu floating half zone by using the method of three-dimensional and unsteady numerical simulation. The convection transits directly from steady and axisymmetric state to oscillatory flow for slender liquid bridge, and transits first from steady and axisymmetric convection to the steady and non-axisymmetric convection, then, secondly to the oscillatory convection for the fatter liquid bridge. This result implies that the volume of liquid bridge is not only a sensitive critical parameter for the onset of oscillation, but also relates to the new mechanism for the onset of instability in the floating half zone convection even in case of large Prandtl number fluid.

受水压振荡时饱和砂床的液化与变形特性

为考察饱和砂土地基受到水压振荡作用时的动力响应,根据两相连续介质理论,分析了液化发展、液化区厚度、颗粒速度等.研究结果表明,随着土强度和渗透性的降低与载荷强度的增加,液化发展速率加快;随着土强度的降低和渗透性、载荷强度的增加,变形量增加;如果土的初始分布不均匀,则在渗透性小的位置,应变和孔隙水压力将剧烈变化,这可以解释为什么在一定条件下土中会出现断裂和水层.

Instability from steady and axisymmetric to steady and asymmetric floating half zone convection in a fat liquid bridge of larger Prandtl number

The linear instability analysis of the present paper shows that the thermocapillary convection in a half floating zone of larger Prandtl number has a steady instability mode w(i) = 0 and m = 1 for a fat liquid bridge V = 1.2 with small geometrical aspect ratio A = 0.6. This conclusion is different from the usual idea of hydrothermal instability, and implies that the instability of the system may excite a steady and axial asymmetric state before the onset of oscillation in the ease of large Prandtl number.

Turbulence modeling for flows around convex features giving rapid eddy distortion

Reynolds averaged Navier-Stokes model performances in the stagnation and wake regions for turbulent flows with relatively large Lagrangian length scales (generally larger than the scale of geometrical features) approaching small cylinders (both square and circular) is explored. The effective cylinder (or wire) diameter based Reynolds number, ReW ≤ 2.5 × 103. The following turbulence models are considered: a mixing-length; standard Spalart and Allmaras (SA) and streamline curvature (and rotation) corrected SA (SARC); Secundov's νt-92; Secundov et al.'s two equation νt-L; Wolfshtein's k-l model; the Explicit Algebraic Stress Model (EASM) of Abid et al.; the cubic model of Craft et al.; various linear k-ε models including those with wall distance based damping functions; Menter SST, k-ω and Spalding's LVEL model. The use of differential equation distance functions (Poisson and Hamilton-Jacobi equation based) for palliative turbulence modeling purposes is explored. The performance of SA with these distance functions is also considered in the sharp convex geometry region of an airfoil trailing edge. For the cylinder, with ReW ≈ 2.5 × 103 the mixing length and k-l models give strong turbulence production in the wake region. However, in agreement with eddy viscosity estimates, the LVEL and Secundov νt-92 models show relatively little cylinder influence on turbulence. On the other hand, two equation models (as does the one equation SA) suggest the cylinder gives a strong turbulence deficit in the wake region. Also, for SA, an order or magnitude cylinder diameter decrease from ReW = 2500 to 250 surprisingly strengthens the cylinder's disruptive influence. Importantly, results for ReW ≪ 250 are virtually identical to those for ReW = 250 i.e. no matter how small the cylinder/wire its influence does not, as it should, vanish. Similar tests for the Launder-Sharma k-ε, Menter SST and k-ω show, in accordance with physical reality, the cylinder's influence diminishing albeit slowly with size. Results suggest distance functions palliate the SA model's erroneous trait and improve its predictive performance in wire wake regions. Also, results suggest that, along the stagnation line, such functions improve the SA, mixing length, k-l and LVEL results. For the airfoil, with SA, the larger Poisson distance function increases the wake region turbulence levels by just under 5%. © 2007 Elsevier Inc. All rights reserved.