Multi-field Coupling Research Articles

Experimental study on the influence of roadway shape on the evolution of outburst fluid static pressure

To explore the static pressure dynamic disaster mechanism of coal-and-gas outburst (CGO) fluid, the self-developed multi-field coupling large-scale physical simulation test system of coal mine dynamic disaster was used to carry out gas outburst and CGO physical simulation tests in straight, L-shaped and T-shaped roadways. The influence of roadway shape on the evolution of static pressure was explored, and the role of pulverized coal in the process of static pressure dynamic disaster was clarified. The results indicated that the static pressure showed a fluctuating downward trend during the outburst process. When gas outburst, the middle and front parts of the roadway in the straight section roadway were the most serious areas of static pressure disasters in the three shapes of roadways. The duration and range of high static pressure disaster in L-shaped roadway were larger than those in T-shaped and straight roadways in turn. When CGO, the most serious area of static pressure disaster in L-shaped and T-shaped roadways moved backward to the middle of the straight section roadway, and there was a rebound phenomenon in the process of static pressure fluctuation decline, which showed the pulse characteristics of CGO. During the outburst, the static pressure dynamic disaster hazard of L-shaped roadway was higher than that of T-shaped roadway, and the static pressure at the bifurcation structure decayed faster than that at the turning structure, which indicated that T-shaped roadway was more conducive to the release of static pressure in roadway, thus reduced the risk of static pressure disaster. When gas outburst, the static pressure attenuation of the fluid in the roadway before and after the turning and bifurcation structure was greater than that of CGO. The peak static pressure and impulse of the fluid during gas outburst were 2 times and 4–5 times that of CGO respectively. The presence of pulverized coal reduced the attenuation of static pressure and the hazard of dynamic disaster, prolonged the release time of energy, and led to the change of the maximum static pressure disaster area.

Quantitative assessment method of thermal cycle accumulation mechanism during the pulsed laser single-channel multilayer cladding process

Pulsed laser cladding provides a periodical thermal accumulation, enabling a greater temperature gradient and cooling rate inside the substrate, which can efficiently ameliorate the microstructure of the cladding layer. It is crucial to quantitatively investigate the mechanism of the interaction for the cladding process parameters on the transient evolution of the multi-field coupling to improve the cladding quality. In this paper, a 3-D model during the single-channel multilayer pulsed laser cladding process was constructed which is grounded on the temperature-variable physical properties parameters of the material. The integrated consideration in the masking impact of the molten powder waist beam on the pulsed laser beam, the molten pool liquid metal surface tension, the influence of buoyancy on the flow, and the transient evolution of the molten pool edge morphology. The temperature, flow, and stress fields of the single-channel multilayer cladding process were computed and resolved. The mechanism in multi-field coupling impact of laser pulse frequency, laser power, and pulse duty cycle on the single-channel multilayer cladding process was emphasized. The response surface equations between each influencing factor and the cladding temperature, temperature change gradient, and thermal stress were determined according to the response surface method, and the sensitivity of each parameter was attained from the Monte Carlo method. Using SEM to observe the profile and microstructure of the cladding layer, the material hardness distribution was measured experimentally, and the effectiveness of the model was confirmed. This research supplies an essential rationale for upgrading the quality of pulsed laser cladding layers.

The effect of friction characteristics on electromagnetic impact deformation mechanism of 2A10 aluminum alloy bars

