A method and device for solving liquid hydrogen leakage and diffusion based on atmospheric boundary layer inlet velocity

Abstract

The present invention discloses a method and device for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity, which belongs to the field of numerical solution, including: based on the atmospheric boundary layer inlet velocity profile characteristics, combining the atmospheric boundary layer characteristics and the liquid hydrogen leakage and diffusion model, using conjugate heat transfer to solve the density, and solving the thermal diffusion coefficient at the same time, so that it is closer to the actual physical process, so that the simulation result is more realistic. PBiCGStab is used to introduce residuals, virtual residuals and virtual search directions, which has a faster convergence speed and higher calculation efficiency. The PIMPLE algorithm predicts and corrects the optimization strategy, which has the advantages of high stability and precision. The present invention can accurately describe the diffusion laws in complex environments such as the airflow state, diffusion behavior characteristics, and hydrogen concentration distribution of hydrogen clouds after liquid hydrogen leakage.

Description

Technical Field

[0001] The invention belongs to the field of numerical solution, and in particular relates to a method and a device for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity. Background Art

[0002] Liquid hydrogen has the characteristics of high energy density, cleanness, environmental protection, storage and transportation, which makes it widely used in aerospace, chemical energy and transportation. However, the safety risks brought by liquid hydrogen leakage have gradually become the focus of attention. ) of low temperature ( ) characteristics and its extremely high volatility make it not only gasify quickly during leakage, but also undergo a complex diffusion process in the atmospheric boundary layer. Especially under different wind speeds, temperatures and atmospheric turbulence conditions, its diffusion behavior is difficult to accurately predict. Therefore, accurate numerical simulation is crucial. At present, the simulation of liquid hydrogen leakage and diffusion mainly faces problems such as large amount of calculation, complex model and difficulty in handling boundary conditions.

[0003] The traditional liquid hydrogen leakage diffusion model has a slow iterative convergence speed and has limitations in describing the velocity profile of the atmospheric boundary layer. It cannot fully consider the complexity of the atmospheric boundary layer and cannot truly and accurately reflect the actual situation. This leads to insufficient prediction accuracy and safety risks. Summary of the invention

[0004] In order to overcome the shortcomings of the prior art and to address the serious safety issues that may be caused by the leakage and diffusion of liquid hydrogen, the present invention proposes a method and device for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity, so as to provide a scientific basis and technical support for the safe use of liquid hydrogen.

[0005] The liquid hydrogen leakage and diffusion numerical simulation solver of the present invention is based on the principle of computational fluid dynamics (CFD), combined with the physical properties of liquid hydrogen and the complex factors in the leakage and diffusion process, to construct a high-precision numerical solution model. The model can accurately describe the diffusion process of liquid hydrogen from the leakage source to the surrounding environment, including multiple stages such as liquid hydrogen phase change, gas jet, turbulent diffusion, etc. Through this numerical simulation solution, the diffusion law after liquid hydrogen leakage can be deeply analyzed, and the diffusion range and impact degree under different conditions (such as wind speed, ambient temperature, leakage amount, etc.) can be evaluated, so as to formulate corresponding safety precautions.

[0006] In order to achieve the above object, the present invention adopts the following technical solution:

[0007] A method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity comprises the following steps:

[0008] Step 1: According to the physical properties of liquid hydrogen and the leakage scenario, a liquid hydrogen leakage diffusion model is established. The liquid hydrogen leakage diffusion model includes the mass conservation equation, momentum conservation equation and energy conservation equation;

[0009] Step 2: Determine the characteristics of the atmospheric boundary layer, and obtain relevant physical parameters of the atmospheric boundary layer inlet velocity through meteorological data or field measurement data, wherein the relevant physical parameters include reference height, average stream wind speed, and gas dynamic roughness;

[0010] Step 3: Combining the atmospheric boundary layer characteristics with the liquid hydrogen leakage and diffusion model, the atmospheric boundary layer characteristics mathematical model is introduced into the liquid hydrogen leakage and diffusion model, and the boundary conditions and initial conditions in the liquid hydrogen leakage and diffusion model are corrected;

