Molecular dynamics simulation of the flow mechanism of shear-thinning fluids in a microchannel

Shear-thinning fluids have been widely used in microfluidic systems,but their internal flow mechanism is still unclear.Therefore,in this paper,molecular dynamics simulations are used to study the laminar flow of shear-thinning fluid in a microchannel.We validated the feasibility of our simulation method by evaluating the mean square displacement and Reynolds number of the solution layers.The results show that the change rule of the fluid system's velocity profile and interaction energy can reflect the shear-thinning characteristics of the fluids.The velocity profile resembles a top-hat shape,intensifying as the fluid's power law index decreases.The interaction energy between the wall and the fluid decreases gradually with increasing velocity,and a high concentration of non-Newtonian fluid reaches a plateau sooner.Moreover,the velocity profile of the fluid is related to the molecule number density distribution and their values are inversely proportional.By analyzing the radial distribution function,we found that the hydrogen bonds between solute and water molecules weaken with the increase in velocity.This observation offers an explanation for the shear-thinning phenomenon of the non-Newtonian flow from a micro perspective.

Explosion resistance performance of reinforced concrete box girder coated with polyurea:Model test and numerical simulation

To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.

Desulfurization characteristics of slaked lime and regulation optimization of circulating fluidized bed flue gas desulfurization process--A combined experimental and numerical simulation study

Circulating fluidized bed flue gas desulfurization(CFB-FGD) process has been widely applied in recent years. However, high cost caused by the use of high-quality slaked lime and difficult operation due to the complex flow field are two issues which have received great attention. Accordingly, a laboratory-scale fluidized bed reactor was constructed to investigate the effects of physical properties and external conditions on desulfurization performance of slaked lime, and the conclusions were tried out in an industrial-scale CFB-FGD tower. After that, a numerical model of the tower was established based on computational particle fluid dynamics(CPFD) and two-film theory. After comparison and validation with actual operation data, the effects of operating parameters on gas-solid distribution and desulfurization characteristics were investigated. The results of experiments and industrial trials showed that the use of slaked lime with a calcium hydroxide content of approximately 80% and particle size greater than 40 μm could significantly reduce the cost of desulfurizer. Simulation results showed that the flow field in the desulfurization tower was skewed under the influence of circulating ash. We obtained optimal operating conditions of 7.5 kg·s^(-1)for the atomized water flow, 70 kg·s^(-1)for circulating ash flow, and 0.56 kg·s^(-1)for slaked lime flow, with desulfurization efficiency reaching 98.19% and the exit flue gas meeting the ultraclean emission and safety requirements. All parameters selected in the simulation were based on engineering examples and had certain application reference significance.

Assessing cutter-rock interaction during TBM tunnelling in granite:Large-scale standing rotary cutting tests and 3D DEM simulations

The widespread utilisation of tunnel boring machines(TBMs)in underground construction engineering requires a detailed investigation of the cutter-rock interaction.In this paper,we conduct a series of largescale standing rotary cutting tests on granite in conjunction with high-fidelity numerical simulations based on a particle-type discrete element method(DEM)to explore the effects of key cutting parameters on the TBM cutter performance and the distribution of cutter-rock contact stresses.The assessment results of cutter performance obtained from the cutting tests and numerical simulations reveal similar dependencies on the key cutting parameters.More specifically,the normal and rolling forces exhibit a positive correlation with penetration but are slightly influenced by the cutting radius.In contrast,the side force decreases as the cutting radius increases.Additionally,the side force shows a positive relationship with the penetration for smaller cutting radii but tends to become negative as the cutting radius increases.The cutter's relative effectiveness in rock breaking is significantly impacted by the penetration but shows little dependency on the cutting radius.Consequently,an optimal penetration is identified,leading to a low boreability index and specific energy.A combined Hertz-Weibull function is developed to fit the cutter-rock contact stress distribution obtained in DEM simulations,whereby an improved CSM(Colorado School of Mines)model is proposed by replacing the original monotonic cutting force distribution with this combined Hertz-Weibull model.The proposed model outperforms the original CSM model as demonstrated by a comparison of the estimated cutting forces with those from the tests/simulations.The findings from this work that advance our understanding of TBM cutter performance have important implications for improving the efficiency and reliability of TBM tunnelling in granite.

