JP2005083452A - Evaluation device for hydrogen supply system, and its evaluation method and evaluation program - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide an evaluation device for a hydrogen supply system using a highly reliable MH alloy. <P>SOLUTION: The evaluation device comprises a time step setting part 27 for dividing an operation time for a container filled with a hydrogen storage alloy into small time intervals for analysis, a heat analysis part 23 for performing heat conduction analysis in each time step divided by the time step setting part assuming zero hydrogen migration and calculating parameters regarding the temperatures of inner portions of the MH container, a hydrogen relative parameter calculation part 24 for calculating parameters regarding hydrogen migration in each time step identical with that of the heat conduction analysis by a heat conduction analysis means assuming the zero hydrogen migration and calculating hydrogen reaction heat, a temperature parameter correction part 25 for correcting the parameters regarding the temperatures using the reaction heat, and a hydrogen parameter correction part 26 for correcting the parameters regarding the hydrogen migration using the parameters regarding the corrected temperatures. <P>COPYRIGHT: (C)2005,JPO&NCIPI

Description

  The present invention relates to an evaluation apparatus, an evaluation method, and an evaluation program for evaluating a hydrogen supply system using a hydrogen storage alloy (hereinafter referred to as MH alloy) for supplying hydrogen to a fuel cell.

  Conventionally, when evaluating the performance of hydrogen absorption / release in a hydrogen supply system using an MH alloy, heat transfer calculation in the alloy layer is performed, and the release / absorption amount of hydrogen is calculated based on the heat and temperature obtained as a result. Calculated as rate limiting by heat transfer. For the alloy physical property values at this time, results obtained separately for the alloy were employed as constants.

  However, in the heat transfer / hydrogen transfer of MH alloy, hydrogen transfer and heat transfer occur at the same time and affect each other. Therefore, as in the past, the effect of hydrogen transfer is ignored and only the heat transfer is analyzed. However, it is inherently undesirable to obtain the amount of hydrogen transfer based on the analysis result, and a highly accurate performance evaluation cannot be performed.

  That is, as another example where hydrogen and heat transfer occur at the same time, there is known an example in which a change in the concentration of hydrogen absorbed in a steel material is obtained as a problem of hydrogen diffusion in the steel material when welding the steel material. However, in this case, the movement of hydrogen is affected by the movement of heat, but the movement of heat can be treated as not affected by the movement of hydrogen. Therefore, in this case, it is possible to analyze the heat conduction first and solve the hydrogen diffusion problem considering the influence of temperature change. Furthermore, since the basic equation approximates heat conduction, hydrogen diffusion can be handled by slightly modifying the heat conduction program.

  On the other hand, in the case of hydrogen in the MH alloy, as described above, the movement of hydrogen and the movement of heat occur at the same time and affect each other. Analyzing the amount of hydrogen transferred from the vessel can accurately determine the hydrogen absorption / release performance of a vessel in which a large pressure difference occurs in the hydrogen pressure due to the site in the MH vessel, particularly in an MH alloy vessel. Even if such an analysis method is applied to an evaluation device or the like, a reliable evaluation result cannot be obtained.

  The present invention has been made to solve such problems, and an object thereof is to provide an evaluation apparatus, an evaluation method, and an evaluation program for a hydrogen supply system using a highly reliable MH alloy. .

  In order to solve the above-described problems, the present invention is an evaluation apparatus for a hydrogen supply system using a hydrogen storage alloy, and in order to perform an operation simulation of a container filled with the hydrogen storage alloy, a time corresponding to the operation time is set. A time step dividing unit that divides into minute time intervals for performing analysis, and for each time step divided by the time step dividing unit, a heat transfer analysis is performed assuming that hydrogen transfer is zero, and each part in the MH container is analyzed. A heat conduction analysis means for calculating a parameter relating to temperature, and a hydrogen transfer amount analysis means for calculating a parameter relating to hydrogen transfer assuming that heat conduction is zero at every same time step in which heat conduction analysis is performed by the heat conduction analysis means, Reaction heat calculation means for calculating heat of hydrogen reaction based on parameters relating to hydrogen movement analyzed by the hydrogen movement amount analysis means; The temperature parameter correcting means for correcting the temperature-related parameter analyzed by the heat conduction analyzing means using the heat of reaction calculated by the heat calculating means, and the temperature-related parameter corrected by the temperature correcting means, Hydrogen parameter correcting means for correcting a parameter relating to hydrogen movement obtained by the movement amount analyzing means.

  In the hydrogen supply system evaluation apparatus according to the present invention, the hydrogen transfer parameters analyzed by the hydrogen transfer amount analysis means include a hydrogen pressure and a hydrogen transfer amount, and are corrected by the hydrogen parameter correction means. The parameter relating to movement includes hydrogen pressure. Here, a pressure loss coefficient can be introduced into the relationship between the hydrogen pressure and the amount of hydrogen transfer.

  Further, the present invention is a method for evaluating a hydrogen supply system using a hydrogen storage alloy, and is a micro method for analyzing a time corresponding to an operation time in order to perform an operation simulation of a container filled with the hydrogen storage alloy. A time step dividing step divided into various time intervals, and for each time step divided by the time step dividing step, a heat conduction analysis is performed on the assumption that hydrogen transfer is zero, and a parameter relating to the temperature of each part in the MH container is calculated. A heat transfer analysis step, a hydrogen transfer amount analysis step for calculating a parameter related to hydrogen transfer on the assumption that heat transfer is zero at every same time step in which heat transfer analysis is performed by the heat transfer analysis step, and the hydrogen transfer amount analysis step A reaction heat calculation step for calculating the heat of hydrogen reaction based on the parameters relating to hydrogen movement analyzed by The temperature parameter correction step of correcting the temperature-related parameter analyzed by the heat conduction analysis step using the reaction heat calculated by the reaction heat calculation step, and the temperature-related parameter corrected by the temperature correction step, A hydrogen parameter correcting step for correcting a parameter relating to hydrogen transfer obtained by the hydrogen transfer amount analyzing step.

  Further, the present invention is a hydrogen supply system evaluation program for causing a computer to execute a hydrogen supply system evaluation method using a hydrogen storage alloy, and for operating simulation of a container filled with the hydrogen storage alloy. A time step division step that divides the time corresponding to the time into minute time intervals for analysis, and for each time step divided by the time step division step, heat conduction analysis is performed assuming that hydrogen transfer is zero, and MH A heat transfer analysis step for calculating a parameter related to the temperature of each part in the container, and a hydrogen transfer for calculating a parameter for hydrogen transfer assuming that heat transfer is zero at each same time step in which heat transfer analysis is performed by the heat transfer analysis step. And a parameter relating to hydrogen transfer analyzed in the hydrogen transfer amount analysis step. A reaction heat calculation step for calculating the heat of hydrogen reaction based on the data, and a temperature parameter correction step for correcting the parameter relating to the temperature analyzed by the heat conduction analysis step using the reaction heat calculated by the reaction heat calculation step And a hydrogen parameter correcting step of correcting a parameter relating to hydrogen movement obtained in the hydrogen movement amount analyzing step using a parameter relating to the temperature corrected in the temperature correcting step.

  As described above, according to the present invention, it is possible to obtain an evaluation apparatus, an evaluation method, and an evaluation program for a hydrogen supply system using a highly reliable MH alloy.

  Hereinafter, embodiments of the present invention will be described with reference to the drawings.

  FIG. 1 is a block diagram showing a hydrogen supply system according to an embodiment of the present invention.

  This hydrogen supply system includes an MH alloy container 1 that is filled with an MH alloy and is an object of evaluation, a hot water supply device 2 that heats the MH alloy that is filled in the MH alloy container 1, and a cooling water supply device that cools the MH alloy. 3, a fuel cell 4 to which hydrogen gas released from the MH alloy container 1 is supplied, and a hydrogen supply device 5 for supplying and occluding hydrogen to the MH alloy. Controllers 6A and 6B for measuring and controlling the flow rate and the like are provided.

  2 and 3, the MH alloy container 1 has a cylindrical shape in which a plurality of cylindrical capsule containers 11 are joined in multiple stages, and each capsule container 11 is partitioned by a plurality of disc-shaped fins 12, MH alloy 13 is filled between these fins. In addition, a ventilation member 14 is provided so as to pass through the center of each disk-like fin 12. The ventilation material 14 extends to the outer surface of the capsule container 11 along the fins of the adjacent capsule containers in each capsule container. A heat exchange unit 15 connected to the hot water supply device 2 or the cooling water supply device 3 is provided on the outer surface of the capsule container 11.

  In such a configuration, when the MH alloy is cooled, the hydrogen supplied from the hydrogen supply device 5 is stored in the MH alloy through the ventilation member 14. On the other hand, when the MH alloy is heated, hydrogen stored in the MH alloy is released and supplied to the fuel cell 4 through the ventilation member 14.

  The present invention performs a simulation for evaluating the storage and release of hydrogen in the MH alloy container. The virtual operation time of the MH container is finely divided into time steps, and (1) Zero hydrogen transfer at each time. (2) Assuming zero heat conduction and analyzing only hydrogen transfer, the actual phenomenon is analyzed as the sum of the solutions of (1) and (2) above. (This method is known as the superposition principle in the physical problem, and is considered to approach the actual phenomenon as much as possible by making the time step finer).

  Known numerical analysis methods for heat conduction problems include the finite element method, the finite difference method (difference method), the boundary element method, and the like. Each of these methods has characteristics, but as a method of analyzing a complex physical phenomenon, it is easy to apply a difference method close to a physical model.

  FIG. 4 is a block diagram showing an evaluation apparatus for a hydrogen supply system using a hydrogen storage alloy in the embodiment, and FIG. 5 is a flowchart showing the operation of this apparatus.

  The evaluation apparatus for the hydrogen supply system using the MH alloy shown in FIG. 4 is based on a data input unit 21 for inputting data such as a predetermined analysis model and analysis conditions and various initial values, and an input initial value and the like. An initial hydrogen amount calculation unit 22 that calculates an initial hydrogen amount, a thermal analysis unit 23 that calculates a temperature of each part by ignoring parameters related to hydrogen transfer, and a hydrogen transfer parameter that ignores parameters related to heat transfer As parameters, the hydrogen pressure and the amount of movement of each part are analyzed, and the hydrogen-related parameter analysis part 24 for calculating the reaction heat of hydrogen based on the analysis result, and each part calculated by the thermal analysis part using this reaction heat A temperature correction unit 25 that corrects the temperature, a hydrogen pressure correction unit 26 that corrects the hydrogen pressure of each part using the corrected temperature of each part, and an operation time step of each part 23 to 26. And a time step setting unit 27 for setting up.

  The calculation of the thermal analysis unit 23 to the hydrogen pressure correction unit 26 is performed by the time step setting unit 27 for each time step obtained by dividing the operation time of the hydrogen supply system (MH alloy container) into minute times. The result is output as an evaluation value.

  According to the flowchart of FIG. 5, first, an analysis model is designated by the data input unit 21 (S1), shape dimension data is input (S2), and PCT (hydrogen pressure P, hydrogen concentration C, temperature T) curve data. Input (S3) of (MH storage characteristics of MH alloy), input of analysis conditions (S4), and input of initial values (MH alloy temperature Ti, hydrogen pressure Pi) (S5).

  Next, the initial hydrogen amount is calculated (S6). This initial hydrogen amount is calculated using the subroutine (1) for calculating the hydrogen amount from the pressure using the PCT diagram shown in FIG. 10, and using the PCT curve from the MH alloy temperature and the hydrogen pressure Pi. Determine the amount of hydrogen.

  Next, a time step is set in step S7, and steps S8 to S11 are processed at every minute time interval to obtain parameters relating to hydrogen transfer and temperature parameters.

  First, in step S8, thermal analysis is performed by the thermal analysis unit, and the amount of heat transfer and the temperature of each unit are obtained as parameters relating to temperature. This thermal analysis is performed by a subroutine (2) for obtaining the thermal conductivity of the MH alloy from the hydrogen pressure Pi using the experimental values shown in FIG.

  Once the thermal analysis is performed in step S8, the hydrogen pressure and transfer amount of each part are analyzed as parameters relating to hydrogen transfer, and the heat of hydrogen reaction is calculated (S9). In this analysis, a subroutine (3) described later is used, and the hydrogen flow velocity is obtained from the hydrogen pressure of each part.

