Numerical analysis of electromagnetic force field through the use of two examples

In modern processing of conducting materials, such as steel,the time-varying electromagnetic field plays a key role in obtaining the desired microstructure or eliminating solidification defects such as porosity and segregation in cast billets. Up to now,few studies on the induced electromagnetic force （also called the Lorentz force） field in liquid metal have been reported. Compared with the magnetic field, the induced force field is the real and only direct cause for flow control. The electromagnetic force is comprised of two components. One is time-independent and the other is time- dependent. The time-dependent component varies with time in both amplitude and direction. When it reaches the extreme value,it can be one dozen times larger than the time-independent component. In this paper, a new method to quantitatively describe the induced electromagnetic force in liquid metal under a harmonic electromagnetic field,including both its time- independent and dependent components, was proposed based on the formula derivation from the data of amplitude and phase angle. Through this method ,the features of the time-dependent component were discussed, including the directions of rotation and the long axis. As a result, the force pattern was described. With two example calculations, the method was explained in detail. The results of both examples show that the force field in liquid metal can be divided into several regions with different force features. Example 1 shows the effect of coil position on the evolution of the force field pattern in liquid metal. Example 2 is a kind of stirring by the travelling magnetic field ,whose results present the sub-structures in metal and show that most of them have almost the same rotating direction.

Hugoniot curve calculation of nitromethane decomposition mixtures:A reactive force field molecular dynamics approach

We investigate the Hugoniot curve, shock-particle velocity relations, and Chapman-Jouguet conditions of the hot dense system through molecular dynamics （MD） simulations. The detailed pathways from crystal nitromethane to reacted state by shock compression are simulated. The phase transition of N2 and CO mixture is found at about 10 GPa, and the main reason is that the dissociation of the C-O bond and the formation of C-C bond start at 10.0-11.0 GPa. The unreacted state simulations of nitromethane are consistent with shock Hugoniot data. The complete pathway from unreacted to reacted state is discussed. Through chemical species analysis, we find that the C-N bond breaking is the main event of the shock-induced nitromethane decomposition.

Modification of short-range repulsive interactions in ReaxFF reactive force field for Fe–Ni–Al alloy

The short-range repulsive interactions of any force field must be modified to be applicable for high energy atomic collisions because of extremely far from equilibrium state when used in molecular dynamics(MD)simulations.In this work,the short-range repulsive interaction of a reactive force field(ReaxFF),describing Fe-Ni-Al alloy system,is well modified by adding a tabulated function form based on Ziegler-Biersack-Littmark(ZBL)potential.The modified interaction covers three ranges,including short range,smooth range,and primordial range.The short range is totally predominated by ZBL potential.The primordial range means the interactions in this range is the as-is ReaxFF with no changes.The smooth range links the short-range ZBL and primordial-range ReaxFF potentials with a taper function.Both energies and forces are guaranteed to be continuous,and qualified to the consistent requirement in LAMMPS.This modified force field is applicable for simulations of energetic particle bombardments and reproducing point defects'booming and recombination effectively.

Potentials of classical force fields for interactions between Na^+ and carbon nanotubes

Carbon nanotubes （CNTs） have long been expected to be excellent nanochannels for use in desalination membranes and other bio-inspired human-made channels owing to their experimentally confirmed ultrafast water flow and theoretically predicted ion rejection. The correct classical force field potential for the interactions between cations and CNTs plays a cru- cial role in understanding the transport behaviors of ions near and inside the CNT, which is key to these expectations. Here, using density functional theory calculations, we provide classical force field potentials for the interactions of Na＋/hydrated Na＋ with （7,7）, （8,8）, （9,9）, and （10,10）-type CNTs. These potentials can be directly used in current popular classical soft- ware such as nanoscale molecular dynamics （NAMD） by employing the tclBC interface. By incorporating the potential of hydrated cation-g interactions to classical all-atom force fields, we show that the ions will move inside the CNT and accu- mulate, which will block the water flow in wide CNTs. This blockage of water flow in wide CNTs is consistent with recent experimental observations. These results will be helpful for the understanding and design of desalination membranes, new types of nanofluidic channels, nanosensors, and nanoreactors based on CNT platforms.

Spatial information recognizing of ocean eddies based on virtual force field and its application