Purpose The purpose objective of this study is to identify the friction coefficient and friction effect in electromagnetic upsetting (EMU) high-speed forming process. Design/methodology/approach Based on numerical simulation and upsetting experiment of 2A10 aluminum alloy bar, the friction coefficient between contact surfaces is obtained by combining the fitting displacement distribution function and the electromagnetic-mechanical coupling numerical model, and the influence of friction effect is analyzed. Findings The maximum impact velocity and acceleration during EMU are 13.9 m/s and −3.3 × 106 m/s2, respectively, and the maximum strain rate is 7700 s−1. The functional distribution relationship between friction coefficient combination (FS, FD) and characteristic parameters [upper diameter (D1) and middle diameter (D2)] is established. The values of FS and FD are 0.1402 and 0.0931, respectively, and the maximum relative error is 2.39%. By analyzing the distribution of equivalent stress and strain, it is found that plastic deformation has obvious zoning characteristics and there is serious failure concentration in the strong shear zone. Originality/value Friction coefficient significantly affects stress or strain distributions in material forming process, but it is difficult to obtain friction coefficients through experimental tests in the high-speed forming process. In this paper, a multi-field coupling numerical model is proposed to determine friction coefficients and applied to the electromagnetic impact loading process (a high-speed forming process). Peer review The peer review history for this article is available at: https://publons.com/publon/10.1108/ILT-05-2024-0154/

Analytical solution on magneto-electro-elasto-hygrothermal coupling in rotating functionally graded piezoelectric/piezomagnetic cylinders on Winkler elastic foundation

This study focuses on the multi-field coupling phenomena in a rotating functionally graded piezoelectric/piezomagnetic (FGPEPM) hollow cylinder under magneto-electro-elasto-hygrothermal loading. The cylinder is assumed to be infinitely long and on a Winkler elastic foundation. The hygrothermo parameter, magneto-electro-elastic properties, hygrothermal expansion, and density within the FGPEPM cylinder follow power laws with different gradient indexes along its thickness. Analytical solutions for the distributions of radial displacement, stresses, electric potential, magnetic potential, temperature, and moisture are derived by solving magneto-electro-elasto-hygrothermal coupling differential equations with appropriate boundary conditions. Validation against the existing analytical results for simplified models confirms the effectiveness of the present analytical solution. Numerical analyses further explore the effects of gradient indexes, rotational velocity, and elastic foundation stiffness on the multi-field coupling behavior of the FGPEPM cylinder. The findings provide valuable insights for designing multi-field coupling characteristics for FGPEPM materials.

Ultra-low stress-hysteresis and huge superelasticity in NiMn-based shape memory microwire

Hysteresis related to first-order phase transformation in shape memory alloys, which is the macroscopic manifestation of energy dissipation, is detrimental to the precise control of actuation and causes structural and functional fatigue of components. It is of vital importance to explore high-performance shape memory alloys with low stress-hysteresis, large superelasticity, and wide temperature range operation in practical applications. Here, we have developed a Ni-Mn-Fe-In shape memory microwire, exhibiting an ultra-low stress-hysteresis and huge tensile superelasticity in a wide temperature range. The microwire shows a smooth surface and a single crystal structure (with ⟨001⟩A-oriented along the axial direction of microwire), and the microstructure of the microwire contains austenite matrix and sparsely distributed precipitates with an average size of 20–80 nm, all of which may be beneficial to obtain low hysteresis and large strains in the microwire. As a result, the microwire exhibits a minimum stress-hysteresis of as low as 8.5 MPa (with overall strain of 15.3%) and corresponding energy dissipation as low as 1.44 MJ/m3. The microwire always shows a low stress-hysteresis (less than 24 MPa) and low energy dissipation (less than 2.86 MJ/m3) above room temperature. The microwire shows a huge superelasticity with recoverable strains higher than 15% in the wide temperature range from 218 to 418 K. Together with the advantages of easy fabrication and no post-processing required, this microwire shows a tremendous potential for cyclic actuators and energy conversion devices under multi-field coupling.

Electricity-driven dealkalization of bauxite residue based on thermodynamics, kinetics, and mineral transformation.