When the hydrogen reaction heat is calculated in step S9, the temperature of each part obtained in step S8 is corrected by adding / subtracting the hydrogen reaction heat to the heat amount fluctuation of each part obtained by the thermal analysis in step S8. When the temperature of each part is corrected in this way, the hydrogen pressure of each part is corrected using the temperature, the hydrogen amount of each part obtained in step S6, and the PCT curve (S11), and the hydrogen generation amount is obtained. By causing the computer to execute the operations in steps S1 to S11 described above, the evaluation program for the hydrogen supply system according to the present invention can be provided.
(Thermal analysis: Step S8)
A model of the MH alloy container in the present embodiment is shown in FIG. 6 (corresponding to FIG. 3). The MH alloy container 1 according to the embodiment divides the inside of a cylindrical capsule container 11 with fins 12 in multiple stages with respect to the longitudinal direction, and follows the surface of each fin from the center of the cylindrical axis to the outer surface of the capsule. A ventilation member 14 is provided, and a region surrounded by the ventilation member 14 and the fins 12 is filled with the MH alloy 13.

  In the embodiment, the portion filled with the MH alloy is divided into, for example, five parts (parameter m) in the radial direction of the container, and the thickness direction (layer direction) is divided into four parts (parameter n), for example. I want to ask. In FIG. 6, PH (M, N) and C (M, N) indicate the hydrogen pressure and hydrogen storage amount at the (M, N) point, respectively, Tm and T (M, N) indicate the temperature of each part, h1 to h3 indicate the thermal conductivity of each part.

  In the difference method, the difference between the values of each point is obtained sequentially. In this difference method, the temperature difference of each part is used for heat transfer calculation, the pressure difference is used for calculation of the hydrogen flow rate, and the change in the hydrogen storage amount is calculated from the hydrogen flow rate obtained from the pressure difference.

  In order to facilitate the understanding of the analysis, consider a model in which heat conduction and hydrogen transfer in an MH alloy are close to a one-dimensional flow as shown in FIG.

