Molecular mechanics and dynamics simulation of hydrogen diffusion in aluminum melt

The main impurities in aluminum melt are hydrogen and Al_2O_3,which can deteriorate melt quality and materials performance.However,the diffusion process of H atoms in aluminum melt and the interactions among Al atoms,Al_2O_3 and hydrogen have been studied rarely.Molecular mechanics and dynamics simulations are employed to study the diffusion behaviors of different types of hydrogen,such as free H atoms,H atoms in H_2 and H^+ions in H_2O using COMPASS force field.Correspondingly,force field types h,h1h and h1o are used to describe different types of hydrogen which are labeled as H_h,H_(h1h) and H_(h1o).The results show that the adsorption areas are maximum for H_(h1o),followed by H_(h1h) and H_h.The diffusion ability of H_(h1o) is the strongest whereas H_h is hard to diffuse in aluminum melt because of the differences in radius and potential well depth of various types of hydrogen.Al_2O_3 cluster makes the Al atoms array disordered,creating the energy conditions for hydrogen diffusion in aluminum melt.Al_2O_3 improves the diffusion of H_h and H_(h1o),and constrains H_(h1h) which accumulates around it and forms gas porosities in aluminum.H_(h1o) is the most dispersive in aluminum melt,moreover,the distance of Al-H_(h1o) is shorter than that of Al-H_(h1h),both of which are detrimental to the removal of H_(h1o).The simulation results indicate that the gas porosities can be eliminated by the removal of Al_2O_3 inclusions,and the dispersive hydrogen can be removed by adsorption function of gas bubbles or molten fluxes.

Numerical Simulation and Dynamical Analysis for Low Salinity Water Lens in the Expansion Area of the Changjiang Diluted Water

The low salinity water lenses（LSWLes） in the expansion area of the Changjiang diluted water（CDW） exist in a certain period of time in some years. The impact of realistic river runoff, ocean currents and weather conditions need to be taken into account in the dynamical analysis of LSWL, which is in need of research. In this paper, the POM-σ-z model is used to set up the numerical model for the expansion of the CDW. Then LSWL in summer 1977 is simulated, and its dynamic mechanism driven by wind, tide, river runoff and the Taiwan Warm Current is also analyzed. The simulated results indicate that the isolated LSWL detaches itself from the CDW near the river mouth, and then moves towards the northeast region outside the Changjiang Estuary. Its maintaining period is from July 26 to August 11. Its formation and development is mainly driven by two factors. One is the strong southeasterly wind lasting for ten days. The other is the vertical tidal mixing during the transition from neap tide to spring tide.

Factors influencing physical property evolution in sandstone mechanical compaction:the evidence from diagenetic simulation experiments

In order to analyze the factors influencing sandstone mechanical compaction and its physical property evolution during compaction processes, simulation exper- iments on sandstone mechanical compaction were carried out with a self-designed diagenetic simulation system. The experimental materials were modem sediments from dif- ferent sources, and the experiments were conducted under high temperature and high pressure. Results of the exper- iments show a binary function relation between primary porosity and mean size as well as sorting. With increasing overburden pressure during mechanical compaction, the evolution of porosity and permeability can be divided into rapid compaction at an early stage and slow compaction at a late stage, and the dividing pressure value of the two stages is about 12 MPa and the corresponding depth is about 600 m. In the slow compaction stage, there is a good exponential relationship between porosity and overburden pressure, while a good power function relationship exists between permeability and overburden pressure. There is also a good exponential relationship between porosity and permeability. The influence of particle size on sandstone mechanical compaction is mainly reflected in the slowcompaction stage, and the influence of sorting is mainly reflected in the rapid compaction stage. Abnormally high pressure effectively inhibits sandstone mechanical com- paction, and its control on sandstone mechanical com- paction is stronger than that of particle size and sorting. The influence of burial time on sandstone mechanical compaction is mainly in the slow compaction stage, and the porosity reduction caused by compaction is mainly con- trolled by average particle size.