[0011] Step 4: Based on the modified liquid hydrogen leakage diffusion model in step 3, numerical algorithms and optimization strategies are used to improve computational efficiency and convergence, and to handle large-scale computational domains and complex boundary conditions;

[0012] Step 5: simulate the leakage and diffusion of liquid hydrogen to obtain simulation results, and verify them through experimental data or known results to ensure the accuracy and reliability of the simulation;

[0013] Step 6: Analyze the simulation results and extract key parameters and features to provide a scientific basis for liquid hydrogen storage, transportation and safety assessment; the key parameters and features include diffusion distance, hydrogen concentration, lower flammable limit, hydrogen cloud volume, and diffusion profile.

[0014] The present invention also provides a liquid hydrogen leakage and diffusion solving device based on the atmospheric boundary layer inlet velocity, comprising the following modules:

[0015] The model building module builds a liquid hydrogen leakage and diffusion model based on the physical properties of liquid hydrogen and the leakage scenario. The liquid hydrogen leakage and diffusion model includes the mass conservation equation, momentum conservation equation and energy conservation equation.

[0016] A parameter acquisition module determines the characteristics of the atmospheric boundary layer and obtains relevant physical parameters of the atmospheric boundary layer inlet velocity through meteorological data or field measurement data, wherein the relevant physical parameters include reference height, average stream wind speed, and gas dynamic roughness;

[0017] The correction module combines the atmospheric boundary layer characteristics with the liquid hydrogen leakage and diffusion model, introduces the atmospheric boundary layer characteristics mathematical model into the liquid hydrogen leakage and diffusion model, and corrects the boundary conditions and initial conditions in the liquid hydrogen leakage and diffusion model;

[0018] The boundary processing module, based on the modified liquid hydrogen leakage diffusion model, adopts numerical algorithms and optimization strategies to improve computational efficiency and convergence, and handle large-scale computational domains and complex boundary conditions;

[0019] The simulation verification module simulates the leakage and diffusion of liquid hydrogen, obtains simulation results, and verifies them through experimental data or known results to ensure the accuracy and reliability of the simulation;

[0020] The result analysis module analyzes the simulation results and extracts key parameters and features to provide a scientific basis for liquid hydrogen storage, transportation and safety assessment; the key parameters and features include dispersion distance, hydrogen concentration, lower flammable limit, hydrogen cloud volume, and diffusion profile.

[0021] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the above-mentioned method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity are implemented.

[0022] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-mentioned method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity.

[0023] The beneficial effects of the present invention compared with the prior art are:

[0024] 1. The present invention is based on the atmospheric boundary layer inlet velocity. By combining the atmospheric boundary layer characteristics with the liquid hydrogen leakage diffusion model, the present invention can more accurately simulate the diffusion process of liquid hydrogen in the atmosphere and improve the accuracy of the simulation.

[0025] 2. Density is solved using conjugate heat transfer to obtain the components Density ,density It is obtained from the mass equation and the advection equation, which can simulate the density change more accurately. At the same time, the thermal diffusion coefficient is solved, and the conjugate heat transfer is closer to the actual physical process, making the simulation results more realistic.

[0026] 3. PBiCGStab (preconditioned stable biconjugate gradient method) is used to iteratively solve the linear equations, which has a faster convergence speed and higher computational efficiency.

[0027] 4. The PIMPLE algorithm prediction-correction strategy is used to solve the liquid hydrogen leakage diffusion problem, which has the advantages of high stability and accuracy.

[0028] 5. This invention is suitable for simulation of liquid hydrogen leakage and diffusion of different scales, including small-scale leakage and large-scale leakage scenarios, and has wide applicability. It evaluates the diffusion area and hydrogen concentration after the leakage, and provides safety warnings in combination with the fire and explosion risks of hydrogen concentration, providing strong support for the safe management of liquid hydrogen.

