Loading...

Archive

    25 May 2024, Volume 41 Issue 3 Previous Issue    Next Issue

    For Selected: Toggle Thumbnails
    A Boundary Variation Diminishing Reconstructed Solver for Mie-Grüneisen Mixture Model
    Zongduo WU, Jin YAN, Jianhua PANG, Yifang SUN, Qingyang MAN
    2024, 41(3): 277-286.  DOI: 10.19596/j.cnki.1001-246x.8712
    Abstract ( )   HTML ( )   PDF (2168KB) ( )  

    A numerical reconstruction is implemented according to a combination of MUSCL(Monotonic Upstream-Centered Scheme for Conservation Laws) and THINC(Tangent of Hyperbola for Interface Capturing) methods. And a new constructed solver is derived by the combination. The main principle of this reconstruction solver is to keep the variation diminishing of every cell boundary to a low level. And a BVD(Boundary Variation Diminishing) principle is set according to the left and right boundary states in MUSCL method, as well as which in THINC methods. To obey the BVD principle, an alternative choice between the MUSCL and THINC is needed here. Thanks to the BVD solver, the discontinuous section becomes smooth and the oscillation is significantly controlled. The numerical performance of BVD reconstruction is then tested by 1D and 2D numerical examples.

    Finite Deformation Theory of Poroviscoelasticity Based on Logarithmic Strain
    Xiong TANG, Pei ZHENG
    2024, 41(3): 287-297.  DOI: 10.19596/j.cnki.1001-246x.8723
    Abstract ( )   HTML ( )   PDF (3120KB) ( )  

    In the framework of finite deformation theory, a theoretical model of poroviscoelasticity is proposed which is based on logarithmic strain and Kelvin rheological model. The model is obtained by assuming a linear relationship between Kirchhoff stress and pore pressure and logarithmic strain and the variation of Lagrangian porosity, and then directly replacing the infinitesimal strain in the linear pore viscoelastic model with the logarithmic strain. As a verification, the theoretical model is used to study the classic Terzaghi's one-dimensional consolidation problem. By comparing with the numerical results of the poroelastic finite deformation model, the results show that the viscoelastic response and elastic response curves of the pore solid skeleton are almost identical in the early stage of consolidation, but with the passage of time, the viscous response of the pore solid skeleton gradually dominates the deformation of the skeleton and affects the final result of consolidation. In addition, the viscous response of the skeleton delays the diffusion of pore pressure. In addition, by setting the viscosity contribution coefficient ζ=0.001, the poroviscoelastic response is numerically "degraded" to the poroelastic response, which verifies the correctness of the model to a certain extent.

    Simulation of Droplet Impact on Cable Surfaces with Different Contact Angles
    Jiru HUAI, Peng WANG, Mengyu YANG
    2024, 41(3): 298-307.  DOI: 10.19596/j.cnki.1001-246x.8717
    Abstract ( )   HTML ( )   PDF (13751KB) ( )  

    Icing of transmission lines will have a negative impact on the safety of line operation and human production and life. The research on anti-icing has gradually developed from the early active anti-icing to the present passive anti-icing, but the passive anti-icing research mostly stays on the improvement of hydrophobicity without in-depth exploration of its anti-icing principle. For common droplet transmission cables, the corresponding calculation model is built. The relationship between the contact time, spreading radius and contact angle between droplet and cable is discussed by using the built-in VOF of Fluent to calculate different hydrophobic surfaces, and the pressure nephogram and velocity vector diagram of droplet under different contact angles are analyzed. The simulation results show that the motion state of the droplet contact with the wall is affected by the wall contact angle, and the spread area and contact time of the droplet are inversely proportional to the wall contact angle. Increasing the wall contact angle will reduce the spread area, reduce the contact time, accelerate the droplet retraction and rebound, and the moving droplet will eventually slide off the wall. The spread area of the wall surface with small contact angle increases and the contact time increases. The droplet finally spreads completely and stays on the wall surface. These results have certain guidance and reference significance for exploring the principle and characteristics of droplet icing on cable wall under different hydrophobicity.

    An Approximate Shock Wave Formula for Real Gases
    Shesheng XUE, Xirui ZHU
    2024, 41(3): 308-315.  DOI: 10.19596/j.cnki.1001-246x.8718
    Abstract ( )   HTML ( )   PDF (1080KB) ( )  

    On the basis of the expanded multi-dimensional Virial equation, the gas equation of state (EOS) is expanded the series to the second term, so that the density effect of real gas can be considered. Using the basic formula of the shock wave and small parameter contained in the Virial equation, the explicit expressions of the gas density, pressure and velocity behind the shock wave changing with the shock wave velocity and the shock wave Mach number are obtained by the perturbation method. The results show that for the same Mach number of the shock waves, compared with the results by the ideal gas model, the pressure, velocity and density behind the waves are all lower, especially the density. With the increase of the Mach number and the density of the gas, the difference in results will be even greater. These relations are the proper corrections to the ideal gas, which reflects the influence of the volume and repulsion effects of the gas molecule, and are fit for the gases whose density is lower than 100 kg·m-3. The formula of the shock wave for the ideal gas can be regarded as its zero-order approximation. It is very convenient to use these relations to analyze the properties of shock waves.

    Lattice Boltzmann Model for Simulating Heat-fluid-solid Interaction in Wellbore
    Chunyu HU
    2024, 41(3): 316-324.  DOI: 10.19596/j.cnki.1001-246x.8715
    Abstract ( )   HTML ( )   PDF (4161KB) ( )  

    A model based on lattice Boltzmann method (LBM) is developed to solve the heat transfer problems in the wellbore, which can simultaneously solve the controlling equations of temperature-pressure coupled flow of fluid in the wellbore, forced thermal convection in fluid flow and fluid-solid heat exchange between the wellbore and the formation, so it can realize the coupling solution of fluid flow field, fluid temperature field and wellbore temperature field, and compared with the traditional model, it not only overcomes the defect that the velocity is constant in wellbore, but also has a wider scope of application and higher computational efficiency. The reliability and accuracy of the model have been verified by comparative analysis with previous studies. The results show that the change of fluid velocity in the wellbore affects the temperature distribution of the fluid in the wellbore. Under the production condition, the temperature distribution at the center of the fluid is high and the temperature at the boundary of the fluid is low at the same depth. The downward trend of fluid temperature along the shaft axis will go through three stages: slow, steady and slow, and these three stages are affected by flow Reynolds number or fluid Prandtl number.

    Study on Influence of Displacement Pressure Gradient on Water Flooding Oil Recovery
    Wenfu CUI
    2024, 41(3): 325-333.  DOI: 10.19596/j.cnki.1001-246x.8726
    Abstract ( )   HTML ( )   PDF (12535KB) ( )  

    For sandstone reservoirs with extremely high water content, the Volume of Fluid (VOF) method is able to trace the dynamic changes of multiphase interface and reproduce the physical process of microscopic seepage. The influence of different displacement pressure gradients on the microscopic occurrence characteristics of remaining oil and water flooding recovery is studied. The influence of displacement pressure gradient on water flooding recovery is revealed by analyzing the seepage characteristics and remaining oil occurrence state in microscopic pore structure. The increase of displacement pressure gradient can break through capillary force of small throat, form a new seepage channel, and increase the water flooding recovery. The viscosity and wettability of crude oil affect the recovery degree under different pressure gradients.

    Study on Radiation Characteristics of Cylindrical Hydrogen and Argon Plasma Antennas under Inhomogeneous Magnetic Field
    Guangzu PAN, Dan DU, Hua ZHOU, Kaijian YANG, Guanjin QIAO, Shaoxiong HU, Weibo YAO, Xueyu GONG
    2024, 41(3): 334-344.  DOI: 10.19596/j.cnki.1001-246x.8722
    Abstract ( )   HTML ( )   PDF (12181KB) ( )  

    In this paper, two-dimensional axisymmetric model of hydrogen and argon plasma antenna is established under non-uniform magnetic field conditions, and the effects of magnetic field distribution on electron density distribution, radiation direction map and gain performance of monopole columnar hydrogen and argon plasma antenna are compared and analyzed under low pressure conditions based on COMSOL calculation results. The results show that under certain conditions: 1) The electron distribution of argon plasma antenna is more affected by the magnetic field than that of hydrogen plasma antenna. The electron density distribution of argon plasma antenna is inhomogeneous, while that of hydrogen plasma antenna is relatively uniform. 2) The radiation direction of the hydrogen and argon plasma antenna can be controlled by adjusting the magnetic field and wave frequency. At some specific wave frequencies, the hydrogen and argon plasma antenna have good directivity with small side lobe. Moreover, gain of hydrogen and argon plasma antenna can be improved by changing the magnetic field distribution while maintaining their directivity. 3) The influence of the magnetic field distribution on the gain of the argon plasma antenna is greater than that of the hydrogen plasma antenna. When the magnetic field is mainly near the anode, the gain of the hydrogen plasma antenna is greater than that of the argon plasma antenna. As the magnetic field strength near the cathode increases gradually, the gain of argon plasma antenna is much larger than that of hydrogen plasma antenna.

    Molecular Dynamics Simulation of Physical Properties of Silicon Modified Phenolic Resin
    Bili XU, Zhao JING, Xiao LIU, Bo DAI, Guangfu JI, Kuibao ZHANG, Nina GE
    2024, 41(3): 345-356.  DOI: 10.19596/j.cnki.1001-246x.8710
    Abstract ( )   HTML ( )   PDF (13497KB) ( )  

    The physical properties of modified nano-SiO2 and methyl-phenyl-dimethoxy-silane modified phenolic resin are studied by molecular dynamics simulation. The results show that the glass transition temperature of unmodified phenolic resin at 300 K is 362 K, the elastic modulus and shear modulus are 5.45 GPa and 2.19 GPa, the thermal conductivity and thermal expansion coefficients are 0.37 W·(m·k)-1 and 3.8×10-5K-1, respectively. The addition of nano-SiO2 increases the glass transition temperature by 1.6%, the elastic modulus and shear modulus by 34.9% and 28.8%, and the thermal conductivity and thermal expansion by 11% and 31.6%, respectively. The thermal conductivity and thermal expansion are reduced by 11% and 31.6%, respectively. SiO2 surface grafting 3%, 5%, 7% and 10% silane coupling agent and methyl-phenyl- dimethoxy-silane modified phenolic resin, the glass transition temperature increased by 10.5%, 15.2%, 16.8%, 19.3% and 1.5% respectively, the elastic modulus increased by 44.4%, 53.2%, 53.8%, 63.5% and 13.4% respectively, and the thermal conductivity decreased by 12.4%, 13.5%, 11.2%, 7% and 10% respectively. Moreover, the thermal expansion coefficient of phenol formaldehyde resin modified by methyl phenyl dimethoxy silane increased by 51.8% compared with the unmodified phenol formaldehyde resin. The study show that the doping of nano-SiO2, the grafting of silane coupling agent on the SiO2 surface, and the modification of methyl-phenyl-dimethoxy-silane can improve the glass transition temperature, and mechanical properties and reduce the thermal conductivity of phenolic resin. Only nano-SiO2 doping can reduce the thermal expansion coefficient, whereas the modification of methyl-phenyl-dimethoxy-silane will increase substantially.

    Application of Genetic Algorithm to Optimal Design of Shielding Materials for Neutron-γ Mixed Radiation Fields
    Wenmin HAN, Yaodong DAI, Chuqing YAO, Jiaxiang TIAN, Danfeng JIANG, Yifan ZHOU
    2024, 41(3): 357-366.  DOI: 10.19596/j.cnki.1001-246x.8714
    Abstract ( )   HTML ( )   PDF (8011KB) ( )  

    Based on the neutron-γ mixed radiation field, the metal oxide filler components in the material are optimized, and the comprehensive shielding performance of WO3/Bi2O3/Gd2O3/B4C mixed filler against low energy neutrons and different energy γ rays is obtained by Monte Carlo simulation. The optimal ratio of filler components is found by using genetic algorithm and neural network. Through the calculation and optimization of the total dose equivalent, it is found that the optimal ratio is different under different radiation environments. And the comprehensive shielding performance can be optimized by using Bi2O3 and B4C (9:1) mixed fillers when the neutron (thermal neutron Maxwell distribution spectrum) flux is equal to γ ray (0.5-3 MeV) flux and the total mass of the shielding filler is constant. The results of Monte Carlo program show that the error is within an acceptable range, which indicates that the optimal design of the shielding filler is feasible. It can save a lot of calculation time and provide a theoretical basis for the design and preparation of shielding materials.

    Calculation Method of 3D Transient Temperature Fields in Functionally Graded Materials Subjected to Annular Heat Source
    Xuefeng SUO, Denghui HE, Huadong WANG, Wei CAO
    2024, 41(3): 367-379.  DOI: 10.19596/j.cnki.1001-246x.8764
    Abstract ( )   HTML ( )   PDF (15289KB) ( )  

    In order to solve the complex and time-consuming problem of 3D transient heat conduction analysis of functionally gradient materials subjected to annular heat source, a semi-analytical method for the heat conduction of exponentially graded materials with an annular heat source is developed through Fourier-Laplace frequency response functions, numerical Laplace inversion and two-dimensional discrete Fourier inversion, which is compared with 3D Fourier transform method in literature and FEM to verify the reliability of the proposed method. Based on the semi-analytical method, the temperature and temperature gradient response, as well as the sensitivity of their peaks under annular and Gaussian heat sources are compared. Case studies show that the temperature and temperature gradient distribution have different responses to the inner and outer diameter of the annular heat source and the material gradient indices. Compared with Gaussian heat sources, annular heat sources have the relatively uniform temperature and temperature gradient fields. Changing the inner and outer diameter of the annular heat source can adjust the temperature peak and temperature gradient peak.

    P-SV Wave Prestack Inversion Based on Hybrid Algorithm
    Qin LI, Min ZHANG, Ying XU
    2024, 41(3): 380-391.  DOI: 10.19596/j.cnki.1001-246x.8730
    Abstract ( )   HTML ( )   PDF (9784KB) ( )  

    The simulated annealing method is introduced into genetic algorithm to form a hybrid algorithm, which can overcome the defects of genetic algorithm in optimization and avoid the premature phenomenon. According to the approximate formula of reflection coefficient of P-SV wave, the objective function of anisotropy parameter inversion is established. The anisotropy parameter inversion of TTI medium is achieved by using single algorithm and mixed algorithm respectively. The results show that the mixed algorithm inversion accuracy is higher and the process is more stable. Furthermore, Gaussian white noise with a signal-to-noise ratio of 5 is added to the data of the theoretical model to test the anti-noise property. The error of the inversion results is small and the anti-noise ability is good. Finally, the anisotropic parameter inversion of the modified Hess model is carried out, and the inversion profile is obtained, which is in good agreement with the original profile, and the effectiveness of the hybrid algorithm is verified. In addition, the research area is selected to perform inversion test on the actual data. The inversion results are good and the accuracy is high, which verifies the effectiveness of the hybrid algorithm in practical application.

    A Novel Random Walk Algorithm for Optimal Configuration of Micro-grid
    Lei PAN, Guomin CUI, Ruifang ZHANG, Hongbin LIU, Yuan XIAO, Zhikang YI
    2024, 41(3): 392-402.  DOI: 10.19596/j.cnki.1001-246x.8703
    Abstract ( )   HTML ( )   PDF (3818KB) ( )  

    In order to solve the configuration optimization problem of isolated micro-grid, wind driven generator, photovoltaic, diesel generator and energy storage battery optimization model described in the form of energy flow matching is established, which can flexibly form the node connection relation representing the output of equipment at each time. At the same time, in view of the precocious convergence of swarm intelligence algorithm applied to optimal configuration of micro-grid, a random walk optimization algorithm suitable for optimal configuration of micro-grid is proposed. Guided by reducing the annual comprehensive cost of the system, the algorithm realizes synchronous optimization of continuous variable (equipment output) and integer variable (equipment quantity) by randomly increasing or decreasing the hourly output of equipment. By accepting the differential solution mechanism, the algorithm has the ability to jump out of the local optimal solution and better take into account the global search and local search in the capacity optimization process of micro-grid. Applying the random walk algorithm to the simulation example, the annual comprehensive cost is 552 826.39 yuan. Compared with particle swarm optimization algorithm, a better result is obtained. The superiority of the algorithm in optimization accuracy is verified.

Share: