Simulation method, device, equipment and medium for hydraulic characteristics of tongue flap gate

Abstract

The invention discloses an analogue simulation method, device and equipment for hydraulic characteristics of a tongue flap gate and a storage medium. The analogue simulation method comprises the steps that a tongue flap gate three-dimensional geometric model containing key parameters is constructed, a computational domain is divided, and a proper grid division technology is adopted to complete grid independence verification; establishing a control equation system of the tongue flap gate based on the three-dimensional geometric model of the tongue flap gate, setting boundary conditions, and carrying out numerical solution on hydraulic characteristics of the tongue flap gate; carrying out postprocessing analysis on the tongue flap gate according to the numerical value of the hydraulic characteristics of the tongue flap gate to obtain the hydraulic characteristics, and carrying out multi-working-condition comparative analysis; the simulation model is verified through a physical model test, structural parameter optimization design is performed based on a verification result, a mapping relation model is established to realize intelligent operation regulation and control, and meanwhile, simulation calculation is accelerated by adopting parallel calculation and a prediction model. According to the method, high-precision simulation of the hydraulic characteristics of the tongue flap gate is achieved, the calculation efficiency is remarkably improved, intelligent regulation and control and multi-objective optimization are supported, and the engineering safety and economy are remarkably enhanced.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of water conservancy engineering, and in particular to a simulation method, device, equipment and medium for the hydraulics of a flap gate. BACKGROUND

[0002] In the field of water conservancy engineering, the flap gate is the core equipment of the layered water intake system, and its hydraulics directly determines the water intake efficiency and equipment safety. The traditional design and performance analysis of the flap gate mainly rely on empirical formulas, physical model tests and simple computational fluid dynamics (CFD) methods. The empirical formulas are usually based on specific working conditions and are difficult to adapt to complex and variable actual operating conditions, resulting in a large deviation between the design results and the actual requirements; although the physical model test can intuitively reflect the working state of the flap gate, the test cost is high and the test period is long, and there is a similarity error between the scaled model and the actual engineering structure, which cannot accurately simulate the real working conditions.

[0003] The existing CFD method also has significant shortcomings, which is mostly limited to single-phase flow field analysis and fails to fully consider the coupling effects of door body movement, cavitation effect and free surface flow and other multi-physical fields, resulting in a simulation error of more than 15% in the flow characteristics, a lag in the identification of cavitation risk, and an inability to provide reliable basis for engineering design. At the same time, in terms of structural parameter design, there is a lack of systematic multi-parameter coupling analysis means, and the design period is as long as 3-6 months, which cannot meet the needs of modern water conservancy engineering for rapid design and optimization. In addition, the application of digital twin model in the operation and regulation of the flap gate is not mature, and due to the lack of high-precision hydraulics data support, it is unable to realize real-time and accurate matching of opening, flow and pressure, and the response time of regulation is more than 30 minutes, which seriously affects the operation efficiency and safety of water conservancy engineering. Therefore, there is an urgent need for a high-precision, parameterized and intelligent simulation method for the hydraulics of the flap gate to solve the many problems existing in the prior art. SUMMARY

[0004] The main purpose of the present application is to provide a simulation method, device, equipment and medium for the hydraulics of a flap gate, which can solve the technical problems of incomplete monitoring, fixed threshold and lack of intelligence in emergency handling existing in the prior art.

[0005] To achieve the above-mentioned purpose, the present application provides the following technical solution: a simulation method for the hydraulics of a flap gate, comprising: constructing a three-dimensional geometric model of the flap gate containing key parameters and dividing the calculation domain, adopting appropriate grid division technology and completing grid independence verification; establishing a control equation system of the flap gate based on the three-dimensional geometric model of the flap gate and setting boundary conditions, and performing numerical solution of the hydraulics of the flap gate; Post-processing analysis of the tongue valve valve is carried out based on the numerical values ​​of its hydraulic characteristics to obtain the hydraulic properties, and multi-condition comparative analysis is performed. The simulation model is verified through physical model experiments. Based on the verification results, the structural parameters are optimized and a mapping relationship model is established to achieve intelligent operation control. At the same time, parallel computing and prediction models are used to accelerate the simulation calculation.

[0006] As a further improvement to this application, the steps of constructing a three-dimensional geometric model containing key parameters and dividing the computational domain are as follows: A 3D modeling software was used to create a model including the opening of the lingual valve. Door body curvature Doorway width A three-dimensional parametric geometric model, wherein, This indicates the minimum opening of the lingual valve. Indicates the maximum opening of the lingual valve; This represents the minimum radius of curvature of the door body surface. Indicates the maximum curvature of the door body surface; Indicates the minimum width of the door slot. Indicates the maximum width of the door slot; Divide the computational domain and set the length of the upstream inlet section as follows: The gate height is set to a multiple of μ, and the downstream outlet section length is set to μ times the gate height. Boolean operations are used to integrate all components into a full flow channel geometric model to ensure the continuity and sealing of the flow channel. The coefficients for the multiple relationship between the upstream inlet section length and the gate height, as well as the multiple relationship coefficient between the downstream outlet section length and the gate height, are also set.