Numerical simulation of dust distribution at a fully mechanized face under the isolation effect of an air curtain

At a fully mechanized working face of a coal mine as prototype,we investigated,by simulation,the flow field and dust distribution during the process of its isolation by a curtain of air,using the CFD software, Fluent.The results show that the air curtain installed on the shearer can effectively prevent the dust （especially the respirable dust）from diffusing into the work area of the operator,reducing the dust concentration on the side of the operator and greatly improving his working environment.The field application of the air curtain shows that the dust-isolation effect of an air curtain is quite noticeable.The isolation efficiency for respiratory dust is over 70%and,as well,it has good dust-isolation effect for nonrespiratory dust.The air curtain is a useful way to resolve the problem of dust-isolation at a fully mechanized working face.It has a practical background elsewhere with more extensive applications.

Numerical simulation of spatial-temporal evolution characteristics of subsurface fluid based on strong body seismogenic model$$$$

Using finite element technique of the plane-strain problem in solid-liquid two-phase medium, we Studied the char acteristics of 'field precursors' and 'focus precursors' of subsurface fluid and their spatial-temporal evolution in case of dip-slip earthquake. The results show that: ① the change of ground fluid is slow and the anomaly is not prominent in the early period which is of elastic accumulation and non-linear; ② dilatancy emerges and anomalyfocus mainly in the source region in the moderate period which is hardening and of local dilatancy. In the period the focus precursors emerge earlier than the field precursors; ③ anomalies spreed continuously in the source area and new regions with big anomaly emerge out of the source region in the middle-short period which is of large scale dilatancy.

Allosteric Mechanism of Calmodulin Revealed by Targeted Molecular Dynamics Simulation

Calmodulin （CAM） is involved in the regulation of a variety of cellular signaling pathways. To accomplish its physiological functions, CaM binds with Ca2＋ at its EF-hand Ca2＋ binding sites which induce the conformational switching of CaM. However, the molecular mechanism by which Ca2＋ binds with CaM and induces conformational switching is still obscure. Here we combine molecular dynamics with targeted molecular dynamics simulation and achieve the state-transition pathway of CaM. Our data show that Ca2＋ binding speeds up the conformational transition of CaM by weakening the interactions which stabilize the closed state. It spends about 6.5 ns and 5.25 ns for transition from closed state to open state for apo and holo CaM, respectively. Regarding the contribution of two EF-hands, our data indicate that the first EF-hand triggers the conformational transition and is followed by the second one. We determine that there are two interaction networks which contribute to stabilize the closed and open states, respectively.

Identification of Forcing Mechanisms of Convective Initiation over Mountains through High-Resolution Numerical Simulations

Convection and its ensuing severe weather, such as heavy rainfall, hail, tornado, and high wind, have significant im- pacts on our society and economy （e.g., Cao et al., 2004; Fritsch and Carbone, 2004; Verbout et al., 2006; Ashley and Black, 2008; Cao, 2008; Cao and Ma, 2009; Zhang et al., 2014）. Due to its localized and transient nature, the initiation of convection or convective initiation remains one of the least understood aspects of convection in the scientific communi- ties, and it is a significant challenge to accurately predict the exact timing and location of convective initiation （e.g., Cai et al., 2006; Wilson and Roberts, 2006; Xue and Martin, 2006; Cao and Zhang, 2016）.

Investigation and Improvement of the Staggered Labyrinth Seal

