Fluid Model Research Articles

1D modeling of plasma streamers at ammonia-air flame conditions

Abstract This work is a self-consistent 1D modeling of streamers in ammonia-oxygen-nitrogen-water mixtures. A fluid model that includes species transport, electrostatic potential, and detailed chemistry was developed and verified. This model is then used to simulate the avalanche, streamer formation and propagation phases, driven by a nanosecond voltage pulse, at different thermochemical conditions derived from a 1D laminar premixed ammonia-air flame. The applicability of the Meek’s criterion in predicting the streamer inception location was successfully confirmed. Streamer formation and propagation duration were found to vary significantly with different thermochemical conditions, due to the difference in ionization rates. The thermochemical state also affected the breakdown characteristics which was tested by maintaining the background reduced electric field constant. Detailed kinetic analyses revealed the importance of O(1D) in the production of key radicals, such as O, OH, and NH2. Furthermore, the contributions of the dissociative electronic excitation of NH3 towards the production of H and NH2 radicals have also been reported. Spatial and temporal evolution of the electron energy loss fractions for various inelastic collision processes at different thermochemical states uncovered the input plasma energy spent of fuel dissociation and the large variability in the dominant processes during the avalanche and streamer propagation phases. The methodology and analyses reported in this work are key towards developing effective strategies for controlled nanosecond-pulsed non-equilibrium plasma sources used for ammonia ignition and flame stabilization.

Aerodynamic breakup of emulsion droplets in airflow

AbstractAerodynamic breakup refers to the process where large droplets are fragmented into small droplets by the aerodynamic force in airflow, which plays a vital role in fluid atomization and spray applications. Previous research has primarily concentrated on the aerodynamic breakup of single‐component droplets, but investigations into the breakup of emulsion droplets are limited. This study experimentally investigated the aerodynamic breakup of water‐in‐oil emulsions in airflow, utilizing high‐speed photography to observe the breakup process and digital in‐line holography to measure fragment sizes. Comparative analyses between emulsion droplets and single‐component droplets are conducted to examine the breakup morphology, breakup regime, deformation characteristics, and fragment size distributions. The emulsion droplets exhibit higher apparent viscosity and shorter stretching lengths of the bag film and peripheral rim due to the presence of a dispersed phase. The breakup regime transitions of emulsions are modeled by integrating the viscosity model of emulsions and the transition model of the pure fluid. The fragment sizes of emulsion droplets are larger due to the shorter lengths of the bag film and peripheral rim.

Numerical Simulation Study on Optimization of Development Parameters of Condensate Gas Reservoirs

Due to the retrograde condensation phenomenon in the development process, the fluid phase change is complex, and it becomes particularly difficult to accurately describe the fluid flow characteristics and residual oil and gas distribution characteristics during the development of condensate gas reservoirs. It is difficult to select the development program and subsequent dynamic adjustment for the efficient, reasonable, and sustainable development of condensate gas reservoirs. In this paper, the phase characteristics of condensate gas reservoirs are clarified; the basic fluid model is created by using computer modeling, using Win-Prop; and in view of the characteristics of the target condensate gas reservoirs, the CMG (Computer Modeling Group) numerical simulation method is applied to study the effects of six factors, the thickness of the reservoir, permeability, porosity, rock compression coefficient, the ratio of the vertical permeability to the horizontal permeability, and the injection of different media, on the development effect through the study of different development parameters of gas condensate reservoirs. The purpose of this study is to provide guidance for the rational development of condensate gas reservoirs in practical production.

Mathematical Modeling on the Multiphase Flow during Top–Bottom Combined Blowing Basic Oxygen Furnace Steelmaking Process