[0007] As a further improvement to this application, the steps of employing a suitable mesh generation technique and completing the mesh independence verification are as follows: A tetrahedral-hexahedral hybrid mesh generation technique is used to generate meshes in the door body and door slot areas. Tetrahedral unstructured mesh and refined boundary layer to Upstream and downstream channels adopt The hexahedral structured mesh, This represents the minimum size of the tetrahedral unstructured mesh for the door body and door slot area. Indicates the maximum size; The dimensionless upper limit of the distance for boundary layer encryption control; This represents the minimum size of the hexahedral structured mesh for the upstream and downstream flow channels. Indicates the maximum size; Control the total number of grids in Preliminary calculations were performed on key hydraulic parameters under different grid densities. This indicates the lower limit of the total number of grid cells. Indicates the maximum number of grid cells; Grid independence verification is carried out, when the number of grids increases by a certain percentage, the change of key hydraulic parameters is less than the set threshold, then it is considered that the grid division scheme meets the requirements.

[0008] As a further improvement of the application, the step of establishing control equation system and setting boundary conditions is as follows: The continuity equation is established as u=0, momentum equation p u t+u u p u u)+F and RNGk- The core of the control equation system is the turbulence model, and the standard wall function is used in the near-wall region; The VOF model is used to track the gas-liquid two-phase interface, the surface tension coefficient is set as κ, the gravitational acceleration is set as q, and the vaporization pressure calculation formula based on water temperature T (°C) is constructed pv T T)) to realize the simulation of cavitation effect, and к is the surface tension coefficient; ζ, θ are coefficients in the vaporization pressure calculation formula; The inlet is set as a velocity inlet , represents the lower limit of the inlet velocity, represents the upper limit of the inlet velocity

[0009] The grid angular velocity ω , represents the lower limit of the angular velocity of the tongue valve around the axis, represents the upper limit of the angular velocity.

[0010] As a further improvement of the application, the step of numerical solution is as follows: The ANSYS Fluent solver is selected, the SIMPLE algorithm is used for pressure-velocity coupling calculation, the second-order upwind format is used for the discretization of the convection term, and the central difference format is used for the diffusion term; The residual convergence criterion is set as the residual curve of each equation being lower than ν, while monitoring the parameters such as flow rate and force on the door body, when the parameter fluctuation amplitude is ν<φ, it is determined that the solution converges, ν is the upper limit value of the residual convergence criterion, and φ is the upper limit value of the parameter fluctuation amplitude; Steady-state and transient-state solutions are carried out respectively, and key data in the solution process are recorded to provide a basis for subsequent analysis.

[0011] As a further improvement of the application, the analysis verification and scheme optimization: Analysis and verification of the scheme optimization: characteristic calculation Flow coefficient wherein, is the flow rate, is the gate height, is the gravitational acceleration, and the over-flow capacity is analyzed; Cavitation index is the local pressure, is the flow rate, and the high-risk area of is identified; Hydraulic torque is the distance from the surface element of the gate body to the rotating shaft, is the surface area of the gate body, and the parameters of the hoistway gate are optimized; Multi-working condition comparison: set water level , in the range, is the lowest water level, is the highest water level, opening , water temperature in the range, is the lowest water temperature, is the highest water temperature, and the orthogonal working conditions are considered for temperature correction, and a group of working condition databases are constructed to analyze the parameter influence law; The model verification includes comparison of simulation and physical test data, and the over-flow flow error is the upper limit of the flow error, and the pressure coincidence degree is the lower limit of the pressure coincidence degree, and the cavitation position deviation is the upper limit of the cavitation position deviation, the parameters are adjusted to meet the water conservancy standards, and the NSGA-II algorithm is used to optimize the modeling parameters of the flap gate, and the objective function of the NSGA-II algorithm can be the over-flow efficiency, cavitation index, and hoist torque. After verification, the optimal parameter combination is determined.

[0012] As a further improvement of the present application, the intelligent hole control and accelerated calculation include: Intelligent hole control: based on a group of simulation data, a mapping model is established through a BP neural network and embedded in a digital twin platform; multiple working condition logics are set to realize the opening regulation error , is the upper limit of the opening regulation error; Accelerated calculation: set CUDA and improve the calculation efficiency by more than 8 times, introduce PINN to fuse CFD data and control equations, and the prediction speed is improved times, is the lower limit of the speed improvement times, and the precision error is the upper limit of the precision error, and a parameterized script is developed to support batch processing.

[0013] To achieve the above object, the application further provides the following technical scheme: A tongue flap gate hydraulics characteristic simulation device is applied to a tongue flap gate hydraulics characteristic simulation method, and the tongue flap gate hydraulics characteristic simulation device comprises: A modeling unit is configured to construct a tongue flap gate three-dimensional geometric model containing key parameters and divide a calculation domain, adopt a suitable mesh division technology, and complete mesh independence verification; A solving unit is configured to establish a control equation system of the tongue flap gate based on the tongue flap gate three-dimensional geometric model, set a boundary condition, and perform numerical solving of the tongue flap gate hydraulics characteristic; An analysis unit is configured to perform post-processing analysis of the tongue flap gate based on the numerical tongue flap gate hydraulics characteristic to obtain the hydraulics characteristic, and perform multi-working condition comparative analysis; A simulation unit is configured to verify a simulation model through a physical model test, perform structure parameter optimization design based on a verification result, establish a mapping relationship model to realize intelligent operation regulation and control, and simultaneously adopt parallel computing and a prediction model to accelerate simulation calculation.