Wavefield simulation in porous media saturated with two immiscible fluids

Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation（capillary pressure）.So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.

ISOLATION DESIGN AND SIMULATIONAL EXPERIMENT OF NANOMETER MICRO MOVEMENT TEST-BED

The relation of mass, stiffness and rate of damping is obtained by using the mechanical analysis of the obstructive vibration system of two dimensions for the design of the obstructive vibration system of more freedom and the micro vibration test bed. The result of stimulational experiment indicates that the isolation of vibration of this system is satisfactory. The design method of vibration can be used as the reference to ultra precision machine tool, super micro orientation machanism and so on.

Model test vs virtual simulation of a VLCC FPSO hookup

This paper describes the model test and the virtual simulation respectively for the VLCC class FPSO hookup, as well as addresses their different applications to the mating operation between the FPSO and the soft yoke mooring system （SYMS） in extremely shallow water. The scope of the model test and the virtual simulation covers various installation stages including a series of positioning trials, positioning keeping and temporary mooring to the pre-installed SYMS mooring tower, pendulum mating, and yoke ballasting to storm-safe. The model test is to accurately verify bollard pull capacity to keep the FPSO in position and assess motion responses and mooring loads for the FPSO and installation vessels during various installation stages. The virtual simulation is to provide a virtual-reality environment, thus realistically replicating the hookup operation at the Simulation Test Center （STC） facility and identifying any deficiencies in key installation personnel, execution plan, or operation procedures. The methodologies of the model test and the virtual simulation addressed here can be easily extended to the deepwater applications such as positioning and installation operations of various floating systems.

Digital Simulation of Full Scale Static Test of Aircraft

Full scale aircraft static test is a very important process of aircraft design, it is costly and time consuming. The testing accuracy and validity mainly depend on the rationality of the test scheme design. When the aircraft is being tested, the specimen＇s safety mainly depends on monitoring and understanding the testing data by way of evaluating the coherence with the digital simulation data synchyononsly. The test digital simulation can aid realizing above requirements and improving the test efficiency significantly during test scheme design stage or testing stage respectively. The key technologies and the solving methods of test digital simulation are presented and the application example is given.

On modeling approach for embedded real-time software simulation testing

Modeling technology has been introduced into software testing field. However, how to carry through the testing modeling effectively is still a difficulty. Based on combination of simulation modeling technology and embedded real-time software testing method, the process of simulation testing modeling is studied first. And then, the supporting environment of simulation testing modeling is put forward. Furthermore, an approach of embedded real-time software simulation testing modeling including modeling of cross-linked equipments of system under testing （SUT）, test case, testing scheduling, and testing system service is brought forward. Finally, the formalized description and execution system of testing models are given, with which we can realize real-time, closed loop, mad automated system testing for embedded real-time software.

Centrifugal model tests and numerical simulations for barrier dam break due to overtopping

The present study focuses on the breaching process and failure of barrier dams due to overtopping. In this work, a series of centrifugal model tests is presented to examine the failure mechanisms of landslide dams. Based on the experimental results, failure process and mechanism of barrier dam due to overtopping are analyzed and further verified by simulating the experimental overtopping failure process. The results indicate that the barrier dam will develop during the entire process of overtopping in the width direction, whereas the breach will cease to develop at an early stage in the depth direction because of the large particles that accumulate on the downstream slope. Moreover, headcut erosion can be clearly observed in the first two stages of overtopping, and coarsening on the downstream slope occurs in the last stage of overtopping. Thus, the bottom part of the barrier dam can survive after dam breaching and full dam failure becomes relatively rare for a barrier dam. Furthermore, the remaining breach would be smaller than that of a homogeneous cohesive dam under the same conditions.

A laboratory acoustic emission experiment and numerical simulation of rock fracture driven by a high-pressure fluid source

