Top Downloaded

    Published in last 1 year| In last 2 years| In last 3 years| All| Most Downloaded in Recent Month | Most Downloaded in Recent Year|

    Published in last 1 year
    Please wait a minute...
    For Selected: Toggle Thumbnails
    Effects of X Direction KSEA Interaction on Geometric Quantum Discord
    Batur AYGUL, Jinfeng ZHANG, Hamutjan AKBAR, Abliz AHMAD
    Chinese Journal of Computational Physics    2024, 41 (5): 663-669.   DOI: 10.19596/j.cnki.1001-246x.8798
    Abstract180)   HTML1)    PDF (3549KB)(598)      

    Based on NMQSD (non-Markovian quantum state diffusion) method, this paper studies containing dynamical evolution of geometric quantum discord in a system with X direction Kaplan-Shekhtman-Entin-Wohlman-Aharony (KSEA) interaction in a hybrid non-Markovian/Markovian environment and a single non-Markovian/Markovian environment. The results show that the geometric quantum discord in the hybrid non-Markovian environment is higher than that in the single non-Markovian environment. When interaction with X direction KSEA in the same initial state and different directions of intensity, the geometric quantum discord of the system achieves better recovery in a hybrid non-Markov environment. And we find that the single bath in the Markovian environment is better than the hybrid bath.

    Table and Figures | Reference | Related Articles | Metrics
    Liquid Heavy Metal Reactor Fuel Rod and Control Rod Design
    Changheng XU, Hui PAN, Mingtao HE, Changyou ZHAO, Dechang CAI, Huaijin XU
    Chinese Journal of Computational Physics    2024, 41 (5): 582-588.   DOI: 10.19596/j.cnki.1001-246x.8863
    Abstract179)   HTML2)    PDF (5070KB)(567)      

    The aim of the paper is to analyse the physical properties of liquid heavy metal fast reactor fuel rods and control rods from a neutronics point of view to support their structural design and material selection.In this paper, the liquid-uranium ratio selection, radial power distribution, cladding material selection, reflective layer material selection for fuel rods, radial structural design of control rods, absorber material selection, and moderated structural design are computationally analysed using the fast neutron reactor analysis program SARAX. Calculation results show that: the Keff of the fuel rod grid element is linearly related to the liquid-uranium ratio, and the selection of the liquid-uranium ratio needs to consider the heat-carrying factor more. The fuel rod cladding is preferred to be made of ferritic/martensitic steel materials. The materials for the reflective layer at the axial ends of the fuel may be chosen to be made of steel materials.Based on the results of the analysis of the Pu migration phenomenon, none of the current physical calculation procedures are coupled with the Pu migration phenomenon, and the power distribution obtained from the calculation is relatively homogeneous, which leads to a bias of unconservative results in the calculation of the fuel centre temperature.The use of thin rod design for control rods can effectively reduce the self-screening effect, but the use of thin rod structure will lead to lower space utilisation. B4C is commonly used as the absorber material in fast reactors.The wrapping of the control rod absorber with a moderator softens the neutron energy spectrum and increases the B-10 neutron absorption cross section significantly, increasing the value of the control rod. At the same outer diameter, the rod value of a pure absorber core block is less than that of a core block wrapped with a moderator (ZrH2).

    Table and Figures | Reference | Related Articles | Metrics
    A Boundary Variation Diminishing Reconstructed Solver for Mie-Grüneisen Mixture Model
    Zongduo WU, Jin YAN, Jianhua PANG, Yifang SUN, Qingyang MAN
    Chinese Journal of Computational Physics    2024, 41 (3): 277-286.   DOI: 10.19596/j.cnki.1001-246x.8712
    Abstract161)   HTML13)    PDF (2168KB)(564)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    Lattice Boltzmann Model for Simulating Heat-fluid-solid Interaction in Wellbore
    Chunyu HU
    Chinese Journal of Computational Physics    2024, 41 (3): 316-324.   DOI: 10.19596/j.cnki.1001-246x.8715
    Abstract110)   HTML4)    PDF (4161KB)(550)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    Finite Deformation Theory of Poroviscoelasticity Based on Logarithmic Strain
    Xiong TANG, Pei ZHENG
    Chinese Journal of Computational Physics    2024, 41 (3): 287-297.   DOI: 10.19596/j.cnki.1001-246x.8723
    Abstract153)   HTML7)    PDF (3120KB)(543)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    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
    Chinese Journal of Computational Physics    2024, 41 (3): 357-366.   DOI: 10.19596/j.cnki.1001-246x.8714
    Abstract153)   HTML3)    PDF (8011KB)(532)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    Quantum Monte Carlo and Its Applications in Condensed and Warm Dense Matters
    Tianxing MA, Ting GUO, Zhongbing HUANG, Haiqing LIN
    Chinese Journal of Computational Physics    2024, 41 (6): 701-716.   DOI: 10.19596/j.cnki.1001-246x.8974
    Abstract209)   HTML13)    PDF (908KB)(529)      

    This paper firstly introduces the development of quantum Monte Carlo and several types of typical quantum Monte Carlo methods in detail, and then summarizes the recent researches of condensed and warm dense matter systems, including transport, magnetism, superconductivity and thermodynamic properties. Finally, the development prospect of quantum Monte Carlo methods is discussed.

    Reference | Related Articles | Metrics
    A Novel Random Walk Algorithm for Optimal Configuration of Micro-grid
    Lei PAN, Guomin CUI, Ruifang ZHANG, Hongbin LIU, Yuan XIAO, Zhikang YI
    Chinese Journal of Computational Physics    2024, 41 (3): 392-402.   DOI: 10.19596/j.cnki.1001-246x.8703
    Abstract131)   HTML4)    PDF (3818KB)(528)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    Multi-block Local Mesh Refinement Method for LBM with Inversion of Distribution Function
    Rupu WEI, Peng DING
    Chinese Journal of Computational Physics    2024, 41 (5): 643-650.   DOI: 10.19596/j.cnki.1001-246x.8797
    Abstract153)   HTML3)    PDF (1685KB)(512)      

    In response to the complex coupling calculation of LBM (lattice Boltzmann method) local mesh refinement, a method of inverting macroscopic physical quantities to obtain particle distribution functions is proposed to avoid conversion of distribution functions. By rewriting the Boltzmann equation, the relationship (depends on Knudsen number) inside the distribution function is derived, and the inversion calculation equation of the distribution function is given. Based on this, a new multi-block mesh refinement algorithm is designed, and its performance is verified by applying the algorithm to an example. The results show that the algorithm can display the numerical characteristics of the flow field under different Reynolds numbers, and the streamline graph is continuous and complete. The position of the vortex center in the flow field obtained is consistent with literature data. This algorithm can be used for flow field simulations and provides guidance for the development of local mesh refinement.

    Table and Figures | Reference | Related Articles | Metrics
    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
    Chinese Journal of Computational Physics    2024, 41 (3): 345-356.   DOI: 10.19596/j.cnki.1001-246x.8710
    Abstract169)   HTML7)    PDF (13497KB)(509)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    Deflated Preconditioned Conjugate Gradient Solvers for Linear Elastic Crack Problems
    Xingkang LIU, Xingding CHEN, Yunlong YU
    Chinese Journal of Computational Physics    2024, 41 (5): 619-629.   DOI: 10.19596/j.cnki.1001-246x.8793
    Abstract106)   HTML3)    PDF (4564KB)(503)      

    This paper focuses on some efficient deflated preconditioners for static elastic crack problems modelled by the geometrical extended finite element method. We not only construct the deflation subspace matrix which is suitable for linear elastic crack problems, but also give the principle for selecting the deflated mesh nodes. To further accelerate the convergence, we combine the deflation technique with the "crack tip" domain decomposition preconditioners through multiplicative way, and propose efficient adapted deflated preconditioned conjugate gradient solvers which can eliminate the high-frequency and low-frequency errors simultaneously in the iterations. Numerical experiments demonstrate the effectiveness of our algorithm.

    Table and Figures | Reference | Related Articles | Metrics
    Simulation of Droplet Impact on Cable Surfaces with Different Contact Angles
    Jiru HUAI, Peng WANG, Mengyu YANG
    Chinese Journal of Computational Physics    2024, 41 (3): 298-307.   DOI: 10.19596/j.cnki.1001-246x.8717
    Abstract155)   HTML8)    PDF (13751KB)(497)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    An Approximate Shock Wave Formula for Real Gases
    Shesheng XUE, Xirui ZHU
    Chinese Journal of Computational Physics    2024, 41 (3): 308-315.   DOI: 10.19596/j.cnki.1001-246x.8718
    Abstract181)   HTML3)    PDF (1080KB)(497)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    Study on HEMP Environment of Ground Based on G-TF/SF Technique
    Chao YANG, Haiyan XIE, Xinyang ZHAI, Hailiang QIAO, Zaigao CHEN, Yinjun GAO
    Chinese Journal of Computational Physics    2024, 41 (5): 589-595.   DOI: 10.19596/j.cnki.1001-246x.8857
    Abstract163)   HTML2)    PDF (4854KB)(470)      

    For the finite-difference time-domain (FDTD) solution of the half-space problem, the conventional total-field/scattered-field (TF/SF) technique cannot be directly applied to introduce excitation fields for electromagnetic scattering from the ground. In this paper, we extend the generalized total-field/scattered-field (G-TF/SF) method to accurately calculate the high altitude electromagnetic pulse (HEMP) environment in 3D ground scenarios and adopt the convolution perfectly matched layer (CPML) absorption boundary condition around it. The G-TF/SF method effectively eliminates edge effects in numerical calculations for finite ground, thereby improving accuracy in the calculation of the ground electromagnetic environment. Moreover, only the incident field needs to be introduced at the G-TF/SF boundary, avoiding the computation of reflected and transmitted fields.

    Table and Figures | Reference | Related Articles | Metrics
    Phase-Field Fracture Method Based on Eshelby Theory for Heterogeneous PBX
    Yun XU, Yao LONG, Meizhen XIANG, Jun CHEN
    Chinese Journal of Computational Physics    2024, 41 (5): 559-568.   DOI: 10.19596/j.cnki.1001-246x.8859
    Abstract211)   HTML17)    PDF (10193KB)(464)      

    Taking the advantages of tracking discontinuous material surfaces explicitly within the continuous mechanics framework, the phase field method has been successfully applied to study crack propagation and damage in brittle materials. Considering that the complex micro-structural inclusion-matrix interaction dominates the damage nucleation for heterogeneous PBX, we develop a phase field inclusion model based on the Eshelby's inclusion theory, and carry out the numerical simulation study of damage initiation and evolution. The phase field energy consists of the elastic energy and the inclusion-matrix interaction energy. Combining with the Mori-Tanaka method, the effective elastic moduli of heterogeneous PBXs with different volume fractions are derived. For the proposed phase field inclusion model, the nonlinear debonding effects and damage distribution can be characterized by the phase field order parameter directly. It owns explicity in physical mechanism, and completeness in mathematical theory. We apply this model to compute the high volume fraction heterogeneous PBXs with typical circular and polygonal inclusions, and investigate the influences of loading, inclusion shape, volume fraction and computational parameters on the debonding mechanism. Numerical results indicate that the inclusion-matrix micro-structural evolution promotes interface debonding and the formation of macroscopic material failure, which coincides with experimental observations.

    Table and Figures | Reference | Related Articles | Metrics
    Study on Reentrant Arrhythmia Caused by Weak Coupling
    Zhijie WEI, Yuxiang MO, Rongmei LIN, Xiaoke LAN, Guoning TANG
    Chinese Journal of Computational Physics    2024, 41 (5): 670-679.   DOI: 10.19596/j.cnki.1001-246x.8761
    Abstract115)   HTML2)    PDF (4973KB)(455)      

    The Luo Rudy phase I heart model is used in this paper to study the propagation of waves generated by point wave sources in one- and two-dimensional weakly coupled myocardial tissues. It is observed that the weak coupling can lead to the reduction of the wave propagation speed. When there is a large coupling intensity gradient in the myocardial tissue, weak coupling can lead to the frustration of wave propagation and early afterdepolarization showed by cardiomyocytes. When the coupling intensity between cardiomyocytes is low enough, circular motion of wave and complex wave propagation can be observed in myocardial tissues. In addition, we also study the propagation of planar wave in two-dimensional weakly coupled myocardial tissues by constructing two different two-dimensional myocardial tissues with different coupling structures. We find that spiral waves can be spontaneously generated when plane wave propagates in the two myocardial tissues. Increasing the coupling intensity between cells can effectively prevent the circus movement of wave and the spontaneous generation of spiral wave. The physical mechanism of the above phenomenon is analyzed in this paper.

    Table and Figures | Reference | Related Articles | Metrics
    Fast Complex-amplitude Expanded Phase Field Crystal Model for Different Crystals through a Ginzburg-Landau Approach
    Kun WANG, Jun CHEN, Pei WANG, Wenjun HU, Zheng ZHONG
    Chinese Journal of Computational Physics    2024, 41 (5): 547-558.   DOI: 10.19596/j.cnki.1001-246x.8855
    Abstract218)   HTML27)    PDF (13492KB)(440)      

    This work extends the idea of the traditional complex-amplitude expanded phase field crystal (APFC) model using the Ginzburg-Landau approach. A fast structural APFC model is proposed as a quick and effective method for describing different crystal structures. Taking square and rectangular phases as examples, we systematically determine the structure-dependent parameters in the fast structure APFC model and validates its effectiveness through numerical simulations. In particular, when dealing with rectangular phases, it is found that this method not only solves the stability problem of the rectangular phase but also describes the structural phase transition between rectangular and orthorhombic layered phases, demonstrating the capability of the model in describing multiple structural phase transitions. Finally, through simulating the classic rotation-shrinking of a circular grain, we confirm the ability of the model for correctly predicting physical laws and reveal the roles of different crystal symmetries on the rotation-shrinking behavior of the grain. The proposed method in this paper can effectively promote the application of APFC models in the simulation research of more and larger material systems.

    Table and Figures | Reference | Related Articles | Metrics
    Two-level Grad-div Stabilized Finite Element Methods for Steady Incompressible Navier-Stokes Equations
    Yali WANG, Bo ZHENG, Yueqiang SHANG
    Chinese Journal of Computational Physics    2024, 41 (4): 418-425.   DOI: 10.19596/j.cnki.1001-246x.8731
    Abstract122)   HTML4)    PDF (5605KB)(430)      

    Accuracy of the approximate velocity of the steady incompressible Navier-Stokes equations computed by the standard mixed finite element methods is often affected by the pressure. In order to circumvent or weaken the influence of pressure on the accuracy of the computed velocity, by combining grad-div stabilized method with two-level finite element method, this paper presents a kind of two-level grad-div stabilized finite element methods for solving the steady incompressible Navier-Stokes equations numerically. The basic idea of the methods is to first solve a grad-div stabilized nonlinear Navier-Stokes problems on a coarse grid, and then solve, respectively, Stokes-linearized, Newton-linearized and Oseen-linearized Navier-Stokes problem with grad-div stabilization on a fine grid. Numerical examples are given to verify the high efficiency of the two-level grad-div stabilized finite element methods.

    Table and Figures | Reference | Related Articles | Metrics
    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
    Chinese Journal of Computational Physics    2024, 41 (3): 334-344.   DOI: 10.19596/j.cnki.1001-246x.8722
    Abstract94)   HTML3)    PDF (12181KB)(426)      

    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.

    Table and Figures | Reference | Related Articles | Metrics
    Active Learning Algorithm Using Neural Operator Networks and Bayesian Neural Networks: Learning Macroscale Models for Collective Behavior from Microscale Data
    Zhengya GAO, Zhiping MAO
    Chinese Journal of Computational Physics    2024, 41 (6): 783-796.   DOI: 10.19596/j.cnki.1001-246x.8979
    Abstract134)   HTML1)    PDF (5030KB)(419)      

    With the development of artificial intelligence and scientific computing, deep learning plays a significant role in mathematical modeling. In this work we develop an active learning algorithm that uses microscopic data to establish a macroscopic model for collective behavior. Specifically, we take the Cucker-Smale model in this work and develop the corresponding active learning algorithm that integrates neural operator networks and Bayesian neural networks by utilizing microscopic particle data and partial physics. This algorithm is used to efficiently establishes the corresponding macroscopic Euler model through microscopic data. Finally, the effectiveness of the active learning algorithm is validated through one-dimensional and two-dimensional numerical simulations.

    Table and Figures | Reference | Related Articles | Metrics
Share: