Flow Field Characteristics of Multi-Trophic Artificial Reef Based on Computation Fluid Dynamics

On the basis of computational fluid dynamics,the flow field characteristics of multi-trophic artificial reefs,including the flow field distribution features of a single reef under three different velocities and the effect of spacing between reefs on flow scale and the flow state,were analyzed.Results indicate upwelling,slow flow,and eddy around a single reef.Maximum velocity,height,and volume of upwelling in front of a single reef were positively correlated with inflow velocity.The length and volume of slow flow increased with the increase in inflow velocity.Eddies were present both inside and backward,and vorticity was positively correlated with inflow velocity.Space between reefs had a minor influence on the maximum velocity and height of upwelling.With the increase in space from 0.5 L to 1.5 L(L is the reef lehgth),the length of slow flow in the front and back of the combined reefs increased slightly.When the space was 2.0 L,the length of the slow flow decreased.In four different spaces,eddies were present inside and at the back of each reef.The maximum vorticity was negatively correlated with space from 0.5 L to 1.5 L,but under 2.0 L space,the maximum vorticity was close to the vorticity of a single reef under the same inflow velocity.

Manipulating the electron dynamics in the non-sequential double ionization process of Ar atoms by an orthogonal two-color laser field

Electron dynamics during non-sequential double ionization(NSDI) is one of the most attractive areas of research in the field of laser–atom or laser–molecule interaction. Based on the classic two-dimensional model, we study the process of NSDI of argon atoms driven by a few-cycle orthogonal two-color laser field composed of 800 nm and 400 nm laser pulses. By changing the relative phase of the two laser pulses, a localized enhancement of NSDI yield is observed at 0.5πand 1.5π, which could be attributed to a rapid and substantial increase in the number of electrons returning to the parent ion within extremely short time intervals at these specific phases. Through the analysis of the electron–electron momentum correlations within different time windows of NSDI events and the angular distributions of emitted electrons in different channels, we observe a more pronounced electron–electron correlation phenomenon in the recollision-induced ionization(RII) channel. This is attributed to the shorter delay time in the RII channel.

Charge self-consistent dynamical mean field theory calculations incombination with linear combination of numerical atomic orbitalsframework based density functional theory

We present a formalism of charge self-consistent dynamical mean field theory(DMFT)in combination with densityfunctional theory(DFT)within the linear combination of numerical atomic orbitals(LCNAO)framework.We implementedthe charge self-consistent DFT+DMFT formalism by interfacing a full-potential all-electron DFT code with threehybridization expansion-based continuous-time quantum Monte Carlo impurity solvers.The benchmarks on several 3d,4fand 5f strongly correlated electron systems validated our formalism and implementation.Furthermore,within the LCANOframework,our formalism is general and the code architecture is extensible,so it can work as a bridge merging differentLCNAO DFT packages and impurity solvers to do charge self-consistent DFT+DMFT calculations.

Dynamic Phase Field Description on Ductile Fracture Process

In comparison to discrete descriptions of fracture process,the recently proposed phase field methodology averts the numerical tracking strategy of discontinuities in solids,which enables the numerical implement simplification.An implicit finite element formulation based on the diffuse phase field is extended for stable and efficient analysis of complex dynamic fracture process in ductile solids.This exhibited formulation is shown to capture entire range of the characteristics of ductile material presenting J2-plasticity,embracing plasticization,cracks initiation,propagation,branching and merging while fulfilling the basic principle of thermodynamics.Herein,we implement a staggered time integration scheme of the dynamic elasto-plastic phase field method into the commercial finite element code.The numerical performance of the present advanced phase field model has been examined through several classic dynamic fracture benchmarks,and in all cases simulation results are in good agreement with the associated experimental data and other numerical results in previous literature.

Charging Properties and Particle Dynamics of Chang’e-5 Lunar Sample in an External Electric Field

Facing the challenges of in-situ utilization of lunar regolith resources,applying an external electric field to manipulate lunar particles has become a promising method for space particle control,which mainly depends on the particle charging properties in the applied electric field.Using the surficial lunar regolith samples brought back from the Moon by the Chang’e-5 mission(CE5 LS),this work successively studied their charging properties,particle dynamics,and their collision damages to aerospace materials under the action of an external electric field in high-vacuum conditions.The results indicated that the charging pro-cess and electrostatic projection of lunar regolith particles under high-vacuum conditions were different from those under atmosphere conditions.The particle diameter range of CE5 LS used in the experiment is 27.7-139.0 lm.For electric field strength of 3-12 kV·cm^(-1),the charge obtained by CE5 LS is 4.8×10^(-15)-4.7×10^(-13) C and the charge-to-mass ratio is 1.2×10^(-5)-6.8×10^(-4) C·kg^(-1).The CE5 LS is easier to be negatively charged in an external electric field.Furthermore,significant damages were observed on the target impact surfaces,indicating severe influences of lunar regolith particles on aerospace materials.Our work contributes to a more comprehensive understanding of physical mechanisms controlling the lunar regolith shielding and utilization,and will inspire broad efforts to develop the lunar in-situ engi-neering solutions.

