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.

Investigation of Projectile Impact Behaviors of Graphene Aerogel Using Molecular Dynamics Simulations

Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.

Study of the Relationship Between New Ionic Interaction Parameters and Salt Solubility in Electrolyte Solutions Based on Molecular Dynamics Simulation

Studying the relationship between ionic interactions and salt solubility in seawater has implications for seawater desalination and mineral extraction.In this paper,a new method of expressing ion-to-ion interaction is proposed by using molecular dynamics simulation,and the relationship between ion-to-ion interaction and salt solubility in a simulated seawater water-salt system is investigated.By analyzing the variation of distance and contact time between ions in an electrolyte solution,from both spatial and temporal perspectives,new parameters were proposed to describe the interaction between ions:interaction distance(ID),and interaction time ratio(ITR).The best correlation between characteristic time ratio and solubility was found for a molar ratio of salt-to-water of 10:100 with a correlation coefficient of 0.96.For the same salt,a positive correlation was found between CTR and the molar ratio of salt and water.For type 1-1,type 2-1,type 1-2,and type 2-2 salts,the correlation coefficients between CTR and solubility were 0.93,0.96,0.92,and 0.98 for a salt-to-water molar ratio of 10:100,respectively.The solubility of multiple salts was predicted by simulations and compared with experimental values,yielding an average relative deviation of 12.4%.The new ion-interaction parameters offer significant advantages in describing strongly correlated and strongly hydrated electrolyte solutions.

Molecular Dynamics Simulation of Shock Response of CL-20 Co-crystals Containing Void Defects

To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.

Flow field, sedimentation, and erosion characteristics around folded linear HDPE sheet sand fence: Numerical simulation study

Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.

Exploring the molecular mechanism of action of curcumin for the treatment of diabetic retinopathy,using network pharmacology,molecular docking,and molecular dynamics simulation

Background:Based on network pharmacology and molecular docking,the present study investigated the mechanism of curcumin(CUR)in diabetic retinopathy treatment.Methods:Based on the DisGeNET,Swiss TargetPrediction,GeneCards,Online Mendelian Inheritance in Man,Gene Expression Omnibus,and Comparative Toxicogenomics Database,the intersection core targets of CUR and diabetic retinopathy were identified.The intersection target was imported into the STRING database to obtain the protein-protein interaction map.According to the Database for Annotation,Visualization and Integrated Discovery database,the intersected targets were enriched in Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes pathways.Then Cytoscape 3.9.1 is used to make the drug-target-disease-pathway network.The mechanism of CUR and diabetic retinopathy was further verified by molecular docking and molecular dynamics simulation.Results:There were 203 intersecting targets of CUR and diabetic retinopathy identified.1320 GO entries were enriched for GO functions,which were primarily involved in the composition of cells such as identical protein binding,protein binding,enzyme binding,etc.It was found that 175 pathways were enriched using Kyoto Encyclopedia of Genes and Genomes pathway enrichment methods,which were mainly included in the lipid and atherosclerosis,AGE-RAGE signaling pathway in diabetic complications,pathways in cancer,etc.In the molecular docking analysis,CUR was found to have a good ability to bind to the core targets of albumin,IL-1B,and IL-6.The binding of albumin to CUR was further verified by molecular dynamics simulation.Conclusion:As a result of this study,CUR may exert a role in the treatment of diabetic retinopathy through multi-target and multi-pathway regulation,which indicates a possible direction of future research.

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.

In situ TEM observations and molecular dynamics simulations of deformation defect activities in Mg via nanoindentation

In this work, we performed in situ nanoindentation in TEM to capture the real-time dislocation and twinning activities in pure Mg during loading and unloading. We demonstrated that the screw component of dislocations glides continuously, while the edge components rapidly become sessile during loading. The twin tip propagation is intermittent, whereas the twin boundary migration is more continuous. During unloading, we observed the elastic strain relaxation causes both dislocation retraction and detwinning. Moreover,we note that the plastic zone comprised of dislocations in Mg is well-defined, which contrasts with the diffused plastic zones observed in face-centered cubic metals under the nanoindentation impressions. Additionally, molecular dynamics simulations were performed to study the formation and evolution of deformation-induced crystallographic defects at the early stages of indentation. We observed that,in addition to dislocations, the I1stacking fault bounded with a <1/2c+p> Frank loop can be generated from the plastic zone ahead of the indenter, and potentially serve as a nucleation source for abundant dislocations observed experimentally. These new findings are anticipated to provide new knowledge on the deformation mechanisms of Mg, which are difficult to obtain through conventional ex situ approaches. These observations may serve as a baseline for simulation work that investigate the dynamics of dislocation slip and twinning in Mg and alloys.