[0014] To achieve the above object, the application further provides the following technical scheme: An electronic device comprises a processor and a memory coupled with the processor, the memory storing program instructions executable by the processor; the processor executes the program instructions stored in the memory to implement the tongue flap gate hydraulics characteristic simulation method.

[0015] To achieve the above object, the application further provides the following technical scheme: A storage medium stores program instructions, and the program instructions are executed by a processor to implement the tongue flap gate hydraulics characteristic simulation method.

[0016] Beneficial effects: Compared with the prior art, the application adopts a distributed optical fiber temperature sensor to collect vertical water temperature data, can comprehensively and accurately obtain water temperature information at different depths of the reservoir, and overcomes the limitations of traditional single-point monitoring. The high-precision and high-resolution characteristics of the sensor enable the timely capture of slight water temperature changes, providing a reliable data basis for subsequent anomaly identification and regulation. Combined with the data of the water quality sensor and the weather station, the fusion of multi-source data is realized, various factors affecting the water temperature are comprehensively considered, and the overall grasp of the water temperature condition is further improved.

[0017] The application adopts three-dimensional parameterized modeling technology, constructs a full-flow passage model of the flap gate in three-dimensional software, accurately sets and flexibly adjusts key parameters such as the flap gate opening degree, the gate body curvature, and the gate slot width, cooperates with scientific calculation domain division and hybrid mesh optimization technology, and ensures that the model can accurately reflect the actual working condition. At the same time, the multi-physical field coupling control equation is used, combined with the turbulence model, the VOF model and the cavitation model, the factors such as water flow, gate body movement, free surface and cavitation effect are comprehensively considered, and high-precision simulation of the hydraulic characteristics of the flap gate is realized. The comprehensive use of such technical means enables in-depth analysis of the influence of different parameters on the performance of the equipment in the design stage, avoids structural defects caused by rough simulation in traditional design, significantly improves the flow efficiency and structural safety of the flap gate, and provides a solid technical support for the optimized design of the core equipment of the water conservancy project.

[0018] The application also uses the Unity to develop a digital twin platform, embeds the opening-flow-pressure relationship atlas obtained through simulation in the platform, and establishes an accurate mapping model. When the system obtains working condition change information through sensors such as water temperature monitors, it can quickly trigger the control logic, automatically calculate and output the optimal flap gate opening degree, and link the gate hoist to accurately adjust. At the same time, CUDA parallel computing and physical information neural network (PINN) technology are introduced, which greatly improves the simulation calculation speed, significantly shortens the single working condition simulation time, and the PINN model can realize fast prediction on the premise of ensuring accuracy. The integration of these technical means enables the flap gate to respond to working condition changes in real time and intelligently during the operation of the water conservancy project, significantly shortens the control response time, reduces manual intervention, improves the intelligentization and automation level of the water conservancy project operation, and ensures stable and efficient operation of the equipment.

[0019] The application also uses the non-dominated sorting genetic algorithm (NSGA-II) to perform multi-objective optimization on the structure parameters of the flap gate. In the embodiment, by adjusting parameters such as the gate body curvature and the gate slot chamfer, the over-flow efficiency, cavitation risk, and gate hoist energy consumption are effectively balanced, the engineering design and operation cost are reduced, and the service life of the equipment is prolonged. In addition, with the help of high-precision simulation and intelligent control technology, the ecological water demand can be accurately met, such as accurately controlling the water temperature during the fish spawning period, effectively protecting the downstream aquatic habitat, and realizing the win-win of economic benefit and ecological benefit of the water conservancy project. At the same time, multi-working condition comparative analysis and model verification technology ensure the reliability and applicability of the method under different environmental conditions, and provide strong technical support for the scientific management and sustainable utilization of water resources. BRIEF DESCRIPTION OF DRAWINGS

[0020] Figure 1 It is a step flowchart for an embodiment of the flap gate hydraulic characteristics simulation method of the application; Figure 2 A functional module schematic diagram of an embodiment of the tongue flap gate hydrodynamic characteristic simulation method of the application; Figure 3 A structural schematic diagram of an embodiment of the electronic device of the application. DETAILED DESCRIPTION

[0021] The technical solutions in the embodiments of the application will be clearly and completely described below with reference to the drawings in the embodiments of the application. Obviously, the described embodiments are only a part of the embodiments of the application, but not all the embodiments of the application. Based on the embodiments in the application, all other embodiments obtained by a person of ordinary skill in the art without creative work are within the protection scope of the application.

[0022] As Figure 1 shown, to solve the above technical problems, the application provides a tongue flap gate hydrodynamic characteristic simulation method, comprising: S1, constructing a three-dimensional geometric model of the tongue flap gate containing key parameters and dividing a calculation domain, adopting a suitable mesh division technology and completing mesh independence verification; S2, establishing a control equation system of the tongue flap gate based on the three-dimensional geometric model of the tongue flap gate and setting boundary conditions, and performing numerical solution of the tongue flap gate hydrodynamic characteristics; S3, performing post-processing analysis of the tongue flap gate according to the numerical solution of the tongue flap gate hydrodynamic characteristics to obtain the hydrodynamic characteristics, and performing multi-working condition comparative analysis; S4, verifying the simulation model through a physical model test, performing structural parameter optimization design based on the verification result and establishing a mapping relationship model to realize intelligent operation regulation and control, and simultaneously adopting parallel computing and a prediction model to accelerate the simulation calculation.