Bioturbation Effects of Benthic Fish on the Phosphorus Dynamic in Overlying Water of Paddy Field

[Objective]The research aimed to investigate the bioturbation effects of benthic fish Misgurnus anguillicaudatus on phosphorus dynamic in overlying water of paddy field,as well as to explore the bioturbation mechanism.[Method]Based on simulation experiment,the phosphorus contents in overlying water were analyzed comparatively with and without Misgurnus anguillicaudatus by the using of ion chromatography and spectrophotometry.[Result]The concentrations of total phosphorus(TP),dissolved total phosphorus(DTP)and particular phosphorus(PP)in bioturbation group had no significant differences with those in control group in initial stage of experiment,and became significantly higher than control group in middle and late stages of experiment(P<0.05). The PP/TP ratios in bioturbation group were bigger than those in control group,the increase of TP concentration in bioturbation group was mainly due to the increase of PP. The ratios of dissolved inorganic phosphorus(DIP) to DTP (DIP/DTP) were significantly bigger than those in control group in middle and late stages of experiment(P<0.05).[Conclusion]The benthic fish had bioturbation effects on phosphorus in overlying water of paddy field,which increased the available phosphorus for rice growth.

Dynamic Field Division of Hydrocarbon Migration,Accumulation and Hydrocarbon Enrichment Rules in Sedimentary Basins

Hydrocarbon distribution rules in the deep and shallow parts of sedimentary basins are considerably different, particularly in the following four aspects. First, the critical porosity for hydrocarbon migration is much lower in the deep parts of basins： at a depth of 7000 m, hydrocarbons can accumulate only in rocks with porosity less than 5%. However, in the shallow parts of basins （i.e., depths of around 1000 m）, hydrocarbon can accumulate in rocks only when porosity is over 20%. Second, hydrocarbon reservoirs tend to exhibit negative pressures after hydrocarbon accumulation at depth, with a pressure coefficient less than 0.7. However, hydrocarbon reservoirs at shallow depths tend to exhibit high pressure after hydrocarbon accumulation. Third, deep reservoirs tend to exhibit characteristics of oil （-gas）-water inversion, indicating that the oil （gas） accumulated under the water. However, the oil （gas） tends to accumulate over water in shallow reservoirs. Fourth, continuous unconventional tight hydrocarbon reservoirs are distributed widely in deep reservoirs, where the buoyancy force is not the primary dynamic force and the caprock is not involved during the process of hydrocarbon accumulation. Conversely, the majority of hydrocarbons in shallow regions accumulate in traps with complex structures. The results of this study indicate that two dynamic boundary conditions are primarily responsible for the above phenomena： a lower limit to the buoyancy force and the lower limit of hydrocarbon accumulation overall, corresponding to about 10%-12% porosity and irreducible water saturation of 100%, respectively. These two dynamic boundary conditions were used to divide sedimentary basins into three different dynamic fields of hydrocarbon accumulation： the free fluid dynamic field, limit fluid dynamic field, and restrain fluid dynamic field. The free fluid dynamic field is located between the surface and the lower limit of the buoyancy force, such that hydrocarbons in this field migrate and accumulate under the influence of, for example, the buoyancy force, pressure, hydrodynamic force, and capillary force. The hydrocarbon reservoirs formed are characterized as ＂four high,＂ indicating that they accumulate in high structures, are sealed in high locations, migrate into areas of high porosity, and are stored in reservoirs at high pressure. The basic features of distribution and accumulation in this case include hydrocarbon migration as a result of the buoyancy force and formation of a reservoir by a caprock. The limit fluid dynamic field is located between the lower limit of the buoyancy force and the lower limit of hydrocarbon accumulation overall; the hydrocarbon migrates and accumulates as a result of, for example, the molecular expansion force and the capillary force. The hydrocarbon reservoirs formed are characterized as ＂four low,＂ indicating that hydrocarbons accumulate in low structures, migrate into areas of low porosity, and accumulate in reservoirs with low pressure, and that oil（-gas）-water inversion occurs at low locations. Continuous hydrocarbon accumulation over a large area is a basic feature of this field. The restrain fluid dynamic field is located under the bottom of hydrocarbon accumulation, such that the entire pore space is filled with water. Hydrocarbons migrate as a result of the molecular diffusion force only. This field lacks many of the basic conditions required for formation of hydrocarbon reservoirs： there is no effective porosity, movable fluid, or hydrocarbon accumulation, and potential for hydrocarbon exploration is low. Many conventional hydrocarbon resources have been discovered and exploited in the free fluid dynamic field of shallow reservoirs, where exploration potential was previously considered to be low. Continuous unconventional tight hydrocarbon resources have been discovered in the limit fluid dynamic field of deep reservoirs; the exploration potential of this setting is thought to be tremendous, indicating that future exploration should be focused primarily in this direction.

