Frontogenesis and Frontolysis of a Cold Filament Driven by the Cross-Filament Wind and Wave Fields Simulated by a Large Eddy Simulation

The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and wave fields are studied.The results show that the intense wind and wave fields further break the symmetries of submesoscale flow fields and suppress the levels of filament frontogenesis. The changes of secondary circulation directions—that is, the conversion between the convergence and divergence of the surface cross-filament currents with the downwelling and upwelling jets in the filament center—are associated with the inertial oscillation. The filament frontogenesis and frontolysis caused by the changes of secondary circulation directions may periodically sharpen and smooth the gradient of submesoscale flow fields.The lifecycle of the cold filament may include multiple stages of filament frontogenesis and frontolysis.

Mathematical modeling and simulations of stress mitigation by coating polycrystalline particles in lithium-ion batteries

A chemo-mechanical model is developed to investigate the effects on the stress development of the coating of polycrystalline Ni-rich LiNixMnyCo_(z)O_(2)(x≥0.8)(NMC)particles with poly(3,4-ethylenedioxythiophene)(PEDOT).The simulation results show that the coating of primary NMC particles significantly reduces the stress generation by efficiently accommodating the volume change associated with the lithium diffusion,and the coating layer plays roles both as a cushion against the volume change and a channel for the lithium transport,promoting the lithium distribution across the secondary particles more homogeneously.Besides,the lower stiffness,higher ionic conductivity,and larger thickness of the coating layer improve the stress mitigation.This paper provides a mathematical framework for calculating the chemo-mechanical responses of anisotropic electrode materials and fundamental insights into how the coating of NMC active particles mitigates stress levels.

B-COPNA resin formation from ethylene tar light fractions:Process development and mechanical exploration by molecular simulation

An efficient utilization strategy of ethylene tar(ET),the main by-product of the ethylene cracking unit,is urgently required to meet demands for modern petrochemical industry.On the other hand,condensed polynuclear aromatic resin of moderate condensation degree(B-COPNA)is a widely used carbon material due to its superb processability,the production of which is,however,seriously limited by the high cost of raw materials.Under such context,an interesting strategy was proposed in this study for producing B-COPNA resin using crosslinked light fractions of ethylene tar(ETLF,boiling point<260℃)facilitated by molecular simulation.1,4-Benzenedimethanol(PXG)was first selected as the crosslinking agent according to the findings of molecular simulation.The effects of operating conditions,including reactions temperature,crosslinking agent,and catalyst content on the softening point and yield of B-COPNA resin products were then investigated to optimize the process.The reaction mechanism of resin production was studied by analyzing the molecular structure and transition state of ETLF and crosslinking agents.It was shown that PXG exhibited a superior capacity of withdrawing electrons and a higher electrophilic reactivity than other crosslinking agents.In addition to the highest yield and greatest heat properties,PXG-prepared resin contained the most condensed aromatics.The corresponding optimized conditions of resin preparation were 180℃,1:1.9(PXG:ETLF),and 3%(mass)of catalyst content with a resin yield of 78.57%.It was the electrophilic substitution reaction that occurred between the ETLF and crosslinking agent molecules that were responsible for the resin formation,according to the experimental characterization and molecular simulation.Hence,it was confirmed that the proposed strategy and demonstrated process can achieve a clean and high value-added utilization of ETLF via B-COPNA resin preparation,bringing huge economic value to the current petrochemical industry.

Facilitated Prediction of Micropollutant Degradation via UV-AOPs in Various Waters by Combining Model Simulation and Portable Measurement

The degradation of micropollutants in water via ultraviolet(UV)-based advanced oxidation processes(AOPs)is strongly dependent on the water matrix.Various reactive radicals(RRs)formed in UV-AOPs have different reaction selectivities toward water matrices and degradation efficiencies for target micropollutants.Hence,process selection and optimization are crucial.This study developed a facilitated prediction method for the photon fluence-based rate constant for micropollutant degradation(K′_(p,MP))in various UV-AOPs by combining model simulation with portable measurement.Portable methods for measuring the scavenging capacities of the principal RRs(RRSCs)involved in UV-AOPs(i.e.,HO^(·),SO_(4)^(·-),and Cl^(·))using a mini-fluidic photoreaction system were proposed.The simulation models consisted of photochemical,quantitative structure–activity relationship,and radical concentration steady-state approximation models.The RRSCs were determined in eight test waters,and a higher RRSC was found to be associated with a more complex water matrix.Then,by taking sulfamethazine,caffeine,and carbamazepine as model micropollutants,the k′_(p,MP) values in various UV-AOPs were predicted and further verified experimentally.A lower k′_(p,MP) was found to be associated with a higher RRSC for a stronger RR competition;for example,k′_(p,MP) values of 130.9 and 332.5 m^(2) einstein^(–1),respectively,were obtained for carbamazepine degradation by UV/H_(2)O_(2) in the raw water(RRSC=9.47×10^(4) s^(-1))and sand-filtered effluent(RRSC=2.87×10^(4) s^(-1))of a drinking water treatment plant.The developed method facilitates process selection and optimization for UV-AOPs,which is essential for increasing the efficiency and cost-effectiveness of water treatment.