[0029] In summary, the present invention has the advantages of high accuracy, high computational efficiency, and easy solution, and has broad application prospects and research value in liquid hydrogen safety assessment. The present invention establishes a liquid hydrogen leakage and diffusion solution method based on the atmospheric boundary layer inlet velocity, and uses the atmospheric boundary layer inlet velocity as an important parameter to describe the characteristics of the atmospheric boundary layer, which plays a vital role in improving the accuracy and precision of liquid hydrogen leakage and diffusion simulation. The PIMPLE algorithm prediction-correction strategy is used for solution to more accurately describe the airflow state, diffusion behavior characteristics, hydrogen concentration distribution, etc. of the hydrogen cloud in the area after the leakage, especially the diffusion law of the hydrogen cloud in a complex environment, which helps to more deeply understand the physical mechanism of liquid hydrogen leakage and diffusion. This makes the liquid hydrogen leakage and diffusion simulation have important application value in the safety assessment and accident emergency disposal in the field of hydrogen energy, and provides a scientific basis for preventing and avoiding safety accidents such as fire and explosion. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Figure 1 This is a flow chart of a method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity of the present invention;

[0031] Figure 2 It is the structural hierarchy diagram of the coupled phase change model;

[0032] Figure 3 It is a geometric structure grid distribution diagram in the embodiment;

[0033] Figure 4 Flowchart for parameter setting of embodiment;

[0034] Figure 5 Graph showing the concentration distribution of hydrogen cloud in the embodiment. DETAILED DESCRIPTION

[0035] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0036] like Figure 1 As shown, a method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity of the present invention comprises the following steps:

[0037] Step 1: Based on the physical properties of liquid hydrogen and the leakage scenario, a liquid hydrogen leakage and diffusion model is established, including the mass conservation equation, momentum conservation equation, and energy conservation equation.

[0038] ;

[0039] in, is the divergence operator, is the density, is the velocity vector, is the net mass transfer rate.

[0040] Components The interface advection equation is:

[0041] ;

[0042] in, It is a component The volume fraction of It is a component The density of It is a component The mass transfer rate, is the mass source term. t represents time.

[0043] ;

[0044] in, is the correction pressure, is the dynamic viscosity, is the acceleration due to gravity, is the surface tension. The superscript T indicates the transpose of the matrix.

[0045] ;

[0046] in, is a variable constructed from the intermediate velocity field, are the diagonal terms of the momentum matrix equation, is the density of the liquid phase, is the density of the gas phase, is about the pressure Poisson source term.

[0047] ;

[0048] in, is the specific heat at constant pressure, is the temperature, is the effective thermal conductivity, is the energy source term, It is latent heat.

[0049] Step 2: Determine the characteristics of the atmospheric boundary layer. Obtain relevant physical parameters of the atmospheric boundary layer inlet velocity through meteorological data or field measurements, including reference altitude, average streamwise wind speed, gas dynamic roughness, etc.

[0050] ;

[0051] ;

[0052] in, is the surface normal streamwise velocity distribution, is the surface normal turbulent kinetic energy profile, is the ground normal coordinate component, is the ground normal displacement height, is the aerodynamic roughness, is the friction speed, is the reference altitude, It is a reference The average streamwise wind speed at .

[0053] Step 3: Combine the atmospheric boundary layer characteristics with the liquid hydrogen leakage and diffusion model, introduce the mathematical model of atmospheric boundary layer characteristics into the liquid hydrogen leakage and diffusion model, and modify the boundary conditions and initial conditions in the model to improve the accuracy of the simulation.

[0054] Step 4: Based on the modified liquid hydrogen leakage and diffusion model in step 3, efficient numerical algorithms and optimization strategies are used to improve computational efficiency and convergence, and to handle large-scale computational domains and complex boundary conditions.

[0055] Step 5: Conduct simulation and verification. Simulate the leakage and diffusion of liquid hydrogen and verify it through experimental data or known results to ensure the accuracy and reliability of the simulation.

