Liquid-gas Phase Research Articles

Insights into hydrophobicity‐tuned catalyst microenvironment in gas diffusion electrode for boosting CO2 electrolysis

AbstractCO2 electrolysis to value‐added chemicals and fuels shows tremendous potential to store renewable electricity. It features with H2O as the hydrogen source of CO2 electroreduction. Here, we aim at the effect of the hydrophobicity‐tuned catalyst microenvironment in gas diffusion electrode (GDE) on boosting CO2 electrolysis. The hydrophobic polytetrafluoroethylene (PTFE) nanoparticles are introduced in the Sn catalyst layer of GDE. The PTFE‐modified GDE compared to the non‐PTFE electrode delivers a higher HCOOH Faraday efficiency and lower overpotential at current densities of 50–250 mA cm−2. A two‐dimensional axisymmetric model jointly with electrochemical impedance spectroscopy is developed to provide mechanistic insights into the hydrophobicity regulation. The simulated polarization current densities are in keeping with experimental results. The boosted CO2 electrolysis is attributed to the decrease of the distance from the gas–liquid phase boundary to the electrocatalyst. The understanding of the catalyst microenvironment established here can be generalized to other gas‐involved electrocatalytic reactions.

Continuous ultrasonic ozone coupling technology-assisted control of ceramic membrane fouling coupled enhanced multiphase mixing to treat dye wastewater and CFD flow field simulation

In this study, ozone catalysts (hydrogenation-modified red mud, HM-RM) successfully prepared by hydrogenation-modification of industrial hazardous solid waste red mud (RM) as a raw material in accordance with the viewpoint of treating waste with waste and using waste. Meanwhile, as for the common phenomenon of membrane fouling, uneven distribution of multiphase solid catalysts and ozone in liquids, the addition of ultrasound can not only disperse materials, but also play a role in online cleaning of ceramic membranes and catalysts. The optimum treatment conditions for Rhodamine B (RhB) solution with volume of 2 L and concentration of 40 mg/L were catalyst concentration of 0.4 mg/L, reaction temperature of 45 °C, ultrasonic time of 1 h, ultrasonic intensity of 600 W, removal rate of RhB was up to 90 %. In addition, the computational fluid dynamics (CFD) simulation method was used to investigate the fluid flow between the two gas-liquid phases and the effect of the negative pressure of the membrane pump on the fluid by the analysis of flow, pressure and ozone flux of the ceramic membrane(CM) reaction apparatus. The CFD simulation results showed that at the inlet gas-liquid flow rate of 3 m/s and the negative pressure of 20,000 Pa, the maximum flow rates of CM-1 were 3 m/s, 0.752 m/s for CM-2, and 0.228 m/s for CM-3, respectively. Vortices, which are beneficial to solid-liquid mixing and gas-liquid mass transfer, formed between the suction port CM-1 of CM-1 and the inlets of CM-2 and CM-3. This discovery is consistent with relevant experimental research results. Significantly higher concentrations of both •OH and dissolved ozone were observed in the US/HM-RM/O3 system compared to other systems, indicating the significant improvement in ozone utilization rate through the application of ultrasound. The superiority of the US/HM-RM/O3 device was demonstrated. The real dye effluent was tested under optimum operating conditions and the results showed that COD and TOC were reduced by 81.34 % and 60.23 % respectively after 180 min of treatment. The above research can provide technical support for the treatment of dye wastewater using Ultrasound-enhanced ozone oxidation ceramic membranes.

New fundamental equations of state for binary hydrogen mixtures containing argon, helium, and neon

In this work, new equations of state for the binary mixtures H2 + Ar, H2 + He, and H2 + Ne are presented. The equations are formulated in terms of the reduced Helmholtz energy and allow for the calculation of all thermodynamic properties over the entire fluid surface including the gas phase, liquid phase, supercritical region, and equilibrium states. The models are validated by comparisons with experimental data and their physical behavior is analyzed. Furthermore, the new equations of state are compared to other state-of-the-art models from the literature.

Optimization of methyl ester synthesis using gas/liquid phase pulsed discharge plasma in a novel oscillatory slug flow reactor