In general, the equation for heat conduction is as follows: T is temperature, t is time, y is a coordinate axis,
^{∂T / ∂t = a · ∂ 2} T / ∂x 2
a = λ / cγ
Here, λ is thermal conductivity, c is specific heat, and γ is specific gravity.
The difference formula in this case is
T (x, t + Δt) = T (x, t) + a / Δt ² [T (x + Δx, t) −2T (x, t) + T (x−Δx, t)]
Therefore, when there is little change in the physical property value due to the one-dimensional problem, it can be analyzed using this equation. However, when there are more complicated shapes and changes in physical property values, it is easier to replace them with more primitive models as shown in FIG.

  Based on FIG. 7, the relationship (unsteady change) of heat transfer in each component is shown as a temperature change due to heat conduction in the heat exchange unit material as follows.

Qin, the heat flux passing through the x-plane for Δt
ΔQ _in = λ _h × (Δθ / Δx _h ) _x × A × Δt
Qout, the heat flux passing through the x + Δx plane during Δt
ΔQ _out = λ _h × (Δθ / Δx _h ) _{x +} Δ _x × A × Δt
The difference between the two is equal to the amount of heat required for the temperature rise of the element. Ρ _h C _ph × A × (Δx) × Δθ = (ΔQ _in −ΔQ _out )
= Λ _h × ((Δθ / Δx _h ) _x − (Δθ / Δx _h ) _{x +} Δ _x ) x × A × Δt
Change in heat per unit time ρ _h C _ph (Δθ / Δt) = (ΔQ _in −ΔQ _out )
= Λ _h × ((Δθ / Δx _h ) _x − (Δθ / Δx _h ) _{x +} Δ _x ) / Δx [kW]
Where ρ is the specific heat, λ is the thermal conductivity, λ _h is the thermal conductivity of the heat exchange unit material, (Δθ / Δx _h ) _x is the temperature gradient at x, and (Δθ / Δx _h ) _{x +} Δ _x is at x + Δx. The temperature gradient, Δx, is the unit length in the x direction of the element.

  A point to be noted in the thermal conductivity analysis of the MH alloy is that the thermal conductivity of the MH alloy changes greatly depending on the pressure of hydrogen and the like.

  The hydrogen pressure varies greatly with the absorption and release of hydrogen. Therefore, a thermal conductivity calculation subroutine for the MH alloy is prepared, the thermal conductivity is calculated at each time step, and the temperature of each part is calculated using the obtained thermal conductivity.

