A wave simulation method, apparatus, equipment, and medium based on a triangular numerical wave tank.

Abstract

This invention relates to the field of wave simulation technology, and discloses a wave simulation method, apparatus, equipment, and medium based on a triangular numerical wave flume. The method includes constructing a right-angled triangular computational domain, meshing the right-angled triangular computational domain, setting the top boundary of the right-angled triangular computational domain as a pressure outlet, and setting the horizontal right-angled side, vertical right-angled side, and inclined side as velocity inlets; loading the theoretical velocity vector of wave water particles onto the velocity inlets corresponding to the horizontal right-angled side, vertical right-angled side, and inclined side; constructing a wave simulation model and superimposing a momentum source term within the wave simulation model; adjusting the difference between the theoretical velocity solution of the wave and the current numerical solution of the flow field through a weighting function of the momentum source term; iteratively solving the wave simulation model, and outputting and updating the wave surface elevation and flow field velocity distribution within the right-angled triangular computational domain. This invention can reduce the problem of severe reflection interference at the lateral boundary and improve computational efficiency.

Description

Technical Field

[0001] This invention relates to the field of wave simulation technology, and more specifically, to a wave simulation method, apparatus, equipment, and medium based on a triangular numerical wave tank. Background Technology

[0002] Numerical wave flumes are indispensable research tools in marine engineering, coastal structure design, and nearshore renewable energy development. Traditional numerical wave flumes often use rectangular or two-dimensional cross-sectional models, which make it difficult to completely eliminate boundary reflection problems. This typically requires complex wave-damping layers, which not only expands the computational domain but also significantly increases computational costs. Furthermore, when simulating multi-directional irregular waves or oblique waves, the sidewalls of rectangular flumes generate severe reflection interference, which seriously affects the realism of the flow field and the accuracy of the calculations.

[0003] In existing technologies, numerical wave pools based on commercial software (such as STAR-CCM+ and Fluent) are mostly constructed by combining velocity inlets or wave-generating plate boundaries with damping wave-absorbing zones. However, these methods still face the problem of incomplete wave absorption when dealing with complex wave conditions such as focused waves and oblique waves, or when analyzing the three-dimensional coupled response of structures in waves. Therefore, existing technologies suffer from severe lateral boundary reflection interference and low computational efficiency. Summary of the Invention

[0004] To overcome the shortcomings of existing wave simulation techniques, such as severe lateral boundary reflection interference and low computational efficiency, the present invention proposes the following technical solution:

[0005] In the first aspect, this invention proposes a wave simulation method based on a triangular numerical wave tank, comprising:

[0006] A right-angled triangle computational domain is constructed; the right-angled triangle computational domain includes a top boundary and three lateral boundaries consisting of a horizontal right-angled side, a vertical right-angled side, and an inclined side;

[0007] The right-angled triangular computational domain is meshed, and the top boundary of the right-angled triangular computational domain is set as the pressure outlet, while the horizontal right-angled side, the vertical right-angled side, and the inclined side are all set as velocity inlets.

[0008] Calculate the theoretical velocity vector of the wave-water particle, and load the theoretical velocity vector onto the velocity inlets corresponding to the horizontal right-angled side, the vertical right-angled side, and the inclined side;

[0009] A wave simulation model is constructed within the right-angled triangular computational domain, and a momentum source term is superimposed on the wave simulation model; the momentum source term is adjusted by a weighting function to control the difference between the theoretical velocity solution of the wave and the current numerical solution of the flow field.

[0010] The wave simulation model is iteratively solved to output and update the wave surface elevation and flow velocity distribution within the right-angled triangular computational domain.

[0011] As a preferred technical solution, the construction of the right-angled triangle computational domain includes:

[0012] Construct an isosceles right triangle as the basic computational domain;

[0013] The normal direction of the inclined side of the isosceles right triangle is configured as the main propagation direction of the wave;

[0014] Configure the horizontal right-angled side of the isosceles right triangle as the bottom boundary of the basic computational domain;

[0015] Configure the perpendicular right-angled side of the isosceles right triangle as the sidewall boundary of the basic computational domain;

[0016] The inclined side of the isosceles right triangle is configured as a wave-shaped incident boundary.

[0017] As a preferred technical solution, a multiphase flow model and a viscous incompressible fluid model are selected and activated;

[0018] The Reynolds-averaged Navier-Stokes equations are set as the governing equations;

[0019] The governing equations are closed using the RNGk-ε model;

[0020] The wave simulation model is obtained by coupling the multiphase flow model, the viscous incompressible fluid model, and the closed-loop governing equations.

[0021] As a preferred technical solution, the calculation of the theoretical velocity vector of wave-water particles includes:

[0022] Based on the Stokes fifth-order wave theory, calculate the theoretical horizontal velocity component and theoretical vertical velocity component of the wave water particle at different times.

[0023] The theoretical velocity vector is constructed using the theoretical horizontal velocity component and the theoretical vertical velocity component.