[0023] As a further improvement of the application, the step of constructing a three-dimensional geometric model containing key parameters and dividing a calculation domain is as follows: A three-dimensional parameterized geometric model containing tongue flap gate opening , door body surface radian , gate slot width is established by adopting a three-dimensional modeling software, wherein, represents the minimum opening of the tongue flap gate, represents the maximum opening of the tongue flap gate; represents the minimum radian of the door body surface, represents the maximum radian of the door body surface; represents the minimum width of the gate slot, represents the maximum width of the gate slot; the calculation domain is divided, and the length of the upstream inlet section is set as The gate height is set to a multiple of μ, and the downstream outlet section length is set to μ times the gate height. Boolean operations are used to integrate all components into a full flow channel geometric model to ensure the continuity and sealing of the flow channel. The coefficients for the multiple relationship between the upstream inlet section length and the gate height, as well as the multiple relationship coefficient between the downstream outlet section length and the gate height, are also set.

[0024] A dynamically adjustable parametric model is established to enable flexible control over core structural parameters such as the flap opening, door curvature, and door slot width. Parametric design (e.g., SolidWorks modeling) supports the rapid generation of different structural schemes, avoiding redundant modeling and improving design iteration efficiency. Clearly defined assembly relationships and dimensional constraints ensure consistency between the geometric model and the actual engineering structure, providing accurate physical boundaries for subsequent flow field calculations.

[0025] As a further improvement to this application, the following steps are taken to adopt a suitable mesh generation technique and complete the mesh independence verification: A tetrahedral-hexahedral hybrid mesh generation technique is used to generate meshes in the door body and door slot areas. Tetrahedral unstructured mesh and refined boundary layer to Upstream and downstream channels adopt The hexahedral structured mesh, This represents the minimum size of the tetrahedral unstructured mesh for the door body and door slot area. Indicates the maximum size; The dimensionless upper limit of the distance for boundary layer encryption control; This represents the minimum size of the hexahedral structured mesh for the upstream and downstream flow channels. Indicates the maximum size; controls the total number of grid cells. Preliminary calculations were performed on key hydraulic parameters under different grid densities. This indicates the lower limit of the total number of grid cells. This indicates the upper limit of the total number of grids; grid independence verification is performed. If the change in key hydraulic parameters is less than the set threshold when the number of grids increases by a certain percentage, then the grid division scheme is considered to meet the requirements.

[0026] For complex flow fields (such as boundary layers and vortex regions) and regular flow channels near the flap gate, a differentiated meshing strategy is adopted to balance computational accuracy and efficiency. Unstructured meshes (tetrahedral) in the gate region adapt to complex curved surfaces, and denser boundary layers accurately capture wall flow characteristics; structured meshes (hexahedral) in the flow channel reduce the number of meshes and improve computational speed; mesh independence verification ensures that the calculation results are independent of mesh density, avoiding simulation deviations caused by mesh errors.

[0027] As a further improvement to this application, the steps for establishing the governing equation system and setting boundary conditions are as follows: Establish a continuity equation system... u=0, momentum equation p u t+u u p u u)+F and RNGk- The core of the control equation system is the turbulence model, and the standard wall function is used to process the near-wall region. The VOF model is used to track the gas-liquid interface, the surface tension coefficient is set as κ, the gravitational acceleration is set as q, and the vaporization pressure calculation formula based on water temperature T(°C) is constructed pv T T)) to realize the simulation of cavitation effect, and к is the surface tension coefficient; ζ, θ are the coefficients in the vaporization pressure calculation formula; the inlet is set as a velocity inlet , , which represents the lower limit of the inlet velocity, , which represents the upper limit of the inlet velocity, and the outlet is a pressure outlet , , which represents the lower limit of the outlet pressure, , which represents the upper limit of the outlet pressure; the flap valve is defined as the grid angular velocity ω , , which represents the lower limit of the angular velocity of the flap valve rotating around the axis, , which represents the upper limit of the angular velocity.

[0028] A multi-physics coupled mathematical model is constructed to simulate the real flow environment (such as turbulence, cavitation, and gas-liquid two-phase flow). The turbulence model accurately simulates complex turbulent flow, and the standard wall function is used to process the near-wall region flow. The VOF model tracks the gas-liquid interface, and the cavitation effect is quantitatively analyzed by combining the vaporization pressure formula. The dynamic mesh technology is used to simulate the rotating motion of the flap valve, and supports dynamic condition (such as opening and closing process) simulation.

[0029] As a further improvement of the present application, the numerical solution steps are as follows: the ANSYS Fluent solver is selected, the SIMPLE algorithm is used for pressure-velocity coupling calculation, the second-order upwind format is used for discretization of the convection term, and the central difference format is used for the diffusion term; the residual convergence criterion is set as the residual curve of each equation being lower than ν, while monitoring the flow rate, the force on the door body, and other parameters; when the parameter fluctuation amplitude is ν<φ, it is determined that the solution converges, ν is the upper limit value of the residual convergence criterion, and φ is the upper limit value of the parameter fluctuation amplitude; steady-state and transient-state solutions are performed respectively, and key data in the solution process are recorded to provide a basis for subsequent analysis.