FIG. 8 shows an example of measurement of thermal conductivity. If this measurement result is used, the thermal conductivity can be set from the pressure for each time. The thermal conductivity of the MH alloy varies depending on the alloy packing density. The figure shows an example of theoretical calculation results, but the thermal conductivity calculation program can also reflect this effect. In subroutine (2) shown in the thermal analysis (step S8), the thermal conductivity of the MH alloy is calculated for each unit time from the hydrogen pressure in FIG. Regarding the calculation of the thermal conductivity, in the initial calculation step, the input initial value is used as the pressure, and in the subsequent steps, the pressure corrected for each step is used.
(Hydrogen transfer analysis: Step S9)
The hydrogen transfer analysis is performed using the segmentation elements used in the heat conduction analysis. FIG. 9 is an explanatory diagram showing a basic model for performing this hydrogen transfer analysis. Subroutine (3) is used for this analysis. Specifically, 1) The hydrogen pressure is calculated based on the PCT diagram from the temperature and the hydrogen amount at each point subjected to the thermal analysis. 2) The pressure difference between the elements and the average pressure are obtained from the hydrogen pressure at each point, and the hydrogen flow velocity flowing into the element is obtained from these values. In the initial step of the calculation, the initial hydrogen amount can be input as an initial value. However, using the initial hydrogen amount calculated using the subroutine (1) is convenient because it saves input. is there.

  The relationship between the pressure loss and the hydrogen flow rate can be experimentally obtained and converted into a subroutine by formulating it. In general, the pressure loss is often approximated as being proportional to the square of the flow velocity, but in the case of a particle-packed layer such as an MH alloy layer, the pressure loss is considered to be proportional to the flow velocity because it is a laminar flow.

  The increase / decrease in the hydrogen amount of each element is obtained by solving Fi-Fc for each element shown in FIG. That is, the differential equation (FIGS. 9C and 9D) for the pressure at each point is solved, and based on the pressure difference at each point, the hydrogen flow rate (Fi: amount of hydrogen gas introduced into the element, Fc: out of the element). The amount of released hydrogen gas) can be obtained, and the increase or decrease can be calculated from the flow rate difference between the points. Note that the pressure at each point is obtained by using the input initial value in the initial stage of the calculation step and solving the difference equation for the value obtained in the immediately preceding step in the subsequent steps. Further, the amount of hydrogen is obtained by adding or subtracting the obtained increase or decrease of the amount of hydrogen to the initial value to obtain the amount of hydrogen in the next step. Furthermore, in this embodiment, a pressure loss coefficient (FIG. 9B) indicating the relationship between the amount of hydrogen introduced (released) into each element and the pressure loss in the alloy layer is incorporated to further improve the evaluation accuracy. .