A 3D mathematical model on the multiphase flow and splashing phenomenon in a 210 t Basic oxygen furnace converter is established. The turbulent multiphase flow of the molten steel, slag, and oxygen is simulated by the combined realizable k−ε model and volume of fluid model. The trajectory of argon bubbles and the interaction between the molten steel and argon is simulated by the two‐way coupled Lagrangian discrete phase model. Effects of the lance height and bottom‐blowing gas flow rate on the multiphase flow and splashing phenomenon are investigated. The results indicated that the cavity depth is decreased from 0.45 to 0.36 m when the lance height increases from 1700 to 2300 mm. The cavity depth is not significantly affected by the bottom‐blowing argon flow rate with the 44 000 m3 h−1 top‐blowing oxygen flow rate and 1700 mm lance height. The increased lance height resulted in a decrease in splashing. When the oxygen lance height increases from 1700 to 2300 mm, the amount of the splashed molten steel decreases from 3774.3 to 1953.0 kg.

Thermodynamic and Kinetic Behavior Study of Gas Bubbles in High‐Temperature Steel Melts

This study provides an exhaustive exploration of the thermal expansion of low‐temperature argon bubbles in high‐temperature molten steel and the resulting dynamic effects. The energy equation of compressible fluid simultaneously with the volume of fluid model under the Euler framework is solved. The generation and rising behavior of bubbles are validated through rigorous physical water model experiments. Findings reveal that, within a heating time of 0.2–0.45 ms, gas bubbles already complete 65–90% of their internal energy increase. The total heat transfer time and heat transfer coefficient of bubbles exhibit distinctive patterns influenced by the initial diameter and temperature of the bubbles. Specifically, under the conditions studied, bubbles exhibit a total heat transfer time of ≈20–60 ms, accompanied by a heat transfer coefficient ranging from 500 to 2600 W (m2 * K)−1. During the thermal expansion process, gas bubbles experience energy transfer and conversion within the two‐phase flow. This article establishes a dynamic model describing bubble thermal expansion in steel melts based on underwater bubble dynamics’ first principles. The model provides a quantitative depiction of instantaneous changes in bubble size, velocity, and kinetic energy under specific conditions.

Construction of Kinetic Model of Vacuum Arc Magnetic Fluid Under External Magnetic Field

The external magnetic field is an important control method of vacuum arc at present. Because of its high efficiency and economy, it has become an important method to study the kinetic model of vacuum arc magnetic fluid. Therefore, the kinetic model of vacuum arc magnetic fluid under the external magnetic field is constructed. Firstly, the mass conservation equation, radial momentum conservation equation, axial momentum conservation equation, energy conservation equation, electromagnetic equations are analyzed. The boundary conditions of the dynamic model are cathode boundary, anode boundary and fluid wall boundary. The fluid dynamics calculation software FLUENT is used as the calculation platform, and its secondary development is carried out. Combining VC++self-programming and self-defined macro functions, according to the boundary conditions, the kinetic model of vacuum arc magnetic fluid under external magnetic field is solved, and the kinetic model analysis of vacuum arc magnetic fluid under external magnetic field is realized. Experimental analysis shows under the action of longitudinal magnetic field, the intensity of longitudinal magnetic field increases, the maximum point of ion rotation speed in the magnetic fluid of vacuum arc moves outward along the edge of vacuum arc, the number density of electrons gradually moves from the central axis to the radial direction, and the high density area increases first and then decreases; Under the action of transverse magnetic field, the intensity of transverse magnetic field increases, the ion pressure in vacuum arc magnetic fluid gradually shifts and increases, and the ion number density also increases.

Thermal management for the prismatic lithium-ion battery pack by immersion cooling with Fluorinated liquid