High alkalinity content of bauxite residue is a major factor that hinders resource reutilization and pollutes the environment. Although acid neutralization is a direct and effective method, the amount of acid and secondary waste of sodium salt are still difficult problems to solve. Herein, we innovatively integrated an electric field into the acid neutralization dealkalization of bauxite residue and analyzed the dealkalization behavior by thermodynamics, kinetics, and mineral transformation. The results show that the pH of the anode chamber was maintained at the acidic levels of 3-6 after 30min of galvanostatic electrolysis, and bauxite residue can realize dealkalization by acid neutralization. In the anode chamber, Na+ was released into the leachate via the reactions of Na3Al3Si3O12 and the removal of encapsulated soluble alkali. The stainless steel wire mesh anode exhibited its superiority and decreased the Na2O content in bauxite residue from 9.48 to 3.13% through convective mass transfer driven by the electric field and steady-state diffusion under stirring. This research provides a promising method for the electricity-driven dealkalization of bauxite residue, thus facilitating the development of multifield coupling theory and the application of electric fields in the alumina industry.

Free vibration of functionally-graded piezoelectric semiconductor rectangular beam under thermal load

In this paper, the free vibration of a functionally-graded piezoelectric semiconductor (FGPS) rectangular beam under thermal load is studied within the framework of the theory of Timoshenko beam and the theory of piezoelectric semiconductor. It is assumed that the material properties in the functional gradient layer change smoothly. The multi-field coupling in the beam is deduced by introducing the first-order shear deformation theory in this paper, and the governing equations including motion equation, Gauss law and current continuity condition are obtained. The vibration of a simply supported FGPS rectangular beam is realized by solving the governing equations. In numerical examples, the natural frequencies of the FGPS rectangular beam under thermal load are analyzed through different perspectives. The results show that for different length scales, the natural frequency of the FGPS beam under thermal load is reduced by increasing the initial electron concentration in a certain range. The damping characteristics are insensitive to temperature change, and the damping characteristics are decreased by increasing length to thickness ratio and the thickness of the FGPS layer.

Shrinkage separation prediction of concrete-filled steel tube arch bridge: A coupling model concerning multiple factors

The shrinkage separation is one of the key factors restricting the development of concrete-filled steel tube (CFST) arch bridge. In this paper, its mechanism and evaluation method concerning multiple factors, including structure, environment, material and construction were investigated. Firstly, the hydration process of in-tube concrete was simulated using the hydration degree as a basic state parameter. Then with the consideration of energy and mass conservation, as well as the moisture and heat transmission effects, the temperature and relative humidity fields, early-age mechanical properties, basic creep, and total deformation of in-tube concrete were further calculated. On the basis, in accordance with the stress criterion and displacement coordination, a multi-field (hydro-thermo-hygro-constraint) coupling model was finally established, taking the ratio of normal stress to bonding strength at the steel-concrete interface as the shrinkage separation risk index. The simulation results obtained from the model were in good agreements with the full-scale CFST model test.

Thermal dynamics in deep shale reservoirs: Influences of the kerogen microstructural behavior on the gas adsorption/desorption capacity

In shale gas extraction projects, an investigation into the mechanisms of energy/mass transfer associated with shale gas adsorption/desorption in organic matter (kerogen) microstructure under high temperature and stress condition is crucial for improving the efficiency of shale gas production. This study presents a coupling thermo-hydro-mechanical model based on an improved fractal method that could explain the microstructural evolution of the kerogen system and the resultant alterations during the gas adsorption/desorption process under varying thermal conduction, gas seepage, and stress conditions. The influence of porosity, diameter, and tortuosity on the abundance, length, and complexity of kerogen networks under coupled multi-field effects is evaluated. The significance of this study is it could address the following aspects quantitively: (1) the spatiotemporal evolution of kerogen fractal dimensions following various extraction timelines; (2) the influence of shale temperatures on kerogen structures; (3) the influence of the kerogen fractal dimension on the shale gas desorption capacity and production efficiency; and (4) under different temperatures, when the fractal dimension/tortuosity fractal dimension of kerogen changes due to extraction disturbances, the volumetric deformation induced by gas adsorption increases by a maximum of 26.1%/decreases by 28.1% and in the later stages of extraction, the maximum gas pressure decreases by 44.7%/increases by 47.1%. The proposed fractal method adeptly reveals shale gas desorption behaviors under multi-field coupling conditions from a microscopic perspective, which cannot be found in the literature.