[0024] As a preferred technical solution, a momentum source term is superimposed within the wave simulation model, including:

[0025] Calculate the velocity difference between the wave theory velocity and the current flow field numerical velocity, and set the velocity difference as a momentum source term within the right-angled triangular computational domain. The expression is as follows:

[0026]

[0027] In the formula, Represents the momentum source term. Indicates the forcing coefficient. Indicates fluid density, This represents the current numerical solution of the flow field. This represents the velocity solution in wave theory.

[0028] As a preferred technical solution, the method further includes obtaining the forcing coefficient. ,include:

[0029] Within the wave-generating and wave-dissipating regions of the right-angled triangular computational domain, a spatial weight distribution based on trigonometric functions is established.

[0030] The forcing coefficient values ​​at different spatial locations are calculated based on the spatial weight distribution, and the expression is as follows:

[0031]

[0032] In the formula, This represents the magnitude of the forcing coefficient. This represents the dimensionless spatial coordinates of the wave-generating region or the wave-dissipating region.

[0033] As a preferred technical solution, the step of loading the theoretical velocity vector onto the velocity inlets corresponding to the horizontal right-angled side, the vertical right-angled side, and the inclined side includes:

[0034] When an obliquely incident wave is generated, the phase delay parameter is calculated based on the spatial position of each grid node on the oblique edge. The theoretical velocity vector is then corrected using the phase delay parameter, and the corrected velocity vector is loaded into the velocity inlet.

[0035] When generating multi-directional irregular waves, based on the principle of linear superposition, the theoretical velocity vectors of multiple wave components with different frequencies and different incident angles are calculated. The theoretical velocity vectors of the multiple wave components are synthesized, and the synthesized velocity vector is loaded into the velocity inlet.

[0036] As a preferred technical solution, after updating the wavefront elevation and flow field velocity distribution within the right-angled triangular computational domain, the method further includes:

[0037] A structural model and an overlapping mesh enclosing the structural model are established within the right-angled triangular computational domain;

[0038] Based on the iterative solution results of the wave simulation control equations, the flow field pressure distribution data within the right-angled triangular computational domain are obtained;

[0039] Perform data interpolation between the overlapping mesh and the base mesh of the right-angled triangular computational domain to transfer the flow field pressure distribution data to the structural model;

[0040] Integrate the flow field pressure distribution data on the surface of the structure model to output the wave load on the structure model.

[0041] In a second aspect, the present invention also proposes an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to perform the operations performed by the wave simulation method based on a triangular numerical wave tank as described in any of the embodiments of the first aspect.

[0042] Thirdly, the present invention also proposes a computer-readable storage medium storing a program that is executed by a processor as described in any of the embodiments of the first aspect, using the wave simulation method based on a triangular numerical wave tank.

[0043] The beneficial effects of this invention include at least the following: By constructing a right-angled triangular computational domain consisting of a top boundary and lateral boundaries composed of horizontal, vertical, and inclined sides, the top boundary is set as the pressure outlet, and the three right-angled or inclined sides are set as velocity inlets, directly loading the theoretical velocity vector of the wave water particles. This specific geometric configuration and boundary condition combination provides precise physical quantity input and forced constraint on the wave flow field from the boundary source. Furthermore, momentum source terms are superimposed within the wave simulation model, and a weighting function is used to dynamically adjust the difference between the theoretical velocity solution of the wave and the current numerical solution of the flow field. By closely combining the theoretical velocity constraint at the boundary with the error correction of the momentum source terms within the computational domain, this invention enables the internal flow field numerical solution to be actively and accurately guided to approach the theoretical solution state during the iterative solution process, effectively suppressing the generation of non-physical reflection waves during wave propagation, thereby reducing the problem of severe reflection interference at the lateral boundary and improving computational efficiency. Attached Figure Description

[0044] Figure 1 This is a flowchart illustrating the wave simulation method based on a triangular numerical wave tank provided in an embodiment of the present invention.

[0045] Figure 2 A schematic diagram of the three-dimensional structure of the triangular numerical wave tank provided in an embodiment of the present invention.

[0046] Figure 3 Provided for embodiments of the present invention Figure 2 Cross-sectional view at point AA.

[0047] Figure 4 This is a schematic diagram of the three-dimensional mesh division of the triangular computational domain provided in an embodiment of the present invention.

[0048] Figure 5 Provided for embodiments of the present invention Figure 4 A magnified view of a section at point B.

[0049] Figure 6 Three-dimensional morphological cloud maps of wave evolution at different times are provided for embodiments of the present invention.

[0050] Figure 7 This is a long-time history comparison and verification diagram of the numerical simulation wavefront elevation with theoretical and experimental values ​​provided in the embodiments of the present invention.

[0051] Figure 8 This is a magnified comparison verification diagram of the local wavefront elevation provided in an embodiment of the present invention.

