Method, device, processor and storage medium for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench
By calculating the dimensionless number of fluids and the average diameter of bubbles, combined with the gas phase humidification performance prediction model, the problem of humidification range design of bubble humidification devices under large flow rates is solved, and rapid performance prediction and accuracy improvement in the early stage of the design is achieved.
Patent Information
- Application Number
- CN202311816222.6
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-12-27
- Publication Date
- 2025-02-07
- Estimated Expiration
- 2043-12-27
AI Technical Summary
The existing bubble humidification device is difficult to meet the design requirements of the gas humidification range under large flow rates, resulting in a large amount of experimental verification required during the design stage, which increases material and time costs.
By calculating the dimensionless number of fluids, predicting the average diameter of the bubbles, and combining the gas phase humidification performance prediction model, the gas phase outlet humidity is calculated, and the design parameters are adjusted to meet the design goals.
It realizes rapid estimate of the performance of humidification device in the early stage of design, shortens debugging time, and improves design accuracy and efficiency.
Smart Images

Figure CN117790846B_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the field of hydrogen fuel cell technology, and in particular to the field of a bubbling humidification device for a fuel cell test bench, and specifically refers to a method, device, processor and computer-readable storage medium thereof for implementing performance prediction processing for a bubbling humidification device for a fuel cell test bench. Background Art
[0002] Proton exchange membrane fuel cells are generally composed of bipolar plates, sealing gaskets and membrane electrode assemblies, in which the membrane electrode assembly is composed of a proton exchange membrane, a catalyst layer and a gas diffusion layer. During the operation of the proton exchange membrane fuel cell, in order to ensure good conductivity of protons, it is generally necessary to humidify the reaction gas to ensure the appropriate humidity of the proton exchange membrane. If the proton exchange membrane is too dry, the activity of the catalyst layer covering the surface of the proton exchange membrane will decrease, and it may even cause the membrane electrode assembly to rupture, causing the anode hydrogen and cathode oxygen to mix with each other, resulting in combustion and explosion. If the proton exchange membrane is too moist, it will cause serious water accumulation inside the fuel cell, which will hinder the flow of gas and weaken the activity of the catalyst.
[0003] Proton exchange membrane fuel cells with different membrane electrode materials, different flow channel structures and different output powers require different intake air flow rates, intake air temperatures and intake air humidity to ensure their normal operation. Especially for high-power and high-stacked proton exchange membrane fuel cells, it is necessary to ensure the widest possible humidification range at a large gas flow rate.
[0004] At present, the gas humidification devices of fuel cell test benches are generally divided into bubbling type and spray type. The basic principle of the bubbling humidification device is that the dry gas passes through the small holes distributed at the bottom of the liquid pool under a certain gas pressure to obtain a large number of bubbles dispersed in the liquid pool, so that the gas phase and the liquid phase can be in direct contact, thereby obtaining better mass transfer and heat transfer effects. At present, for the bubbling humidification device, the gas humidification range under large flow needs to be estimated through experience and ultimately verified by a large number of experiments. If the humidification device cannot meet the design requirements, it is necessary to re-adjust the existing structure based on experience and test it again, and the required material cost and time cost are high. Therefore, for the bubbling humidification device, a humidification performance prediction method is urgently needed. This method can be used to quickly estimate the performance of the humidification equipment and optimize the design parameters during the design stage, so as to minimize the repeated debugging process after the humidification device is manufactured. Summary of the invention
[0005] The purpose of the present invention is to overcome the shortcomings of the above-mentioned prior art and provide a method, device, processor and computer-readable storage medium for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench, which has the requirements of short debugging time, high accuracy and wide application range.
[0006] In order to achieve the above-mentioned purpose, the method, device, processor and computer-readable storage medium thereof for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench of the present invention are as follows:
[0007] The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench is mainly characterized in that the method comprises the following steps:
[0008] (1) Calculate the dimensionless number of the fluid according to the structural parameters and physical properties of the bubbling device;
[0009] (2) According to the dimensionless number, the average bubble diameter d at the corresponding airflow velocity is calculated through the bubble average diameter prediction model. b ;
[0010] (3) According to the average bubble diameter d b , and input the gas phase inlet humidity H1, liquid pool liquid phase temperature T l , gas saturation humidity H corresponding to liquid phase temperature s and area mass transfer coefficient k c , the gas phase outlet humidity H2 is obtained according to the gas phase humidification performance prediction model;
[0011] (4) Determine whether the gas phase outlet humidity H2 meets the design target. If so, the design is completed and the step ends; otherwise, readjust the design parameters and continue with step (1).
[0012] Preferably, the dimensionless number of the fluid calculated in step (1) includes the Bond number Bo, and the Bond number Bo is calculated in step (1) as follows:
[0013] The Bond number Bo is calculated according to the following formula:
[0014]
[0015] Where g is the acceleration of gravity, N is the number of pores, d0 is the pore diameter, ρ l is the liquid density, σ l is the surface tension of the liquid.
[0016] Preferably, the dimensionless number of the fluid calculated in step (1) includes the Galileo number Ga. The Galileo number Ga calculated in step (1) is specifically:
[0017] The Galilean number Ga is calculated according to the following formula:
[0018]
[0019] Where g is the acceleration of gravity, N is the number of pores, d0 is the pore diameter, and v lis the dynamic viscosity of the liquid.
[0020] Preferably, the dimensionless number of the fluid calculated in step (1) includes the Froude number Fr. The Froude number Fr calculated in step (1) is specifically:
[0021] The Froude number Fr is calculated according to the following formula:
[0022]
[0023] Among them, u0 is the pore flow velocity, g is the gravitational acceleration, N is the number of pores, and d0 is the pore diameter.
[0024] Preferably, the bubble average diameter prediction model in step (2) is specifically:
[0025] The bubble average diameter prediction model is expressed according to the following formula:
[0026]
[0027] Among them, k, a, b, c are coefficients and exponents that can be corrected, Bo is the Bond number, Ga is the Galileo number, Fr is the Froude number, N is the number of pores, and d0 is the pore diameter.
[0028] Preferably, the method further comprises the step of adjusting the coefficients and exponents k, a, b, c according to computational fluid dynamics simulation results, specifically comprising the following steps:
[0029] (1-1) Establish a corresponding humidification equipment CAD model according to the required tank cavity;
[0030] (1-2) Discretize the fluid region inside the tank;
[0031] (1-3) Select a suitable calculation model and input the boundary conditions of the humidification device.
[0032] (1-4) Delineate the liquid area and gas area inside the tank according to the liquid distribution of the tank, and perform initialization and solution;
[0033] (1-5) Post-processing is performed on the calculation results, including extracting the size of each bubble in the liquid pool to obtain the average diameter and diameter distribution of the bubbles.
[0034] Preferably, the gas phase humidification performance prediction model in step (3) is specifically:
[0035] The gas phase humidification performance prediction model is expressed according to the following formula:
[0036]
[0037] Among them, kc is the mass transfer coefficient, d b is the average bubble diameter, H1 is the gas phase inlet humidity, H s is the saturated humidity of the gas at the liquid phase temperature.
[0038] Preferably, the method further comprises the step of determining the mass transfer coefficient k c The steps of fitting by bubbling experiment specifically include the following steps:
[0039] (2-1) delivering liquid at a specified temperature T to the bubbling body through a liquid delivery system, so that a liquid pool with a height of h is formed in the bubbling area;
[0040] (2-2) Opening the gas transmission channel, delivering dry air to the gas distribution plate at the bottom of the bubbling body through the air compressor, and forming bubbles through the gas distribution plate; recording the flow rate Q and air pressure of the input gas in the bubbling area through the flow sensor, air pressure sensor and humidity sensor arranged at the outlet of the air compressor;
[0041] (2-3) Turn on the humidity sensor and air pressure sensor at the top air outlet. After the humidity sensor reading stabilizes, record the value;
[0042] (2-4) Change different input parameters, repeat (2-1) to (2-3), and obtain the mass transfer coefficient k under different input parameters c , and fit multiple sets of recorded data.
[0043] Preferably, the step (2-2) further includes calculating the pore flow rate u0, specifically:
[0044] The pore flow velocity u0 is calculated according to the following formula:
[0045]
[0046] Among them, Q is the flow rate of gas input into the bubbling zone, N is the number of pores, is the pore diameter.
[0047] Preferably, the mass transfer coefficient k obtained in step (2-4) is c , specifically:
[0048] The mass transfer coefficient k is obtained according to the following formula c :
[0049]
[0050] Among them, d b is the average bubble diameter, H1 is the gas phase inlet humidity, H2 is the gas phase outlet humidity, H s is the saturated humidity of the gas at the liquid phase temperature.
[0051] The main feature of the device for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench is that the device comprises:
[0052] a processor configured to execute computer executable instructions;
[0053] The memory stores one or more computer executable instructions. When the computer executable instructions are executed by the processor, the steps of the method for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench are realized.
[0054] The processor for implementing performance prediction processing for a bubbling humidification device of a fuel cell test bench has the main feature that the processor is configured to execute computer executable instructions. When the computer executable instructions are executed by the processor, the various steps of the above-mentioned method for implementing performance prediction processing for a bubbling humidification device of a fuel cell test bench are implemented.
[0055] The main feature of the computer-readable storage medium is that a computer program is stored thereon, and the computer program can be executed by a processor to implement the various steps of the above-mentioned method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench.
[0056] The method, device, processor and computer-readable storage medium for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench of the present invention are adopted. Compared with the existing design stage that completely relies on experience to design the structure of the bubbling humidification device, the method uses a semi-theoretical-semi-empirical model to quickly design the key parameters of the bubbling humidification device in the early stage of design, which greatly reduces the debugging time after the humidification device is manufactured. At the same time, the method can combine the operation data of simulation, experiment or actual device to make subtle corrections to the model, so that more accurate data can be obtained in the subsequent design of other bubbling humidification equipment. BRIEF DESCRIPTION OF THE DRAWINGS
[0057] Figure 1 It is a flow chart of a method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to the present invention.
[0058] Figure 2 The present invention provides a computational fluid dynamics simulation result for correcting the bubble average diameter prediction model coefficient of the method for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench.
[0059] Figure 3 A schematic diagram of a humidifying device for a fuel cell test bench according to the method for implementing performance prediction processing for a bubbling humidifying device for a fuel cell test bench of the present invention.
[0060] Figure 4 It is a fitted surface diagram of the average diameter of bubbles in the method for realizing performance prediction processing for the bubbling humidification device of a fuel cell test bench according to the present invention.
[0061] Figure 5 It is a fitted surface diagram of the area mass transfer coefficient k_c of the method for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench of the present invention.
[0062] Figure 6 The present invention is a relationship curve between the pore flow rate, the bubble diameter and the outlet gas water content calculated by the method for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench.
[0063] Figure 7 The graph is a relationship curve between the pore size, bubble diameter and outlet gas water content calculated by the method for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench of the present invention.
[0064] Figure 8 The graph is a relationship curve between the number of pores, bubble diameter and outlet gas water content calculated by the method for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench of the present invention. DETAILED DESCRIPTION
[0065] In order to more clearly describe the technical content of the present invention, further description is given below in conjunction with specific embodiments.
[0066] The method for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench of the present invention comprises the following steps:
[0067] (1) Calculate the dimensionless number of the fluid according to the structural parameters and physical properties of the bubbling device;
[0068] (2) According to the dimensionless number, the average bubble diameter d at the corresponding airflow velocity is calculated through the bubble average diameter prediction model. b ;
[0069] (3) According to the average bubble diameter d b , and input the gas phase inlet humidity H1, liquid pool liquid phase temperature T l , gas saturation humidity H corresponding to liquid phase temperature s and area mass transfer coefficient k c , the gas phase outlet humidity H2 is obtained according to the gas phase humidification performance prediction model;
[0070] (4) Determine whether the gas phase outlet humidity H2 meets the design target. If so, the design is completed and the step ends; otherwise, readjust the design parameters and continue with step (1).
[0071] As a preferred embodiment of the present invention, the dimensionless number of the fluid calculated in step (1) includes the Bond number Bo. The Bond number Bo is calculated in step (1) as follows:
[0072] The Bond number Bo is calculated according to the following formula:
[0073]
[0074] Where g is the acceleration of gravity, N is the number of pores, d0 is the pore diameter, ρ l is the liquid density, σ l is the surface tension of the liquid.
[0075] As a preferred embodiment of the present invention, the dimensionless number of the fluid calculated in step (1) includes the Galileo number Ga. The Galileo number Ga calculated in step (1) is specifically:
[0076] The Galilean number Ga is calculated according to the following formula:
[0077]
[0078] Where g is the acceleration of gravity, N is the number of pores, d0 is the pore diameter, and v l is the dynamic viscosity of the liquid.
[0079] As a preferred embodiment of the present invention, the dimensionless number of the fluid calculated in step (1) includes the Froude number Fr. The Froude number Fr calculated in step (1) is specifically:
[0080] The Froude number Fr is calculated according to the following formula:
[0081]
[0082] Among them, u0 is the pore flow velocity, g is the gravitational acceleration, N is the number of pores, and d0 is the pore diameter.
[0083] As a preferred embodiment of the present invention, the bubble average diameter prediction model in step (2) is specifically:
[0084] The bubble average diameter prediction model is expressed according to the following formula:
[0085]
[0086] Among them, k, a, b, c are coefficients and exponents that can be corrected, Bo is the Bond number, Ga is the Galileo number, Fr is the Froude number, N is the number of pores, and d0 is the pore diameter.
[0087] As a preferred embodiment of the present invention, the method further includes the step of adjusting the coefficients and exponents k, a, b, c according to the computational fluid dynamics simulation results, specifically including the following steps:
[0088] (1-1) Establish a corresponding humidification equipment CAD model according to the required tank cavity;
[0089] (1-2) Discretize the fluid region inside the tank;
[0090] (1-3) Select a suitable calculation model and input the boundary conditions of the humidification device.
[0091] (1-4) Delineate the liquid area and gas area inside the tank according to the liquid distribution of the tank, and perform initialization and solution;
[0092] (1-5) Post-processing is performed on the calculation results, including extracting the size of each bubble in the liquid pool to obtain the average diameter and diameter distribution of the bubbles.
[0093] As a preferred embodiment of the present invention, the gas phase humidification performance prediction model in step (3) is specifically:
[0094] The gas phase humidification performance prediction model is expressed according to the following formula:
[0095]
[0096] Among them, k c is the mass transfer coefficient, d b is the average bubble diameter, H1 is the gas phase inlet humidity, H s is the saturated humidity of the gas at the liquid phase temperature.
[0097] As a preferred embodiment of the present invention, the method further comprises the step of determining the mass transfer coefficient k c The steps of fitting by bubbling experiment specifically include the following steps:
[0098] (2-1) delivering liquid at a specified temperature T to the bubbling body through a liquid delivery system, so that a liquid pool with a height of h is formed in the bubbling area;
[0099] (2-2) Opening the gas transmission channel, delivering dry air to the gas distribution plate at the bottom of the bubbling body through the air compressor, and forming bubbles through the gas distribution plate; recording the flow rate Q and air pressure of the input gas in the bubbling area through the flow sensor, air pressure sensor and humidity sensor arranged at the outlet of the air compressor;
[0100] (2-3) Turn on the humidity sensor and air pressure sensor at the top air outlet. After the humidity sensor reading stabilizes, record the value;
[0101] (2-4) Change different input parameters, repeat (2-1) to (2-3), and obtain the mass transfer coefficient k under different input parameters c , and fit multiple sets of recorded data.
[0102] As a preferred embodiment of the present invention, the step (2-2) further includes calculating the pore flow rate u0, specifically:
[0103] The pore flow velocity u0 is calculated according to the following formula:
[0104]
[0105] Among them, Q is the flow rate of gas input into the bubbling zone, N is the number of pores, is the pore diameter.
[0106] As a preferred embodiment of the present invention, the mass transfer coefficient k obtained in step (2-4) is c , specifically:
[0107] The mass transfer coefficient k is obtained according to the following formula c :
[0108]
[0109] Among them, d b is the average bubble diameter, H1 is the gas phase inlet humidity, H2 is the gas phase outlet humidity, H s is the saturated humidity of the gas at the liquid phase temperature.
[0110] The device for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench of the present invention comprises:
[0111] a processor configured to execute computer executable instructions;
[0112] The memory stores one or more computer executable instructions. When the computer executable instructions are executed by the processor, the steps of the method for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench are realized.
[0113] The processor of the present invention is used to implement performance prediction processing for the bubbling humidification device of the fuel cell test bench, wherein the processor is configured to execute computer-executable instructions. When the computer-executable instructions are executed by the processor, the various steps of the above-mentioned method for implementing performance prediction processing for the bubbling humidification device of the fuel cell test bench are implemented.
[0114] The computer-readable storage medium of the present invention stores a computer program thereon, and the computer program can be executed by a processor to implement the various steps of the above-mentioned method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench.
[0115] In a specific embodiment of the present invention, a method for predicting the performance of a bubbling humidification device of a fuel cell test bench is proposed.
[0116] The structure of the bubbling humidification equipment is generally divided into: tank, gas distribution plate, gas inlet pipe, liquid inlet pipe, gas outlet pipe and liquid outlet pipe. The tank contains a certain height of liquid, and the gas distribution plate is distributed with a large number of pores of uniform size. Among them, the number and size of pores are the key parameters that determine the humidification effect.
[0117] In view of the above structure and key parameters, the purpose of the present invention can be achieved through the following technical solutions:
[0118] In the first step, dimensionless numbers are calculated based on the structural parameters and physical properties of the bubbling equipment.
[0119] The structural parameters of the bubbling device include pore diameter d0, pore number N, and pore flow rate u0;
[0120] The physical properties include gravitational acceleration g, liquid surface tension σ l , liquid density ρ l and the dynamic viscosity v of the liquid l ;
[0121] Specifically, the acceleration due to gravity g = 9.8 m / s 2 , the surface tension of water is 0.07N / m, and the density of water is 1000kg / m 3 , the dynamic viscosity of water is 1×10 -6 m 2 / s.
[0122] The dimensionless numbers include Bond number Bo, Galileo number Ga and Froude number Fr.
[0123] The Bond number can be calculated by the following formula:
[0124]
[0125] The Galileo number can be calculated by the following formula:
[0126]
[0127] The Froude number can be calculated by the following formula:
[0128]
[0129] In the second step, based on the dimensionless number obtained in the first step, the average bubble diameter d at the corresponding airflow velocity can be calculated through formula (4): b .
[0130]
[0131] The third step is to calculate the average bubble diameter d obtained in the second step. b , and input the gas phase inlet humidity H1, liquid pool liquid phase temperature T l , gas saturation humidity H corresponding to liquid phase temperature s and area mass transfer coefficient k c According to formula (5), the gas phase outlet humidity H2 is obtained.
[0132]
[0133] Step 4: If the gas phase outlet humidity H2 obtained in step 3 does not meet the design target, the design parameters need to be readjusted and the process returns to step 1 to start over. If it meets the target, the design is completed.
[0134] Specifically, the inner cavity of the tank is a cylindrical tank, k = 8507, a = 0.40, b = 0.88, c = 2.23. The coefficients are obtained by computational fluid dynamics simulation, including:
[0135] S1: Establish a corresponding humidification equipment CAD model according to the required tank inner cavity. The humidification equipment CAD model must include the tank wall, the tank internal fluid area, the gas phase inlet pore and the gas phase outlet.
[0136] S2: Discretize the fluid region inside the tank. Preferably, the discretization method is a finite volume method. Preferably, the discretization length is less than 1 mm.
[0137] S3: Select a suitable calculation model and input the boundary conditions of the humidification device. The boundary conditions include the inlet gas phase flow rate u0, the roughness of the tank wall and the gas phase outlet pressure. Preferably, the multiphase flow model uses the heterogeneous Euler model; the viscosity model uses the SSTk-ω model.
[0138] S4: Delineate the liquid area and gas area inside the tank according to the liquid distribution in the tank, and perform initialization and solution.
[0139] S5: Figure 2 This is the calculation result of a single 1mm hole bubbling device at a speed of 0.7m / s. The calculation result is post-processed to extract the size of each bubble in the liquid pool and obtain the average diameter and diameter distribution of the bubbles, such as Figure 4 shown.
[0140] Furthermore, the mass transfer coefficient k c It will be affected by the inlet pore flow rate u0, pore diameter d0 and liquid temperature T, and needs to be fitted through bubbling experiments. The experimental device generally consists of two parts: the bubbling body and the gas-liquid transmission system. Figure 3 As shown, 1 is the bubbling body, 2 is the gas flow control device, 3 is the temperature sensor, 4 is the humidity sensor, 5 is the liquid pump, and 6 is the liquid level observation window.
[0141] The experimental steps include:
[0142] T1: Liquid at a specified temperature T is delivered to the bubbling body through the liquid delivery system, and the height h of the liquid pool in the bubbling area can be observed through 6.
[0143] T2: Open the gas transmission channel, and deliver dry air to the gas distribution plate at the bottom of the bubbling body through the air compressor, and form bubbles through the gas distribution plate. At this time, the flow rate Q and air pressure of the input gas in the bubbling area are recorded by the flow sensor, air pressure sensor and humidity sensor arranged at the outlet of the air compressor. The pore flow rate is calculated by the following formula:
[0144]
[0145] T3: Turn on the humidity sensor and pressure sensor at the top air outlet. After the humidity sensor reading stabilizes, record the value.
[0146] T4: Change different input parameters and repeat T1 to T3. The mass transfer coefficient k under different input parameters is obtained by formula (7): c , and fit multiple sets of recorded data.
[0147]
[0148] When the absolute pressure of the gas is 3 atm and the water temperature is 95°C, the saturated humidity of the gas is 0.19 kg of water per kg of dry air. At this time, according to a performance prediction method for a bubbling humidification device of a fuel cell test bench, when the single hole diameter is 1 mm, the relationship between the pore flow rate, the bubble diameter and the outlet gas water content calculated by the example of the present invention is as follows: Figure 6 As shown; According to a fuel cell test bench bubbling humidification device performance prediction method, when the single hole flow rate is 0.5m / s, the relationship between the pore size, bubble diameter and outlet gas water content calculated by the example of the present invention is as follows Figure 7 As shown; According to a fuel cell test bench bubbling humidification device performance prediction method, when the pore flow rate is 0.5m / s and the pore diameter is 1mm, the relationship between the number of pores, the bubble diameter and the outlet gas water content calculated by the example of the present invention is as follows Figure 8 shown.
[0149] The specific implementation scheme of this embodiment can refer to the relevant description in the above embodiment, which will not be repeated here.
[0150] It can be understood that the same or similar parts of the above embodiments can be referenced to each other, and the contents not described in detail in some embodiments can refer to the same or similar contents in other embodiments.
[0151] It should be noted that, in the description of the present invention, the terms "first", "second", etc. are only used for descriptive purposes and cannot be understood as indicating or implying relative importance. In addition, in the description of the present invention, unless otherwise specified, the meaning of "plurality" refers to at least two.
[0152] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code that includes one or more executable instructions for implementing the steps of a specific logical function or process, and the scope of the preferred embodiments of the present invention includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present invention belong.
[0153] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution device. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.
[0154] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the corresponding program may be stored in a computer-readable storage medium, which, when executed, includes one of the steps of the method embodiment or a combination thereof.
[0155] In addition, each functional unit in each embodiment of the present invention may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.
[0156] The storage medium mentioned above can be a read-only memory, a magnetic disk or an optical disk, etc.
[0157] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "examples", "specific examples", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.
[0158] The method, device, processor and computer-readable storage medium for realizing performance prediction processing for the bubbling humidification device of the fuel cell test bench of the present invention are adopted. Compared with the existing design stage that completely relies on experience to design the structure of the bubbling humidification device, the method uses a semi-theoretical-semi-empirical model to quickly design the key parameters of the bubbling humidification device in the early stage of design, which greatly reduces the debugging time after the humidification device is manufactured. At the same time, the method can combine the operation data of simulation, experiment or actual device to make subtle corrections to the model, so that more accurate data can be obtained in the subsequent design of other bubbling humidification equipment.
[0159] In this specification, the present invention has been described with reference to specific embodiments thereof. However, it is apparent that various modifications and variations may be made without departing from the spirit and scope of the present invention. Therefore, the specification and drawings should be regarded as illustrative rather than restrictive.
Claims
1. A method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench, characterized in that: The method comprises the following steps: (1) Calculate the dimensionless number of the fluid according to the structural parameters and physical properties of the bubbling device; (2) According to the dimensionless number, the average bubble diameter d at the corresponding airflow velocity is calculated through the bubble average diameter prediction model. b ; (3) According to the average bubble diameter d b , and input the gas phase inlet humidity H1, liquid pool liquid phase temperature T l , gas saturation humidity H corresponding to liquid phase temperature s and area mass transfer coefficient k c , the gas phase outlet humidity H2 is obtained according to the gas phase humidification performance prediction model; (4) Determine whether the gas phase outlet humidity H2 meets the design target. If so, the design is completed and the step ends; otherwise, readjust the design parameters and continue with step (1); The gas phase humidification performance prediction model in step (3) is specifically: The gas phase humidification performance prediction model is expressed according to the following formula: Among them, k c is the area mass transfer coefficient, d b is the average bubble diameter, H1 is the gas phase inlet humidity, H s is the saturated humidity of the gas at the liquid phase temperature.
2. The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to claim 1, characterized in that: The dimensionless number of the fluid calculated in step (1) includes the Bond number Bo. The Bond number Bo is calculated in step (1) as follows: The Bond number Bo is calculated according to the following formula: Where g is the acceleration of gravity, N is the number of pores, d0 is the pore diameter, ρ l is the liquid density, σ l is the surface tension of the liquid.
3. The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to claim 1, characterized in that: The dimensionless number of the fluid calculated in step (1) includes the Galileo number Ga. The Galileo number Ga calculated in step (1) is specifically: The Galilean number Ga is calculated according to the following formula: Where g is the acceleration of gravity, N is the number of pores, d0 is the pore diameter, and v l is the dynamic viscosity of the liquid.
4. The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to claim 1, characterized in that: The dimensionless number of the fluid calculated in step (1) includes the Froude number Fr. The Froude number Fr calculated in step (1) is specifically: The Froude number Fr is calculated according to the following formula: Among them, u0 is the pore flow velocity, g is the gravitational acceleration, N is the number of pores, and d0 is the pore diameter.
5. The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to claim 1, characterized in that: The bubble average diameter prediction model in step (2) is specifically: The bubble average diameter prediction model is expressed according to the following formula: Among them, k, a, b, c are coefficients and exponents that can be corrected, Bo is the Bond number, Ga is the Galileo number, Fr is the Froude number, N is the number of pores, and d0 is the pore diameter.
6. The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to claim 5, characterized in that: The method further comprises the step of adjusting the coefficients and exponents k, a, b, c according to the computational fluid dynamics simulation results, specifically comprising the following steps: (1-1) Establish a corresponding humidification equipment CAD model according to the required tank cavity; (1-2) Discretize the fluid region inside the tank; (1-3) Select a suitable calculation model and input the boundary conditions of the humidification device; (1-4) Delineate the liquid area and gas area inside the tank according to the liquid distribution of the tank, and perform initialization and solution; (1-5) Post-processing is performed on the calculation results, including extracting the size of each bubble in the liquid pool to obtain the average diameter and diameter distribution of the bubbles.
7. The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to claim 1, characterized in that: The method further comprises determining the area mass transfer coefficient k c The steps of fitting by bubbling experiment specifically include the following steps: (2-1) delivering liquid at a specified temperature T to the bubbling body through a liquid delivery system, so that a liquid pool with a height of h is formed in the bubbling area; (2-2) Opening the gas transmission channel, delivering dry air to the gas distribution plate at the bottom of the bubbling body through the air compressor, and forming bubbles through the gas distribution plate; recording the flow rate Q and air pressure of the input gas in the bubbling area through the flow sensor, air pressure sensor and humidity sensor arranged at the outlet of the air compressor; (2-3) Turn on the humidity sensor and air pressure sensor at the top air outlet, and record the value after the humidity sensor reading stabilizes; (2-4) Change different input parameters, repeat (2-1) to (2-3), and obtain the mass transfer coefficient k under different input parameters c , and fit multiple sets of recorded data.
8. The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to claim 7, characterized in that: The step (2-2) further includes calculating the pore flow rate u0, specifically: The pore flow velocity u0 is calculated according to the following formula: Among them, Q is the flow rate of gas input into the bubbling zone, N is the number of pores, is the pore diameter.
9. The method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench according to claim 1, characterized in that: The area mass transfer coefficient k obtained in the steps (2-4) is c , specifically: The area mass transfer coefficient k is obtained according to the following formula c : Among them, d b is the average bubble diameter, H1 is the gas phase inlet humidity, H2 is the gas phase outlet humidity, H s is the saturated humidity of the gas at the liquid temperature.
10. A device for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench, characterized in that: The device comprises: a processor configured to execute computer executable instructions; A memory storing one or more computer executable instructions, wherein when the computer executable instructions are executed by the processor, the steps of the method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench as described in any one of claims 1 to 9 are implemented.
11. A processor for implementing performance prediction processing for a bubbling humidification device of a fuel cell test bench, characterized in that: The processor is configured to execute computer executable instructions. When the computer executable instructions are executed by the processor, the various steps of the method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench as described in any one of claims 1 to 9 are implemented.
12. A computer-readable storage medium, characterized in that: A computer program is stored thereon, and the computer program can be executed by a processor to implement the various steps of the method for realizing performance prediction processing for a bubbling humidification device of a fuel cell test bench as described in any one of claims 1 to 9.