Temperature effect on nanotwinned Ni under nanoindentation using molecular dynamic simulation

Temperature effect on atomic deformation of nanotwinned Ni (nt-Ni) under localized nanoindentation is investigated in comparison with nanocrystalline Ni (nc-Ni) through molecular simulation.The nt-Ni exhibits enhanced critical load and hardness compared to nc-Ni,where perfect,stair-rod and Shockley dislocations are activated at (111),(111) and (111) slip planes in nt-Ni compared to only SSockley dislocation nucleation at (111) and (111) slip planes of nc-Ni.The nt-Ni exhibits a less significant indentation size effect in comparison with nc-Ni due to the dislocation slips hindrance of the twin boundary.The atomic deformation associated with the indentation size effect is investigated during dislocation transmission.Different from the decreasing partial slips parallel to the indenter surface in nc-Ni with increasing temperature,the temperaturedependent atomic deformation of nt-Ni is closely related to the twin boundary:from the partial slips parallel to the twin boundary (~10 K),to increased confined layer slips and decreased twin migration(300 K–600 K),to decreased confined layer slips and increased dislocation interaction of dislocation pinning and dissociation (900 K–1200 K).Dislocation density and atomic structure types through quantitative analysis are implemented to further reveal the above-mentioned dislocation motion and atomic structure alteration.Our study is helpful for understanding the temperature-dependent plasticity of twin boundary in nanotwinned materials.

Molecular Dynamic Simulation on the Absorbing Process of Isolating and Coating of α-olefin Drag Reducing Polymer

The absorbing process in isolating and coating process of α-olefin drag reducing polymer was studied by molecular dynamic simulation method, on basis of coating theory of α-olefin drag reducing polymer particles with polyurethane as coating material. The distributions of sodium laurate, sodium dodeeyl sulfate, and sodium dodeeyl benzene sulfonate on the surface of α-olefin drag reducing polymer particles were almost the same, but the bending degrees of them were obviously different. The bending degree of SLA molecules was greater than those of the other two surfactant molecules. Simulation results of absorbing and accumulating structure showed that, though hydrophobie properties of surfactant molecules were almost the same, water density around long chain sulfonate sodium was bigger than that around alkyl sulfate sodium. This property goes against useful absorbing and accumulating on the surface of α-olefin drag reducing polymer particles; simulation results of interactions of different surfactant and multiple hydroxyl compounds on surface of particles showed that, interactions of different surfaetant and one kind of multiple hydroxyl compound were similar to those of one kind of surfaetant and different multiple hydroxyl compounds. These two contrast types of interactions also exhibited the differences of absorbing distribution and closing degrees to surface of particles. The sequence of closing degrees was derived from simulation; control step of addition polymerization interaction in coating process was absorbing mass transfer process, so the more closed to surface of particle the multiple hydroxyl compounds were, the easier interactions With isoeyanate were. Simulation results represented the compatibility relationship between surfactant and multiple hydroxyl compounds. The isolating and coating processes of α-olefin drag reducing polymer were further understood on molecule and atom level through above simulation research, and based on the simulation, a referenced theoretical basis was provided for practical optimal selection and experimental preparation of α-olefin drag reducing polymer particles suspension isolation agent.

Molecular Dynamic Simulation of Melting Points of Trans-1,4,5,8- tetranitro-1,4,5,8-tetraazadacalin （TNAD） with Some Propellants

Molecular dynamic simulation was employed to predict the melting points Tm of TNAD/HMX, TNAD/RDX, TNAD/DINA, and TNAD/DNP systems （tans-1,4,5,8- tetranitro-1,4,5,8-tetraazadacalin （TNAD）, dinitropiperazine （DNP）, cyclotetramethylenetetranitroamine （HMX）, cyclotrimethylenetrinitramine （RDX）, and N-nitrodihydroxyethylaminedinitrate （DINA））. Tm was determined from the inflexion point on the curve of mean specific volume vs. temperature. The result shows that the Tm values of TNAD/HMX, TNAD/RDX, and TNAD/DINA systems are 500, 536, and 488 K, respectively. The TNAD/DNP system has no obvious Tm value, which shows the system is insoluble. Using Tm, the solubility of the four systems was analyzed. The radial distribution functions of the four systems were analyzed and the main intermolecular forces between TNAD and other energetic components are short-range interactions. The better the solubility is, the stronger the intermoleenlar interaction is. In addition, the force field energy at different temperature was also analyzed to predict Tm of the four systems.