[0052] Figure 9 This is a time history diagram of wave surface elevation under different wave incident angles provided in an embodiment of the present invention, wherein, Figure 9 (a) Provided by an embodiment of the present invention, where the wave incident angle is... A time history diagram of wavefront elevation. Figure 9 (b) Provided by the embodiment of the present invention, at a wave incident angle of... A time history diagram of wavefront elevation. Figure 9 (c) Provided by the embodiment of the present invention at a wave incident angle of... A time history diagram of wavefront elevation. Figure 9 (d) is provided in the embodiment of the present invention when the wave incident angle is... A time history diagram of wavefront elevation.

[0053] Figure 10 The three-dimensional wavefront morphology cloud maps under different wave incident angles provided in the embodiments of the present invention, wherein, Figure 10 (a) Provided by an embodiment of the present invention, where the wave incident angle is... Three-dimensional wavefront morphology cloud map, Figure 10 (b) Provided by the embodiment of the present invention at a wave incident angle of... Three-dimensional wavefront morphology cloud map, Figure 10 (c) Provided by the embodiment of the present invention at a wave incident angle of... Three-dimensional wavefront morphology cloud map, Figure 10 (d) is provided in the embodiment of the present invention when the wave incident angle is... Three-dimensional wavefront morphology cloud map.

[0054] Figure 11 This is a comparison and verification diagram of the oblique incident wavefront elevation and the analytical solution provided in an embodiment of the present invention.

[0055] Figure 12 The three-dimensional wavefront morphology cloud map of multi-directional irregular waves provided in the embodiments of the present invention.

[0056] Figure 13 This is a time history comparison verification diagram of the horizontal wave force on the structure provided in the embodiment of the present invention.

[0057] Figure 14 This is a time history comparison and verification diagram of the pitching moment acting on the structure provided in the embodiment of the present invention.

[0058] Figure 15 The time history curves of horizontal wave forces on a structure under different wave incident directions are provided for embodiments of the present invention.

[0059] Figure 16 The time history curves of pitching moment under different wave incident directions are provided for embodiments of the present invention.

[0060] Figure 17 The above diagrams show the time history of the float's rise and fall under different conditions, as provided in the embodiments of the present invention.

[0061] Figure 18 The float rise and fall time history curves under different incident directions are provided for embodiments of the present invention.

[0062] Figure 19 This is a schematic diagram of the structure of the electronic device provided in an embodiment of the present invention. Detailed Implementation

[0063] The embodiments of the present invention will be described below with reference to the accompanying drawings and preferred technical solutions. Those skilled in the art can easily understand other advantages and effects of the present invention from the content disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention. It should be understood that the preferred technical solutions are only for illustrating the present invention and not for limiting the scope of protection of the present invention.

[0064] It should be noted that the illustrations provided in the following embodiments are only schematic representations of the basic concept of the present invention. Therefore, the drawings only show the components related to the present invention and are not drawn according to the actual number, shape and size of the components in the actual implementation. In the actual implementation, the form, quantity and proportion of each component can be arbitrarily changed, and the layout of the components may also be more complex.

[0065] In the following description, numerous details are explored to provide a more thorough explanation of embodiments of the invention. However, it will be apparent to those skilled in the art that embodiments of the invention may be practiced without these specific details. In other embodiments, well-known structures and devices are shown in block diagram form rather than in detail to avoid obscuring embodiments of the invention.

[0066] Example 1

[0067] This embodiment proposes a wave simulation method based on a triangular numerical wave tank, such as... Figure 1 As shown, Figure 1 This is a flowchart illustrating a wave simulation method based on a triangular numerical wave tank provided in this embodiment. The method includes the following steps:

[0068] S1: Construct a right-angled triangle computational domain; the right-angled triangle computational domain includes a top boundary and three lateral boundaries consisting of a horizontal right-angled side, a vertical right-angled side, and an inclined side;

[0069] S2: Mesh the right-angled triangle computational domain, set the top boundary of the right-angled triangle computational domain as the pressure outlet, and set the horizontal right-angled side, the vertical right-angled side, and the inclined side as the velocity inlet;

[0070] S3: Calculate the theoretical velocity vector of the wave water particles, and load the theoretical velocity vector onto the velocity inlets corresponding to the horizontal right-angled side, the vertical right-angled side, and the inclined side;

[0071] S4: Construct a wave simulation model within the right-angled triangular computational domain, and superimpose a momentum source term within the wave simulation model; the momentum source term adjusts the difference between the wave theoretical velocity solution and the current flow field numerical solution through a weighting function;

[0072] S5: Iteratively solve the wave simulation model, output and update the wave surface elevation and flow velocity distribution within the right-angled triangle computational domain.