In the present research, an innovative oscillatory slug flow reactor (OSFR) under the treatment of gas–liquid phase pulsed discharge plasma was developed for biodiesel production. The main goal was continuous production of high quality biodiesel at low temperature and pressure. Experimental tests were carried out under the influence of four main operating parameters including applied voltage, molar ratio, reactant flow rate, and catalyst concentration. The response surface method was employed to optimize experimental tests. The results showed that the proposed technology provided 94% production efficiency under the optimal conditions of voltage 19.4 kV, molar ratio 6.4, flow rate 2.7 ml/s, and catalyst 0.9 wt. %. According to the statistical analysis, increasing the applied voltage and reducing the flow rate have a strong effect on the Fatty Acid Methyl Ester yield, while the concentration of potassium hydroxide and methanol have less effect on the overall efficiency. In addition, the characteristics of the produced biodiesel were in accordance with ASTM D6751 standards. Surprisingly, the optimal energy consumption in this system was 95 kJ/l, which is more economically viable. In general, this study showed that the combined system of gas/liquid phase plasma in the OSFR reactor has a high synergistic potential for the transesterification reaction.

Advancing immunotherapy in triple negative breast Cancer: A novel multimodal theranostic nanoplatform integrating synergetic ferroptosis and photothermal therapy

Triple negative breast cancer (TNBC) is a subtype of breast cancer with high invasiveness and the worst prognosis. Immunotherapy has demonstrated significant potential in its treatment. However, the response rate remains unsatisfactory. In this work, we prepared multifunctional theranostic nanoplatforms (P-R@P/U-V) with synergetic ferroptosis and photothermal therapy (PTT) in improving TNBC tumor microenvironment for improving local treatment in situ, and synergistic anti-programmed cell death-ligand 1 (PD-L1) immunotherapy to prevent distant metastasis. With P-R@P/U-V, real-time monitoring of the specific distribution in tumors is possible due to the greatly improved ultrasound imaging (US), photoacoustic imaging (PAI), and magnetic resonance imaging (MRI) capabilities. Moreover, P-R@P/U-V has excellent photothermal conversion ability. Whereas PTT killed TNBC cells after laser irradiation, the liquid–gas phase transition not only improved US but also triggered controllable release of RSL3 in situ, enhancing ferroptosis of TNBC cells. Besides, the release of tumor-associated antigens (TAAs) increased via PTT and ferroptosis, strongly stimulating T cell activation, dendritic cell maturation, and anti-tumor cytokine secretion. The lung metastases of TNBC were significantly inhibited following the combined use of P-R@P/U-V and PD-L1 inhibitors. The multimodal and multifunctional nanoplatforms demonstrated excellent effects in early diagnosis and synergistic treatment of TNBC. This represents a novel strategy based on local treatment for additional enhanced anti-PD-L1 immunotherapy, with high clinical application prospects.

Mechanocatalytic Synthesis of Ammonia by Titanium Dioxide with Bridge-Oxygen Vacancies: Investigating Mechanism from the Experimental and First-Principle Approach.

Mechanochemical ammonia (NH3) synthesis is an emerging mild approach derived from nitrogen (N2) gas and hydrogen (H) source. The gas-liquid phase mechanochemical process utilizes water (H2O), rather than conventional hydrogen (H2) gas, as Hsources, thus avoiding carbon dioxide (CO2) emission duringH2 production. However, ammonia yield is relatively low to meet practical demand due to huge energy barriers of N2 activation and H2O dissociation. Here, six transition metal oxides (TMO) such as titanium dioxide (TiO2), iron(III) oxide (Fe2O3), copper(II) oxide (CuO),niobium(V) oxide(Nb2O5), zinc oxide (ZnO), and copper(I) oxide (Cu2O) are investigated as catalysts in mechanochemical N2 fixation. Among them, TiO2 shows the best mechanocatalytic effect and the optimum reaction rate constant is 3.6-fold higher than the TMO-free process. The theoretical calculations show that N2 molecules prefer to side-on chemisorb on the mechano-induced bridge-oxygen vacancies in the (101) crystal plane of TiO2 catalyst, while H2O molecules can dissociate on the same sites more easily to provide free H atoms, enabling an alternative-way hydrogeneration process of activated N2 molecules to release NH3 eventually. This work highlights the cost-effective TiO2 mechanocatalyst for ammonia synthesis under mild conditions and proposes a defect-engineering-induced mechanocatalytic mechanism to promote N2 activation and H2O dissociation.

Generalized Maxwell equal area law and black holes in complex free energy

Maxwell equal area law is an important and traditional analytical tool in thermodynamic phase transition, especially in the calculation of gas-liquid phase transition, which reconciles the theoretical calculation with the experimental results. Undoubtedly, its importance is also self-evident for the black hole thermodynamic system. In this study, we construct a generalized Maxwell equal area law, which allows different states of thermodynamic systems to be within the generalized free energy. The black hole thermodynamic characteristics are spontaneously emerged in the free energy landscape. Furthermore, by analytic continuation, we utilize the properties of analytical functions to investigate some universal characteristics of thermodynamic phase transitions in black holes, and preliminarily establish the counterpart of thermodynamic phase transitions in the complex domain.

Preparation of monocyclic aromatic hydrocarbons from industrial lignin residue and polyethylene co-pyrolysis by microwave-assisted in fluidized bed based on bimetal-loaded HZSM-5/MCM-41 core-shell catalyst

Enzymatic hydrolysis lignin (EHL), lignosulfonate (LS), alkali lignin (AL), and hydrolysis lignin (HL) residues are mostly burned for heating in industrial production, which greatly wastes resources and pollutes the environment. In this study, we investigated the effects of different temperatures and catalyst-to-biomass mass (C/B) ratios on the relative content of hydrocarbons in the products by microwave-assisted pyrolysis in a fluidized bed. Co and Ni bimetal-loaded molecular sieve with core-shell structure (Co-Ni@H/M) was used as the catalyst. The main component of waste plastic, polyethylene (PE), was pyrolyzed with these four lignin residues to prepare monocyclic aromatic hydrocarbons (MAHs), and the mechanism of synergistic interaction in co-pyrolysis was investigated. It was shown that the maximum liquid phase yield and hydrocarbons relative content were obtained at 550 °C when lignin residues were pyrolyzed alone. Compared with the pyrolysis of lignin residues alone, the liquid phase and gas phase yields of co-pyrolysis were improved, while the solid phase yield decreased. Moreover, the temperature for maximum liquid phase yield and hydrocarbons relative content were reduced to 500 °C. Further exploration of the liquid phase products revealed that the relative contents of hydrocarbons in the four lignin residues and PE (EHL-PE, LS-PE, AL-PE, and HL-PE) co-pyrolysis products at 500 °C were 49.6 %, 33.1 %, 37.7 %, and 43.2 %, respectively, with MAHs as the main components. The results of the effect of different C/B ratios on the distribution of pyrolysis products showed that Co-Ni@H/M promoted the deoxygenation of oxygenated compounds, especially phenolic compounds, and the conversion of aliphatic hydrocarbons into MAHs by cyclization and aromatization.

A novel synergistic intensification method for CO2 absorption by metal salt and task-specific ionic liquids hybrid solvent in microchannel reactor

A novel metal salt-promoting mass transfer intensification method for CO2 absorption into task-specific ionic liquids aqueous solution was proposed. The mass transfer intensification and intrinsic mechanisms of different concentrations of metal salts (NaBF4, KBF4, KCl, and CuCl2) with task-specific ionic liquids (1-aminopropyl-3-methylimidazolium tetrafluoroborate [Apmim] [BF4]) aqueous solutions for CO2 capture process were visually investigated using a microchannel reactor. The results showed that adding metal salts could remarkably enhance the interaction between metal cations and ionic liquid anions, thereby decreasing the interaction between anions and cations within the ionic liquid molecules and promoting the chemical reaction between ionic liquid cations and CO2. Under the condition of the same salt concentration, the addition of KBF4 has more obvious enhancement for mass transfer in CO2 absorption process. The mass transfer coefficient increases as the gas-liquid phase flow rate rise under slug flow condition. In addition, based on the effects of metal salt concentration, absorbent density, and gas-liquid phase flow rate on the CO2 volumetric mass transfer coefficient kLa, a new correlation for predicting the volumetric mass transfer coefficient of CO2 absorption into the hybrid solvent in microchannel reactor was proposed and validated.

Exploration of the coupled lattice Boltzmann model based on a multiphase field model: A study of the solid–liquid–gas interaction mechanism in the solidification process

A multiphase field model coupled with a lattice Boltzmann (PF-LBM) model is proposed to simulate the distribution mechanism of bubbles and solutes at the solid-liquid interface, the interaction between dendrites and bubbles, and the effects of different temperatures, anisotropic strengths and tilting angles on the solidified organization of the SCN-0.24wt.% butanedinitrile alloy during the solidification process. The model adopts a multiphase field model to simulate the growth of dendrites, calculates the growth motions of dendrites based on the interfacial solute equilibrium; and adopts a lattice Boltzmann model (LBM) based on the Shan–Chen multiphase flow to simulate the growth and motions of bubbles in the liquid phase, which includes the interaction between solid–liquid–gas phases. The simulation results show that during the directional growth of columnar dendrites, bubbles first precipitate out slowly at the very bottom of the dendrites, and then rise up due to the different solid–liquid densities and pressure differences. The bubbles will interact with the dendrite in the process of flow migration, such as extrusion, overflow, fusion and disappearance. In the case of wide gaps in the dendrite channels, bubbles will fuse to form larger irregular bubbles, and in the case of dense channels, bubbles will deform due to the extrusion of dendrites. In the simulated region, as the dendrites converge and diverge, the bubbles precipitate out of the dendrites by compression and diffusion, which also causes physical phenomena such as fusion and spillage of the bubbles. These results reveal the physical mechanisms of bubble nucleation, growth and kinematic evolution during solidification and interaction with dendrite growth.

Numerical Study of the Influence of the Type of Gas on Drag Reduction by Microbubble Injection

In this work, a novel numerical method for studying the influence of gas types on drag reduction by microbubble injection is presented. Aimed at the microbubble drag reduction (MBDR) process for different types of gases, the mass transfer velocity of different types of gases in the gas–liquid phase is defined by writing a user-defined function (UDF), which reflected the influence of gas solubility on the drag reduction rate. An Eulerian multiphase flow model and the Realizable k−ε turbulence model are used for numerical calculation. The population balance model is used to describe the coalescence and breakup phenomena of the microbubble groups. Henry’s theorem is used to calculate the equilibrium concentration of the microbubble mixed flow. The interphase mass transfer rate of the microbubble injection process for different types of gases is studied by using permeation theory. The local mass fraction of the mixed flow is solved by the component transport equation. It is found that the larger the solubility of the gas, the lower the efficiency of MBDR. When the volume flow rate of the same type of gas is the same but the injection speed is different, the larger the solubility of the gas is, the greater the difference in the drag reduction ratio.

DIRECT NUMERICAL SIMULATIONS OF TWO-PHASE FLUIDS INTERFACE IN 2D ROTATING DRUMS USING A COUPLED VOF–IBM NUMERICAL APPROACH

Rotating drums play important roles in numerous industrial applications, such as mineral processing. This work is focused on the numerical study of the interface evolution in liquid-liquid and liquid-gas phase rotating drums. A new coupling strategy between volume of fluid (VOF) and immersed boundary method (IBM) approaches is developed. Relevant dimensionless numbers, including Reynolds, Froude, and Bond numbers, alongside viscosity and density ratios, are considered for the flow pattern characterization. Direct numerical simulations are performed in order to explore flow regimes within the rotating drum, addressing a gap in the literature concerning less-explored flow patterns, particularly in the rotating drum containing liquid-liquid phases. The flow pattern families characterizing rotating drums carrying liquid-liquid phases found in this study are (i) gravity stratified, (ii) mixing, (iii) annular, and (iv) rotation stratified flows. Additionally, the characteristic flow pattern families, (i) gravity stratified, (ii) pool, (iii) annular with pool, and (iv) annular flows, are identified in rotating drums carrying liquid-gas phases. The difference in the transitory responses between the rotating drum featuring liquid-liquid and liquid-gas phases is also shown and discussed. The main results highlight significant contributions for understanding the dynamics of rotating drums, particularly concerning the transitional interface development. By identifying new flow patterns and exploring transitional phenomena, this study enriches the understanding of complex fluid behavior within rotating drum configurations.

Ionic liquid–solvent systems in the gas–liquid phase for acetylene hydrochlorination

Ionic liquid–solvent system is an excellent homogeneous catalyst for acetylene hydrochlorination due to its low viscosity, low cost and high activity.

Structure-dependent spin-polarized electron transport in twin-crystal Cu1-xEuxO semiconductors.