This study constructs a novel FS49-based battery thermal management system (BTMS), proposing an optimization method for the system energy density and an indirect control method for the system cooling capacity. The boiling of dielectric refrigerant occurred at the battery surface, which provided strong and uniform cooling for each battery cell. The results show that the peak temperature difference of liquid immersion cooling (LIC) module during 1C rate discharging and charging was reduced by 91.3% and 94.44%, respectively, compared to the natural convection (NC) module. The reduction of temperature nonuniformity greatly reduced the state of charge (SOC) inhomogeneity of different cells within the module. Moreover, the energy density of LIC module can be optimized by reducing the cell spacing and liquid filling ratio. Comparing with 100% filling rate, the module maximum temperature corresponding to 25% filling rate (with wick) during 2C rate discharging is only increased by 1.6℃, however, the module peak temperature difference is reduced by 53.3%, and the energy density is increased by 13.14%. In addition, Equivalent circuit model (ECM) and volume of fluid (VOF) model were used to simulate the LIC module, and the numerical results are in good agreement with the experimental data. This study provides a systematic design plan and numerical method for the engineering application of LIC in the field of BTMS.

Numerical simulation of He atmospheric pressure plasma jet impinging on the tilted dielectric surface

The target surface to be treated in reality is often not smooth and horizontal and may also be in different tilting angles. The treatment of the tilted dielectric surface by the atmospheric pressure plasma jet (APPJ) undoubtedly increases the complexity of surface modification. Therefore, a two-dimensional fluid model is established to reveal the internal mechanism of the interaction between the He APPJ and the tilted dielectric surface by means of numerical simulation. The distribution of the gas flow in a small angular range (0°, 3°, 5°, 8°, 10°, and 15°) is studied. In addition, the effects of the tilt angle on the jet morphology, discharge dynamic properties, and species distribution of the He APPJ are emphatically discussed. It is found that the jet morphology and parameters are no longer symmetrical under the tilted surface. With the increase in the tilt angle, the enhanced electric field in the upper surface region leads to the increase in the ionization rate and electron density here, and also accelerates the propagation speed of the jet to the dielectric surface in the atmospheric environment. Driven by the electric field force, the jet is closer to the dielectric surface, resulting in a decrease in the thickness of the cathode sheath and an increase in the surface charge density in the area to the right of the central axis. The influence of the gas flow structure leads to the shortening of the jet development distance and a decrease in the jet velocity on the upper surface. N and O also form higher fluxes on the upper surface due to the increase in the electron density.

A meshless model for three-dimensional direct numerical simulation of local scour around cylinder bridge foundation

A meshless local scour model for cylinder bridge foundation based on the smoothed particle hydrodynamics (SPH) method is proposed in this paper. Differing from the traditional Eulerian mesh model used for monopile local scour, it overcomes the high requirements of mesh quality and density for capturing large deformations at the fluid interface. Unlike the transient sediment erosion SPH model used for dam-break, sediment flushing, and water jet, the meshless approach is innovatively applied to the continuous simulation of bridge local scour under steady flow conditions. The present SPH model is established based on the SPH multi-phase fluid model, integrating a classical fluid model combined with numerical stabilization algorithms to govern water particle motion and a sediment model from Zhang et al. (2024) to govern sediment evolution, and tests a new scheme for boundary treatment in the model based on DualSPHysics. Furthermore, a novel scour visualization method that aims at extracting both visible and actual bed surfaces from discrete particles is proposed. The model is validated through comparison with rigid bed and live bed experiments, demonstrating favorable agreement in terms of flow and scour simulation. Importantly, the model results reveal a discrepancy between the elevation of the post-scour visible bed surface derived from the conventional Eulerian mesh model and experimental data, compared to the actual bed surface unaffected by water flow erosion. This indicates the necessity for incorporating a safety margin in foundation design when utilizing experimental results to determine the maximum scour depth.

Reliability of two-dimensional steady magnetized Jeffery fluid over shrinking sheet with chemical effect