Experimental and numerical simulation study of the effect for the anode positions on the discharge characteristics of 300W class low power Hall thrusters

Low-power Hall thruster(LHT) generally has poor discharge efficiency characteristics due to the large surface-to-volume ratio.Aiming to further refine and improve the performance of 300 W class LHT in terms of thrust and efficiency,and to obtain the most optimal operating point,the experimental study of the discharge characteristics for three different anode positions was conducted under the operation of various discharge voltages(100-400 V) and anode mass flow rates(0.65 mg·s-1and 0.95 mg·s-1).The experimental results indicated that the thruster has the most excellent performance in terms of thrust and efficiency etc at a channel length of 27 mm for identical operating conditions.In addition,particle in cell simulations,employed to reveal the underlying physical mechanisms,show that the ionization and acceleration zone is pushed downwards towards the channel exit as the anode moves towards the exit.At the identical operating point,when the channel length is reduced from 32 to 27 mm,the ionization and acceleration zone moves towards the exit,and the parameters such as thrust and efficiency increase due to the high ionization rate,ion number density,and axial electric field.When the channel length is further moved to 24 mm,the parameters in terms of thrust(F) and efficiency(ηa) are reduced as a result of the decreasing ionization efficiency(ηm) and the larger plume divergence angle(α).In this paper,the results indicated that an optimum anode position(ΔL=27 mm) exists for the optimum performance.

Understanding the interfacial behaviors of benzene alkylation with butene using chloroaluminate ionic liquid catalyst: A molecular dynamics simulation

To better understand the benzene alkylation with chloroaluminate ionic liquids(ILs) as catalyst, the interfacial properties between the benzene/butene binary reactants and chloroaluminate ILs with varying cation alkyl chain length and different anions were investigated using molecular dynamics(MD) simulations. The results indicate that ILs can obviously improve the interfacial width, solubility and diffusion of reactants compared to H_(2)SO_(4). The longer alkyl chains of cations present a density enrichment at the interface and protrude into the binary reactants phase. Furthermore, the ILs consisting of 1-octyl-3-methylimidazolium cations([Omim]^(+)) and the stronger acidity heptachlorodialuminate anions([Al_(2)Cl_(7)]^(-)) are more beneficial to promote the interfacial width and facilitate the dissolution and diffusion of benzene in both the IL bulk and the interfacial region in comparison to the ones with shorter alkyl chains cations and weaker acidity anions. The information gives us a better guideline for the design of ILs for benzene alkylation.

Molecular Dynamics Simulation of Interface Properties between Water-Based Inorganic Zinc Silicate Coating Modified by Organosilicone and Iron Substrate

The interface properties of Fe(101)/zinc silicate modified by organo-siloxane(KH-570)was studied by using the method of molecular dynamics simulation.By calculating the temperature and energy fluctuation of equilibrium state,equilibrium concentration distribution,MSD of layer and different groups,and interaction energy of two interface models,the influencing mechanism on the interface properties of adding organosiloxane into coating system was studied at the atomic scale.It shows that the temperature and energy of interface oscillated in a small range and it was exited in a state of dynamic equilibrium within the initial simulation stage(t<20 ps).It can be seen from the multiple peak states of concentration distribution that the iron substrate,organo-siloxane and zinc silicate are distributed in the form of a concentration gradient in the real environment.The rapid diffusion of free zinc powder in zinc silicate coating was the essential reason that affected the comprehensive properties of coating.The interface thickness decreased from 7.45 to 6.82Å,the MSD of free zinc powder was effectively reduced,and the interfacial energy was increased from 104.667 to 347.158 kcal/mol after being modified by organo-siloxane.

