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.

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.

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.

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.

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.

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.

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.

Lattice Boltzmann method for simulating particle-fluid interactions

The lattice Boltzmann method （LBM） has gained increasing popularity in the last two decades as an alternative numerical approach for solving fluid flow problems. One of the most active research areas in the LBM is its application in particle-fluid systems, where the advantage of the LBM in efficiency and parallel scalability has made it superior to many other direct numerical simulation （DNS） techniques. This article intends to provide a brief review of the application of the LBM in particle-fluid systems. The numerical techniques in the LBM pertaining to simulations of particles are discussed, with emphasis on the advanced treatment for boundary conditions on the particle-fluid interface. Other numerical issues, such as the effect of the internal fluid, are also briefly described. Additionally, recent efforts in using the LBM to obtain closures for particle-fluid drag force are also reviewed.

Discrete element method for high-temperature spread in compacted powder systems

The discrete element method is applied to investigate high-temperature spread in compacted metallic particle systems formed by high-velocity compaction. Assuming that heat transfer only occurs at contact zone between particles, a discrete equation based on continuum mechanics is proposed to investigate the heat flux. Heat generated internally by friction between moving particles is determined by kinetic equations. For the proposed model, numerical results are obtained by a particle-flow-code-based program. Temperature profiles are determined at different locations and times. At a fixed location, the increase in temperature shows a logarithmic relationship with time. Investigation of three different systems indicates that the geometric distribution of the particulate material is one of the main influencing factors for the heat conduction process. Higher temperature is generated for denser packing, and vice versa. For smaller uniform particles, heat transfers more rapidly.

Combination of direct-forcing fictitious domain method and sharp interface method for dielectrophoresis of particles

In this paper, we combine the direct-forcing fictitious domain （DF/FD） method and the sharp interface method to resolve the problem of particle dielectrophoresis in two dimensions. The flow field and the motion of particles are solved with the DF/FD method, the electric field is solved with the sharp inter- face method, and the electrostatic force on the particles is computed using the Maxwell stress tensor method. The proposed method is validated via three problems： effective conductivity of particle compos- ite between two planar plates, cell trapping in a channel, and motion of particles due to both conventional and traveling wave dielectrophoretic forces.

Time-dependent surrounding soil pressure and mechanical response of tunnel lining induced by surrounding soil viscosity

A numerical simulation method of shield tunnel excavation is developed to capture the time-dependent deformation behaviour of surrounding soil. The simulation method consists of four parts:(i) an elastic-plastic-viscous constitutive model that can not only reasonably describe the viscous deformation behaviour of soil, but also appropriately calculate the plastic deformation under typical stress paths of excavation;(ii) simulation of main factors related to shield tunnel excavation, including the shield machine, face pressure, lining, grout behavior, and contacts between multiple media;(iii) a simulation procedure for excavation to reflect the process of shield tunnel excavation and achieve reasonable stress and strain fields at the end of the construction stage;(iv) a creep process that is used to investigate the long-term mechanical behaviours of the surrounding soil and tunnel lining. Taking the CK570H tunnel project in Taipei as the background, a numerical simulation is conducted by adopting the developed simulation method. Based on the simulation results, the radial and circumferential stresses acting on the lining, which are induced by the surrounding soil viscosity, are analysed. The rule of the mechanical response of lining, including its deformation, bending moment, and axial force, with time is revealed. On this basis, the long-term safety of the lining is evaluated.