[0056] Step 6. Analyze the simulation results and extract key parameters and characteristics to provide a scientific basis for liquid hydrogen storage, transportation and safety assessment; key parameters and characteristics include diffusion distance, hydrogen concentration, lower flammable limit, hydrogen cloud volume, diffusion profile, etc.

[0057] Specifically, the net mass transfer rate in step 1 is solved as follows:

[0058] ;

[0059] in, is the mass transfer rate from liquid to gas, is the mass transfer rate from the gas phase to the liquid phase.

[0060] Specifically, the density in step 1 is solved by conjugate heat transfer to obtain the component Density ,density It is obtained from the mass equation and the advection equation, and the thermal diffusion coefficient is solved at the same time. The solution formula is as follows:

[0061] ;

[0062] ;

[0063] in, is the thermal diffusivity, is the thermal conductivity, is the turbulent thermal conductivity.

[0064] Specifically, in step 4, the residual is introduced based on the conjugate gradient method , virtual residual and virtual search directions , using the preconditioned stable biconjugate gradient method (PBiCGStab), satisfying:

[0065] ;

[0066] ;

[0067] ;

[0068] in, and is an integer not less than 0, and , is a symmetric positive definite matrix, and the superscript T indicates the transpose of the matrix.

[0069] Specifically, the step four includes: solving the coupled equations of the velocity field and the pressure field in step one, and solving them using the PIMPLE (a combination of an implicit operator partitioning algorithm for pressure and a semi-implicit method for a set of pressure coupled equations) algorithm prediction-correction strategy.

[0070] Figure 2 The Lee model is used to accurately describe and predict phase change phenomena, including the transition between liquid hydrogen and gaseous hydrogen. It is used to calculate the different source terms in the control equation of step 1 above, including the mass source term. Pressure Poisson source term and energy source term . By solving the phase transition Lee model to calculate the values ​​of different source terms, it provides a powerful tool for in-depth understanding and prediction of phase transition processes.

[0071] Perform mesh division, Figure 3This is a distribution diagram of the geometric structure grid in the embodiment. The uniform distribution of grid units in different areas is achieved by adjusting the segmentation ratio. Specifically, the grids in some areas are denser, while other areas are relatively sparse. This flexible grid division strategy is essential for capturing key flow characteristics in fluid dynamics, especially in areas that require high-precision simulation. The accuracy of the simulation can be improved by increasing the grid density. On the contrary, in areas where the flow changes more slowly, the number of grids can be appropriately reduced to improve computational efficiency and reduce the demand for computing resources. Before solving the calculation, the "geometric model" stage is first established. Then, the key "parameter setting" stage is entered, such as Figure 4 As shown, it includes the following sub-steps:

[0072] Step (1) Model parameter setting: According to the geometric model and simulation requirements, set the mesh parameters such as mesh size, shape, and number, and perform mesh division and verification.

[0073] Step (2) Physical property parameter setting: define the physical properties of the fluid, such as density, viscosity, temperature, etc.

[0074] Step (3) Boundary condition setting: According to the simulation objectives and actual conditions, set reasonable boundary conditions, such as inlet velocity, outlet pressure, etc.

[0075] Step (4) Numerical solution settings: Select appropriate solution algorithm, iteration steps, time step, convergence residual and other solution parameters. These settings together form the basis of numerical solution and affect the accuracy and computational efficiency of simulation results.

[0076] Step (5) Initialization settings: Initialization requires assumptions about parameters such as the initial velocity distribution and pressure distribution, which helps accelerate the convergence of the flow. If a better initial solution is used at the beginning of the calculation, the convergence state can be reached more effectively.

[0077] The calculation results are post-processed. Figure 5 shown. Figure 5 The figure shows the distribution of hydrogen cloud concentration. It can be clearly observed that driven by atmospheric wind, hydrogen clouds tend to diffuse to the upper right. This diffusion phenomenon is the result of the synergistic effect of wind and buoyancy. Wind, as the main driving force, guides hydrogen clouds to migrate along the wind direction; while the buoyancy effect, because the density of hydrogen is less than that of air, causes the clouds to gradually diffuse during the rising process. The combined effect of these two mechanisms shapes the diffusion pattern and path of hydrogen clouds in the atmosphere.