Molecular Dynamics Simulation of Hydrate Decomposition Under Action of Alcohol Inhibitors

The molecular dynamics method is used to investigate decomposition of methane hydrate at different temperatures,pressures and concentrations of inhibitor.By analyzing the parameters of system conformation,mean square displacement and radial distribution function,the decomposition of hydrate in the presence of alcohol inhibitors ethylene glycol and glycerol is explored.The results show that the hydroxyl groups in alcohol molecules can destroy the cage structure of hydrate,and form hydrogen bonds with nearby water molecules to effectively prevent the reformation of hydrate.Therefore,ethylene glycol and glycerol serve as inhibitors of methane hydrate,furthermore,in terms of inhibition effect,glycerol is better than ethylene glycol by comparing rate of hydrate decomposition.

Numerical Simulation of Contamination Accumulation Characteristics of Composite Insulators in Salt Fog Environment

To investigate the fouling characteristics of the composite insulator surface under the salt fog environment,the FXBW-110/120-2 composite insulator was taken as the research object.Based on the field-induced charge mechanism,the multi-physical field coupling software COMSOL was used to numerically simulate the fouling characteristics,explored the calculation method of ESDD,and demonstrated its rationality.Based on this method,the pollution characteristics of the composite insulator under the pollution fog environment were studied,and the influence of wind speed,droplet size,and voltage type on the pollution characteristics of the composite insulator was analyzed.The results showed that:with the increase in wind speed,the amount of accumulated pollution of insulator increases in the range of droplet size,and the relationship between wind speed and accumulated pollution is approximately linear;at the same wind speed,the amount of accumulated pollution increases with the increase of droplet size under the action of DC voltage;when there is no voltage,the amount of dirt on the upper surface of the insulator is more than that on the lower surface,while it is the opposite under DC voltage.

基于Plant Simulation的化纤自动落丝系统仿真实验分析

基于Plant Simulation软件,构建一种化纤自动落丝系统的3D数字化模型,对系统效率关键输入因子进行参数化设置,将系统各工位三维数字模型导入软件中建立层次化运动机构图形结构,并运用Simltalk语言实现三维动作仿真,通过基础物理参数设置,保证了仿真数字化模型与现实系统更具一致性。通过多级实验设计分析了系统影响因子对于系统效率的影响特性曲线,并进一步通过动态参数化实验方法,计算出双输入因子对系统效能影响的敏感度。该实验结果可为化纤自动落丝系统的建设成本控制与可行性分析提供指导,具有较好的工程应用价值。

Elucidating the molecular targets of Curcuma longa for breast cancer treatment using network pharmacology,molecular docking and molecular dynamics simulation

Background:To elucidate the molecular mechanisms of Curcuma longa(C.longa)in breast cancer treatment.Methods:Phytocompounds of C.longa were obtained from Dr.Duke’s Phytochemical and Ethnobotanical Database.Potential active targets were retrieved from Bindingdb,SEA and Swiss Target Prediction databases.Breast cancer targets were retrieved from the Therapeutic Target Database.Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were done using DAVID and KOBAS3.0 databases respectively.The Cytoscape software was used to construct the phytocompound-target-pathway network.The PyRx and Desmond software were utilized for molecular docking and molecular dynamics simulation respectively.Results:Out of one hundred and fifty-six phytocompounds,fifty-four modulated proteins involved in breast cancer.Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated C.longa exerts its therapeutic effect through regulating several key pathways.Molecular docking analysis revealed that most phytocompounds of C.longa had a good affinity with the key targets.Molecular dynamics simulation showed that ethinylestradiol formed stable ligand-protein complexes.Conclusion:The results of this study will enhance our understanding of the potential molecular mechanisms by which C.longa inhibits breast cancer and lay a foundation for future experimental studies.

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.

Coupled CFD-DEM Numerical Simulation of the Interaction of a Flow-Transported Rag with a Solid Cylinder