A new approach to detecting ocean eddies automatically from remote sensing imageries based on the ocean eddy＇s eigen-pattern in remote sensing imagery and ＂force field-based shape extracting method＂ is proposed. First, the analysis on extracting eddies＇ edges from remote sensing imagery using conventional edge detection arithmetic operators is performed and returns digitized vector edge data as a result. Second, attraction forces and fusion forces between edge curves were analyzed and calculated based on the vector eddy edges. Thirdly, the virtual significant spatial patterns of eddy were detected automatically using iterative repetition followed by optimized rule. Finally, the spatial form auto-detection of different types of ocean eddies was done using satellite images. The study verified that this is an effective way to identify and detect the ocean eddy with a complex form.

PARTITION RATIO AT SOLID-LIQUID INTERFACE UNDER INFLUENCE OF TEMPERATURE GRADIENT AND EXTERNAL FORCE FIELD

The relationship between the partition ratio at a solid-liquid interface and the temperature gradient or the external force field has been theoretically analysed.It is shown that under the influence of a temperature gradient or an external force field,the partition ratio at a solid-liquid interface will deviate from the equilibrium value.

MM CALCULATIONS OF THE METAL-PROTEIN MODEL COMPLEXES Ⅰ.MOLECULAR FORCE FIELD FOR COBALT COMPEXES

Cobalt-protein complexes play an important role in biochemical processes.The structure of the model molecule,Co(H_2O)_3SO_4(phen) has been studied by molecular mechanics.The molecular force field (MM2) parameters have been developed for the particular class of the complexes.

Ab Initio Calculations on Nitramide and Its Methyl Derivatives( Ⅱ )——The Harmonic Force Fields and Vibrational Spectra of Nitramide and Its Isotopic Derivatives

The harmonic force field and the vibrational spectrum of nitramide were calculated by using the ab initio gradient program TEXAS at the Hartree-Fock level with a 4-21G basis set. The directly computed theoretical harmonic force field was scaled by using empirical scale factors which are transferred from other molecules and provided an a priori prediction of fundamental frequencies and intensities. The average deviations between predicted vibrational frequencies of nitramide and experimental IR spectrum in an argon matrix are 63 cm-1 for symmetric vibrations and 41 cm-1 for antisymmetric modes. A new set of scale factors was optimized in this paper. These scale factors reduced the average deviations to 2. 3 cm-1 for symmetric modes and 0. 8 cm-1 for antisymmetric ones. The vibrational spectra of three isotopic derivatives of nitramide were predicted by using the force field resulted from the optimized set of scale factors, which are in good agreement with their experimental data in an argon matrix.

High-Fidelity Hugoniots of α Phase RDX Solid from High-Quality Force Field with Thermal,Zero-Point Vibration,and Anharmonic Effects

It is shown that the introduction of thermal effect, zero-point vibration, and phonon anharmonicity to a high quality and first-principle-Sased force field （atomic potential） results in a significant improvement in predict- ing the densities for the α phase crystalline hexahydro-1,3,5-trinitro-l,3,5-triazine （RDX）, and derivation of its high-fidelity Hugoniot locus and Mie-Grfineisen equation of state covering a very wide range of pressures and temperatures. This work can be used to efficiently and accurately predict the thermophysical properties of solid explosives over the pressures and temperatures to which they are subjected, which is a long-standing issue in the field of energetic materials.

Theoretical Model for Particle Behavio at Solidifying Front in Electromagnetic Force Field

The particle migrating behavior at solidifying front is discussed in theory with the application of electromagnetic force field(EMFF), on the basis of foregone analysis of force upon particle ahead of solidifying front without electromagnetic force field. The critical solidification velocities of particle pushing/engulfment transition ahead of horizontal and vertical interface are derived respectively when a certain EMFF is applied. And the critical electromagnetic forces of particle pushing/engulfment transition ahead of horizontal and vertical interface are also derived separately when a certain solidifying velocity is given.

Image Deformation by User-Defined Force Fields

We present a new method for image deformation. The warping technique provides smooth distortion with intuitive and easy manipulation. Driven by a restrained force field, the input image is deformed gradually and continuously. The method allows us to customize the force fields and the region of interest manually through some simple steps. Experimental results demonstrate the effectiveness and convenience of the approach.

The action mechanism of the work done by the electric field force on moving charges to stimulate the emergence of carrier generation/recombination in a PN junction