[0030] The control equation is solved by numerical method to obtain the key parameters of the flow field (pressure, flow rate, gas holdup, etc.). Technical value: the SIMPLE algorithm efficiently handles pressure-velocity coupling, the second-order upwind scheme improves the accuracy of the convection term discretization and reduces numerical diffusion; the residual and parameter fluctuation double convergence criteria ensure the stability of the solution and avoid false convergence. Covering static conditions (steady state) and dynamic processes (transient state) to meet different analysis needs. Steady-state solution is used to obtain the characteristics of the stable flow field (such as the flow capacity under rated opening); transient solution captures dynamic processes (such as pressure pulsation and vortex evolution during opening and closing), providing data for structural dynamics analysis.

[0031] As a further improvement of the present application, analysis verification and scheme optimization: Analysis verification and scheme optimization: characteristic calculation: flow coefficient , wherein, is the flow rate, is the gate height, is the gravitational acceleration, analysis of flow capacity; cavitation index is the local pressure, is the flow rate, identify high-risk areas; hydrodynamic moment is the distance from the surface element of the gate body to the rotating shaft, is the surface area of the gate body, optimize the parameters of the opening and closing machine; Multi-condition comparison: set water level , in the range, is the lowest water level, is the highest water level, opening , water temperature in the range, is the lowest water temperature, is the highest water temperature, orthogonal working condition, considering temperature correction, constructing a database of groups of working conditions, analyzing the influence of parameters; Model verification includes comparing simulation and physical test data, calculating the flow error is the upper limit of the flow error, the pressure coincidence degree is the lower limit of the pressure coincidence degree, the cavitation position deviation is the upper limit of the cavitation position deviation, adjust the parameters to meet the water conservancy standards, and use the NSGA-II algorithm to optimize the modeling parameters of the flap gate. The objective function of the NSGA-II algorithm can be flow efficiency, cavitation index, and opening and closing moment. After verification, the optimal parameter combination is determined.

[0032] Quantify the key performance indicators of the flap gate, and reveal the relationship between the flow field characteristics and the structure parameters. The flow coefficient evaluates the flow efficiency, and provides pressure data for cavitation risk assessment; the hydrodynamic moment M provides key mechanical parameters for the selection of the opening and closing machine and the design of the driving system. Further, in the comparison of multiple working conditions and the construction of a database, the system analyzes the influence of different working conditions (water level, opening degree, water temperature) on the performance of the flap gate, and constructs a standardized data asset. The orthogonal test design (such as water level and opening degree) covers the actual operation range of the project; the database supports parameter sensitivity analysis, and identifies key influencing factors, such as the influence of water temperature on the cavitation threshold. The reliability of the simulation model is verified by physical test data, ensuring the engineering credibility of the numerical simulation results. The flow error, pressure coincidence degree, and other indicators quantify the precision, which meets the water conservancy industry standards (such as SL155-2020); the corrected model can be used as a reliable basis for design optimization, avoiding the problem of "garbage in, garbage out". Based on the simulation data, the optimal structure parameter combination is found through intelligent algorithms, balancing multiple performance indicators. Technical value: NSGA-II algorithm handles the conflict of multiple objectives such as flow efficiency, cavitation risk, and opening and closing torque, generating a Pareto solution set; the optimized structure can significantly improve the engineering performance (among them, the flow efficiency is improved by 12%, and the cavitation index is reduced by 18%).

[0033] As a further improvement of the present application, intelligent hole control and accelerated calculation include: Intelligent hole control: based on Simulation data, a mapping model is established through a BP neural network and embedded in a digital twin platform; multiple working condition logics are set to achieve the upper limit of the opening degree control error , for the opening degree control error; accelerated calculation: set CUDA and improve the calculation efficiency by more than 8 times, introduce PINN to integrate CFD data and control equations, and improve the prediction speed times, is the lower limit of the speed improvement, and the precision error is the upper limit of the precision error, and a parameterized script is developed to support batch processing.

[0034] Establish a real-time mapping between "simulation model → actual device" to realize the automatic adjustment of the flap gate running state. The mapping relationship model (such as BP neural network) quickly matches the optimal opening degree, and the response time is shortened from 30 minutes of manual intervention to seconds; the working condition trigger logic (such as water temperature < T3) automatically adjusts the opening degree, which improves the engineering safety and avoids cavitation damage. Break through the traditional CFD calculation efficiency bottleneck, support large-scale parameterized simulation and real-time prediction. Technical value: CUDA parallel calculation utilizes GPU computing power, single machine efficiency is improved by more than 8 times, and the solution time of complex working conditions is compressed from 48 hours to 6 hours; PINN model integrates physical constraints and data-driven, prediction speed is improved by 5 times, precision error is less than 5%, suitable for real-time control scenarios.

[0035] The following is an example of a tongue flap gate in a certain water conservancy project, which will explain in detail the implementation process of the simulation method of the tongue flap gate hydraulics characteristics: The water conservancy project is located on a river, mainly for flood control, irrigation and power generation. The tongue flap gate as the key water discharge equipment of the project, its hydraulics characteristics directly affect the operation efficiency and safety of the project. In order to optimize the design and operation of the tongue flap gate, it is decided to use the simulation method of the invention for research.

[0036] A three-dimensional parametric geometric model of the tongue flap gate is constructed using professional three-dimensional modeling software SolidWorks. According to the design requirements of the project, the opening range of the tongue flap gate is set to (ie , , ), the range of the door body curvature is (ie , , ), and the range of the door slot width is (ie , , ). The assembly relationship and size constraints of each component are determined, and the full flow channel model is generated by Boolean operation. The length of the upstream inlet section is set to 4 times the height of the door , and the length of the downstream outlet section is set to 8 times the height of the door .