Critical anomaly and finite size scaling of the self-diffusion coefficient for Lennard Jones fluids by non-equilibrium molecular dynamic simulation

We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with temperature decreases by increasing density. For density ρ＊ = ρσ3 = 0.84 we observe a peak at the value of the self-diffusion coefficient and the critical temperature T＊ = kT/ε = 1.25. The value of the self-diffusion coefficient strongly depends on system size. The data of the self-diffusion coefficient are fitted to a simple analytic relation based on hydrodynamic arguments. This correction scales as N-α, where α is an adjustable parameter and N is the number of particles. It is observed that the values of a 〈 1 provide quite a good correction to the simulation data. The system size dependence is very strong for lower densities, but it is not as strong for higher densities. The self-diffusion coefficient calculated with non-equilibrium molecular dynamic simulations at different temperatures and densities is in good agreement with other calculations fronl the literature.

3-D fracture network dynamic simulation based on error analysis in rock mass of dam foundation

Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.

Research on the Dynamic Simulation System of Maize Growth

Dynamic simulation system of maize growth is developed by the physiological and ecological model,morphological structure model,computer science and virtual reality technology,to improve the level of precise operation of maize production.The computer graphics algorithms,virtual reality technology,animation design and information integration technology are applied to maize production by this system.establishment of dynamic simulation system of maize growth is conducive to raise level of precise operation in maize production.The system also can assist the relevant production research and testing,to reduce cost and improve efficiency.

Molecular dynamic simulation of lubricant spreading： effect from the substrate and endbead

Molecular dynamic simulations based on a coarse-gralned, bead-spring model are adopted to investigate the spreading of both nonfunctional and functional perfluoropolyether （PFPE） on solid substrates. For nonfunctional PFPE, the spreading generally exhibits a smooth profile with a precursor film. The spreading profiles on different substrates are compared, which indicate that the bead-substrate interaction has a significant effect on the spreading behaviour, especially on the formation of the precursor film. For functional PFPE, the spreading generally exhibits a complicated terraced profile. The spreading profiles with different endbeads are compared, which indicate that the endbead-substrate interaction and the endbead-endbead interaction, especially the latter, have a significant effect on the spreading behaviour.

Molecular dynamic simulation of the thermodynamic and kinetic properties of nucleotide base pair

A nucleotide base pair is the basic unit of RNA structures. Understanding the thermodynamic and kinetic properties of the closing and opening of a base pair is vital for quantitative understanding the biological functions of many RNA molecules. Due to the fast transition rate, it is difficult to directly observe opening and closing of single nucleic acid base pair in experiments. This review will provide a brief summary of the studies about the thermodynamic and kinetic properties of a base pair opening and closing by using molecular dynamic simulation methods.

Molecular Dynamic Simulation of Lattice Distortion Region Produced by Rounded Grain Boundary in Nanocrystalline Materials

The distortion structure in nanocrystalline NiAl is studied using molecular dynamics simulation. The rounded grain boundaries in these nanograins are a direct source for the observed lattice distortion. The change of grain size affects directly the volume fraction of the distorted lattice in the nanograin.

Study of structural and magnetic properties of Fe(80)P-9B(11) amorphous alloy by ab initio molecular dynamic simulation

The structural and magnetic properties of Fe80P9B11 amorphous alloy are investigated through ab initio molecular dynamic simulation. The structure evolution of Fe（80）P9B（11） amorphous alloy can be described in the framework of topological fluctuation theory, and the fluctuation of atomic hydrostatic stress gradually decreases upon cooling. The left sub peak of the second peak of Fe–B partial pair distribution functions（PDFs） becomes pronounced below the glass transition temperature, which may be the major reason why B promotes the glass formation ability significantly. The magnetization mainly originates from Fe 3d states, while small contribution results from metalloid elements P and B. This work may be helpful for developing Fe-based metallic glasses with both high saturation flux density and glass formation ability.

Atomic scale imaging of monocrystalline Si (001) surface by molecular dynamic simulation

Non-contact atomic force microscopy（nc-AFM） atomic-scale imaging process of monocrystalline silicon surface using capped single-wall carbon nanotube tip is simulated by molecular dynamic method. The simulation resuhs show that the nc-AFM imaging force mainly comes from the C-Si and C-C chemical covalent bonding forces, especially the former, the nonbonding Van der Waals force change is small during the range of stable imaging height. When the tip-surface distance is smaller than the stable imaging height, several neighboring carbon atoms at the tip apex are attracted, and some of them jump onto the sample surface. Finally the tip apex configuration is destroyed with the tip indenting further.