It is discovered that the product of the current and the electric field in a PN junction should be regarded as the rate of work(power)done by the electric field force on moving charges(hole current and electron current),which was previously misinterpreted as solely a Joule heating effect.We clarify that it is exactly the work done by the electric field force on the moving charges to stimulate the emergence of non-equilibrium carriers,which triggers the novel physical phenomena.As regards to Joule heat,we point out that it should be calculated from Ohm’s law,rather than simply from the product of the current and the electric field.Based on this understanding,we conduct thorough discussion on the role of the electric field force in the process of carrier recombination and carrier generation.The thermal effects of carrier recombination and carrier generation followed are incorporated into the thermal equation of energy.The present study shows that the exothermic effect of carrier recombination leads to a temperature rise at the PN interface,while the endothermic effect of carrier generation causes a temperature reduction at the interface.These two opposite effects cause opposite heat flow directions in the PN junction under forward and backward bias voltages,highlighting the significance of managing device heating phenomena in design considerations.Therefore,this study possesses referential significance for the design and tuning on the performance of piezotronic devices.

Vortex tuning of a submarine by Liutex force field model

Vortical structures of a submarine with appendages are fully turbulent and complex.Thus,flow control and vortex manipulation are of great importance for the hydrodynamic performance and acoustic characteristics.Take the generic submarine model DARPA Suboff as the test case,a vortex tuning method based on the Liutex force field is proposed to manipulate the vorticity field.Viscous flow past the submarine model in straight-line motion at a Reynolds number of 1.2×107 is achieved by solving the Reynolds averaged Navier-Stokes(RANS)equations.Multi-block structured mesh topology is used to discretize the computational domain,and the shear stress transport(SST)k-ωturbulence model is implemented to close the equations.The control of vortex is achieved by introducing additional source terms based on Liutex vortex definition and identification system to the RANS equations.The resistance acting on the submarine,flow field as well as the vortical structures are compared and analyzed.Results show that Liutex force model can effectively reduce the resistance by 9.31%and change the vortical structures apparently.

Impact of chosen force fields and applied load on thin film lubrication

The rapid development of molecular dynamics(MD)simulations,as well as classical and reactive atomic potentials,has enabled tribologists to gain new insights into lubrication performance at the fundamental level.However,the impact of adopted potentials on the rheological properties and tribological performance of hydrocarbons has not been researched adequately.This extensive study analyzed the effects of surface structure,applied load,and force field(FF)on the thin film lubrication of hexadecane.The lubricant film became more solid‐like as the applied load increased.In particular,with increasing applied load,there was an increase in the velocity slip,shear viscosity,and friction.The degree of ordering structure also changed with the applied load but rather insignificantly.It was also significantly dependent on the surface structure.The chosen FFs significantly influenced the lubrication performance,rheological properties,and molecular structure.The adaptive intermolecular reactive empirical bond order(AIREBO)potential resulted in more significant liquid‐like behaviors,and the smallest velocity slip,degree of ordering structure,and shear stress were compared using the optimized potential for liquid simulations of united atoms(OPLS‐UAs),condensed‐phase optimized molecular potential for atomic simulation studies(COMPASS),and ReaxFF.Generally,classical potentials,such as OPLS‐UA and COMPASS,exhibit more solid‐like behavior than reactive potentials do.Furthermore,owing to the solid‐like behavior,the lubricant temperatures obtained from OPLS‐UA and COMPASS were much lower than those obtained from AIREBO and ReaxFF.The increase in shear stress,as well as the decrease in velocity slip with an increase in the surface potential parameterζ,remained conserved for all chosen FFs,thus indicating that the proposed surface potential parameterζfor the COMPASS FF can be verified for a wide range of atomic models.

Validation of polarizable force field parameters for nucleic acids by inter-molecular interactions

Modeling structural and thermodynamic prop- erties of nucleic acids has long been a challenge in the development of force fields. Polarizable force fields are a new generation of potential functions to take charge redistribution and induced dipole into account, and have been proved to be reliable to model small molecules, polypeptides and proteins, but their use on nucleic acids is still rather limited. In this article, the interactions between nucleic acids and a small molecule or ion were modeled by AMOEBAbio09, a modern polarizable force field, and conventional non-polarizable AMBER99sb and CHARMM36 force fields. The resulting intermolecular interaction energies were compared with those calculated by ab initio quantum mechanics methods. Although the test is not sufficient to prove the reliability of the polarizable force field, the results at least validate its capability in modeling energetics of static configurations, which is one basic component in force field parameterization.