[0078] The present invention also provides a liquid hydrogen leakage and diffusion solving device based on the atmospheric boundary layer inlet velocity, comprising the following modules:

[0079] Model building module, which builds a liquid hydrogen leakage and diffusion model based on the physical properties and leakage scenarios of liquid hydrogen, including the mass conservation equation, momentum conservation equation, and energy conservation equation;

[0080] The parameter acquisition module determines the characteristics of the atmospheric boundary layer and obtains the relevant physical parameters of the atmospheric boundary layer inlet velocity, including reference altitude, average streamwise wind speed, and gas dynamic roughness, through meteorological data or field measurements;

[0081] The correction module combines the atmospheric boundary layer characteristics with the liquid hydrogen leakage and diffusion model, introduces the atmospheric boundary layer characteristics mathematical model into the liquid hydrogen leakage and diffusion model, and corrects the boundary conditions and initial conditions in the liquid hydrogen leakage and diffusion model to improve the accuracy of the simulation;

[0082] The boundary processing module, based on the modified liquid hydrogen leakage diffusion model, adopts numerical algorithms and optimization strategies to improve computational efficiency and convergence, and handle large-scale computational domains and complex boundary conditions;

[0083] The simulation verification module performs simulation and verification, simulates the leakage and diffusion of liquid hydrogen, and verifies it through experimental data or known results to ensure the accuracy and reliability of the simulation;

[0084] The result analysis module analyzes the simulation results and extracts key parameters and features to provide a scientific basis for liquid hydrogen storage, transportation and safety assessment; key parameters and features include dispersion distance, hydrogen concentration, lower flammable limit, hydrogen cloud volume, and diffusion profile.

[0085] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the above-mentioned method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity are implemented.

[0086] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-mentioned method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity.

[0087] It should be understood by those skilled in the art that the embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes. The solutions in the embodiments of the present invention may be implemented in various computer languages, for example, object-oriented programming language Java and interpreted scripting language JavaScript, etc.

[0088] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0089] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0090] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0091] Although the preferred embodiments of the present invention have been described, those skilled in the art may make other changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present invention.

[0092] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

Claims

1. A method for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity, characterized in that: The steps include: Step 1: According to the physical properties of liquid hydrogen and the leakage scenario, a liquid hydrogen leakage diffusion model is established. The liquid hydrogen leakage diffusion model includes the mass conservation equation, momentum conservation equation and energy conservation equation; Step 2: Determine the characteristics of the atmospheric boundary layer, and obtain relevant physical parameters of the atmospheric boundary layer inlet velocity through meteorological data or field measurement data, wherein the relevant physical parameters include reference height, average stream wind speed, and gas dynamic roughness; Step 3: Combining the atmospheric boundary layer characteristics with the liquid hydrogen leakage and diffusion model, the atmospheric boundary layer characteristics mathematical model is introduced into the liquid hydrogen leakage and diffusion model, and the boundary conditions and initial conditions in the liquid hydrogen leakage and diffusion model are corrected; Step 4: Based on the modified liquid hydrogen leakage diffusion model in step 3, numerical algorithms and optimization strategies are used to improve computational efficiency and convergence, and to handle large-scale computational domains and complex boundary conditions; Step 5: simulate the leakage and diffusion of liquid hydrogen to obtain simulation results, and verify them through experimental data or known results to ensure the accuracy and reliability of the simulation; Step 6: Analyze the simulation results and extract key parameters and features to provide a scientific basis for liquid hydrogen storage, transportation and safety assessment; the key parameters and features include diffusion distance, hydrogen concentration, lower flammable limit, hydrogen cloud volume, and diffusion profile.