In this work, Eu-doped twin copper oxide (twin Cu1-xEuxO) was synthesized using the gas-liquid phase chemical deposition method in combination with high-temperature oxidation. The incorporation of Eu3+ ions was affected by their diffusivity and the related charge trapping mechanisms. The twin Cu1-xEuxO configuration exhibited significant room-temperature ferromagnetism. From our analysis, it was demonstrated that as the Eu3+ doping concentration increased, the saturation magnetization first increased and then gradually decreased, reaching a peak at 0.82 at%. A p-type to an n-type semiconducting transition was also recorded as the doping concentration increased. A significant anomalous Hall effect characterized by a maximum anomalous Hall coefficient of 1.65, and a maximum Hall conductivity mobility of 16.50 Ohm-1 cm-1 and 250.59 cm2 v-1 s-1, respectively, were derived for the twin Cu1-xEuxO, doped with 0.82 at% at room temperature. First-principles computational simulations were also conducted to elucidate the underlying mechanisms of the magnetic properties, the p-type to n-type transition, and the interplay between the spin-polarized states associated with 4f and carriers. In twin Cu1-xEuxO, the anomalous Hall effect originated from the contribution of the edge-to-jump scattering mechanism. The latter can be significantly enhanced by doping with Eu atoms, which yields the manifestation of the oblique scattering mechanism. Our work paves the way for the development of twin Cu1-xEuxO material structures, which emerge as an ideal candidate for future spintronic applications.

Odd Response-Induced Phase Separation of Active Spinners.

Due to the breaking of time-reversal and parity symmetries and the presence of non-conservative microscopic interactions, active spinner fluids and solids respectively exhibit nondissipative odd viscosity and nonstorage odd elasticity, engendering phenomena unattainable in traditional passive or active systems. Here, we study the effects of odd viscosity and elasticity on phase behaviors of active spinner systems. We find the spinner fluid under a simple shear experiences an anisotropic gas-liquid phase separation driven by the odd-viscosity stress. This phase separation exhibits equilibrium-like behavior, with both binodal-like and spinodal curves and critical point. However, the formed dense liquid phase is unstable, since the odd elasticity instantly takes over the odd viscosity to condense the liquid into a solid-like phase. The unusual phase behavior essentially arises from the competition between thermal fluctuations and the odd response-induced effective attraction. Our results demonstrate that the cooperation of odd viscosity and elasticity can lead to exotic phase behavior, revealing their fundamental roles in phase transition.

Model analysis of flow-induced vibration (LIV) in submarine oil and gas pipeline

Abstract In order to realistically simulate the flow-induced vibration of submarine oil and gas pipelines, this paper assembles a gas-liquid two-phase flow vibration device, develops a detailed and comprehensive experimental procedure, and sets the flow aperture, flow pressure, flow depth, flow direction, and gas-liquid phase apparent flow rate and other variables. According to the different types of pipelines, the multiphase flow-induced vibration of submarine oil and gas pipelines can be categorized into straight pipe vibration and elbow pipe vibration, and combined with the theoretical knowledge of fluid dynamics and the equations of motion to construct a multiphase flow-induced vibration model of submarine oil and gas pipelines and discuss the vibration mechanism of submarine oil and gas pipelines in two-phase liquid flow. Combined with the corresponding experimental parameters, the vibration phenomenon induced by gas-liquid two-phase flow in submarine oil and gas pipelines is experimentally analyzed. It is analyzed that the first-order intrinsic frequency of the submarine oil and gas pipeline is kept at 23.2 Hz when the value of the Reynolds number of the in-pipe flow ranges from 2.53 × 104 to 1.08 × 107. The first-order intrinsic frequency of the submarine oil and gas pipeline monotonically decreases with the increase of the Reynolds number of the in-pipe flow. In addition, when Ul=0.4m/s, the average standard deviation of the data is only 1.114m/s² with the rise of Ug, indicating that the vibration intensity of the fluid on the subsea oil and gas pipeline increases gradually with the increase of the apparent flow velocity of the liquid phase. By analyzing the multiphase flow problem of oil and gas pipelines, this study is of great significance in improving the efficiency of oil and gas transportation.

Growth mechanism under the supply-limited regime in mist chemical vapor deposition: presumption of mist droplet state in high-temperature field

In mist CVD, investigation of the dependence of the growth rate on the solution concentration and precursor supply amount in the growth of ZnO films using a methanol-based solution of Zn(acac)2 revealed that the reaction is a rate-determination process of the supply-limited regime. Based on this result, the diffusion coefficient of gasified Zn(acac)2 at standard ambient temperature and pressure (SATP) was calculated to be approximately 10−6m2/s, where the diffusion rate was derived from the relation between the precursor supply amount and the film formation amount, assuming that all precursors were gasified in the reactor. The derived values were approximately one order of magnitude lower than the typical diffusion coefficient of 10−5m2/s in the gas phase at SATP. Several factors considered were unable to explain the much lower diffusion coefficient. Therefore, this result contradicts the assumption that all precursors are gasified in the reactor, and it is presumed that droplets might remain in the liquid–gas phase.