Simulation and Analysis on the Two-Phase Flow Fields in a Rotating-Stream-Tray Absorber by Using Computational Fluid Dynamics

The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, and a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall; the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.

Effect of an external electric field on liquid water using molecular dynamics simulation with a flexible potential

Molecular dynamics simulations of liquid water were performed at 258 K and density of 1.0 g/cm^3 under different strengths of an external electric field, ranging from 0 to 8.0×10^9V/m, to investigate the influence of an external field on structural and dynamic properties of water. The flexible simple point charge model is used for water molecules. An enhancement of the water hydrogen bond structure with increasing strength of the electric field has been deduced from the radial distribution functions and the analysis of hydrogen bond structure. With increasing field strength, water system has a more perfect structure, which is shnilar to ice structure. However, the electrofreezing phenomenon of liquid water has not been detected because of a too large self-diffusion coefficient. The self-diffusion coefficient decreases remarkably with increasing strength of electric field, and the self-diffusion coefficient is anisotropic.

Evaluation of dynamic characteristics of silt in Yellow River Flood Field after freeze-thaw cycles

Frothing is a main disease of highways in Yellow River Flood Field, due to the loss of dynamic strength of roadbed soils under the couple effects of temperature, salt, and vehicle traffic load. This is strongly linked to the dynamic characteristics of silt in this region. To analyze these couple effects on the dynamic characteristics of silt, a series of tests(i.e., freeze-thaw cycling tests, vibration triaxial tests and ultrasonic wave velocity tests) were conducted and two kinds of silt(i.e., salt-free and 3%-salt silt) were designed. The results indicate that the dynamic shear strength and dynamic modulus decrease with increasing freeze-thaw cycles, while the damping ratio simultaneously increases. Furthermore, compared to salt-free silt, the decrement of dynamic shear strength and dynamic modulus of silt with 3% salt is more significant, but the damping ratio of 3%-salt silt is larger. In ultrasonic wave velocity tests, ultrasonic wave velocity of frozen soil specimens decreases as the number of freeze-thaw cycles increases. Based on the results of ultrasonic wave velocity tests, a preliminary model is proposed to evaluate damage of silt through field measurement ultrasonic data. The study could provide a theoretical basis for the treatment of silty soil highway.

Size dependency and potential field influence on deriving mechanical properties of carbon nanotubes using molecular dynamics

A thorough understanding on the mechanical properties of carbon nanotube （CNT） is essential in extending the advanced applications of CNT based systems. However, conducting experiments to estimate mechanical properties at this scale is extremely challenging. Therefore, development of mechanistic models to estimate the mechanical properties of CNTs along with the integration of existing continuum mechanics concepts is critically important. This paper presents a comprehensive molecular dynamics simulation study on the size dependency and potential function influence of mechanical properties of CNT. Commonly used reactive bond order （REBO） and adaptive intermolecular reactive bond order （A1REBO） potential functions were considered in this regard. Young＇s modulus and shear modulus of CNTs are derived by integrating classical continuum mechanics concepts with molecular dynamics simulations. The results indicate that the potential function has a significant influence on the estimated mechanical properties of CNTs, and the influence of potential field is much higher when studying the torsional behaviour of CNTs than the tensile behaviour.

Perfect elastic-viscoplastic field at mode I dynamic propagating crack-tip

The viscosity of material is considered at propagating crack-tip. Under the assumption that the artificial viscosity coefficient is in inverse proportion to power law of the plastic strain rate, an elastic-viscoplastic asymptotic analysis is carried out for moving crack-tip fields in power-hardening materials under plane-strain condition. A continuous solution is obtained containing no discontinuities. The variations of numerical solution are discussed for mode Ⅰ crack according to each parameter. It is shown that stress and strain both possess exponential singularity. The elasticity, plasticity and viscosity of material at crack-tip only can be matched reasonably under linear-hardening condition. And the tip field contains no elastic unloading zone for mode I crack. It approaches the limiting case, crack-tip is under ultra-viscose situation and energy accumulates, crack-tip begins to propagate under different compression situations.