2. The method for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity according to claim 1 is characterized in that: The equations of the liquid hydrogen leakage diffusion model in step 1 include: The mass equation is: ； in, is the divergence operator, is the density, is the velocity vector, is the net mass transfer rate; Components The interface advection equation is: ； in, It is a component The volume fraction of It is a component The density of It is a component The mass transfer rate, is the mass source term, t represents time; ； in, is the correction pressure, is the dynamic viscosity, the superscript T indicates the transpose of the matrix, is the acceleration due to gravity, is the surface tension; ； in, is a variable constructed from the intermediate velocity field, are the diagonal terms of the momentum matrix equation, is the density of the liquid phase, is the density of the gas phase, is about the pressure Poisson source term; ； in, is the specific heat at constant pressure, is the temperature, is the effective thermal conductivity, is the energy source term, It is latent heat.

3. The method for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity according to claim 1 is characterized in that: The calculation formula of step 2 includes: ； ； in, is the surface normal streamwise velocity distribution, is the surface normal turbulent kinetic energy profile, is the ground normal coordinate component, is the ground normal displacement height, is the aerodynamic roughness, is the friction speed, is the reference altitude, is the reference height The average streamwise wind speed at .

4. The method for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity according to claim 2 is characterized in that: The net mass transfer rate in step 1 The calculation formula is: ； in, is the mass transfer rate from liquid to gas, is the mass transfer rate from the gas phase to the liquid phase.

5. The method for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity according to claim 2 is characterized in that: In step 1, the density is solved by conjugate heat transfer to obtain the composition Density ,density Obtained from the mass equation and the interface advection equation, the thermal diffusion coefficient is solved at the same time , the solution formula is as follows: ； ； in, is the thermal diffusivity, is the thermal conductivity, is the turbulent thermal conductivity.

6. The method for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity according to claim 1 is characterized in that: The fourth step comprises: Introducing residuals based on the conjugate gradient method , virtual residual and virtual search directions , using the pre-conditioned stable biconjugate gradient method, satisfying: ； ； ； in, and is an integer not less than 0, and , is a symmetric positive definite matrix.

7. The method for solving the leakage and diffusion of liquid hydrogen based on the atmospheric boundary layer inlet velocity according to claim 2 is characterized in that: The step 4 uses the PIMPLE algorithm prediction-correction strategy to solve the equation of the liquid hydrogen leakage diffusion model in step 1.

8. A liquid hydrogen leakage and diffusion solution device based on the atmospheric boundary layer inlet velocity, characterized in that: Includes the following modules: The model building module builds a liquid hydrogen leakage and diffusion model based on the physical properties of liquid hydrogen and the leakage scenario. The liquid hydrogen leakage and diffusion model includes the mass conservation equation, momentum conservation equation and energy conservation equation. A parameter acquisition module determines the characteristics of the atmospheric boundary layer and obtains relevant physical parameters of the atmospheric boundary layer inlet velocity through meteorological data or field measurement data, wherein the relevant physical parameters include reference height, average stream wind speed, and gas dynamic roughness; The correction module combines the atmospheric boundary layer characteristics with the liquid hydrogen leakage and diffusion model, introduces the atmospheric boundary layer characteristics mathematical model into the liquid hydrogen leakage and diffusion model, and corrects the boundary conditions and initial conditions in the liquid hydrogen leakage and diffusion model; The boundary processing module, based on the modified liquid hydrogen leakage diffusion model, adopts numerical algorithms and optimization strategies to improve computational efficiency and convergence, and handle large-scale computational domains and complex boundary conditions; The simulation verification module simulates the leakage and diffusion of liquid hydrogen, obtains simulation results, and verifies them through experimental data or known results to ensure the accuracy and reliability of the simulation; The result analysis module analyzes the simulation results and extracts key parameters and features to provide a scientific basis for liquid hydrogen storage, transportation and safety assessment; the key parameters and features include dispersion distance, hydrogen concentration, lower flammable limit, hydrogen cloud volume, and diffusion profile.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the steps of a method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity according to any one of claims 1 to 7 are implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of a method for solving liquid hydrogen leakage and diffusion based on the atmospheric boundary layer inlet velocity as described in any one of claims 1 to 7 are implemented.

Images