Glutathione Peroxidase Revisited—Simulation of the Catalytic Cycle by Computer-Assisted Molecular Modelling

Glutathione peroxidase, the first example of selenoproteins identified in mammals, was subjected to force field calculations and molecular dynamics in order to enable a clearer comprehension of enzymatic selenium catalysis. Starting from the established X-ray structure of bovine GPX, all kinetically defined intermediates and enzyme substrate complexes were modelled. The models thus obtained support the hypothesis that the essential steps of the catalysis are three distinct redox changes of the active site selenium which, in the ground state, presents itself at the surface of selenoperoxidases as the center of a characteristic triad built by selenocysteine, glutarnine and tryptophan. In GPX, four arginine residues and a lysine residue provide an electrostatic architecture which, in each reductive step, directs the donor substrate GSH towards the catalytic center in such a way that 1ts sulfhydryl group must react with the selenium moiety. To this end, different equally efficient modes of substrate binding appear possible. The models are consistent with substrate specificity data, kinetic pattern and other functional characteristics of the enzyme. Comparison of molecular models of GPX with those of other members of the GPX superfamily reveals that the cosubstrate binding mechanisrns are unique for the classical type of cytosolic glutathione peroxidases but cannot operate e. g. in plasma GPX and phospholipid hydroperoxide GPX. The structural differences between the selenoperoxidases, shown to be relevant to their specificities, are discussed in terms of functional diversification within the GPX superfamily

Cold filament frontogenesis and frontolysis induced by thermal convection turbulence using large eddy simulation

The frontogenetic processes of a submesoscale cold filament driven by the thermal convection turbulence are studied by a non-hydrostatic large eddy simulation.The results show that the periodic changes in the direction of the cross-filament secondary circulations are induced by the inertial oscillation.The change in the direction of the secondary circulations induces the enhancement and reduction of the horizontal temperature gradient during the former and later inertial period,which indicates that the frontogenetical processes of the cold filament include both of frontogenesis and frontolysis.The structure of the cold filament may be broken and restored by frontogenesis and frontolysis,respectively.The magnitude of the down-filament currents has a periodic variation,while its direction is unchanged with time.The coupling effect of the turbulent mixing and the frontogenesis and frontolysis gradually weakens the temperature gradient of the cold filament with time,which reduces frontogenetical intensity and enlarges the width of cold filament.

Recent advances in protein conformation sampling by combining machine learning with molecular simulation

The rapid advancement and broad application of machine learning(ML)have driven a groundbreaking revolution in computational biology.One of the most cutting-edge and important applications of ML is its integration with molecular simulations to improve the sampling efficiency of the vast conformational space of large biomolecules.This review focuses on recent studies that utilize ML-based techniques in the exploration of protein conformational landscape.We first highlight the recent development of ML-aided enhanced sampling methods,including heuristic algorithms and neural networks that are designed to refine the selection of reaction coordinates for the construction of bias potential,or facilitate the exploration of the unsampled region of the energy landscape.Further,we review the development of autoencoder based methods that combine molecular simulations and deep learning to expand the search for protein conformations.Lastly,we discuss the cutting-edge methodologies for the one-shot generation of protein conformations with precise Boltzmann weights.Collectively,this review demonstrates the promising potential of machine learning in revolutionizing our insight into the complex conformational ensembles of proteins.

Process Optimization of Cu-en/AP Composite Microspheres Preparation by Electrostatic Spray Method Based on ANSYS Simulation