[0073] In the specific implementation process, a right-angled triangular computational domain was first constructed, consisting only of the top boundary and lateral boundaries formed by horizontal, vertical, and inclined sides. This significantly reduced the required computational space in terms of physical topology. Subsequently, based on the spatial meshing of this computational domain, its boundary conditions were specifically configured. The top was set as a pressure outlet to simulate an open atmospheric connection environment, and the horizontal, vertical, and inclined sides were unconventionally set as velocity inlets. By calculating the theoretical velocity vector of wave-water particles, it was directly forced onto the velocity inlets of the aforementioned three sides, thereby forming a precise, enclosed physical quantity input constraint from the boundary source. To further precisely control the evolution of the flow field within the computational domain, a momentum source term was superimposed on the basic wave simulation model. This momentum source term actually acts as a real-time feedback regulator for the flow field error. It uses a preset weighting function to dynamically calculate the difference between the expected theoretical velocity of the wave and the numerical velocity actually calculated in the current flow field in space, and directly applies momentum compensation to the flow field based on this difference. Finally, by continuously iterating the wave simulation model with the above dynamic adjustment mechanism over time, the wave surface elevation fluctuation and flow field velocity distribution inside the right-angled triangular computational domain can be output efficiently and stably and continuously updated, thus completing a high-precision wave numerical simulation.

[0074] Understandably, by constructing a right-angled triangular computational domain consisting of a top boundary and lateral boundaries formed by horizontal, vertical, and inclined sides, setting the top boundary as the pressure outlet, and setting the aforementioned three right-angled or inclined sides as velocity inlets and directly loading the theoretical velocity vectors of wave water particles, this specific geometric configuration and boundary condition combination provides precise physical quantity input and forced constraints to the wave flow field from the boundary source. Furthermore, a momentum source term is superimposed within the wave simulation model, and a weighting function is used to dynamically adjust the difference between the theoretical wave velocity solution and the current numerical solution of the flow field. This invention, by closely combining the theoretical velocity constraints at the boundary with the momentum source term error correction within the computational domain, allows the internal flow field numerical solution to be actively and precisely guided to approximate the theoretical solution state during iterative solving, effectively suppressing the generation of non-physical reflection waves during wave propagation, thereby reducing the problem of severe reflection interference at the lateral boundary and improving computational efficiency.

[0075] Example 2

[0076] This embodiment is an improvement on the wave simulation method based on a triangular numerical wave tank proposed in Embodiment 1.

[0077] In this embodiment, constructing the right-angled triangle computational domain includes:

[0078] An isosceles right triangle is constructed as the basic computational domain on the horizontal projection plane; the basic computational domain is a three-dimensional spatial domain formed by extending vertically, and the upper surface of the three-dimensional spatial domain constitutes the top boundary;

[0079] The normal direction of the inclined side of the isosceles right triangle is configured as the main propagation direction of the wave;

[0080] Configure the horizontal right-angled side of the isosceles right triangle as the bottom boundary of the basic computational domain in the plan view;

[0081] Configure the perpendicular right-angled side of the isosceles right triangle as the sidewall boundary of the basic computational domain in the plan view;

[0082] The inclined side of the isosceles right triangle is configured as a wave-shaped incident boundary.

[0083] It should be noted that by constructing an isosceles right triangle as the basic computational domain, and configuring the normal direction of the inclined side of the isosceles right triangle as the main propagation direction of the wave, configuring the horizontal right-angled side as the bottom boundary, the vertical right-angled side as the sidewall boundary, and the inclined side as the wave incident boundary, the orthogonal symmetry characteristics of isosceles right-angled geometry can be utilized to make the distribution of the flow field components in the horizontal and vertical directions more balanced and symmetrical when the wave is incident along the normal direction of the inclined side.

[0084] By using this specific triangular configuration and boundary condition allocation, non-physical reflection interference generated by oblique waves at conventional boundaries can be effectively avoided. While ensuring the stable evolution of waves, the length of the computational domain is shortened to the maximum extent, thereby significantly improving the computational efficiency of wave simulation.

[0085] In this embodiment, a wave simulation model is constructed within the right-angled triangle computational domain, including:

[0086] Select and activate the multiphase flow model and the viscous incompressible fluid model;

[0087] The Reynolds-averaged Navier-Stokes equations are set as the governing equations;

[0088] The governing equations are closed using the RNGk-ε model;

[0089] The wave simulation model is obtained by coupling the multiphase flow model, the viscous incompressible fluid model, and the closed-loop governing equations.

[0090] It should be noted that the deep coupling of multiple fluid dynamics models achieves high-fidelity reproduction of the real wave dynamics environment. Specifically, the multiphase flow model accurately tracks and reconstructs the complex gas-liquid interface (i.e., the free surface), the viscous incompressible fluid model conforms to the fundamental physical properties of real water, and the Reynolds-averaged governing equations with RNGk-ε turbulence closure effectively handle the complex turbulent characteristics such as high strain rates and streamline curvature that accompany wave evolution. This invention significantly improves the accuracy of calculating flow field density and pressure gradient at the gas-liquid interface, effectively suppressing non-physical numerical dissipation and wave surface ambiguity during long-distance wave propagation. This provides an extremely stable and highly realistic basic wave simulation model for subsequent momentum source term superposition wave generation and the output of high-precision wave surface elevation and flow field velocity distribution.