Abstract A numerical framework is established for a two-dimensional steady flow of the magnetized Jeffery fluid model over elongated/shrinking sheets, with potential applications such as coating sheets, food products, fiber optics, drilling fluids, and the manufacturing processes of thermoplastic polymers. The model also demonstrates the influence of chemical reaction, magnetic field, and stability analysis which provide a novel contribution to this study. To ensure the ease and effectiveness of this analysis, we transform the set of difference equations governing the system into ordinary equations using the similarity transformation. The reliability of the solution is examined by using stability analysis. The Navier–Stokes equations have been transformed into self-similar equations by adopting appropriate similarity transformations and subsequently solved numerically using the bvp4c (three-stage Labatto-three-A formula) approach. The comparison between the derived asymptotic solutions and previously documented numerical results is quite remarkable. The self-similar equations display a duality of solutions within a limited range of the shrinking parameter, as observed from the data. For each stretching scenario, there is a unique solution. Hence, an examination of temporal stability has been conducted through linear analysis to establish the most fundamentally viable solution. The determination of stability in the analysis is based on the sign of the smallest eigenvalue, which indicates whether a solution is unstable or stable. The analysis of stability reveals that the first solution, which describes the primary flow, remains stable. Through the utilization of graphs, we thoroughly examine and discuss the influence of emerging factors. The numerical results obtained from this analysis demonstrate multiple solutions within a certain range of M 1 ≥ M ci {M}_{1}\ge {M}_{{ci}} , i = 1 , 2 , 3 i=1,2,3 , and no solution in the range M 1 < M ci {M}_{1}\lt {M}_{{ci}} . M ci {M}_{{ci}} denotes the critical values, which increase as the quantities of Sc {\rm{Sc}} increase from 0.3 to 0.9. Similarly, multiple solutions exist for λ ≥ λ ci \lambda \ge {\lambda }_{{ci}} , i = 1 , 2 , 3 i=1,2,3 , and no solution in the range λ < λ ci \lambda \lt {\lambda }_{{ci}} is observed.

Genomic analysis and functional properties of Lactobacillus johnsonii GJ231 isolated from healthy beagles.

Probiotics are one of the management tools to improve the host's healthy microbiota. The positive effects of probiotics on host health are species-specific, so probiotics isolated from host's own gut may be most beneficial. Many of the metabolites (e.g., short-chain fatty acids, bacteriocins, and hydrogen peroxide) produced by Lactobacillus johnsonii have specific inhibitory profiles against invading pathogens. In this study, we isolated L. johnsonii GJ231 from the intestinal tract of healthy female beagles. The genome size of 1.763 M encoded a total of 1,691 predicted genes. Many carbohydrate-active enzymes responsible for carbohydrate degradation and the production of short-chain fatty acids were also predicted. The metabolic profile of short-chain fatty acids in L. johnsonii GJ231 was determined using LC-MS/MS. The bacteriocin-producing gene bacteriocin (lactacin F) in L. johnsonii GJ231 was also predicted. In vitro, experiments demonstrated that GJ231 can thrive in weak acids, 0.3% bile salts, and artificial gastrointestinal fluid models. It was tolerant of to high temperatures up to 70°C, was non- hemolytic, inhibited pathogenic bacteria, and had a high antioxidant capacity. In vivo safety experiments conducted in mice revealed that oral administration of GJ231 not only had no toxic side effect but also increased their antioxidant capacity. In conclusion, combining the above test results, which collectively demonstrate that canine-derived L. johnsonii GJ231 was a non-pathogenic, acid-tolerant and bile-salt-tolerant probiotic strain that inhibits pathogenic bacteria and improves host antioxidant function. This may make it a promising candidate for the development of innovative functional foods for pets.

Improved efficiency of continuous conversion of alcohols to chlorides in flow microreactors with DESs/HCl via predictive fluid model