Analysis of methane diffusion on permeability rebound and recovery in coal reservoirs: Implications for deep coalbed methane-enhanced extraction

A proper understanding of the effect of methane diffusion on coal reservoir permeability rebound and recovery is essential, as coal reservoir permeability is the key parameter influencing the efficiency of coalbed methane migration and computational research on it is lacking. In this paper, the multifield coupling model for methane migration was established. Then, two parameters, the influence coefficient of diffusion on permeability rebound (DPRB) and the influence coefficient of diffusion on permeability recovery (DPRC), were proposed to quantify the effect of methane diffusion on rebound and recovery of coal reservoir permeability. Subsequently, we used COMSOL software to study the variation rules of the coal reservoir permeability rebound time, permeability recovery time, and permeability rebound value, DPRB, and DPRC for different geologic parameters. The results shown that the permeability rebound time and recovery time are proportional to the coal seam initial pressure, but inversely proportional to the initial permeability and initial diffusion coefficient. The rebound value decreases with increasing coal seam initial pressure and initial permeability, but ascends with rising initial diffusion coefficient. DPRB declines with increasing coal seam initial pressure, initial permeability, and initial diffusion coefficient, but they are all greater than 0.7, indicating that methane diffusion has a significant effect on permeability rebound. The DPRC values for different coal seam initial pressures, initial permeabilities, and initial diffusion coefficients are above 0.98, which implies that methane diffusion dominates the permeability recovery process. Finally, a conceptual model was presented to research the mechanism of diffusion influence on rebound and recovery of coal reservoir permeability, and the implications for enhanced drainage of deep coalbed methane were discussed. Therefore, the results of this paper can provide a theoretical foundation for deep coalbed methane-enhanced extraction.

Damage characterization of high- and low-temperature performance of porous asphalt mixtures under multi-field coupling

Porous asphalt (PA) have attracted considerable interest due to their functional advantages such as significant water permeability and noise attenuation properties. Nevertheless, its vulnerability to damage significantly hinders its widespread application. Therefore, this paper aims to investigate the damage characteristics of PA's high- and low-temperature performance under the coupled effects of moisture and temperature. Firstly, Hamburg Wheel-Tracking Test (HWTT) and three-point bending tests were conducted on PA under different coupling conditions. Subsequently, its damage characteristics were analyzed using the three-stage permanent deformation model, logistic model, and two-way analysis of variance (ANOVA) method. The findings indicate that under dry conditions, PA-13 exhibits a linear decrease in rutting resistance within the temperature range of 40–60 °C, followed by a significant decrease at 70 °C. During immersion in water, PA-13 undergoes a second and third-stage rutting evolution process, with an accelerated creep rate observed at 70°C. High temperatures exacerbate the viscoelastic damage of PA-13, while multi-field coupling promotes the simultaneous occurrence of viscoelastic damage and stripping damage. Furthemore, the results of the two-way ANOVA indicate that moisture has a more significant impact on high-temperature performance damage. In addition, in a dry environment, the low-temperature crack resistance of PA-13 deteriorates with decreasing temperature, but there is little difference at −20 °C and −30 °C. This is attributed to PA-13 exhibiting flexible cracking at 0 °C and −10 °C, while undergoing brittle cracking at −20 °C and −30 °C. Additionally, the coupled effect of temperature and moisture reduces the damage tolerance of PA-13 and accelerates its cracking process. The results of the two-way ANOVA indicate that both low temperature and treatment methods have a statistically significant impact on low-temperature performance, with treatment methods having a higher degree of influence. The findings of this paper contribute to an in-depth understanding of the damage mechanism of PA-13.

Automated Shape Correction for Wood Composites in Continuous Pressing