In order to improve our understanding of rock fracture and fault instability driven by high-pressure fluid sources, the authors carried out rock fracture tests using granite under a confining pressure of 80 MPa with fluid injection in the laboratory. Furthermore, we tested a number of numerical models using the FLAC;modeling software to find the best model to represent the experimental results. The high-speed multichannel acoustic emission（AE） waveform recording system used in this study made it possible to examine the total fracture process through detailed monitoring of AE hypocenters and seismic velocity.The experimental results show that injecting high-pressure oil into the rock sample can induce AE activity at very low stress levels and can dramatically reduce the strength of the rock. The results of the numerical simulations show that major experimental results, including the strength, the temporal and spatial patterns of the AE events, and the role of the fluid can be represented fairly well by a model involving（1） randomly distributed defect elements to model pre-existing cracks,（2） random modification of rock properties to represent inhomogeneity introduced by different mineral grains, and（3）macroscopic inhomogeneity. Our study, which incorporates laboratory experiments and numerical simulations, indicates that such an approach is helpful in finding a better model not only for simulating experimental results but also for upscaling purposes.

Numerical Simulation of Fluid Migration during Ore Formation of Carboniferous Exhalation-Sedimentary Massive Sulfide Deposits in the Tongling District, Anhui Province

Numerical simulation of fluid migration during the ore-forming process of the Carboniferous exhalation- sedimentary (Sedex) massive sulfide deposits in the Tongling district shows that fluid and thermal activities in lying-wall rocks were limited to a small area around the main draining passage, which led to weak mineralization and alteration in the lying-wall rock. Temperature and fluid fields indicate that mineralization and alteration in the lying-wall rock of the Sedex-type deposits are usually weaker than those of volcano-hosted massive sulfide deposits. Fluid migration involves the following processes: seawater penetrating and leaching the lying-wall rocks, then mixing with ascending hydrothermal fluids in the main draining passage, and finally jetting into seafloor. Although fluid activity-influenced area is rather small, the content of metals leached out from the lying-wall rocks is high enough to form large-scale ore deposits. Tensional contemporaneous faults accompanied with strong heat flows controlled the formation and distribution of Sedex deposits. The tensional tectonic regime on the northern margin of the Yangtze block during the Hercynian provided Sedex deposits with a prerequisite geodynamic condition.

Grassland Evolution Under Soil Degradation: Numerical Simulation and Test

Both theoretical and field observations were examined to study the close relationship between soil degeneration and the evolution of grassland vegetation. A general n-species model of equal competition under different degrees of soil degradation was applied to field data in order to probe the dynamic processes and mechanisms of vegetation evolution due to the effects of the soil's ecological deterioration on grassland vegetation. Comparisons were made between the theoretical results and the practical surveys with satisfactory results.

Simulation of large-scale numerical substructure in real-time dynamic hybrid testing

A solution scheme is proposed in this paper for an existing RTDHT system to simulate large-scale finite element （FE） numerical substructures. The analysis of the FE numerical substructure is split into response analysis and signal generation tasks, and executed in two different target computers in real-time. One target computer implements the response analysis task, wherein a large time-step is used to solve the FE substructure, and another target computer implements the signal generation task, wherein an interpolation program is used to generate control signals in a small time-step to meet the input demand of the controller. By using this strategy, the scale of the FE numerical substructure simulation may be increased significantly. The proposed scheme is initially verified by two FE numerical substructure models with 98 and 1240 degrees of freedom （DOFs）. Thereafter, RTDHTs of a single frame-foundation structure are implemented where the foundation, considered as the numerical substructure, is simulated by the FE model with 1240 DOFs. Good agreements between the results of the RTDHT and those from the FE analysis in ABAQUS are obtained.

Progress in molecular-simulation-based research on the effects of interface-induced fluid microstructures on flow resistance