Recent studies on staggered labyrinth seals have focused on the effects of different parameters,such as the pressure ratio and rotational speed on the leakage flow rate.However,few investigations pay sufficient attention to flow details and the sealing mechanism,which would be of practical importance in designing seals having higher performance.This paper establishes a theoretical model to study the seal mechanism,thus revealing that leakage is determined by the pressure ratio and geometric structure.Numerical simulation is implemented to illustrate details of the flow field within the seal structure.Viscous dissipation is used to quantitatively investigate the contribution that each location makes to the seal performance,revealing that orifices and stagnation points are the most important positions in the seal structure,generating the most dissipation.The orifice is carefully studied by using the theoretical model.Experiments for different pressure ratios are conducted and the results match well with those of the theoretical model and numerical simulation,verifying the theoretical model and analysis of the seal mechanism.Three new designs,based on a good understanding of the seal mechanism,are presented,with one reducing leakage by 24.5%.

The Erosion Effect of Kapton Film in a Ground-based Atomic Oxygen Irradiation Simulator

In order to investigate the effect of space environmental factors on spacecraft materials, a ground-based simulation facility for space atomic oxygen（AO） irradiation was developed in our laboratory. Some Kapton film samples were subjected to AO beam generated by this facility. The Kapton films before and after AO exposure were analyzed comparatively using optical microscopy, scanning electronic microscopy, atomic force microscopy, high-precision microbalance, and X-ray photoelectron spectroscopy. The experimental results indicate that the transmittance of Kapton film will be reduced by AO irradiation notably, and its color deepens from pale yellow to brown. Surface roughness of the AO-treated sample is already increased obviously after AO irradiation for 5 hours, and exhibits a flannel-like appearance after 15 hours’ exposure in AO beam. The imide rings and benzene rings in kapton molecule are partially decomposed, and some new bonds form during AO irradiation. The mass loss of kapton film increases linearly with the increase of AO fluence, which is resulted from the formation of volatile products, such as CO, CO2 and NOx. The breakage in structure and degradation in properties of AO-treated Kapton film can be attributed to the integrated effect ofimpaction and oxidization of AO beam. The test results agree well with the space flight experimental data.

DNS analysis of incipient drop impact dynamics using an accurate level set method

A series of 2D direct numerical simulations were performed with an accurate level set method for single drop impacts.The adopted ACLS method was validated to be efficient with perfect mass conservation in both normal and oblique impacts.A square-root correction for neck bases was modified in accuracy as well as scope of applications.In addition,process of jet formation and evolution was studied to reveal internal dynamics in drop impacts.It's found that pressure gradient and vortex are coexisting and completive reasons for jet topology while the inclined angle has a significant effect on them.Mechanisms of jet formation and evolution are different in the front and back necks.With the help of PDF distribution and correction calculation,a compromise in the competition is observed.This work lays a solid foundation for further studies of dynamics in gas-liquid flows.

Influence of gas transport mechanisms on the productivity of multi-stage fractured horizontal wells in shale gas reservoirs

In order to investigate the influence on shale gas well productivity caused by gas transport in nanometer- size pores, a mathematical model of multi-stage fractured horizontal wells in shale gas reservoirs is built, which considers the influence of viscous flow, Knudsen diffusion, surface diffusion, and adsorption layer thickness. A dis- crete-fracture model is used to simplify the fracture mod- cling, and a finite element method is applied to solve the model. The numerical simulation results indicate that with a decrease in the intrinsic matrix permeability, Knudsen diffusion and surface diffusion contributions to production become large and cannot be ignored. The existence of an adsorption layer on the nanopore surfaces reduces the effective pore radius and the effective porosity, resulting in low production from fractured horizontal wells. With a decrease in the pore radius, considering the adsorption layer, the production reduction rate increases. When the pore radius is less than 10 nm, because of the combined impacts of Knudsen diffusion, surface diffusion, and adsorption layers, the production of multi-stage fractured horizontal wells increases with a decrease in the pore pressure. When the pore pressure is lower than 30 MPa, the rate of production increase becomes larger with a decrease in pore pressure.

Effects of dipolar interactions on magnetic properties of Co nanowire arrays