Thus, when the hydrogen amount V _H2MH flowing into each element is obtained, the reaction heat Q _MH is obtained according to the following equation using the hydrogen amount, and the temperature change due to the reaction heat Q _MH can be obtained.

Q _MH = V _H2MH /22.42×(±ΔH ^*)
Here, ΔH ^* is a value obtained by experiment, + indicates occlusion, and − indicates release.
(Temperature and hydrogen pressure correction)
As described above, when the heat of reaction is obtained, the temperature of each part obtained by the thermal analysis is corrected using the heat, and the hydrogen pressure is obtained from the PCT diagram shown in FIG. 11 using the corrected temperature. To correct the hydrogen pressure used so far. The program for obtaining the hydrogen pressure from the PCT diagram shown in FIG. 11 is used as subroutine (4).

  And it progresses to the calculation in the next time step using this corrected hydrogen pressure and the temperature of each part.

  For the analysis according to the present invention, the difference method is basically used. However, the difference method is likely to diverge if each point division is inappropriate due to its nature. is there.

  If an element with a small heat capacity is connected, the part is calculated as steady heat conduction, and it is necessary to devise measures such as removing unstable elements.

  Further, since there is a possibility that an error is included in the analysis value, it is necessary to verify the accuracy for each step, and to make the time step finer if the error becomes large.

It is a figure which shows the hydrogen supply system in embodiment of this invention. It is sectional drawing which shows MH alloy container. It is a figure which shows the calculation model in MH alloy container. It is a block diagram which shows the evaluation apparatus in embodiment. It is a flowchart which shows operation | movement of embodiment. It is a figure which shows an example of the division | segmentation in MH alloy container for applying the difference method. It is a figure which shows the basic model of heat conduction analysis. It is a figure which shows hydrogen pressure and the heat conductivity of MH alloy. It is a figure which shows the basic model for showing hydrogen movement. It is a figure which shows the example of a calculation of the subroutine (1) using a PCT diagram. It is a figure which shows the example of a calculation of the subroutine (4) using a PCT diagram.

Explanation of symbols

  DESCRIPTION OF SYMBOLS 1 MH alloy container, 2 Hot water supply apparatus, 3 Cooling water supply apparatus, 4 Fuel cell, 5 Hydrogen supply apparatus, 6A, 6B Controller, 11 Capsule container, 12 Fin, 13 MH alloy, 14 Ventilation material, 15 Heat exchange unit , 21 Data input unit, 22 Initial hydrogen amount calculation unit, 23 Thermal analysis unit, 24 Hydrogen related parameter analysis unit, 25 Temperature correction unit, 26 Hydrogen pressure correction unit, 27 Time step setting unit.

