Study of 3-D Numerical Simulation for Gas Transfer in the Goaf of the Coal Mining

In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.

Numerical simulation of hydraulic fracturing and associated microseismicity using finite-discrete element method

Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

3D Numerical Simulation Analysis of Passive Drag near Free Surface in Swimming

The aim of this work is to build a 3D numerical model to study the characteristics of passive drag on competitive swimmers taking into account the impact of the free surface. This model solves the 3D incompressible Navier-Stokes equations using RNG k-ε turbulence closure. The volume of fluid（VOF） method is used to locate the free surface. The 3D virtual model is created by Computer Aided Industrial Design（CAID） software, Rhinoceros. Firstly, a specific posture of swimming is studied. The simulation results are in good agreement with the data from mannequin towing experiments. The effects of a swimmer＇s arms and legs positions on swimming performance are then studied. Finally, it is demonstrated that the present method is capable of simulating gliding near the free surface.

Numerical simulation of influences of the earth medium's lateral heterogeneity on co- and post-seismic deformation

Many studies revealed that the Earth medium＇s lateral heterogeneity can cause considerable effects on the co- and post-seismic deformation field. In this study, the threedimensional finite element numerical method are adopted to quantify the effects of lateral heterogeneity caused by material parameters and fault dip angle on the co- and postseismic deformation in the near- and far-field. Our results show that： 1） the medium＇s lateral heterogeneity does affect the co-seismic deformation, with the effects increasing with the medium＇s lateral heterogeneity caused by material parameters; 2） the Lame parameters play a more dominant role than density in the effects caused by lateral heterogeneity; 3） when a fault＇s dip angle is smaller than 90, the effects of the medium＇s lateral heterogeneity on the hanging wall are greater than on the footwall; 4） the impact of lateral heterogeneity caused by the viscosity coefficient on the post-seismic deformation can affect a large area, including the near- and far-field.

水下航行体湍流数值模拟方法研究进展综述

In this paper,we present an overview of numerical simulation methods for the flow around typical underwater vehicles at high Reynolds numbers,which highlights the dominant flow structures in different regions of interest.This overview covers the forebody,midbody,stern,wake region,and appendages and summarizes flow phenomena,including laminar-to-turbulent transition,turbulent boundary layers,flow under the influence of curvatures,wake interactions,and all associated complex vortex structures.Furthermore,the current issues and challenges of capturing these flow structures are addressed.This overview provides a deep insight into the use of numerical simulation methods,including the Reynolds-averaged Navier–Stokes(RANS)method,large eddy simulation(LES)method,and the hybrid RANS/LES method,and evaluates their applicability in capturing detailed flow features.

DEM simulation of particle flow on a single deck banana screen

A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.

NUMERICAL VALIDATION OF COMPUTATIONAL MODEL FOR SHEET CAVITATING FLOWS

A computational modeling for the sheet cavitating flows is presented. The cavitation model is implemented in a viscous Navier-Stokes solver. The cavity interface and shape are determined using an iterative procedure matching the cavity surface to a constant pressure boundary. The pressure distribution, as well as its gradient on the wall, is taken into account in updating the cavity shape iteratively. Numerical computations are performed for the sheet cavitating flows at a range of cavitation numbers across the hemispheric headform/cylinder body with different grid numbers. The influence of the relaxation factor in the cavity shape updating scheme for the algorithm accuracy and reliability is conducted through comparison with other two cavity shape updating numerical schemes. The results obtained are reasonable and the iterative procedure of cavity shape updating is quite stable, which demonstrate the superiority of the proposed cavitation model and algorithms.

Monte-Carlo simulation of a stochastic differential equation

For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.

A Modified Phase-Field Model for Polymer Crystal Growth