The effective and comprehensive utilization of forest resources has become the theme of the global “dual-carbon strategy”. Forestry restructured wood is a kind of wood-based panel made of wood-based fiber composite material by high-temperature and high-pressure restructuring–molding, and has become an important material in the field of construction, furniture manufacturing, as well as derivative processing for its excellent physical and mechanical properties, decorative properties, and processing performance. Taking Medium Density Fiberboard (MDF) as the recombinant material as the research object, an event-triggered synergetic control mechanism based on interventional three-way decision making is proposed for the viscoelastic multi-field coupling-distributed agile control of the “fixed thickness section” in the MDF continuous flat-pressing process, where some typical quality control problems of complex plate shape deviations including thickness, slope, depression, and bump tend to occur. Firstly, the idea of constructing the industrial event information of continuous hot pressing based on information granulation is proposed, and the information granulation model of the viscoelastic plate shape process mechanism is established by combining the multi-field coupling effect. Secondly, an FMEA-based cyber granular method for diagnosing and controlling the plate thickness diagnosis and control failure information expression of continuous flat pressing is proposed for the problems of plate thickness control failure and plate thickness deviation defect elimination that are prone to occur in the continuous flat-pressing process. The precise control of the plate thickness in the production process is realized based on event-triggered control to achieve the intelligent identification and processing of the various types of faults. The application test is conducted in the international mainstream production line of a certain type of continuous hot-pressing equipment for the production of 18 mm plate thickness; the synergistic effect is basically synchronized after 3 s, the control accuracy reaches 30%, and the average value of the internal bond strength is 1.40, which ensures the integrity of the slab. Practical tests show that the method in the actual production is feasible and effective, with detection and control accuracy of up to ±0.05 mm, indicating that in the production of E0- and E1-level products, the rate of superior products can reach more than 95%.

Shaking table tests and numerical analysis of monopile-supported offshore wind turbines under combined wind, wave and seismic loads

Currently, large-diameter monopile-supported offshore wind turbines (OWTs) are widely used in offshore areas of China. However, considerable numbers of wind farms are located in seismically active areas with liquefiable soil layers. Therefore, shaking table tests and numerical simulation based on a scaled model of Siemens SWT-4.0-130 OWTs were carried out to investigate the failure mechanism of the OWT system in liquefiable sand under combined wind, wave and seismic loads. The scaling law focusing on the governing physical mechanism of the OWT system considering multi-field coupling are presented and verified. The dynamic responses of the OWT and the excess pore pressure fluctuation around the soil-monopile domain in different operating states are discussed and compared. The results show that the proposed scaling laws can well restore the dynamic characteristics of the prototype. The seismic excitation is the dominant load for inducing the structural response under combined wind, wave and seismic loads. Due to the coupling effect, the dynamic responses under combined wind and wave loads were slightly less than the superposition of structural responses under separated wind or wave loads. Liquefaction of the shallow soil layer within the range of four times the pile diameter was detected when the PGA reached 0.4 g (three-dimensional input) or 0.64 g (one-dimensional input), which induced obvious tilt and settlement of the OWT model after the test. The sudden increase in the pore pressure resulted in the liquefaction of soil, which further led to a significant decrease in the natural frequencies and decay of the acceleration responses at the tower top. The above findings are expected to provide valuable experience for the 1 g underwater shaking table tests on OWT systems and insights for the practical design of OWTs in seismically active areas.

Experimental analysis of the thermo-hydro-mechanical (THM) coupling in freezing vertical shafts of unsaturated sandy soil