[0091] In this embodiment, calculating the theoretical velocity vector of wave-water particles includes:

[0092] Based on the Stokes fifth-order wave theory, calculate the theoretical horizontal velocity component and theoretical vertical velocity component of the wave water particle at different times.

[0093] The theoretical velocity vector is constructed using the theoretical horizontal velocity component and the theoretical vertical velocity component.

[0094] It should be noted that, based on the aforementioned Stokes fifth-order wave theory, the theoretical horizontal and vertical velocity components of wave particles at different times are calculated, and a theoretical velocity vector is constructed using these components. The principle is that, compared to linear micro-amplitude wave theory, higher-order Stokes wave theory introduces higher-order nonlinear terms for expansion, which can more accurately describe the nonlinear hydrodynamic characteristics under finite water depth or steep wave conditions. The velocity vector obtained using Stokes fifth-order wave theory has higher physical accuracy and can accurately reflect the asymmetric characteristics of real waves, such as sharper crests and flatter troughs, thus ensuring the authenticity of the generated wave morphology and flow field from the source of the boundary input.

[0095] In this embodiment, a momentum source term is superimposed within the wave simulation model, including:

[0096] Calculate the velocity difference between the wave theory velocity and the current flow field numerical velocity, and set the velocity difference as a momentum source term within the right-angled triangular computational domain. The expression is as follows:

[0097]

[0098] In the formula, Represents the momentum source term. Indicates the forcing coefficient. Indicates fluid density, This represents the current numerical solution of the flow field. This represents the velocity solution in wave theory.

[0099] It should be noted that a penalty feedback mechanism based on velocity difference is introduced into the momentum equation to force the internal flow field to evolve towards the theoretical expectation. By directly adjusting the error through this momentum source term, real-time flexible coupling between the numerical solution and the theoretical solution is achieved, effectively suppressing wave reflection accumulation in the computational domain and accelerating the iterative convergence process of the wave simulation control equation.

[0100] In this embodiment, the method further includes obtaining the forcing coefficient. ,include:

[0101] Within the wave-generating and wave-dissipating regions of the right-angled triangular computational domain, a spatial weight distribution based on trigonometric functions is established.

[0102] The forcing coefficient values ​​at different spatial locations are calculated based on the spatial weight distribution, and the expression is as follows:

[0103]

[0104] In the formula, This represents the magnitude of the forcing coefficient. This represents the dimensionless spatial coordinates of the wave-generating region or the wave-dissipating region.

[0105] It should be noted that this embodiment utilizes the smooth transition characteristics of the cosine square function to construct a weight field that gradually varies with spatial position. This smooth spatial weight distribution eliminates the abrupt changes in momentum source terms during boundary action, resulting in a gradual transition of the forced adjustment force of the source terms. This completely avoids the non-physical pressure fluctuations in the flow field caused by forced wave generation or wave dissipation, greatly enhancing the robustness of the numerical solution.

[0106] In this embodiment, loading the theoretical velocity vector onto the velocity inlets corresponding to the horizontal right-angled side, the vertical right-angled side, and the inclined side includes:

[0107] When an obliquely incident wave is generated, the phase delay parameter is calculated based on the spatial position of each grid node on the oblique edge. The theoretical velocity vector is then corrected using the phase delay parameter, and the corrected velocity vector is loaded into the velocity inlet.

[0108] When generating multi-directional irregular waves, based on the principle of linear superposition, the theoretical velocity vectors of multiple wave components with different frequencies and different incident angles are calculated. The theoretical velocity vectors of the multiple wave components are synthesized, and the synthesized velocity vector is loaded into the velocity inlet.

[0109] It should be noted that when generating obliquely incident waves, the phase delay parameter is calculated based on the spatial position of each grid node on the oblique edge to correct the theoretical velocity vector. When generating multi-directional irregular waves, multiple theoretical velocity vectors of wave components with different frequencies and incident angles are synthesized and loaded based on the principle of linear superposition. The principle lies in the spatial position dependence of wave phase and the additivity of multi-frequency wave components under linear conditions. This embodiment does not require reconstructing complex physical boundaries or modifying the grid topology. By only correcting and synthesizing the boundary data at the velocity inlet, it can flexibly and efficiently generate complex oblique waves and multi-directional irregular waves within the right-angled triangular computational domain, breaking through the bottleneck of traditional flumes in handling multi-directional waves when simulating complex sea conditions.

[0110] In this embodiment, after updating the wavefront elevation and flow field velocity distribution within the right-angled triangle computational domain, the method further includes:

[0111] A structural model and an overlapping mesh enclosing the structural model are established within the right-angled triangular computational domain;

[0112] Based on the iterative solution results of the wave simulation control equations, the flow field pressure distribution data within the right-angled triangular computational domain are obtained;

[0113] Perform data interpolation between the overlapping mesh and the base mesh of the right-angled triangular computational domain to transfer the flow field pressure distribution data to the structural model;