Quasimolecular Dynamic Simulation for Bending Fracture of Laminar Composite Materials

Recently, quasimolecular dynamics has been successfully used to simulate the deformation characteristics of actual size solid materials. In quasimolecular dynamics, which is an attempt to bridge the gap between atomistic and continuum simulations, molecules are aggregated into large units, called quasimolecules, to evaluate large scale material behavior. In this paper, a 2-dimensional numerical simulation using quasimolecular dynamics was performed to investigate laminar composite material fractures and crack propagation behavior in the uniform bending of laminar composite materials. It was verified that under bending deformation laminar composite materials deform quite differently from homogeneous materials

Studies on sensitivity to tension and gating pathway of MscL by molecular dynamic simulation

Mechanosensitive（MS） ion channels play an important role in various physiological processes.Although the determination of the structure of mechanosensitive channel of large conductance（MscL） makes the simulation study possible,it has not so far been possible to directly simulate the gating mechanism of MscL in atomic detail.In this article,MscL has been studied via molecular dynamic（MD） simulations to gain a detailed description of the sensitivity to lateral tension and the gating pathway.MscL undergoes conformational rearrangement in sustaining lateral tension,and the open state is obtained when 2.0 MPa lateral tension is directly applied on the pure protein.During the opening process,Loop region responds to tension first,and the mechanical sensitivity is followed by S1 domain.Transmembrane（TM） bundle is the key position for channel opening,and the motion of TM1 helices finally realizes the significant expansion of the constricted gating pore.C-terminus domain presents expansion later during the TM opening.In our study,return of the whole protein to the initial closed state is achieved only in the early opening stage.During the relaxation from the open state,the TM helices are the most mobile domain,which is different from the opening process.

Design and dynamic simulation of hydraulic system of a new automatic transmission

A new hydraulic system of a novel automatic transmission (AT) was designed. The dimension and structure of valves and cylinders were designed by theoretical calculation. The dynamic simulation model of hydraulic system of AT was established by ITI-SimulationX. Simulation results and theoretical design results were compared to confirm the simulation model. Based on the confirmed simulation model, the simulation results of pressure and flow of the hydraulic system were analyzed. The dynamic simulation method is very helpful for designing and analyzing the performance of hydraulic system and further optimization design. The theoretical design method and dynamic simulation model are feasible for the real industrial applications. The research results can be used in hydraulic system design and optimization.

A dynamic simulation of the wheel–rail impact caused by a wheel flat using a 3D rolling contact model

A three-dimensional （3-D） wheel-rail rolling contact model with a wheel fiat was built using commercial software Hypermesh, and the dynamic finite element simulation was conducted using LS-DYNA 3D/explicit code. Influences of the train speed, flat length and axle load on the vertical wheel-rail impact response were discussed, respectively. The results show that the maximum vertical wheel-rail impact force induced by the wheel flat is higher than that generated by the perfect wheel, and these two dynamic impact forces are much greater than the static axle load. Besides, the maximum von Mises equivalent stress and maximum equivalent plastic strain are observed on the wheel-rail contact surface, and both of them as well as the maximum wheel-rail impact force are sensitive to train speed, fiat length and axle load.

Theoretical design and dynamic simulation of new mining paths of tracked miner on deep seafloor

With comprehensive considerations of the operational safety and collection efficiency for the tracked miner collecting the seafloor poly-metallic nodules, two new improved mining paths for the miner on the deep seafloor were proposed. Compared to the conventional mining path, the design principles and superiorities of the two new paths are that the miner turning with relative long radius should avoid large sinkage and high slip, so as to ensure its operational safety, while the space between its straight-line trajectories before and after the turning is optimum, which is designed as the total width of the miner, and collect nodules as more as possible, so as to ensure its collection efficiency. To realize the new mining paths, theoretical designs and quantitative calculations were carried out to determine the exact positions for the speed controls of the miner during its whole operation process. With the new dynamic model of the miner, and through regulations of the speeds of the left and right tracks of the miner on the exact motion positions according to the theoretical calculations, the two new improved mining paths for the miner on the seafloor were successfully simulated, thus the turning radius of the miner in the simulation is about 21.8 m, while the distance between the straight-line trajectories before and after the turning is about 5.2 m. The dynamic simulation results preliminarily prove the feasibility of these two new mining paths, and further can provide important theoretical guidance and useful technical reference for the practical tracked miner operation and control on the seafloor.