The evolutionary mechanism of frozen soil involves complex dynamic coupling among temperature, moisture and the stress field. However, existing research has struggled to adequately describe the interplay between these factors. To address this, we independently designed and developed a multifunctional loading system appropriate for geotechnical engineering experimentation and a corresponding loading technology. Using a model of vertical shaft freezing, we studied the spatiotemporal evolution of thermo-hydro-mechanical (THM) multi-field coupling. The research findings indicate that, compared to the freezing interface, the main section experiences not only more intense variations in the temperature field but also heightened activity in terms of in-situ freezing and moisture migration. Both the radial and circumferential characteristic faces exhibit wave-like variations in moisture gradient evolution. The circumferential face features a critical gradient at 3.8 m−1, whereas the moisture gradient curve of the radial face undergoes temporal elongation at the peak, resulting in no discernible extremities. Each characteristic position of the frost heave force growth undergoes three distinct phases: incubation, rapid increase and stabilisation. During the same phases, the response times for the growth of frost heave forces on the radial characteristic face are roughly equivalent. However, when moving outward along the equivalent freezing radius, the response time on the circumferential face becomes progressively delayed.

Study on similarity criterion for fluid-thermal coupling responses of metal plate simultaneously exposed to laser irradiation and airflow

Similarity criterion considering the multifield coupling effect is essential for the scaled thermal test under both laser irradiation and airflow. This paper establishes the similarity criterion for the fluid-thermal coupling responses of metal plate simultaneously exposed to the laser irradiation and high-speed airflow by nondimensionalizing the governing equations. The main dimensionless parameters are obtained, and the scaling law when the temperature scaling factor is equal to 1, which means the temperature remains unchanged in the real and scaled model, is firstly established. Moreover, to broaden the applicability, the scaling law when the temperature scaling factor is not equal to 1 is further deduced, and the correction method based on numerical simulation is proposed to eliminate the error caused by the difference of temperature-dependent material properties in the real and scaled model. Numerical cases are conducted and compared to verify the established scaling laws. Results show that errors of the predicted fluid-thermal responses between the real and scaled model are less than 1 %, affirming the reliability of the scaling laws.

Constructing a Microdiffusion-Seepage-Stress Multifield Coupling Model for Nanopore Gas.

Existing research is difficult to fully capture the correlation between gas molecules and pore wall interactions, multiphase flow, and stress distribution in nanopores. Taking gas as an example, a microscopic model was constructed. At the same time, diffusion, seepage, and stress were considered to accurately predict and manage gas transport in nanopores. First, molecular dynamics (MD) simulation methods were adopted to simulate the motion trajectories and interactions of gas molecules in nanopores. Second, a multiscale model was established based on continuum mechanics to consider the interaction between pore walls and gas molecules, and a diffusion equation was established to describe the diffusion process of gas molecules in pores. Then, finite element analysis and porous media models were used to simulate the seepage behavior of gas in the nanopores. Finally, the stress distribution in the pores was analyzed, and the influence of the interaction between the pore wall and gas molecules on stress was considered. The multifield coupling model was experimentally evaluated from three aspects: diffusion coefficient, seepage behavior, and stress distribution. The root-mean-square error (RMSE) and mean absolute error (MAE) of the model in different testing directions were calculated using different simulation tools, such as COMSOL, ANSYS, OpenFOAM, and CFX. The mean values of RMSE and MAE were lower than 0.20 and 0.17, respectively. The constructed model can comprehensively describe gas transmission within nanopores, improving the management accuracy and efficiency.

Investigation on fatigue crack propagation failure mechanism of hydraulic lifting pipe in deep‐ocean natural gas hydrate exploitation

AbstractIn deep‐ocean natural gas hydrate exploitation operation, the fatigue failure mechanism has attracted more and more attention from scholars, but it has not been effectively disclosed. Therefore, in this work, a multifield coupling and multiple‐nonlinear vibration model of lifting pipe is established, which can accurately determine the alternating stress of deep‐ocean lifting pipe. Second, according to Forman theory, the calculation method of crack propagation length and depth on the surface of a deep‐water riser is established, which is verified by the comparison between the experimental and theoretical model calculation results. The results demonstrate that, first, the effect of residual stress in the welded joint of deep‐ocean lifting pipe should be considered in the later parameter influence analysis. Second, the fatigue growth life of deep‐water pipe with small outflow velocity is mainly determined by tensile stress, and that is determined by both tensile stress and bending stress with large outflow velocity. Third, more attention should be paid to the vibration of the lower pipe on‐site to reduce its vibration frequency and vibration stress amplitude, which can effectively reduce the surface crack propagation state of the deep‐ocean pipe and improve the service life of the pipe. Fourth, properly adjusting the tension coefficient of the tensioner during field operation can effectively improve the safety of the pipe, and the optimal tension coefficient is related to the configuration of the deep‐ocean pipe system, which can be analyzed and determined by the model.