To investigate the process optimizationof Cu-en/AP composite microspheres preparation via electrostatic spraying,and to reveal the effects of droplet properties and flow rate variations on the experimental results during the electrostatic spraying process,the prepared process parameters of Cu-en/AP composite microspheres by electrostatic spray method under the orthogonal experimental design simulated by ANSYS(Fluent).The influence of flow rate,solvent ratio,and solid mass on the experimental results is examined using a controlled variable method.The results indicate that under the conditions of a flow rate of 2.67×10^(-3)kg/s an acetone-to-deionized water ratio of 1.5∶1.0,and a solid mass of 200 mg,the theoretical particle size of the composite microspheres can reach e nanoscale.Droplet trajectories in the electric field remain stable without significant deviation.The simulation results show that particle diameter decreases with increasing flow rate,with the trend leveling off around a flow rate of 1×10^(-3)kg/s.As the solvent ratio increases(with higher acetone content),particle diameter initially decreases,reaching a minimum around a ratio of 1.5∶1.0 before gradually increasing.Increasing the solid mass also reduces the particle diameter,with a linear increase in diameter observed at around 220 mg.Cu-en/AP composite microspheres with nanoscale dimensions were confirmed under these conditions by the final SEM images.

THEORETICAL PREDICTION OF TOOL-CHIP CONTACT LENGTH IN ORTHOGONAL METAL MACHINING BY COMPUTER SIMULATION

A method for determination of tool-chip contact length is theoreticallypresented in orthogonal metal machining. By using computer simulation and based on the analyses ofthe elastro-plastic deformation with lagrangian finite element method in the deformation zone, theaccumulated representative length of the low layer, the tool-chip contact length of the chipcontacting the tool rake are calculated, experimental studies are also carried out with 0.2 percentcarbon steel. It is shown that the tool-chip contact lengths obtained from computer simulation havea good agreement with those of measured values.

Numerical simulations of earthquake rupture induced by pressure perturbation

The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.

Molecular dynamics simulation study of nitrogen vacancy color centers prepared by carbon ion implantation into diamond

Nitrogen vacancy(NV)color centers in diamond have useful applications in quantum sensing andfluorescent marking.They can be gen-erated experimentally by ion implantation,femtosecond lasers,and chemical vapor deposition.However,there is a lack of studies of the yield of NV color centers at the atomic scale.In the molecular dynamics simulations described in this paper,NV color centers are pre-pared by ion implantation in diamond with pre-doped nitrogen and subsequent annealing.The differences between the yields of NV color centers produced by implantation of carbon(C)and nitrogen(N)ions,respectively,are investigated.It is found that C-ion implantation gives a greater yield of NV color centers and superior location accuracy.The effects of different pre-doping concentrations(400–1500 ppm)and implantation energies(1.0–3.0 keV)on the NV color center yield are analyzed,and it is shown that a pre-doping concentra-tion of 1000 ppm with 2 keV C-ion implantation can produce a 13%yield of NV color centers after 1600 K annealing for 7.4 ns.Finally,a brief comparison of the NV color center identification methods is presented,and it is found that the error rate of an analysis utiliz-ing the identify diamond structure coordination analysis method is reduced by about 7%compared with conventional identification+methods.

Experimental and Three-Dimensional Numerical Simulation of Phenomena Induced by Submerged Oblique Jet Scouring

Scouring experiments were conducted using a three-dimensional laser scanning technology for angles of the jet spanning the interval from 0°to 30°,and the characteristics of the scour hole in equilibrium conditions were investigated accordingly.The results indicate that the optimal scouring effects occur when the jet angle is in the ranges between 15°and 20°.Moreover,the dimensionless profiles of the scour hole exhibit a high degree of similarity at different jet angles.Numerical simulations conducted using the Flow-3D software to investigate the bed shear stress along the jet impingement surface have shown that this stress is influenced by both the resultant force and the jet impingement surface area.It reaches its maximum value when the jet is vertical,decreases rapidly as the jet starts to tilt,then increases slightly,and decreases again significantly when the angle exceeds 20°.

Study on Effect of Gd (III) Speciation on Ca (II) Speciationin Human Blood Plasma by Computer Simulation

Ca (II) speciation and effect of Gd (III) speciation on Ca (II) speciation in human blood plasma were studied by computer simulation. [CaHCO3](-) is a predominant compound species of Ca (II). Gd (III) can compete with Ca (II) for biological molecules. The presence of Gd (III) results in a increase of concentration of free Ca (II) and a decrease of concentration of Ca (II) compounds.

COMPUTERIZED SIMULATION OF MOLTEN SALT SOLUTION OF Li,KF,Cl SYSTEM BY MOLECULAR DYNAMIC METHOD

The structure and properties of molten salt solution o J Li,K|F,Cl system have been investiged by computerized simulation of molecular dynamic method.The partial RDF,the partial molar energy of mixing and the diffusion coeffients of Li^+,K^+,F^- and Cl^- have been calculated. The results are in agreement with the experimental values.The regularities of the distribution of ions and mieroscopic holes are discussed based on the results of computerized simulation.