The efficient and safe synthesis of chlorides is attractive in pharmaceutical chemicals yet presents challenges. Herein, a novel route has been developed for the conversion of alcohols into the corresponding alkyl chlorides in continuous flow microreactor. Meanwhile, the system utilizes deep eutectic solvents (DESs) hydrochloride solutions (DESs/HCl) as both solvents and catalysts for this reaction. The results suggest that the premixed mode continuous flow microreactor using DESs/HCl can significantly improve the harsh reaction conditions in the synthesis reaction of chlorides. Furthermore, tuning the HCl concentration under 383 K and 30 min conditions significantly increased the yields of n-Butanol chloride and Benzyl chloride to 86 % and 98 %, respectively. This work shows the great application potential of the continuous flow microreactor combined with DESs/HCl to convert alcohols into corresponding alkyl chlorides or aryl chlorides. Besides, the fluid model and HCl absorption kinetics were established based on the calculated diffusion coefficients (DR), Damköhler number (Da), Bodenstein number (Bo), residence time (τ), and experimental results. The flow model of the mixture of DESs/HCl and alcohols in the microreactor was a plug flow, which is beneficial for improving the efficiency of this reaction.

Laminar planar jets of elastoviscoplastic fluids

We perform numerical simulations of planar jets of elastoviscoplastic (EVP) fluid (Saramito (2007) model) at low Reynolds number. Three different configurations are considered: (i) EVP jet in EVP ambient fluid, (ii) EVP jet in Newtonian ambient fluid (miscible), and (iii) EVP jet in Newtonian ambient fluid (immiscible). We investigate the effect of the Bingham number, i.e. of the dimensionless yield stress of the EVP fluid, on the jet dynamics, and find a good agreement with the scaling for laminar, Newtonian jets for the centerline velocity uc/u0∝(x/h)−1/3 and for the jet thickness δmom/h∝(x/h)2/3 at small Bingham number. This is lost once substantial regions of the fluid become unyielded, where we find that the spreading rate of the jet and the decay rate of the centerline velocity increase with the Bingham number, due to the regions of unyielded fluid inducing a blockage effect on the jet. The most striking difference among the three configurations we considered is the extent and position of the regions of unyielded fluid: large portions of ambient and jet fluids (in particular, away from the inlet) are unyielded for the EVP jet in EVP ambient fluid, whereas for the other two configurations, the regions of unyielded fluid are limited to the jet, as expected. We derive a power law scaling for the centerline yield variable and confirm it with the results from our numerical simulations. The yield variable determines the transition from the viscoelastic Oldroyd-B fluid (yielded) to the viscoelastic Kelvin–Voigt material (unyielded).

Solute Dispersion in Casson Blood Flow through a Stenosed Artery with the Effect of Temperature and Electric Field

This study develops a mathematical model to address solute dispersion in arterial blood flow, particularly considering the effects of temperature and electric fields using the Casson fluid model. Given the physiological importance of drug delivery within the human vascular system, understanding these dynamics is crucial, especially in the presence of cardiovascular diseases (CVD). The Generalized Dispersion Model (GDM) is employed to simulate the dispersion function. Results demonstrate that elevated temperatures and electric fields enhance drug uptake and distribution by increasing blood flow. The model accurately predicts drug behavior in narrowed arteries, optimizing dosage regimens and improving therapeutic outcomes. Visual representations and detailed discussions of these findings are presented in the subsequent sections. The results are validated against a previous study that did not consider the effects of stenosis height, temperature and electric field. The findings suggest a strong agreement between the two studies. Additionally, an increase in mean absorption leads to an increase in velocity. When mean absorption increases, the functions of steady dispersion and overall dispersion decrease, while the function of unsteady dispersion increases. The reverse behavior is observed for the electric field. The study concludes that integrating temperature and electric field effects into drug delivery systems can significantly enhance treatment efficacy and reduce side effects, providing a valuable tool for advanced targeted therapies in CVD and cancer management.

A higher dimensional model of geophysical fluid with the complete Coriolis force and vortex structure