Numerical Analysis on Lorentz Force Field in φ160 Billet Under Harmonic Magnetic Field

A method was developed to quantitatively describe the time-varying Lorentz force field induced in billet under harmonic magnetic field.Based on the numerical analysis of the Lorentz force,the numeral features were concluded, which comprised of the mean force,the maximum force,the long axis,short axis and rotation direction of ellipse,the phase difference index,and so on.Using this method,the Lorentz force field in φ160 billet under 3 Hz harmonic magnetic field were analyzed.

Thermodynamic Insights of Base Flipping in TNA Duplex: Force Fields, Salt Concentrations, and Free-Energy Simulation Methods

Threofuranosyl nucleic acid(TNA)is an analogue of DNA with a shift in the internucleotide linkages from the wild-type 5’-to-3’direction to 3’-to-2.’This alteration leads to higher chemical stability,less reactive groups,and lower conformational flexibility.Experimental observations indicate that these characteristic changes are attributable to a minimal perturbation of the interaction network,but the thermodynamic stability of the duplex remains unaltered in the TNA mutation.We applied the equilibrium and nonequilibrium free-energy simulations employing three popular assisted model building with energy refinement(AMBER)force fields for nucleotides to investigate this mutation-dependent behavior in the base flipping from T(DNA)residue to the T-to-TFT mutation(TNA)computationally.The force fields were performed similarly,as described in the base-paired state.However,after exploring the high-energy regions with free-energy simulations,we observed that these three force fields behaved differently.Previous reports conclude that the net-neutral and excess-salt simulations provided similar results.Nonetheless,our free-energy simulation indicated that the presence of excess salt affected the thermodynamic stability.The free-energy barrier along the base-flipping pathway was generally elevated upon the addition of excess salts,but the relative height of the free-energy barriers in DNA and TNA duplexes did not change significantly.This phenomenon emphasizes the importance of adding sufficient salts in the simulation scheme to reproduce the experimental condition.

Studies on nitramide and its methyl derivatives with ab initio calculations——Ⅳ.The harmonic force field and vibrational spectra of dimethylnitramine and its isotopic derivatives

The complete harmonic vibrational force field of dimethylnitramine has been calculated at the Hartree-Fock level using 4-21G basis set.The harmonic force field was then scaled with scale factors previously derived from N-methylnitramine,and the vibrational spectrum of dimethylnitramine was computed.This apriori prediction,made with no reference to observations on dimethylnitramine, agrees with the experimental IR spectrum in gas phase with a mean deviation of 8.4 cm^(-1).Some of the scale factors were reoptimized by fitting of the computed force field to experimental data.The new set of scale factors reduced the mean deviation to 4.5 cm^(-1),and was used to predict the vibrational spectrum of deuterated form of dimethylnitramine(-6D).Dipole moment derivatives were also cal- culated and used to predict infrared intensities which are comparable with experimental values.

Ab initio study of vibrational frequencies, force field and normal coordinate analysis of the cis- and trans-isomers of nitrosomethanol

In this paper, the equilibrium geometries of two isomers of the newly found compound ——nitrosomethanol——have been optimized by ab initio SCF MO method with 3-21G basic set by gradient technique. And the second derivatives of potential energy (i. e. the force constant matrix elements) have been calculated analytically. Hence the entire force fields of the two isomers of nitro- somethanol have been obtained theoretically. The theoretical vibrational frequencies and the corres- ponding normal modes were obtained and compared with the experimental values. and the structures of two isomeric forms of nitrosomethanol are established.

FORCES ON CYLINDERS IN COEXISTING FIELD OF REGULAR WAVES AND CURRENTS

Based on the linear wave theory, the lift force on the cylinder under the action of both regular waves and currents related to inline force and the resultant force has been investigated. The relationship between the hydrodynamic coefficients of resultant force Cf, of drag force Cd, of inertia force Cm and of lift force CL and the redefined KC number is reported. It is indicated that in a certain region of KC number, the influence of lift force on the resultant force can not be ignored.