Dynamic collision avoidance for cooperative fixed-wing UAV swarm based on normalized artificial potential field optimization

Cooperative path planning is an important area in fixed-wing UAV swarm.However,avoiding multiple timevarying obstacles and avoiding local optimum are two challenges for existing approaches in a dynamic environment.Firstly,a normalized artificial potential field optimization is proposed by reconstructing a novel function with anisotropy in each dimension,which can make the flight speed of a fixed UAV swarm independent of the repulsive/attractive gain coefficient and avoid trapping into local optimization and local oscillation.Then,taking into account minimum velocity and turning angular velocity of fixed-wing UAV swarm,a strategy of decomposing target vector to avoid moving obstacles and pop-up threats is proposed.Finally,several simulations are carried out to illustrate superiority and effectiveness.

Machine learning molecular dynamics simulations of liquid methanol

As the simplest hydrogen-bonded alcohol,liquid methanol has attracted intensive experimental and theoretical interest.However,theoretical investigations on this system have primarily relied on empirical intermolecular force fields or ab initio molecular dynamics with semilocal density functionals.Inspired by recent studies on bulk water using increasingly accurate machine learning force fields,we report a new machine learning force field for liquid methanol with a hybrid functional revPBE0 plus dispersion correction.Molecular dynamics simulations on this machine learning force field are orders of magnitude faster than ab initio molecular dynamics simulations,yielding the radial distribution functions,selfdiffusion coefficients,and hydrogen bond network properties with very small statistical errors.The resulting structural and dynamical properties are compared well with the experimental data,demonstrating the superior accuracy of this machine learning force field.This work represents a successful step toward a first-principles description of this benchmark system and showcases the general applicability of the machine learning force field in studying liquid systems.

Structure-preserving algorithms for guiding center dynamics based on the slow manifold of classical Pauli particle

The classical Pauli particle(CPP) serves as a slow manifold, substituting the conventional guiding center dynamics. Based on the CPP, we utilize the averaged vector field(AVF) method in the computations of drift orbits. Demonstrating significantly higher efficiency, this advanced method is capable of accomplishing the simulation in less than one-third of the time of directly computing the guiding center motion. In contrast to the CPP-based Boris algorithm, this approach inherits the advantages of the AVF method, yielding stable trajectories even achieved with a tenfold time step and reducing the energy error by two orders of magnitude. By comparing these two CPP algorithms with the traditional RK4 method, the numerical results indicate a remarkable performance in terms of both the computational efficiency and error elimination. Moreover, we verify the properties of slow manifold integrators and successfully observe the bounce on both sides of the limiting slow manifold with deliberately chosen perturbed initial conditions. To evaluate the practical value of the methods, we conduct simulations in non-axisymmetric perturbation magnetic fields as part of the experiments,demonstrating that our CPP-based AVF method can handle simulations under complex magnetic field configurations with high accuracy, which the CPP-based Boris algorithm lacks. Through numerical experiments, we demonstrate that the CPP can replace guiding center dynamics in using energy-preserving algorithms for computations, providing a new, efficient, as well as stable approach for applying structure-preserving algorithms in plasma simulations.

Simplified Dynamic Models of Grass Field Ecosystem

Some simplified dynamic models of grass field ecosystem are developed and investigated. The maximum simplified one consists of two variables, living grass biomass and soil wetness. The analyses of such models show that there exists only desert regime without grasses if the precipitation p is less than a critical value pc; the grass biomass continuously depends on p if the interaction between grass biomass and the soil wetness is weak, but the strong interaction results in the bifurcation of grass biomass in the vicinity of pc: the grass biomass is rich as p > pc, but it becomes desertification as p<pc. Periodic solutions also exist in the model, if the seasonal cycle of model's parameters is introduced. An improved model consists of three variables, i.e. the living grass biomass x, the nonliving grass biomass accumulated on the ground surface y and the soil wetness z. The behaviours of such three variables model are more complicated. The initial values of y and z play a very important role.

Half-metallic ferromagneticWeyl fermions related to dynamic correlations in the zinc-blende compound VAs

