Rossby Vortex Simulation on a Paraboloidal Coordinate System Using the Lattice Boltzmann Method

In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc,, can be expected.

Direct numerical simulation of transition and turbulence in compressible mixing layer

The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach number Mc = 0.8 and Re = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of Lambda-vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.

Determination of Proper Processing Conditions of Material Surface Heat Treatment Through Finite Element Method

A two-dimensional model has been developed based on the experimental results of stainless steel remelting with the laminar plasma technology to investigate the transient thermo-physical characteristics of the melt pool liquids. The influence of the temperature field, temperature gradient, solidification rate and cooling rate on the processing conditions has been investigated numerically. Not only have the appropriate processing conditions been determined according to the calculations, but also they have been predicted with a criterion established based on the concept of equivalent temperature area density (ETAD) that is actually a function of the processing parameters and material properties. The comparison between the resulting conditions shows that the ETAD method can better predict the optimum condition.

Plot Erosion Model Using Gray Relational Analysis Method

The main factors affecting interrill erosion-including runoff discharge, rainfall intensity, mean flow velocity, and slope gradient-were analyzed by using a gray relational analysis. An equation for interrill erosion was derived by coupling this analysis with dimensional and regression analyses. The values of erosion rates predicted by this equation were in good agreement with experimental observations.

Effective elastic moduli of two-dimensional solids with multiple cracks

An accurate method which directly accounts for the interactions between different microcracks is used for analyzing the elastic problem of multiple cracks solids. The effective elastic moduli for randomly oriented cracks and parallel cracks are evaluated for the representative volume element (RVE) with microcracks in infinite media. The numerical results are compared with those from various micromechanics models and experimental data. These results show that the present method is simple and provides a direct and efficient approach to dealing with elastic solids containing multiple cracks.

Iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact

We present in this paper an iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact. To retain the advantage of integration on an element-by-element basis which is at the heart of modern hydrocodes, we suggest that the first step should be to solve for accelerations at an advanced time step by using the lumped mass approach, then iterate using a consistent mass matrix to improve the estimate. Examples are given to show the improved resolution with the new method.

An Energy-Based Method for Analyzing Instrumented Spherical Indentation Experiments

Using dimensional analysis and finite element calculation, we studied spherical indentation in elastic-plastic solids with work hardening. We report two previously unknown relationships between hardness, reduced modulus, indentation depth, indenter radius, and work of indentation. These relationships, together with the relationship between initial unloading stiffness and reduced modulus, provide an energy-based method for determining contact area, reduced modulus, and hardness of materials from instrumented spherical indentation measurements. This method also provides a means for calibrating the effective radius of imperfectly shaped spherical indenters. Finally, the method is applied to the analysis of instrumented spherical indentation experiments on copper, aluminum, tungsten, and fused silica.

Interactions of Penny-Shaped Cracks in Three-Dimensional Solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is developed in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.

Pressure Controlled Self-Assembly of High Quality Three-Dimensional Colloidal Photonic Crystals

A concise pressure controlled isothermal heating vertical deposition (PCIHVD) method is developed, which provides an optimal growing condition with better stability and reproducibility for fabricating photonic crystals (PCs) without the limitation of colloidal sphere materials and sizes. High quality PCs are fabricated with PCIHVD from polystyrene spheres with diameters ranging from 200 nm to 1 mu m. The deep photonic band gap and steep photonic band edge of the samples are most favorable for realizing ultrafast optical devices, photonic chips, and communications. This method makes a meaningful advance in the quality and diversity of PCs and greatly promotes their wide applications.

Analysis of two-dimensional finite solids with microcracks

A general method is presented for solving the plane elasticity problem of finite plates with multiple microcracks. The method directly accounts for the interactions between different microcracks and the effect of outer boundary of a finite plate. Analysis is based on a superposition scheme and series expansions of the complex potentials. By using the traction-free conditions on each crack surface and resultant forces relations along outer boundary, a set of governing equations is formulated. The governing equations are solved numerically on the basis of a boundary collocation procedure. The effective Young's moduli for randomly oriented cracks and parallel cracks are evaluated for rectangular plates with microcracks. The numerical results are compared with those from various micromechanics models and experimental data. These results show that the present method provides a direct and efficient approach to deal with finite solids containing multiple microcracks.

Study of Growth Mechanism of Lysozyme Crystal by Batch Crystallization Method

The lysozyme crystals were made by batch crystallization method and the distribution of aggregate in solution were measured by dynamic light scattering. The results showed that the dimension of aggregate increased with the increase of the concentration of lysozyme and NaCl, lysozyme molecules aggregated gradually in solution and finally arrived at balance each other. The higher the concentrations of lysozyme and NaCl were, the faster the growth rate of (I 10) face was. The growth rates of lysozyme crystal were obtained by a Zeiss microscope, and the effective surface energy (a) of growing steps were calculated about 4.01 X 10(-8) J.cm(-2) according to the model of multiple two-dimensional nucleation mechanism.

Flows with Inertia in a Three-Dimensional Random Fiber Network

A fiber web is modeled as a three-dimensional random cylindrical fiber network. Nonlinear behavior of fluid flowing through the fiber network is numerically simulated by using the lattice Boltzmann (LB) method. A nonlinear relationship between the friction factor and the modified Reynolds number is clearly observed and analyzed by using the Fochheimer equation, which includes the quadratic term of velocity. We obtain a transition from linear to nonlinear region when the Reynolds numbers are sufficiently high, reflecting the inertial effect of the flows. The simulated permeability of such fiber network has relatively good agreement with the experimental results and finite element simulations.

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.

An Improved Energy Method For Determining Young'S Modulus By Instrumented Indentation Using A Berkovich Tip

We previously proposed a method for estimating Young's modulus from instrumented nanoindentation data based on a model assuming that the indenter had a spherical-capped Berkovich geometry to take account of the bluntness effect. The method is now further improved by releasing the constraint on the tip shape, allowing it to have a much broader arbitrariness to range from a conical-tipped shape to a flat-ended shape, whereas the spherical-capped shape is just a special case in between. This method requires two parameters to specify a tip geometry, namely, a volume bluntness ratio V-r and a height bluntness ratio h(r). A set of functional relationships correlating nominal hardness/reduced elastic modulus ratio (H-n/E-r) and elastic work/total work ratio (W-e/W) were established based on dimensional analysis and finite element simulations, with each relationship specified by a set of V-r and h(r). Young's modulus of an indented material can be estimated from these relationships. The method was shown to be valid when applied to S45C carbon steel and 6061 aluminum alloy.