Magnetic properties and magnetization processes of Co nanowire arrays with various packing densities are investigated by means of object-oriented micromagnetic framework（OOMMF） software package with finite difference micromagnetic simulations. The packing density of nanowires is changed with the diameter, number of nanowires and center-to-center spacing between the wires. The magnetization reversal mechanism and squareness of the hysteresis loops of the nanowire arrays are very sensitive to the packing density of nanowires. Clear steps and plateaux on the demagnetization are visible,which turns out that dipolar interactions among the wires have a significant influence on switching field.

Numerical investigation of unsteady mixing mechanism in plate film cooling

A large-scale large eddy simulation in high performance personal computer clusters is carried out to present unsteady mixing mechanism of film cooling and the development of films. Simulation cases include a single-hole plate with the inclined angle of 30° and blowing ratio of 0.5, and a single-row plate with hole-spacing of 1.5D and 2D （diameters of the hole）. According to the massive simulation results, some new unsteady phenomena of gas films are found. The vortex system is changed in different position with the development of film cooling with the time marching the process of a single-row plate film cooling. Due to the mutual interference effects including mutual exclusion, a certain periodic sloshing and mutual fusion, and the structures of a variety of vortices change between parallel gas films. Macroscopic flow structures and heat transfer behaviors are obtained based on 20 million grids and Reynolds number of 28600.

Molecular Structural and Properties of 3-chloro-4 (dichloromethyl)-5-hydroxy-2[5H]-furanone (MX)

3 - chloro - 4 (dichloromethyl) - 5 - hydroxy - 2 [5H] - furanone (MX) formed during chlorination of water containing natural organic substances, is a very potent bacterial mutgen. Molecular mechanics calculations to evaluate the conformation of structure, and to determine structure relationship properties are put forward. The investigations allow the correlation of molecular structures of MX with its properties, such as mass,partial charges, steric energy, frontier molecular orbital. The VRML molecular models have been investigated using Virtual Reality software. The spectral simulation of MX is illustrated. The principal aim is to develop an efficient method which control of MX.

Clinical Data-Driven Finite Element Analysis of the Kinetics of Chewing Cycles in Order to Optimize Occlusal Reconstructions Dedicated to Professor Karl Stark Pister for his 95th birthday

The occlusal design plays a decisive role in the fabrication of dental restorations.Dentists and dental technicians depend on mechanical simulations of mandibular movement that are as accurate as possible,in particular,to produce interference-free yet chewing-efficient dental restorations.For this,kinetic data must be available,i.e.,movements and deformations under the influence of forces and stresses.In the present study,so-called functional data were collected from healthy volunteers to provide consistent information for proper kinetics.For the latter purpose,biting and chewing forces,electrical muscle activity and jaw movements were registered synchronously,and individual magnetic resonance tomograms(MRI)were prepared.The acquired data were then added to a large complex finite element model of the complete masticatory system using the functional information obtained and individual anatomical geometries so that the kinetics of the chewing process and teeth grinding could be realistically simulated.This allows developing algorithms that optimize computer-aided manufacturing of dental prostheses close to occlusion.In this way,a failure-free function of the dental prosthesis can be guaranteed and its damage during usage can be reduced or prevented even including endosseous implants.

Sulfur poisoning mechanism of three way catalytic converter and its grey relational analysis

Combining a detailed catalytic surface reaction mechanism with noble metal and promoter elementary reactions, a new three-way catalytic converter(TWC) reaction mechanism is established. Based on the new mechanism, steady condition numerical simulation is carried out, and the change of light-off temperatures and conversion efficiency with various SO2 contents is obtained. By grey relational analysis(GRA), the relational grade between conversion efficiency and SO2 content is obtained. And, the result shows that SO2 content has the most important influence on C3H6 and NOX conversion efficiency. This provides an important reference to the improvement of activity design of TWC, and may provide guidance for the condition design and optimization of TWC.

Quantitative Method of the Structural Damage Identification of Gas Explosion Based on Case Study:The Shanxi “11. 23” Explosion Investigation

In order to present a retrospective analysis of exposition accidents using input data from investigation processes,data from a specific accident was examined,in which we analyzed possible involved gas species（ liquefied petroleum gas; nature gas） and computed their concentrations and distributions based on the interactions between the structures and the effects of the explosion. In this study,5 scenarios were created to analyze the impact effect. Moreover,a coupling algorithm was put into practice,with a practical outflow boundary and joint strength are applied. Finally,the damage effects of each scenario were simulated. Our experimental results showed significant differences in the 5 scenarios concerning the damage effects on the building structures. The results from scenario 3 agree with the accident characteristics,demonstrating the effectiveness of our proposed modeling method. Our proposed method reflects gas properties,species and the concentration and distribution,and the simulated results validates the root cause,process,and consequences of accidental explosions. Furthermore,this method describes the evolution process of explosions in different building structures. Significantly,our model demonstrates the quantatative explosion effect of factors like gas species,gas volumes,and distributions of gases on explosion results. In this study,a feasible,effective,and quantitative method for structure safety is defined,which is helpful to accelerate the development of safer site regulations.

A NUMERICAL SIMULATION OF 3-D INNER FLOW IN UP-STREAM PUMPING MECHANICAL SEAL

Numerical simulation of 3-D inner flow between Up-stream Pumping Mechanical Face Seals （UPMFS） faces was initially done by CFD software, which made the flow visualization come true. Simulation results directly discover the action of hydrodynamic lubrication, and by comparison with that of Conventional Mechanic Face Seals （CMFS）, the advantage over bigger bearing capability, less friction and much less leakage are explained clearly. Otherwise there are also some different ideas and results from precedent analysis and computational research results： dynamic and static pressure profiles can be obtained respectively instead of the analytic total pressure distribution only, pressure distribution is nonlinear, while always be solved as linear, lower pressure is observed at the area of inner diameter caused by the grooves, but its possible cavitations effects to the performance of UPMFS still need further study.

Research on mechanics properties of crop stalks:A review

As one of the foundational research areas in agricultural engineering,the study of agricultural materials includes an important branch:mechanics properties of crop stalks.This branch is now widely associated with design and improvement of agricultural equipment,breeding of new species,and development of reinforced composites.This article presents an overview of research on mechanics properties of crop stalks,including the mechanics measurement experiments,mechanics model and theory,microstructure and composition,and mechanics simulation experiments.The current research relies mainly on the mechanics theory of engineering material,neglects the fact that crop stalk is a biological material with complex structure,anisotropy,viscoelasticity and rheology.In addition,some simulation experiments are conducted in idealized conditions,beyond the application.Hence,the paper proposes to create new material models for crop stalks,conduct more practical and effective simulation experiments,and promote cooperation in interdisciplinary studies for the development of mechanics properties of crop stalks in the future research.

Modeling and simulation of temperature control system in plant factory using energy balance

Closed production systems,such as plant factories and vertical farms,have emerged to ensure a sustainable supply of fresh food,to cope with the increasing consumption of natural resource for the growing population.In a plant factory,a microclimate model is one of the direct control components of a whole system.In order to better realize the dynamic regulation for the microclimate model,energy-saving and consumption reduction,it is necessary to optimize the environmental parameters in the plant factory,and thereby to determine the influencing factors of atmosphere control systems.Therefore,this study aims to identify accurate microclimate models,and further to predict temperature change based on the experimental data,using the classification and regression trees(CART)algorithm.A random forest theory was used to represent the temperature control system.A mechanism model of the temperature control system was proposed to improve the performance of the plant factories.In terms of energy efficiency,the main influencing factors on temperature change in the plant factories were obtained,including the temperature and air volume flow of the temperature control device,as well as the internal relative humidity.The generalization error of the prediction model can reach 0.0907.The results demonstrated that the proposed model can present the quantitative relationship and prediction function.This study can provide a reference for the design of high-precision environmental control systems in plant factories.