The realization of 100%polarized topologicalWeyl fermions in half-metallic ferromagnets is of particular importance for fundamental research and spintronic applications.Here,we theoretically investigate the electronic and topological properties of the zinc-blende compound VAs,which was deemed as a half-metallic ferromagnet related to dynamic correlations.Based on the combination of density functional theory and dynamical mean field theory,we uncover that the half-metallic ferromagnet VAs exhibits attractive Weyl semimetallic behaviors which are very close to the Fermi level in the DFT+U regime with effect U values ranging from 1.5 eV to 2.5 eV.Meanwhile,we also investigate the magnetization-dependent topological properties;the results show that the change of magnetization directions only slightly affects the positions of Weyl points,which is attributed to the weak spin–orbital coupling effects.The topological surface states of VAs projected on semi-infinite(001)and(111)surfaces are investigated.The Fermi arcs of all Weyl points are clearly visible on the projected Fermi surfaces.Our findings suggest that VAs is a fully spin-polarized Weyl semimetal with many-body correlated effects in the effective U values range from 1.5 eV to 2.5 eV.

Recent tectonic stress field zoning in Sichuan-Yunnan region and its dynamic interest

In this paper, we have carefully determined the stress zones in the Sichuan-Yunnan region with reference to the in-situ stress data of hydraulic fracturing and the inverted fault slip data by using the step-by-step convergence method for stress zoning based on focal mechanism solutions. The results indicate that the tectonic stress field in the Sichuan-Yunnan region is divided into 3 stress zones by 2 approximately parallel NNW-trending stress transition belts. The area between the 2 belts is the Sichuan-Yunnan stress zone where the maximum principal stress σ1 is just in the NNW direction. The eastern boundary of Sichuan-Yunnan stress zone （the eastern stress transition belt） is basically consistent with the eastern boundary of Sichuan-Yunnan rhombic block. The western boundary of Sichuan-Yunnan stress zone （the western stress transition belt） is not totally consistent with the western boundary of Sichuan-Yunnan rhombic block. The northern segment of the western stress transition belt extends basically along the Jinshajiang fault and accords with the western boundary of Sichuan-Yunnan rhombic block, while its southern segment does not extend along the southwestern boundary of the rhombic block, i.e., Honghe fault and converge with the eastern stress transition belt, but stretches continuously in the NNW direction and accords with the Yingpanshan fault. We therefore consider that under the combined influence from the northward motion of India Plate, the southeastward shift of east Qinghai-Xizang Plateau and the strong obstruction of South China block, the tectonic stress field in the Sichuan-Yunnan region might not be totally controlled by the previous tectonic frame and new stress transition belt may have possibly formed.

A simplified dynamic method for field capacity estimation and its parameter analysis

This paper presents a simplified dynamic method based on the definition of field capacity. Two soil hydraulic characteristics models, the Brooks-Corey （BC） model and the van Genuchten （vG） model, and four soil data groups were used in this study. The relative drainage rate, which is a unique parameter and independent of the soil type in the simplified dynamic method, was analyzed using the pressure-based method with a matric potential of -1/3 bar and the flux-based method with a drainage flux of 0.005 cm/d. As a result, the relative drainage rate of the simplified dynamic method was determined to be 3% per day. This was verified by the similar field capacity results estimated with the three methods for most soils suitable for cultivating plants. In addition, the drainage time calculated with the simplified dynamic method was two to three days, which agrees with the classical definition of field capacity. We recommend the simplified dynamic method with a relative drainage rate of 3% per day due to its simple application and clearly physically-based concept.

Effect of dynamic loading orientation on fracture properties of surrounding rocks in twin tunnels

For expedited transportation,vehicular tunnels are often designed as two adjacent tunnels,which frequently experience dynamic stress waves from various orientations during blasting excavation.To analyze the impact of dynamic loading orientation on the stability of the twin-tunnel,a split Hopkinson pressure bar(SHPB)apparatus was used to conduct a dynamic test on the twin-tunnel specimens.The two tunnels were rotated around the specimen’s center to consider the effect of dynamic loading orientation.LS-DYNA software was used for numerical simulation to reveal the failure properties and stress wave propagation law of the twin-tunnel specimens.The findings indicate that,for a twin-tunnel exposed to a dynamic load from different orientations,the crack initiation position appears most often at the tunnel corner,tunnel spandrel,and tunnel floor.As the impact direction is created by a certain angle(30°,45°,60°,120°,135°,and 150°),the fractures are produced in the middle of the line between the left tunnel corner and the right tunnel spandrel.As the impact loading angle(a)is 90°,the tunnel sustains minimal damage,and only tensile fractures form in the surrounding rocks.The orientation of the impact load could change the stress distribution in the twin-tunnel,and major fractures are more likely to form in areas where the tensile stress is concentrated.