[0114] Integrate the flow field pressure distribution data on the surface of the structure model to output the wave load on the structure model.

[0115] It should be noted that a structural model and an overlapping mesh enclosing the structural model are established within a right-angled triangular computational domain. Flow pressure distribution data is acquired and interpolated before being transferred to the structural model. Finally, the flow pressure distribution data on the surface of the structural model is integrated to output the wave load. The principle lies in using mesh decoupling technology, utilizing contributing and receiving cells in the overlapping region to implicitly transfer physical quantities (such as pressure) between the background wave field and the near-field of the structure. The overlapping mesh mechanism avoids the large mesh deformation or negative volume problems caused by the presence of the structure in traditional methods, ensuring mesh quality and flow field calculation accuracy. Combined with accurate pressure interpolation and surface integration, the wave load on the structure can be output with high fidelity, providing extremely reliable data support for structural hydrodynamic response assessment.

[0116] Example 3

[0117] Based on the wave simulation method based on the triangular numerical wave flume proposed in Example 2, this embodiment conducts a comprehensive simulation verification and result analysis on the effectiveness, accuracy and applicability of the simulation method to marine engineering.

[0118] Figure 2 This is a schematic diagram of the three-dimensional structure of the triangular numerical wave tank provided in an embodiment of the present invention. Figure 3 Provided for embodiments of the present invention Figure 2 Cross-sectional view at point AA. It should be noted that... Figure 2 The invention visually demonstrates the spatial structure of the isosceles right-angled triangle computational domain. By using the inclined side as the wave incident boundary, it cleverly utilizes geometric orthogonal symmetry to effectively shorten the required large-scale damping and wave-dissipating zone, thereby significantly reducing the number of grids.

[0119] like Figure 4 and Figure 5 As shown, Figure 4 This is a schematic diagram of the three-dimensional meshing of the triangular computational domain provided in an embodiment of the present invention. Figure 5 Provided for embodiments of the present invention Figure 4 A magnified view of a portion at point B. It should be noted that by refining the local mesh near the free surface of the wave, it is possible to accurately capture minute changes in wave morphology and the gas-liquid interface while maintaining overall computational efficiency.

[0120] like Figure 6 As shown, Figure 6 This is a three-dimensional morphological cloud map of wave evolution at different times, provided as an embodiment of the present invention. It should be noted that... Figure 6 The text illustrates the wave's behavior over time within the computational domain. :2 to 18; This represents the current actual evolution time. (Representing a complete wave cycle.) The stable propagation process shows that no non-physical wave reflection disturbances occurred at the bottom and sidewall boundaries, directly proving the effectiveness of the combination of this configuration and the momentum source term.

[0121] like Figure 7 and Figure 8 As shown, Figure 7 This is a long-time history comparison and verification diagram of the numerical simulation wavefront elevation with theoretical and experimental values ​​provided in an embodiment of the present invention. Figure 8This is a magnified comparison verification diagram of the local wavefront elevation provided in an embodiment of the present invention. It should be noted that the current simulated numerical wavefront results (Present) maintain a high degree of agreement with the theoretical solution and experimental data of the Stokes 5th Order wave over a long period of evolution, and accurately characterize the nonlinear characteristics of the wave crests becoming sharper and the troughs becoming flatter, verifying the high computational accuracy of the method.

[0122] like Figure 9 and Figure 10 As shown, Figure 9 The wave surface elevation time history diagrams provided in this embodiment of the invention under different wave incident angles specifically include: Figure 9 (a) Provided by an embodiment of the present invention, where the wave incident angle is... A time history diagram of wavefront elevation. Figure 9 (b) Provided by the embodiment of the present invention at a wave incident angle of... A time history diagram of wavefront elevation. Figure 9 (c) Provided by the embodiment of the present invention at a wave incident angle of... A time history diagram of wavefront elevation. Figure 9 (d) is provided in the embodiment of the present invention when the wave incident angle is... A time history diagram of wavefront elevation; Figure 10 The three-dimensional wavefront morphology cloud maps under different wave incident angles provided in the embodiments of the present invention specifically include: Figure 10 (a) Provided by an embodiment of the present invention, where the wave incident angle is... Three-dimensional wavefront morphology cloud map, Figure 10 (b) Provided by the embodiment of the present invention at a wave incident angle of... Three-dimensional wavefront morphology cloud map, Figure 10 (c) Provided by the embodiment of the present invention at a wave incident angle of... Three-dimensional wavefront morphology cloud map, Figure 10 (d) is provided in the embodiment of the present invention when the wave incident angle is... Three-dimensional wavefront morphology cloud map.

[0123] like Figure 11 As shown, Figure 11 This is a comparison and verification diagram of the wavefront elevation of the obliquely incident wave and the analytical solution provided in an embodiment of the present invention. It should be noted that the peaks and troughs of the wave extracted by the numerical simulation (Present) and the analytical solution (Analytic) perfectly coincide, further proving the correctness of the mechanism of correcting the theoretical velocity vector based on the phase delay parameter.