The irrationality of existing phase field model is analyzed and a modified phase-field model is proposed for polymer crystal growth, in which the parameters are obtained from real materials and very simple to use, and most importantly, no paradoxical parameters appeared in the model. Moreover, it can simulate different microstructure patterns owing to the use of a new different free energy function for the simulation of morphologies of polymer. The new free energy function considers both the cases of T〈Tm and T≥Tm, which is more reasonable than that in published literatures that all ignored the T≥Tm case. In order to show the validity of the modified model, the finite difference method is used to solve the model and different crystallization morphologies during the solidification process of isotactic polystyrene are obtained under different conditions. Numerical results show that the growth rate of the initial secondary arms is obviously increased as the anisotropy strength increases. But the anisotropy strength seems to have no apparent effect on the global growth rate. The whole growth process of the dendrite depends mainly upon the latent heat and the latent heat has a direct effect on the tip radius and tip velocity of side branches.

Optimizing the Qusai-static Folding and Deploying of Thin-Walled Tube Flexure Hinges with Double Slots

The thin-walled tube flexure（TWTF） hinges have important potential application value in the deployment mechanisms of satellite and solar array, but the optimal design of the TWTF hinges haven＇t been completely solved, which restricts their applications. An optimal design method for the qusai-static folding and deploying of TWTF hinges with double slots is presented based on the response surface theory. Firstly, the full factorial method is employed to design of the experiments. Then, the finite element models of the TWTF hinges with double slots are constructed to simulate the qusai-static folding and deploying non-linear analysis. What＇s more, the mathematical model of the TWTF flexure hinge quasi-static folding and deploying properties are derived by the response surface method. Considering of small mass and high stability, the peak moment of quasi-static folding and deploying as well as the lightless are set as the objectives to get the optimal performances. The relative errors of the objectives between the optimal design results and the FE analysis results are less than 7%, which demonstrates the precision of the surrogate models. Lastly, the parameter study shows that both the slots length and the slots width both have significant effects to the peak moment of quasi-static folding and deploying of TWTF hinges with double slots. However, the maximum Mises stress of quasi-static folding is more sensitive to the slots length than the slots width. The proposed research can be applied to optimize other thin-walled flexure hinges under quasi-static folding and deploying, which is of great importance to design of flexure hinges with high stability and low stress.

Direct numerical simulation of particle-fluid systems by combining time-driven hard-sphere model and lattice Boltzmann method

A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method （LBM）, Particle-fluid coupling is realized by an immersed boundary method （IBM）, which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle-fluid systems.

Flow Dynamics of a Spiral-groove Dry-gas Seal

The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.

Interaction analyses between tunnel and landslide in mountain area

This paper focuses on the analytical derivation and the numerical simulation analyses to predict the interaction influences between a landslide and a new tunnel in mountain areas. Based on the slip-line theory, the disturbance range induced by tunneling and the minimum safe distance between the tunnel vault and the sliding belt are obtained in consideration of the mechanical analyses of relaxed rocks over the tunnel opening. The influence factors for the minimum safe crossing distance are conducted,including the tunnel radius, the friction angle of surrounding rocks, the inclination angle of sliding belt,and the friction coefficient of surrounding rocks. Secondly, taking account of the compressive zone and relaxed rocks caused by tunneling, the Sarma method is employed to calculate the safety factor of landslide. Finally, the analytical solutions for interaction between the tunnel and the landslide are compared with a series of numerical simulations, considering the cases for different perpendicular distances between the tunnel vault and the sliding belt. Resultsshow that the distance between the tunnel vault and the slip zone has significant influence on the rock stress and strain. For the case of the minimum crossing distance, a plastic zone in the landslide traversed by tunneling would be formed with rather large range, which seriously threatens the stability of landslide. This work demonstrates that the minimum safe crossing distance obtained from numerical simulation is in a good agreement with that calculated by the proposed analytical solutions.

The mass and heat transfer process through the door seal of refrigeration

As one of the main reasons causing leakage heat load in a refrigerator,mass and heat transfer through refrigerator door seal is of great importance to be studied.In this paper,a model is presented for numerical simulation of mass and heat transfer process through refrigerator door seal,and an experiment apparatus is designed and set up as well for comparison.A two-dimensional model and tracer gas method are used in simulation and experiment,respectively.It can be found that the relative deviations of air infiltration rate between the simulated results and experimental results were less than 1%,and the temperature difference errors at two special points of the door seal were less than 2.03℃.In conclusion,the simulated results are in good agreement with the experimental results.This paper initially sets up a model that can accurately simulate the heat and mass transfer through the refrigerator door seal,and the model can be used in refrigerator door seal optimization research in the follow-up study.