Here, we present a higher dimensional model from the vorticity equation, which describes the dynamic characteristics of large scale Rossby waves by utilizing the Gardner-Morikawa coordinate transformation and the perturbation method. To reveal the influence of physical parameters on the higher dimensional model, we first give the dispersion relation of the model and the N-soliton solutions by Hirota method. Subsequently, the lump solutions are derived by using the long wave limit method. It demonstrates that the horizontal component of the Coriolis force acts as a forcing force on the nonlinear Rossby waves, and affects the amplitude of the meridional structure. Moreover, under the background of secondary zonal basic flow, for different lump solutions, the flow field will appear dipole blocking or double vortex structure. It is also indicated that the horizontal Coriolis force only causes the vortex to move in the latitudinal direction.

Innovative molecular descriptors in QSPR modeling: Integrating Carnahan-Starling EoS for predicting diffusion coefficients in hydrocarbons and mixtures

The rapid and accurate prediction of self-diffusion coefficients for pure fluids and mutual diffusion coefficients for binary mixtures remains a focal point of current research. In this work, molecular descriptors are innovatively defined by introducing molecular repulsive and attractive pressures based on the Carnahan-Starling equation of state. This approach aims to improve the quantitative structure–property relationship (QSPR) models for predicting thermophysical properties that vary with thermodynamic conditions. The shapely additive explanation method is employed for a comprehensive analysis of the model mechanism. Among the six molecular descriptors selected, our defined attractive pressure (patt), MATS2e and SpMin3_Bh(v) contribute significantly to the self-diffusion coefficient QSPR model. Based on this approach, QSPR models are constructed for 15 hydrocarbons and five of their binary mixtures. The results demonstrate that the GA-BPNN-based QSPR model for pure fluids achieves high prediction accuracy, with an external Rext2 of 0.978 and an external RMSEext of 5.349 × 10−9·m2·s−1. The optimal mixture descriptor type is x1d1+x2d2, and the GS-SVM-based mixture QSPR model predicts 98.24 % of values within a 5 % error margin, with Rext2 of 0.991 and RMSEext of 0.042 × 10−9·m2·s−1.

Research on the impact of overflow well shut-in on wellbore integrity and safety

If an overflow occurs, the well must be closed immediately. To ensure the safety of the drilling wellbore in ultra-deep wells, it is critical to analyze the transient flow process and pressure fluctuations of the fluid in the annular space of the overflow well shut-in. This paper investigates the pressure fluctuation model of multi-phase fluid in the annular space after overflow occurs during gas well drilling, analyzes the influence factors of multi-phase fluid's pressure wave velocity, and simulates and analyses pressure fluctuation changes in the annular space in light of the impact of overflow well shut-in on wellbore safety. The model's correctness is verified by comparing it to experimental data published in the literature. The model is developed and examined using the example of a genuine gas well in western China. According to the study, as soon as the BOP closes, pressure rapidly increases to a high of 6.09 MPa at 5.89 s before falling to a valley of 4.68 MPa at 35.36 s. When the mild shut-in time lasts 20 s, the wellhead annular pressure reaches a maximum of 5.8 MPa. If the wellhead equipment can withstand the annular space's first pressure wave, a hard shut-in is recommended. To optimize the overflow well shut-in system and ensure wellbore safety, the change in pressure in the wellhead annular space as a result of varied shut-in timings and drilling fluid displacement is studied.

Numerical simulation of the water flow field on subway stairs and evaluation of people safety evacuation

ABSTRACT Stairs in subway stations are vulnerable to floods when rainstorm disasters occur in cities. The stairs, as a critical way for human evacuation, can affect the safe evacuation of people on flood-prone stairs. To evaluate the risk of people evacuating through different slopes and forms of stairs when floods invade subway stations, a numerical model for the water flow on stairs based on the volume of fluid model and the realizable k-ε model was established. The water flow patterns on stairs at the subway station entrance under different slope conditions and with/without rest platforms were simulated. The real-time water flow process on stairs at different inlet depths was obtained, and the escape control index F was used to evaluate the risk of people evacuating through stairs at different slopes and water depths. The results indicate that the presence of a rest platform can cause an increase in water velocity and depth on pedestrian stairs, and people should choose stairs without a rest platform for evacuation during the evacuation process. The research results hope to provide a reference for the people evacuation on stairs, and further improve the theory of safe evacuation of personnel on flood-prone stairs.