Succession Features and Dynamic Simulation of Subalpine Forest in the Gongga Mountain, China

The Gongga Mountain of eastern Tibet Plateau is a representative of the alpine regions with high peaks and deep valleys. Climate change over the last thousand years has controlled the dynamics of glacier and debris flow occurrence, which resulted in substantial changes in the mountainous environment. The authors surveyed the community structure of primary forests in Gongga Mountain and forest succession processes in woodland plots. The changing features in the subalpine environment are discussed in this paper. Tree species and sizes between the glacier shrinking areas and debris flow fans in Hailuogou Valley are compared. The pioneer species that settle in debris flow fans and the glacier shrinking areas are Salix spp. and Populus purdomii. Abies fabri and Picea brachytyla are the climax tree species. The succession process of primary vegetation in Hailuogou (2700 ~ 3200 m) can be divided into four stages: Slash surface (20~ 200 yr) Salix-Populus seeding community (10 ~ 30 yr) Populus-Salix sapling community (30 ~ 100 yr) Populus-Abies mixed community (50 ~ 100 yr) Abies-Picea climax In a natural and undisturbed environment, trees compete for light, water and nutrients. Disaster disturbance in mountains is a very important driving factor for regeneration of woody plants. Repeated destruction of forests by glacier movement or debris flows generated additional forest gaps that allow young plants to grow. In this study the Gongga Forest Succession Model (GFSM) was developed for simulation of forest community succession processes on different scales in Gongga Mountain. A soil succession module was added to the GFSM model to simulate soil formation and chemical element change of woodland. In order to represent major features of forests in Southwestern China, many field works has been done to identify ecological parameters of various trees in the subalpine region. On the basis of simulation of tree life history, the GFSM combines forest succession with soil change in both material components and nutrition content. The Monte-Carlo method was applied to simulate random weather fluctuation and the uncertainty of tree death. These modeled processes agreed with the field investigation results in this region. The elevation distribution of different tree species was also simulated; and the results are consistent with field observations of ecological features of tree species. The modeling approach reflects well the succession dynamics of primary forests in Southwestern China. These results are very useful for improving the management policies and prediction technology for restoration and conservation of primary forests in Southwestern China.

Optimal Design and Dynamic Simulation of Mountain Bike with Rear Suspension

This paper investigates the dynamic design methodology of mountain bikes with rear suspension. Firstly, a multi-rigid body dynamic model of rider and mountain bike coupled system is constructed. The rider model includes 19 skeletons, 18 joints and 118 main muscles. Secondly, to validate the feasibility of the model, an experiment test is designed to reflect the real cycling status. Finally, aiming at enhancing the performance of the rider vibration comfort, the scale parameters of rear suspension are optimized with computer simulation and uniform design. The mathematical model in the vibration performance and the design variables is constructed with regression analysis. The result shows that when the length of side link is 90 mm, the length of connected rod is 336.115 1 mm and the included angle between absorber and side link is 60°, the mountain bike has better vibration comfort. This study and relevant conclusions are of practical importance to the design of the mountain bike＇s rear suspension system.

Phase Transition Temperature of HPT(2,3,6,7,10,11-hexa-n-pentyloxytriphenylene) from Molecule Dynamic Simulation

Molecule dynamics simulation was used on HPT(2,3,6,7,10,11-hexa-n-pentyloxytriphenylene), which is a discotic Liquid crystal. From analyzing the energy and displacement varying with the temperature, the phase transition temperature of PM6MPP can be predicted. The deviations of T-g, T-m and T-i due to the MD time scale are small enough that it should be possibly used to predict the material properties especially when more powerful computers are available.

Process Modeling and Dynamic Simulation for EAST Helium Refrigerator

In this paper,the process modeling and dynamic simulation for the EAST helium refrigerator has been completed.The cryogenic process model is described and the main components are customized in detail.The process model is controlled by the PLC simulator,and the realtime communication between the process model and the controllers is achieved by a customized interface.Validation of the process model has been confirmed based on EAST experimental data during the cool down process of 300-80 K.Simulation results indicate that this process simulator is able to reproduce dynamic behaviors of the EAST helium refrigerator very well for the operation of long pulsed plasma discharge.The cryogenic process simulator based on control architecture is available for operation optimization and control design of EAST cryogenic systems to cope with the long pulsed heat loads in the future.