[0037] Meshing is performed using ANSYS ICEM CFD software. In the door body and door slot area, tetrahedral unstructured mesh with size range of (0.02m-0.06m) (ie a =0.02m, b =0.06m) is used, and the boundary layer is encrypted on the wall surface, so that (y + <25) (c=25). The upstream and downstream flow channels use hexahedral structured mesh with size range of (0.06m-0.25m) (ie (d=0.06m), (e=0.25m)), and the total number of meshes is controlled between 1.5 million and 2.5 million (ie (f=1500000), (g=2500000)). Mesh independence verification is performed by changing the mesh density to ensure that the key hydraulic parameters change less than (3%) under different mesh densities.

[0038] In the ANSYS Fluent software, the control equations and boundary conditions are set. The continuity equation , momentum equation and RNG turbulence model are used, and the standard wall function is used for near-wall region treatment. The VOF model is used to track the gas-liquid interface, and the surface tension coefficient Based on water temperature (assuming the water temperature is in this simulation), the vaporization pressure is calculated by the formula . The inlet is set as a velocity inlet with a velocity range of (i.e. , ), the outlet is set as a pressure outlet with a pressure range of (i.e. , ), and the door body is defined as a moving mesh rotating around an axis with an angular velocity range of (i.e. , ).

[0039] ANSYS Fluent is selected as the solver, the SIMPLE algorithm is used for pressure-velocity coupling calculation, the second-order upwind scheme is used for discretization of the convection term, and the central difference scheme is used for the diffusion term. The residual convergence criterion is set to be small -6 (i.e. h =10 -6 ), real-time monitoring of flow rate, door force and other parameters, and when the parameter fluctuation amplitude is less than i (i.e. =1.5%) is determined to be convergent. Steady and transient solutions are performed respectively to obtain key data such as flow field pressure, flow rate, and gas holdup.

[0040] Flow capacity analysis: calculate the flow coefficient , draw the opening-flow coefficient curve. Through analysis of the curve, it is found that when the opening is about , the flow coefficient reaches the maximum value, indicating that the flow capacity of the flap door is the strongest at this time. Cavitation analysis: extract the door body surface pressure cloud picture, calculate the cavitation index , and identify the high-risk cavitation area near the door slot corner and the bottom of the door body. Hydrodynamic moment calculation: visualize the door slot and the flow field vector diagram behind the door, and calculate the hydrodynamic moment . The calculation results show that the hydrodynamic moment increases with the increase of the opening, and reaches the maximum value when the opening is .

[0041] A number of typical working conditions are set, and the working condition parameters cover water level (i.e. , , ), opening (20°, 40°, , i.e. , , ), water temperature (15° i.e. , ), considering the influence of temperature on water density and viscosity , simulation calculation was carried out for each group of working conditions to construct the working condition database. By comparing and analyzing the differences in hydraulic characteristics under different working conditions, it was found that the increase of water temperature would lead to the increase of vaporization pressure, thereby increasing the risk of cavitation; the increase of water level would increase the flow rate, but would also increase the water power moment.

[0042] By comparing the simulation results with the physical model test data, the error of flow rate was (meeting the error i.e. ), the pressure distribution coincidence degree was (meeting the coincidence degree i.e. ), and the cavitation position deviation was (meeting the deviation i.e. ). According to the comparison results, some parameters of the simulation model were fine-tuned, such as adjusting the parameters of the turbulence model, to further improve the simulation accuracy.

[0043] Based on the verification results, the non-dominated sorting genetic algorithm (NSGA-II) was used to optimize the gate structure parameters (such as gate arc, gate groove chamfer, etc.) for multiple objectives, such as improving the flow efficiency, reducing the cavitation index, and reducing the opening and closing moment. After multiple iterations, a set of optimal structure parameter combination was obtained: the gate arc , and the gate groove chamfer radius was . By comparing the performance indicators before and after optimization, the flow efficiency was improved by , the cavitation index was reduced by 15%, and the opening and closing moment was reduced by 12%.

[0044] Based on 1200 groups of simulation data, a BP neural network was established using Python combined with TensorFlow library to establish the mapping relationship model of the gate opening , flow rate , and pressure . The model was embedded in the digital twin platform, and the working condition triggering logic was set. When the water temperature (i.e. ), the opening was automatically adjusted to , and the opening adjustment error was controlled within (i.e. ).

[0045] CUDA parallel computing technology was employed, distributing computational tasks across an NVIDIA Tesla P100 GPU cluster. Simultaneously, a Physical Information Neural Network (PINN) was introduced to fuse CFD data with physical constraints, enabling rapid flow field prediction. Compared to traditional CFD methods, the simulation speed was improved by 8 times (meeting the improvement requirements). More than twice, here The calculation accuracy error is... (satisfies error) ,Right now A parameterized script interface was developed to automate batch processing of operating conditions, improving simulation efficiency. Through the above implementation process, the hydraulic characteristics of the tongue valve in this water conservancy project were successfully simulated, and design optimization and intelligent control were achieved, providing strong technical support for the actual operation of the project.