Claims (5)

An apparatus for evaluating a hydrogen supply system using a hydrogen storage alloy,
In order to perform an operation simulation of the container filled with the hydrogen storage alloy, a time step dividing means for dividing the time corresponding to the operation time into minute time intervals for analyzing,
For each time step divided by the time step dividing means, heat conduction analysis is performed assuming that hydrogen transfer is zero, and heat conduction analysis means for calculating parameters relating to the temperature of each part in the MH container;
For each same time step in which heat conduction analysis is performed by the heat conduction analysis means, a hydrogen transfer amount analysis means for calculating a parameter relating to hydrogen transfer assuming that heat conduction is zero, and
Reaction heat calculation means for calculating the heat of hydrogen reaction based on the parameters relating to hydrogen movement analyzed by the hydrogen movement amount analysis means;
A temperature parameter correcting means for correcting a parameter relating to the temperature analyzed by the heat conduction analyzing means using the reaction heat calculated by the reaction heat calculating means;
An evaluation apparatus for a hydrogen supply system, comprising: a hydrogen parameter correcting unit that corrects a parameter related to hydrogen transfer obtained by the hydrogen transfer amount analyzing unit using a parameter related to the temperature corrected by the temperature correcting unit. In the hydrogen supply system evaluation apparatus according to claim 1,
The hydrogen transfer parameters analyzed by the hydrogen transfer amount analysis means include hydrogen pressure and hydrogen transfer amount,
The hydrogen supply system evaluation apparatus, wherein the hydrogen parameter is corrected by the hydrogen parameter correcting means, and the hydrogen pressure is included in the parameter. The evaluation apparatus for a hydrogen supply system according to claim 2,
An apparatus for evaluating a hydrogen supply system, wherein a pressure loss coefficient is introduced into a relationship between the hydrogen pressure and a hydrogen transfer amount. A method for evaluating a hydrogen supply system using a hydrogen storage alloy,
In order to perform an operation simulation of the container filled with the hydrogen storage alloy, a time step dividing step for dividing the time corresponding to the operation time into minute time intervals for analyzing,
For each time step divided by the time step dividing step, a heat conduction analysis is performed assuming that hydrogen transfer is zero, and a parameter relating to the temperature of each part in the MH container is calculated.
For each same time step in which heat conduction analysis is performed by the heat conduction analysis step, a hydrogen transfer amount analysis step for calculating a parameter relating to hydrogen transfer on the assumption that heat conduction is zero, and
A reaction heat calculation step for calculating a heat of hydrogen reaction based on a parameter relating to hydrogen transfer analyzed in the hydrogen transfer amount analysis step;
A temperature parameter correction step of correcting a parameter relating to the temperature analyzed by the heat conduction analysis step using the reaction heat calculated by the reaction heat calculation step;
A hydrogen parameter correcting step for correcting a parameter relating to hydrogen transfer obtained by the hydrogen transfer amount analyzing step using a parameter related to temperature corrected by the temperature correcting step. A hydrogen supply system evaluation program for causing a computer to execute a hydrogen supply system evaluation method using a hydrogen storage alloy,
In order to perform an operation simulation of the container filled with the hydrogen storage alloy, a time step dividing step for dividing the time corresponding to the operation time into minute time intervals for analyzing,
For each time step divided by the time step dividing step, a heat conduction analysis is performed assuming that hydrogen transfer is zero, and a parameter relating to the temperature of each part in the MH container is calculated.
For each same time step in which heat conduction analysis is performed by the heat conduction analysis step, a hydrogen transfer amount analysis step for calculating a parameter relating to hydrogen transfer on the assumption that heat conduction is zero, and
A reaction heat calculation step for calculating a heat of hydrogen reaction based on a parameter relating to hydrogen transfer analyzed in the hydrogen transfer amount analysis step;
A temperature parameter correction step of correcting a parameter relating to the temperature analyzed by the heat conduction analysis step using the reaction heat calculated by the reaction heat calculation step;
A hydrogen supply system evaluation program that causes a computer to execute a hydrogen parameter correction step of correcting a parameter related to hydrogen transfer obtained by the hydrogen transfer amount analysis step using a parameter related to temperature corrected by the temperature correction step.

Images