In modern chemical engineering processes, solid interface involvement is the most important component of process intensification techniques, such as nanoporous membrane separation and heterogeneous catalysis. The fundamental mechanism underlying interfacial transport remains incompletely understood given the complexity of heterogeneous interfacial molecular interactions and the high nonideality of the fluid involved. Thus, understanding the effects of interface-induced fluid microstructures on flow resistance is the first step in further understanding interfacial transport. Molecular simulation has become an indispensable method for the investigation of fluid microstructure and flow resistance. Here, we reviewed the recent research progress of our group and the latest relevant works to elucidate the contribution of interface-induced fluid microstructures to flow resistance.We specifically focused on water, ionic aqueous solutions, and alcohol–water mixtures given the ubiquity of these fluid systems in modern chemical engineering processes. We discussed the effects of the interfaceinduced hydrogen bond networks of water molecules, the ionic hydration of ionic aqueous solutions, and the spatial distributions of alcohol and alcohol–water mixtures on flow resistance on the basis of the distinctive characteristics of different fluid systems.

COMPUTATIONAL FLUID DYNAMICS(CFD) SIMULATIONS OF DRAG REDUCTION WITH PERIODIC MICRO-STRUCTURED WALL

Computational fluid dynamics（CFD） simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds number. The purpose of the current study is to numerically find out the effects of periodic micro-structured wall on the flow resistance in rectangular microchannel with the different spacings between microridges ranging from 15 to 60 pm. The simulative results indicate that pressure drop with different spacing between microridges increases linearly with flow velocity and decreases monotonically with slip velocity; Pressure drop reduction also increases with the spacing between microridges at the same condition of slip velocity and flow velocity. The results of numerical simulation are compared with theoretical predictions and experimental results in the literatures. It is found that there is qualitative agreement between them.

Simulation and Analysis on the Two-Phase Flow Fields in a Rotating-Stream-Tray Absorber by Using Computational Fluid Dynamics

The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, and a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall; the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.

Simulation and Field Testing of Multiple Vehicles Collision Avoidance Algorithms

A global planning algorithm for intelligent vehicles is designed based on the A* algorithm, which provides intelligent vehicles with a global path towards their destinations. A distributed real-time multiple vehicle collision avoidance(MVCA)algorithm is proposed by extending the reciprocal n-body collision avoidance method. MVCA enables the intelligent vehicles to choose their destinations and control inputs independently,without needing to negotiate with each other or with the coordinator. Compared to the centralized trajectory-planning algorithm, MVCA reduces computation costs and greatly improves the robustness of the system. Because the destination of each intelligent vehicle can be regarded as private, which can be protected by MVCA, at the same time MVCA can provide a real-time trajectory planning for intelligent vehicles. Therefore,MVCA can better improve the safety of intelligent vehicles. The simulation was conducted in MATLAB, including crossroads scene simulation and circular exchange position simulation. The results show that MVCA behaves safely and reliably. The effects of latency and packet loss on MVCA are also statistically investigated through theoretically formulating broadcasting process based on one-dimensional Markov chain. The results uncover that the tolerant delay should not exceed the half of deciding cycle of trajectory planning, and shortening the sending interval could alleviate the negative effects caused by the packet loss to an extent. The cases of short delay(< 100100 ms) and low packet loss(< 5%) can bring little influence to those trajectory planning algorithms that only depend on V2 V to sense the context, but the unpredictable collision may occur if the delay and packet loss are further worsened. The MVCA was also tested by a real intelligent vehicle, the test results prove the operability of MVCA.

Fast and Stable Surface Feature Simulation for Particle-Based Fluids

In order to efficiently and realistically capture microscopic features of fluid surface,a fast and stable surface feature simulation approach for particle-based fluids is presented in this paper.This method employs a steady tension and adhesion model to construct surface features with the consideration of the adsorption effect of fluid to solid.Molecular cohesion and surface area minimization are appended for surface tension,and adhesion is added to better show the microscopic characteristics of fluid surface.Besides,the model is integrated to an implicit incompressible smoothed particle hydrodynamics(SPH)method to improve the efficiency and stability of simulation.The experimental results demonstrate that the method can better simulates surface features in a variety of scenarios stably and efficiently.

Computational fluid dynamics simulations of respiratory airflow in human nasal cavity and its characteristic dimension study

To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.