Study on the effect of different electrode forms on the surface modification of Ti–6Al–4V alloy by near-dry electrical discharge machining

The surface modification of Ti–6Al–4V alloy was achieved by using an electrical discharge machining (EDM) surface modification technology under the action of pulsed discharge to induce an in-situ synthesis reaction between the molten metallic material and the ionized elements in the atomized medium. In this study, nitrogen and urea solutions are used as the gas-liquid phases of the atomization medium. ANSYS simulation software was used to simulate and analyze the interpolar pressure and flow field of the different electrodes. The effect of different electrode forms on the EDM surface modification layer was investigated by designing comparative trials. Scanning electron microscopy (SEM), X-ray diffraction (XRD), and other detection methods were used to characterize the process indicators, such as surface morphology, roughness, hardness, and wear resistance of the reinforced layer. It was found that a thick and uniform continuous layer can be obtained on the Ti–6Al–4V alloy substrate using multi-channel electrode machining, and the thickness is about 65 μm. XRD analysis indicates that the reinforced layer mainly contains TiN, Ti2N, Ti2AlN, Ti2O, and CuO phases. The samples with reinforced layers were shown to have higher hardness (1578 HV) and better wear resistance than the titanium substrates. The main wear mechanism is spalling wear and light oxidation wear.

Phase separation of a magnetic fluid: Asymptotic states and nonequilibrium kinetics.

We study self-assembly in a colloidal suspension of magnetic particles by performing comprehensive molecular dynamics simulations of the Stockmayer (SM) model, which comprises spherical particles decorated by a magnetic moment. The SM potential incorporates dipole-dipole interactions along with the usual Lennard-Jones interaction and exhibits a gas-liquid phase coexistence observed experimentally in magnetic fluids. When this system is quenched from the high-temperature homogeneous phase to the coexistence region, the nonequilibrium evolution to the condensed phase proceeds with the development of spatial as well as magnetic order. We observe density-dependent coarsening mechanisms-a diffusive growth law ℓ(t)∼t^{1/3} in the nucleation regime and hydrodynamics-driven inertial growth law ℓ(t)∼t^{2/3} in the spinodal regimes. [ℓ(t) is the average size of the condensate at time t after the quench.] While the spatial growth is governed by the expected conserved order parameter dynamics, the growth of magnetic order in the spinodal regime exhibits unexpected nonconserved dynamics. The asymptotic morphologies have density-dependent shapes which typically include the isotropic sphere and spherical bubble morphologies in the nucleation region, and the anisotropic cylinder, planar slab, cylindrical bubble morphologies in the spinodal region. The structures are robust and nonvolatile, and exhibit characteristic magnetic properties. For example, the oppositely magnetized hemispheres in the spherical morphology impart the characteristics of a Janus particle to it. The observed structures have versatile applications in catalysis, drug delivery systems, memory devices, and magnetic photonic crystals, to name a few.

FRACTAL MODEL FOR EFFECTIVE THERMAL CONDUCTIVITY OF COMPOSITE MATERIALS EMBEDDED WITH A DAMAGED TREE-LIKE BIFURCATION NETWORK

Scientists worldwide have always been interested in the study of networks that resemble trees and porous materials. Therefore, based on fractal theory, this paper systematically studies the heat transfer problem of the damaged tree-like networks under different saturations in multi-medium composite materials and derives their dimensionless thermal conductivity (DTC). The research has shown that the dimensionless thermal conductivity (DTC) decreases with an increase in damaged channels. The saturation significantly impacts the heat transfer of the damaged tree networks, and it exhibits linear changes with the ratio of the gas phase, liquid phase, and solid phase. Additionally, it can be observed that the fractal dimension of length distribution and diameter distribution, bifurcation numbers, bifurcation series, porosity, and other parameters in the networks are related to thermal conductivity. The comparison between the data derived from this model and experimental data shows that the proposed model can effectively deepen our understanding of the heat transfer mechanism of the damaged networks in composite materials under different saturation levels. Additionally, the model in this paper does not have an empirical constant, which avoids the influence of potential factors.