[0046] like Figure 2 As shown, an apparatus for simulating the hydraulic characteristics of a tongue-shaped valve gate according to a specific embodiment of the present invention is described. The apparatus includes: Modeling unit 201 is used to construct a three-dimensional geometric model of the tongue flap gate containing key parameters and to divide the computational domain. It adopts an appropriate mesh generation technique and completes the mesh independence verification. Solver 202 is used to establish the governing equation system of the tongue valve based on the three-dimensional geometric model of the tongue valve and set the boundary conditions, and to numerically solve the hydraulic properties of the tongue valve. Analysis unit 203 is used to perform post-processing analysis of the valve based on the numerical values ​​of the valve's hydraulic characteristics to obtain the hydraulic characteristics and to perform multi-condition comparative analysis. Simulation unit 2044 is used to verify the simulation model through physical model experiments, optimize the structural parameters based on the verification results, and establish a mapping relationship model to realize intelligent operation control. At the same time, parallel computing and prediction models are used to accelerate simulation calculation.

[0047] Figure 3 The hardware structure diagram of a computing device 30 for simulating the hydraulic characteristics of a tongue valve is shown according to an embodiment of this specification. The computing device 30 may include at least one processor 301, a memory 302 (e.g., non-volatile memory), a main memory 303, and a communication interface 304, and the at least one processor 301, memory 302, main memory 303, and communication interface 304 are connected together via a bus 305. At least one processor 301 executes at least one computer-readable instruction stored or encoded in the memory 302.

[0048] In several embodiments provided in the present application, it should be understood that the disclosed system, device and method can be implemented in other manners. For example, the division of the units is only a logical function division, and there can be another division manner for the actual implementation, for example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the displayed or discussed mutual couplings or direct couplings or communication connections between different units, or the among different units, can be indirect couplings or communication connections through some interfaces, devices or units, and can be in electrical, mechanical or other forms.

[0049] In addition, each function unit in the various embodiments of the present application can be integrated into a processing unit, or each unit can exist alone physically, or two or more units can be integrated into one unit. The integrated unit can be implemented in the form of hardware, or in the form of a software functional unit. The above is only an embodiment of the present application, and does not limit the patent scope of the present application. Any equivalent structure or equivalent process transformation using the content of the present application specification and drawings, or direct or indirect application in other related technical fields, are also included in the patent protection scope of the present application.

[0050] The specific embodiments of the application are described in detail above, but they are only examples. The present application is not limited to the specific embodiments described above. Any equivalent modification or substitution to the application by those skilled in the art is also within the scope of the present application. Therefore, any equivalent transformation and modification, improvement, etc. made without departing from the spirit and principle range of the present application should be included in the scope of the present application.

Claims

1. A simulation method for the hydraulic characteristics of a tongue-shaped valve, characterized in that, The method includes: A three-dimensional geometric model of the tongue flap gate containing key parameters was constructed and the computational domain was divided. An appropriate mesh generation technique was adopted and the mesh independence was verified. Based on the three-dimensional geometric model of the tongue valve gate, the governing equation system of the tongue valve gate is established and boundary conditions are set to numerically solve the hydraulic properties of the tongue valve gate. Post-processing analysis of the tongue valve valve is carried out based on the numerical values ​​of its hydraulic characteristics to obtain the hydraulic properties, and multi-condition comparative analysis is performed. The simulation model is verified through physical model experiments. Based on the verification results, the structural parameters are optimized and a mapping relationship model is established to achieve intelligent operation control. At the same time, parallel computing and prediction models are used to accelerate the simulation calculation.

2. The simulation method for the hydraulic characteristics of the tongue-shaped valve gate according to claim 1, characterized in that, The steps for constructing a three-dimensional geometric model containing key parameters and dividing the computational domain are as follows: A 3D modeling software was used to create a model including the opening of the lingual valve. Door body curvature Doorway width A three-dimensional parametric geometric model, wherein, This indicates the minimum opening of the lingual valve. Indicates the maximum opening of the lingual valve; This represents the minimum radius of curvature of the door body surface. Indicates the maximum curvature of the door body surface; Indicates the minimum width of the door slot. Indicates the maximum width of the door slot; Divide the computational domain and set the length of the upstream inlet section as follows: The gate height is set to a multiple of μ, and the downstream outlet section length is set to μ times the gate height. Boolean operations are used to integrate all components into a full flow channel geometric model to ensure the continuity and sealing of the flow channel. The coefficients for the multiple relationship between the upstream inlet section length and the gate height, as well as the multiple relationship coefficient between the downstream outlet section length and the gate height, are also set.

3. The simulation method for the hydraulic characteristics of the ancient clairmonte as described in claim 1, characterized in that, The steps for employing a suitable mesh generation technique and completing the mesh independence verification are as follows: A tetrahedral-hexahedral hybrid mesh generation technique is used to generate meshes in the door body and door slot areas. Tetrahedral unstructured mesh and refined boundary layer to Upstream and downstream channels adopt The hexahedral structured mesh, This represents the minimum size of the tetrahedral unstructured mesh for the door body and door slot area. Indicates the maximum size; The dimensionless upper limit of the distance for boundary layer encryption control; This represents the minimum size of the hexahedral structured mesh for the upstream and downstream flow channels. Indicates the maximum size; Control the total number of grids in Preliminary calculations were performed on key hydraulic parameters under different grid densities. This indicates the lower limit of the total number of grid cells. Indicates the maximum number of grid cells; Grid independence verification is performed. If the change in key hydraulic parameters is less than the set threshold when the number of grids increases by a certain proportion, then the grid division scheme is considered to meet the requirements.