Scrutiny of pseudoplastic nanofluid flow under the influence of magnetic hydrodynamics with chemical reaction across the cylinder with slip boundary condition

Buoyantly induced flows possess diversified utilizations in numerous engineering processes for instance, reactor cooling through passive strategy, LED lights, pipes manufacturing, ship funnel and many more. Galvanized sheets sticked to form a hollow cylinder is used for fast cooling in naval ship and chimney of steam power plant. In view of such mesmerizing significance of flow and thermal attributes of fluid over vertical cylinder, current effort is articulated. For this purpose, Williamson fluid model is accounted and nanoparticles are also induced to envision advanced thermal features in the flow over vertically oriented cylinder. Physical effectiveness of magnetic field implications and chemical reactive species are entertained to notice change in hydrothermal and mass fields. Slip boundary constraint along with stagnant flow at the surface of cylinder is considered to inspect behavior in free stream region. The fundamental governing equations of problem are attained in the sense of PDEs by utilizing the concept of boundary layer approximation and converted into coupled nonlinear ODEs by substituting specific similar variables. Numerically the resulting ODEs are resolved by making the use of bvp4c built in MATLAB routine, the outcomes are attained and arranged in graphical and tabular manner. The influence of pertinent flow parameters such as (0.1 ≤ We ≤ 1.0), (1.0 ≤ M ≤ 3.0), (1.0 ≤ Pr ≤ 11), thermophoresis parameter (0.5 ≤ Tp ≤ 4.5), (1.0 ≤ Bp ≤ 5.0), (0.05 ≤ Nr ≤ 0.7), and (0.5 ≤ Kc ≤ 5.0) are examined on the momentum, thermal and concentration distributions. It is manifested that the velocity distribution depreciates by uplifting magnetic field strength. Increment in magnitude of wall drag and convective heat transfer is perceived by enhancing Prandtl number. It is inferred that friction coefficient in absolute sense and heat flux coefficient tends to exceed against elevation in (We). From comprehensive analysis, declination in velocity and thermal field is depicted versus Reynold number whereas, contrary aspects are visualized in concentration. Enhancement in the value of surface drag and thermal flux is revealed versus (Re).

Simulation on the hollow cathode discharge in hydrogen

A rectangular hollow cathode discharge (HCD) in hydrogen with a pressure of 2 Torr is simulated using a 2-D fluid model. The potential, electric field, particle density, and average electron temperature are calculated. The discharge space consists of the cathode sheath region near the cathode electrode and the negative glow (NG) region in the central region of the discharge cell. A high electric field of thousands of V/cm and a low electric field of tens of V/cm appear in the cathode sheath region and NG region, respectively. The average electron temperature in the cathode sheath region is tens of eV, which is significantly higher than that in the NG region. Electrons and H3 + are the main negative particles and positive ions, whose peaks appear in the NG region, and the peak magnitude is on the order of 1010 cm−3. H atom is the highest-density neutral particle other than H2 with a peak density of 1013 cm−3. The reaction kinetics of the generation and consumption of different particles are explored. The results show that each reaction generates certain particles while consuming other particles, ultimately achieving a dynamic equilibrium in the density of various particles. The electrons mainly originate from the ground state ionization between electron and H2 (e+H2 → e+H2 ++e) and are consumed by the dissociative attachment (e+H2 → H−+H). The charge transfer collision reaction (H2 ++H2 → H3 ++H) is the only reaction that produces H3 + ions. Different reactions to the consumption of H3 + ions do not differ significantly. The generation and consumption of H mainly originate from the electron collision dissociation reaction (e+H2 → e+H+H) and the ionization reaction (e+H→H++2e).