[0124] like Figure 12 As shown, Figure 12This is a three-dimensional wavefront morphology cloud map of a multi-directional irregular wave provided in an embodiment of the present invention. It should be noted that the wavefront exhibits a typical three-dimensional short-peak wave staggered distribution characteristic, indicating that the present invention, based on the principle of linear superposition, has successfully overcome the bottleneck of sidewall reflection in traditional rectangular flume simulation of multi-directional waves, and can reproduce complex real sea conditions with high fidelity.

[0125] like Figure 13 and Figure 14 As shown, Figure 13 The horizontal wave force (WW) on the structure provided in the embodiments of the present invention Time-history comparison verification chart, Figure 14 The pitching moment (PMT) of the structure provided in the embodiments of the present invention Time history comparison verification diagram. It should be noted that the numerical solution obtained by this invention based on overlapping grids and pressure integration is compared with the results of existing technologies (SumandYue, Kumar, etc.) and potential flow software (WAMIT). The two are highly consistent, which fully proves the extremely high fidelity of this invention in calculating wave loads on structures.

[0126] like Figure 15 and Figure 16 As shown, Figure 15 Different wave incident directions provided for embodiments of the present invention ( Time history curve of horizontal wave force on the substructure. Figure 16 The diagram shows the time history curves of the pitching moment under different wave incident directions, as provided in this embodiment of the invention. It should be noted that this verifies that the water tank can stably and accurately evaluate the force evolution of a structure under wave load impacts from different directions.

[0127] like Figure 17 As shown, Figure 17 This is a comparison and verification diagram of the float's rise and fall time history provided in an embodiment of the present invention. It should be noted that when the dimensionless response (Z / A) of the float's rise and fall predicted by the method of the present invention is compared with the data in existing literature (Kim and Ferrant), the period and amplitude of the two are highly consistent, which confirms that the present method has extremely high computational accuracy in predicting the motion response of ocean floating bodies.

[0128] Figure 18 This diagram illustrates the heave-sinking time history curves of the float under different incident directions, as provided in this embodiment of the invention. It should be noted that this diagram further evaluates the heave-sinking response of the float under different wave incident directions (e.g., 90°, 75°, 45°). The results show that the evolution curves are smooth and clearly defined, ultimately confirming that the right-angled triangular numerical wave flume of this invention has extremely high engineering practical value in predicting the hydrodynamic and motion response of floating bodies under complex multi-directional sea conditions.

[0129] Example 4

[0130] Figure 19 This is a schematic diagram of the structure of the electronic device 100 provided in this embodiment. The electronic device 100 includes: a memory 101, a processor 102, and a computer program stored in the memory 101 and executable on the processor 102.

[0131] When the processor 102 executes the program, it implements the wave simulation method based on the triangular numerical wave tank provided in the above embodiments.

[0132] Furthermore, the electronic device 100 also includes a communication interface 103 for communication between the memory 101 and the processor 102.

[0133] The memory 101 may include high-speed RAM (Random Access Memory) memory, and may also include non-volatile memory, such as at least one disk storage.

[0134] If the memory 101, processor 102, and communication interface 103 are implemented independently, then the communication interface 103, memory 101, and processor 102 can be interconnected via a bus to complete communication between them. The bus can be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus, or an EISA (Extended Industry Standard Architecture) bus, etc. The bus can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 19 The bus is represented by a single thick line, but this does not mean that there is only one bus or one type of bus.

[0135] Optionally, in a specific implementation, if the memory 101, processor 102, and communication interface 103 are integrated on a single chip, then the memory 101, processor 102, and communication interface 103 can communicate with each other through an internal interface.

[0136] Processor 102 may be a CPU (Central Processing Unit), an ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement embodiments of the present invention.

[0137] This invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the wave simulation method based on a triangular numerical wave tank as described above.

[0138] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of the present invention. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0139] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this invention, "N" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0140] Any process or method description in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or more N executable instructions for implementing custom logic functions or processes, and the scope of preferred embodiments of the invention includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the functions involved, as should be understood by those skilled in the art to which embodiments of the invention pertain.

[0141] It should be understood that various parts of the present invention can be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, the N steps or methods can be implemented in software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware as in another embodiment, it can be implemented using any of the following techniques known in the art, or a combination thereof: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (FPGAs), field-programmable gate arrays (FPGAs), etc.

[0142] Those skilled in the art will understand that all or part of the steps of the methods described in the above embodiments can be implemented by a program instructing related hardware, and the program can be stored in a computer-readable storage medium. When executed, the program includes one or a combination of the steps of the method embodiments.

[0143] Obviously, the above embodiments of the present invention are merely examples for clearly illustrating the present invention, and are not intended to limit the implementation of the present invention. Those skilled in the art can make other variations or modifications based on the above description. It is neither necessary nor possible to exhaustively describe all embodiments here. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention should be included within the scope of protection of the claims of the present invention.

Claims