A coupled Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)approach is used to calculate the interaction of a flexible rag transported by a fluid current with a fixed solid cylinder.More specifically a hybrid Eulerian-Lagrangian approach is used with the rag being modeled as a set of interconnected particles.The influence of various parameters is considered,namely the inlet velocity(1.5,2.0,and 2.5 m/s,respectively),the angle formed by the initially straight rag with the flow direction(45°,60°and 90°,respectively),and the inlet position(90,100,and 110 mm,respectively).The results show that the flow rate has a significant impact on the permeability of the rag.The higher the flow rate,the higher the permeability and the rag speed difference.The angle has a minor effect on rag permeability,with 45°being the most favorable angle for permeability.The inlet position has a small impact on rag permeability,while reducing the initial distance between the rag an the cylinder makes it easier for rags to pass through.

Refinement of Adaptive Dynamical Simulation of Quantum Mechanical Double Slit Interference Phenomenon

We applied adaptive dynamics to double slit interference phenomenon using particle model and obtained partial successful results in our previous report. The patterns qualitatively corresponded well with experiments. Several properties such as concave single slit pattern and large influence of slight displacement of the emission position were different from the experimental results. In this study we tried other slit conditions and obtained consistent patterns with experiments. We do not claim that the adaptive dynamics is the principle of quantum mechanics, but the present results support the probability of adaptive dynamics as the candidate of the basis of quantum mechanics. We discuss the advantages of the adaptive dynamical view for foundations of quantum mechanics.

Exploring unbinding mechanism of drugs from SERT via molecular dynamics simulation and its implication in antidepressants

The human serotonin transporter(SERT)terminates neurotransmission by removing serotonin from the synaptic cleft,which is an essential process that plays an important role in depression.In addition to natural substrate serotonin,SERT is also the target of the abused drug cocaine and,clinically used antidepressants,escitalopram,and paroxetine.To date,few studies have attempted to investigate the unbinding mechanism underlying the orthosteric and allosteric modulation of SERT.In this article,the conserved property of the orthosteric and allosteric sites(S1 and S2)of SERT was revealed by combining the high resolutions of x-ray crystal structures and molecular dynamics(MD)simulations.The residues Tyr95 and Ser438 located within the S1 site,and Arg104 located within the S2 site in SERT illustrate conserved interactions(hydrogen bonds and hydrophobic interactions),as responses to selective serotonin reuptake inhibitors.Van der Waals interactions were keys to designing effective drugs inhibiting SERT and further,electrostatic interactions highlighted escitalopram as a potent antidepressant.We found that cocaine,escitalopram,and paroxetine,whether the S1 site or the S2 site,were more competitive.According to this potential of mean force(PMF)simulations,the new insights reveal the principles of competitive inhibitors that lengths of trails from central SERT to an opening were~18A for serotonin and~22 A for the above-mentioned three drugs.Furthermore,the distance between the natural substrate serotonin and cocaine(or escitalopram)at the allosteric site was~3A.Thus,it can be inferred that the potent antidepressants tended to bind at deeper positions of the S1 or the S2 site of SERT in comparison to the substrate.Continuing exploring the processes of unbinding four ligands against the two target pockets of SERT,this study observed a broad pathway in which serotonin,cocaine,escitalopram(at the S1 site),and paroxetine all were pulled out to an opening between MT1b and MT6a,which may be helpful to understand the dissociation mechanism of antidepressants.

Dynamic simulation analysis of molten salt reactor-coupled air-steam combined cycle power generation system

A nonlinear dynamic simulation model based on coordinated control of speed and flow rate for the molten salt reactor and combined cycle systems is proposed here to ensure the coordination and stability between the molten salt reactor and power system.This model considers the impact of thermal properties of fluid variation on accuracy and has been validated with Simulink.This study reveals the capability of the control system to compensate for anomalous situations and maintain shaft stability in the event of perturbations occurring in high-temperature molten salt tank outlet parameters.Meanwhile,the control system’s impact on the system’s dynamic characteristics under molten salt disturbance is also analyzed.The results reveal that after the disturbance occurs,the controlled system benefits from the action of the control,and the overshoot and disturbance amplitude are positively correlated,while the system power and frequency eventually return to the initial values.This simulation model provides a basis for utilizing molten salt reactors for power generation and maintaining grid stability.