4. The simulation method for the hydraulic characteristics of the tongue-shaped valve according to claim 1, characterized in that, The steps for establishing the governing equation system and setting boundary conditions are as follows: Establish a continuity equation u=0, momentum equation p u t+u u p u u)+F and RNGk- The governing equations are based on the turbulence model, and the standard wall function is used to handle the near-wall region. The VOF model was used to track the gas-liquid interface. With the surface tension coefficient κ and gravitational acceleration q as the constant, a formula for calculating the vaporization pressure pv based on water temperature T was constructed. T T)) is used to simulate cavitation effects, where к is the surface tension coefficient; ζ and θ are coefficients in the formula for calculating vaporization pressure; Set the import as the speed entry point. , This indicates the lower limit of the import speed. Indicates the upper limit of import speed Grid angular velocity ω , This indicates the lower limit of the angular velocity of the tongue flap rotation around its axis. This indicates the upper limit of angular velocity.

5. The simulation method for the hydraulic characteristics of the tongue valve gate according to claim 1, characterized in that, The numerical solution steps are as follows: The ANSYS Fluent solver was selected, and the SIMPLE algorithm was used for pressure-velocity coupled calculations. The convection term was discretized using a second-order upwind scheme, and the diffusion term was discretized using a central difference scheme. The residual convergence criterion is set as the residual curve of each equation being lower than ν. At the same time, parameters such as flow rate and door force are monitored. When the parameter fluctuation amplitude ν < φ, the solution is judged to be converged. ν is the upper limit of the residual convergence criterion, and φ is the upper limit of the parameter fluctuation amplitude. Steady-state and transient solutions were performed separately, and key data during the solution process were recorded to provide a basis for subsequent analysis.

6. The simulation method for the hydraulic characteristics of the tongue-shaped valve according to claim 1, characterized in that, The analysis, verification, and scheme optimization are as follows: Analysis, verification, and solution optimization: Feature calculation: Flow coefficient ,in, For traffic, For the height of the door, Analyze the flow capacity using gravitational acceleration; cavitation index For local pressure, For flow rate, identify High-risk areas; Hydrodynamic torque Let be the distance from the infinitesimal element on the door surface to the axis of rotation. Optimize the opening and closing mechanism parameters based on the gate's surface area; Multi-condition comparison: set water level , Within the range, The lowest water level, For the highest water level, opening degree Water temperature Within the range, The lowest water temperature, For the orthogonal operating condition with the highest water temperature, temperature correction is considered, and the following is constructed: Establish a working condition database to analyze the influence patterns of parameters; The model validation includes comparing simulation and physical test data, and calculating the flow rate error. Upper limit of flow error, pressure matching degree For the lower limit of pressure mating degree and cavitation position deviation To determine the upper limit of cavitation position deviation, the parameters were adjusted to meet the hydraulic standards. The NSGA-II algorithm was then used to optimize the modeling parameters of the tongue valve. The objective function of the NSGA-II algorithm can be the flow efficiency, cavitation index, and opening and closing torque. After verification, the optimal parameter combination was determined.

7. The simulation method for the hydraulic characteristics of the tongue valve gate according to claim 1, characterized in that, The intelligent hole control and accelerated computing include: Intelligent cave control: based on A set of simulation data was used to establish a mapping model through a BP neural network and embedded into a digital twin platform; multiple operating condition logics were set to realize the opening adjustment error. , This is the upper limit of the opening degree control error; Accelerated computation: CUDA is configured and computational efficiency is improved; PINN is introduced to fuse CFD data with control equations, resulting in faster prediction speeds. times, The lower limit of the speed improvement factor, accuracy error To limit the accuracy error, parameterized scripts were developed to support batch processing.

8. A simulation device for the hydraulic characteristics of a tongue-shaped valve gate, which is applied to the simulation method for the hydraulic characteristics of a tongue-shaped valve gate as described in any one of claims 1 to 7, characterized in that, The simulation device for the hydraulic characteristics of the tongue valve includes: The modeling unit is used to construct a three-dimensional geometric model of the tongue flap gate containing key parameters and divide the computational domain, and to use appropriate meshing techniques and complete mesh independence verification. The solution unit is used to establish the governing equation system of the tongue valve based on the three-dimensional geometric model of the tongue valve and set the boundary conditions to numerically solve the hydraulic properties of the tongue valve. The analysis unit is used to perform post-processing analysis of the valve based on the numerical values ​​of the valve's hydraulic characteristics to obtain the hydraulic characteristics and to conduct multi-condition comparative analysis. The simulation unit is used to verify the simulation model through physical model experiments, optimize the structural parameters based on the verification results, and establish a mapping relationship model to achieve intelligent operation control. At the same time, parallel computing and prediction models are used to accelerate simulation calculations.

9. An electronic device, characterized in that, The method includes a processor and a memory coupled to the processor, the memory storing program instructions executable by the processor; when the processor executes the program instructions stored in the memory, it implements the simulation method for the hydraulic characteristics of the tongue valve as described in any one of claims 1 to 7.

10. A storage medium, characterized in that, The storage medium stores program instructions, which, when executed by a processor, enable the simulation method for the hydraulic characteristics of the tongue valve as described in any one of claims 1 to 7.

Images