1. A wave simulation method based on a triangular numerical wave flume, characterized by, include: Constructing a right-angled triangular computational domain includes: An isosceles right triangle is constructed as the basic computational domain on the horizontal projection plane; the basic computational domain is a three-dimensional spatial domain formed by extending vertically, and the upper surface of the three-dimensional spatial domain constitutes the top boundary; The normal direction of the inclined side of the isosceles right triangle is configured as the main propagation direction of the wave; Configure the horizontal right-angled side of the isosceles right triangle as the bottom boundary of the basic computational domain in the plan view; Configure the perpendicular right-angled side of the isosceles right triangle as the sidewall boundary of the basic computational domain in the plan view; Configure the inclined side of the isosceles right triangle as a wavy incident boundary; The right-angled triangular computational domain is meshed, and the top boundary of the right-angled triangular computational domain is set as the pressure outlet, while the horizontal right-angled side, the vertical right-angled side, and the inclined side are all set as velocity inlets. Calculate the theoretical velocity vector of the wave-water particle, and load the theoretical velocity vector onto the velocity inlets corresponding to the horizontal right-angled side, the vertical right-angled side, and the inclined side; A wave simulation model is constructed within the right-angled triangular computational domain, and a momentum source term is superimposed on the wave simulation model. The momentum source term adjusts the difference between the theoretical wave velocity solution and the current flow field numerical solution through a weighting function. The construction of the wave simulation model within the right-angled triangular computational domain includes: Select and activate the multiphase flow model and the viscous incompressible fluid model; The Reynolds-averaged Navier-Stokes equations are set as the governing equations; The governing equations are closed using the RNGk-ε model; The wave simulation model is obtained by coupling the multiphase flow model, the viscous incompressible fluid model, and the closed-loop governing equations. The wave simulation model is iteratively solved to output and update the wave surface elevation and flow velocity distribution within the right-angled triangular computational domain.

2. The wave simulation method based on a triangular numerical wave flume according to claim 1, characterized in that, Calculating the theoretical velocity vector of wave-water particles includes: Based on Stokes' fifth-order wave theory, the theoretical horizontal velocity component and theoretical vertical velocity component of wave water particles at different times are calculated. A theoretical velocity vector is constructed using the theoretical horizontal velocity component and the theoretical vertical velocity component.

3. The triangular numerical wave tank based wave simulation method according to claim 1, wherein, The momentum source term is superimposed within the wave simulation model, including: Calculate the velocity difference between the wave theory velocity and the current flow field numerical velocity, and set the velocity difference as a momentum source term within the right-angled triangular computational domain. The expression is as follows: wherein represents the momentum source term, represents the forcing coefficient, represents the fluid density, represents the current flow field numerical solution, represents the wave theoretical velocity solution.

4. The wave simulation method based on a triangular numerical wave tank according to claim 3, characterized in that, The method also includes obtaining the forcing coefficient. ,include: Within the wave-generating and wave-dissipating regions of the right-angled triangular computational domain, a spatial weight distribution based on trigonometric functions is established. The forcing coefficient values ​​at different spatial locations are calculated based on the spatial weight distribution, and the expression is as follows: In the formula, This represents the magnitude of the forcing coefficient. This represents the dimensionless spatial coordinates of the wave-generating region or the wave-dissipating region.

5. The wave simulation method based on a triangular numerical wave tank according to claim 1, characterized in that, The step of loading the theoretical velocity vector onto the velocity inlets corresponding to the horizontal right-angled side, the vertical right-angled side, and the inclined side includes: When an obliquely incident wave is generated, the phase delay parameter is calculated based on the spatial position of each grid node on the oblique edge. The theoretical velocity vector is then corrected using the phase delay parameter, and the corrected velocity vector is loaded into the velocity inlet. When generating multi-directional irregular waves, based on the principle of linear superposition, the theoretical velocity vectors of multiple wave components with different frequencies and different incident angles are calculated. The theoretical velocity vectors of the multiple wave components are synthesized, and the synthesized velocity vector is loaded into the velocity inlet.

6. The wave simulation method based on a triangular numerical wave tank according to claim 1, characterized in that, After updating the wavefront elevation and flow velocity distribution within the right-angled triangular computational domain, the method further includes: A structural model and an overlapping mesh enclosing the structural model are established within the right-angled triangular computational domain; Based on the iterative solution results of the wave simulation control equations, the flow field pressure distribution data within the right-angled triangular computational domain are obtained; Perform data interpolation between the overlapping mesh and the base mesh of the right-angled triangular computational domain to transfer the flow field pressure distribution data to the structural model; Integrate the flow field pressure distribution data on the surface of the structure model to output the wave load on the structure model.

7. An electronic device, characterized in that, The electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to perform the operations performed by the wave simulation method based on the triangular numerical wave tank as described in any one of claims 1 to 6.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a program that is executed by a processor as described in any one of claims 1 to 6, representing the wave simulation method based on a triangular numerical wave tank.

Images