Performance Simulation of a Double Tube Heat Exchanger Based on Different Nanofluids by Aspen Plus

A heat exchanger’s performance depends heavily on the operating fluid’s transfer of heat capacity and thermal conductivity.Adding nanoparticles of high thermal conductivity materials is a significant way to enhance the heat transfer fluid’s thermal conductivity.This research used engine oil containing alumina(Al_(2)O_(3))nanoparticles and copper oxide(CuO)to test whether or not the heat exchanger’s efficiency could be improved.To establish the most effective elements for heat transfer enhancement,the heat exchangers thermal performance was tested at 0.05%and 0.1%concentrations for Al_(2)O_(3)and CuO nanoparticles.The simulation results showed that the percentage increase in Nusselt number(Nu)for nanofluid at 0.05%particle concentration compared to pure oil was 9.71%for CuO nanofluids and 6.7%for Al_(2)O_(3)nanofluids.At 0.1%concentration,the enhancement percentage in Nu was approximately 23%for CuO and 18.67%for Al_(2)O_(3)nanofluids,respectively.At a concentration of 0.1%,CuO nanofluid increased the LMTD and overall heat transfer coefficient(U)by 7.24 and 5.91%respectively.Both the overall heat transfer coefficient(U)and the heat transfer coefficient(hn)for CuO nanofluid at a concentration of 0.1%increased by 5.91%and 10.68%,respectively.The effectiveness(εn)of a heat exchanger was increased by roughly 4.09%with the use of CuO nanofluid in comparison to Al_(2)O_(3)at a concentration of 0.1%.The amount of exergy destruction in DTHX goes down as Re and volume fractions go up.Moreover,at 0.05%and 0.1%nanoparticle concentrations,the percentage increase in dimensionless exergy is 10.55%and 13.08%,respectively.Finally,adding the CuO and Al_(2)O_(3)nanoparticles improved the thermal conductivity of the main fluid(oil),resulting in a considerable increase in the thermal performance and rate of heat transfer of a heat exchanger.

Dynamic consequences of mutating the typical HPGG motif of apocytochrome b_5 revealed by computer simulation

Apocytochrome b5 with a typical heme-binding motif of HPGG, and its variants with mutated motifs, GPGG, GPGH, HVGG, and HPGP, have been subjected to molecular dynamics simulation. Comparison of the dynamic consequences has revealed the crucial role of HPGG in assembling the heine group of cytochrome b5 and in modulating protein structure, property and function.

Studies on Insoluble Species of Gd(Ⅲ)in Human Blood Plasma by Computer Simulation

The insoluble species of Gd ( Ⅲ ) in human blood plasma were investigated by computer simulation. The distribution of the Gd(Ⅲ ) species was obtained. It was found that most of the Gd( Ⅲ ) ions were bound to phosphate to form precipitate GdPO4 at the concentration of 1. 000 10-7 mol/L and when the concentration of the Gd (Ⅲ ) increased to 3. 750 X 10-4 mol/L, in excess of the concentration of phosphate, the Gd ( Ⅲ ) ions were bound to carbonate to form another kind of precipitate, Gd2(CO3)3.

Molecular simulation study of the microstructures and properties of pyridinium ionic liquid[HPy][BF_(4)]mixed with acetonitrile

The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.

COMPUTERIZED SIMULATION OF THERMODYNAMICS FOR MOLTEN LiF-KF SOLUTION BY MONTE CARLO METHOD

The structure and behaviour of LiF-KF solution,as a typical common-anion system,has been simulated by Monte Carlo method.The calculation of partial radial distribution function of ions,heat of mixing and potential energy distribution shows that the average distance be- tween Li^+ and F^- ions will significantly narrow after mixing of molten LiF and KF.This is very similar to the lean-on-one-side effect in molten LiF-KCl solution.The calculated heat of mixing is in fair agreement with the measured one.The dominant source of the energy of mixing may be that the decrease of the repulsion energy between cations,the decrease of the attraction energy between cations and anions,and the decrease of the repulsion energy be- tween anions.

Computer Simulation of Distribution of Reinforced Particles in Functionally Graded Materials Fabricated by Centrifugal Casting

A mathematic model is developed with viscosity of molten matrix,centrifugalforce,casting temperature,mold temperature and other parameters taken intoconsideration for preditction of the distribution of reinforced particles during centrifugalcasting of FGM,and the simulation of distribution of reinforced particles and thesolidification process during centrifugal casting is performed with the aid of computergraphics.SiC_p/A356 FGM is fabricated by centrifugal casting.The results of computersimulation of distribution of reinforced particles are in good agreement with experimentalobservations.