Radiation heat transfer model for complex superalloy turbine blade in directional solidification process based on finite element method

For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.

A fluctuating lattice-Boltzmann model for direct numerical simulation of particle Brownian motion

A single-relaxation-time fluctuating lattice-Boltzmann （LB） model for direct numerical simulation （DNS） of particle Brownian motion is established by adding a fluctuating component to the lattice-Boltzmann equations （LBEs）. The fluctuating term is proved to be the random stress tensor in fluctuating hydrodynamics by recovering Navier-Stokes equations from LBEs through a Chapman-Enskog expansion. A three-dimensional implementation of the model is also presented, along with simulations of a single spherical particle and 125 spherical particles at short times. Numerical results including the meansquare displacement, velocity autocorrelation function and self-diffusion coefficient of particles compare favorably with theoretical results and previous numerical results.

A Simplified Numerical Simulation Method Considering Active Devices for Opto-Electronic Mixed Modules

A simplified simulation method based on the FDTD technique that can handle active devices is proposed. This method well suits the electrical crosstalk analysis of multi-channel integrated, opto-electronic mixed modules. We apply this method to an 8-channel integrated super-compact high-sensitivity optical module. The results show good agreement between simulations and measurements.

Rheological Simulating Study of Mechanical Conditions for Syn-extensional Basin and Mountain Coupling System

This article gives a mechanical model, in which the layers of lithosphere are assumed to be the creep materials, to study the coupling mechanism of a syn-basin-mountain system quantitatively by using the numerical simulating method. A geological dynamic extensional mode given by some geologists is theoretically discussed and verified. The study shows that lithosphere thickening or thinning is closely related to the thermal activity, or in other words, thermal convection beneath the lithosphere. It is one of the important factors affecting the formation of the basin-mountain coupling system. As an essential condition, only the upward buoyant force and the horizontal dragging force caused by the thermal convection jointly act on the bottom of the lithosphere, the stress and strain states in rock's layers are advantageous to forming the tectonic-landforms of the basin-mountain coupling system. A study on the creep features of the lithosphere shows that the stress and strain in the rock's layers vary with time when the lasting forces act on the boundary. They increase rapidly at initial stage and decrease steadily after reaching the peak value. Phenomena of stress relaxation are significant for studying the tectonic evolution.

Hydrodynamics characterization of a choanoid fluidized bed bioreactor used in the bioartificial liver system： Fully resolved simulation with a fctitious domain method

Choanoid fluidized bed bioreactors （CFBBs） are newly developed core devices used in bioartificial liver- support systems to detoxify blood plasma of patients with microencapsulated liver cells. Direct numerical simulations （DNS） with a direct-forcing/fictitious domain （DF/FD） method were conducted to study the hydrodynamic performance of a CFBB. The effects of particle-fluid density ratio, particle number, and fil- ter screens preventing particles flowing out of the reactor were investigated. Depending on density ratio, two flow patterns are evident： the circulation mode in which the suspension rises along one sidewall and descends along the other sidewall, and the non-circulation mode in which the whole suspension roughly flows upward. The circulation mode takes place under non-neutral-buoyancy where the particle sedimentation dominates, whereas the non-circulation mode occurs under pure or near-neutral buoy- ancy with particle-fluid density ratios of unity or near unity. With particle-fluid density ratio of 1.01, the bioartificial liver reactor performs optimally as the significant particle accumulation existing in the non-circulation mode and the large shear forces on particles in the circulation mode are avoided. At higher particle volume fractions, more particles accumulate at the filter screens and a secondary counter circulation to the primary flow is observed at the top of the bed. Modelled as porous media, the filter screens play a negative role on particle fluidization velocities; without screens, particles are fluidized faster because of the higher fluid velocities in the jet center region. This work extends the DF/FD-based DNS to a fluidized bed and accounts for effects from inclined side walls and porous media, providing some hydrodynamics insight that is important for CFBB design and operation optimization.