Multi-field coupling vibration patterns of the multiphase sink vortex and distortion recognition method

Fluid-induced vibration sensing technology plays an essential role in the safety and stable operation of fluid control equipment. The sensing works oriented to the multiphase sink vortex-induced vibration (MSVIV) face significant challenges in state recognition, such as the aero-engine oil pipe, metallurgical casting ladle, and microfluidic equipment. This paper presents a multi-field coupling vibration modelling and solving strategy to investigate the MSVIV transition patterns. Based on the volume of fluid and level set method, an efficient volume-correction strategy is presented to consider the mass variation and keep velocity fields without divergence. Then, a fluid–structure-acoustic dynamic model is conducted to explore multiphase vortex transport regularities. Considering the solution singularity of the fluid–solid-acoustic coupling system, a residue theorem-based solution strategy is presented to obtain the MSVIV evolution behaviors. Finally, a multichannel sensing MSVIV water-model experiment (MSVIV-WME) platform is conducted to verify the above results, and a wavelet transform-based processing approach can be utilized to obtain mutation energy features. Research results illustrate that the proposed strategy has revealed the MSVIV transition behaviors. The MSVIV process in critical states dominates nonlinear vibration patterns. The self-development sensing platform can monitor the MSVIV evolution process in real-time, and the wavelet transform-based sensing attribute with the dw3 frequency band can identify distortion characteristics. Relevant results can offer useful references for MSVIV state recognition and offer technical solving strategies for microfluidics monitoring equipment.

Analysis and Optimization of the Flow Path within the Protonic Ceramic Fuel Cell Stack

With the development of human society, the demand for energy is increasing, and energy, an essential global necessity, is decreasing significantly and is in danger of being depleted. Fuel cell is the fourth power generation facility after hydroelectric, thermal and nuclear power. Fuel cell flow channel is an important part of fuel cell, which has the functions of transferring reaction gas, discharging products and excess unreacted gas, and balancing fuel cell temperature. This thesis focuses on the model design and optimization of a 10-layer flat-plate proton ceramic fuel cell (PCFC) reactor. First, we design a 10-layer flat-plate PCFC reactor model structure, then we simulate the reactor using Fluent, analyze the multi-field coupling operating characteristics inside the target reactor, optimize the reactor based on the results, and finally, we analyze and discuss the optimization results.

Enabling multi-stage high-temperature strength evolution prediction of ceramizable composites using a novel multi-field coupled model

Strength varies significantly under high-temperature environment, due to the inherent thermomechanical behavior of the ceramizable material and its coupling with possible chemical reactions. The complexity amplifies for composite materials, considering their multi-phase and multi-scale features, and more importantly, their complicated chemical reactions under high-temperature service conditions. This study proposes an innovative multi-field coupling theory framework for predicting the multi-stage evolution behavior of high-temperature mechanical properties of a ceramizable composite, through incorporating an extended chemical kinetics method, coupled deformation, mass diffusion and heat conduction. The developed model enables direct coupling and simultaneous solving of physical, chemical and thermal variables. It captures well the degradation of mechanical properties for the initial stage and the increase of strength for the later stage, along with the increasing of temperature. The validated model also enables well prediction of time-dependent mechanical properties at high service temperature, with an average error of 8.67% against experimental measured results. The developed method can serve as a general method for the prediction of high-temperature mechanical property of thermal protection composites and structures.