Wave absorption method and device for dynamically adjusting relaxation area based on wave characteristics and computer storage medium
By dynamically adjusting the relaxation zone parameters, the problem of energy accumulation from reflection in numerical wave tanks is solved, achieving efficient absorption and stable simulation of waves at different frequencies, and applicable to various wave conditions.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-10-21
- Publication Date
- 2026-03-27
AI Technical Summary
Existing numerical wave flumes are difficult to adapt to different wave frequency components in a fixed-length relaxation region, especially under random waves or multi-peak spectrum wave fields where reflection still occurs. Furthermore, parameter settings rely on experience and cannot be adaptively adjusted, affecting numerical stability.
By monitoring wavefront time series, the length, shape parameters, and intensity of the relaxation zone are dynamically adjusted, and wave damping measures are optimized in conjunction with prediction mechanisms to achieve directional absorption and energy dissipation of waves of different frequencies.
It significantly reduces reflected energy, improves the stability and accuracy of numerical flumes, and is suitable for regular waves, random waves, and multi-peak spectrum wave conditions, making it widely applicable.
Smart Images
Figure CN121744963A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of numerical hydrodynamic simulation, in particular to a wave absorption method based on dynamic adjustment of relaxation zone according to wave characteristics, a device and a computer storage medium. BACKGROUND
[0002] The existing numerical wave tank generally adopts the relaxation wave absorption method, a fixed length relaxation zone is set at both ends of the tank, and the relaxation function is used to smoothly transition the numerical solution and the target solution, thereby absorbing wave energy and reducing boundary reflection. However, the main problems of the existing method are: Fixed length: the relaxation zone is usually set to 1-3 times the dominant wavelength, which is difficult to balance different wave frequency components, especially under random waves or multi-peak spectrum wave field, some frequencies still produce reflection; Parameter-dependent experience: the length, shape function and strength of the relaxation zone need to be set manually, which is difficult to adapt to different wave conditions; Non-adaptive: unable to dynamically adjust according to the real-time characteristics of the wave, resulting in the accumulation of residual reflection energy in long-time simulation, affecting the numerical stability.
[0003] Therefore, a technical solution is needed that can dynamically adjust the parameters of the relaxation zone to achieve efficient absorption of multi-frequency and long-time waves. SUMMARY
[0004] The embodiments of the present application provide a wave absorption method based on dynamic adjustment of relaxation zone according to wave characteristics, a device and a computer storage medium, which extracts reflection information and energy envelope by monitoring the time sequence of the wave surface, dynamically optimizes the parameters of the relaxation zone by combining the prediction mechanism, thereby effectively reducing reflection and improving the stability and accuracy of the numerical tank.
[0005] A wave absorption method based on dynamic adjustment of relaxation zone according to wave characteristics, comprising: S1, setting a numerical wave tank, arranging at least two monitoring points in the working area, and obtaining a time sequence of the wave surface; S2, performing incident / reflection wave decomposition and short-time spectral analysis on the time sequence of the wave surface to obtain reflection coefficients of each frequency band; S3, dynamically adjusting the length, relaxation function shape parameter and relaxation strength of the relaxation zone according to the obtained reflection coefficients; S4, obtaining the wave energy envelope, predicting the wave group propagation position, and increasing the relaxation strength in advance at the position; S5, in the numerical calculation, the flow field variables are mixed and corrected by the dynamic relaxation function to realize the gradual dissipation of wave energy.
[0006] Further, in step S1, the monitoring points are arranged linearly along the wave propagation direction at the upstream of the working area and the middle and lower reaches of the working area.
[0007] Further, in step S2, obtaining the reflection coefficient of each frequency band includes the following steps: S21, obtaining the time sequence of wave surface; S22, segmenting the time sequence of wave surface to form several time periods, and performing frequency analysis in each time period to obtain the wave energy distribution at different time and frequency; S23, establishing an incident / reflected wave decomposition model and solving the incident wave and reflected wave amplitude; S24, calculating the reflection coefficient of each frequency band to obtain the reflection coefficient varying with time and frequency.
[0008] Further, in step S3, the control adjustment strategy is divided into two time scales of closed loop, including fast loop: for making immediate adjustment to short-time reflection surge, small adjustment amount and frequent, mainly changing the length of relaxation zone and the shape parameter of relaxation function; slow loop: for periodic and global optimization, evaluating and adjusting the overall length of relaxation zone and frequency band priority, balancing the wave absorption effect and the cost of calculation / domain occupation.
[0009] Further, the process of dynamically adjusting the relaxation zone includes: increasing the length of the relaxation zone when the low-frequency wave reflection is enhanced, and increasing the relaxation strength or changing the shape of the relaxation function when the high-frequency wave reflection is enhanced.
[0010] Further, the adjustment process of the length of the relaxation zone includes: Step S31a, monitoring the low-frequency reflection coefficient to determine whether it exceeds the preset threshold in consecutive multiple time windows; Step S32a, when the low-frequency reflection continues to increase and the current length of the relaxation zone does not reach 2-3 times the dominant wavelength, triggering length adjustment; Step S33a, lengthening the end position of the relaxation zone to increase the length, and connecting the new area with the original area in a smooth transition manner; Step S34a, continue to monitor the reflection, when the low-frequency reflection decreases, keep the length; if the reflection is still high, lengthen again until the upper limit is reached; Step S35a, when the low-frequency reflection continues to be below the preset threshold for a long time, gradually shorten the length of the relaxation zone to reduce the calculation overhead.
[0011] Further, the adjustment process of the shape parameter of the relaxation function includes: Step S31b, monitoring the reflection distribution of each frequency band to determine whether the high-frequency reflection is enhanced or the low-frequency reflection is enhanced; Step S32b, when the high-frequency reflection coefficient continues to exceed the threshold, triggering the shape parameter adjustment of the relaxation function, increasing the absorption weight at the front end of the relaxation zone to make the wave be quickly dissipated when entering the initial section of the relaxation zone; Step S33b, when the low-frequency reflection coefficient continues to exceed the threshold value, triggering the relaxation function shape parameter adjustment, modifying the relaxation function to a gentle distribution, and gradually expanding the absorption process in the entire relaxation zone; Step S34b, monitoring the adjusted results, and keeping the relaxation function shape parameter unchanged when the reflection drops below the threshold value; if the reflection still fluctuates greatly, continue to fine-tune the relaxation function shape parameter until the reflection is stable.
[0012] Further, the adjustment process of the relaxation strength includes: Step S31c, monitoring the reflection coefficients of the main frequency band and the high frequency band, and judging whether they are significantly increased in a short time; Step S32c, when the high-frequency or main-frequency reflection exceeds the preset threshold value, or when it is identified that a wave group with energy exceeding the threshold value is about to arrive, triggering the relaxation strength adjustment; Step S33c, gradually increasing the relaxation strength coefficient in the predicted wave group propagation area or the reflection enhanced area to enhance the wave energy dissipation; Step S34c, preferentially increasing the relaxation strength at the surface layer position to more effectively consume high-frequency wave energy; Step S35c, when the wave group passes or the reflection returns to the normal range, gradually reducing the relaxation strength to return to the normal level.
[0013] A computer device, comprising: a memory and a processor, which are communicatively connected, the memory stores computer instructions, and the processor executes the computer instructions to perform the above method.
[0014] A computer readable storage medium, which stores computer instructions for causing a computer to perform the above method.
[0015] The above technical solutions provided by the embodiments of the present application have at least the following beneficial effects: 1. Through real-time monitoring and decomposition analysis, directional absorption of different frequency waves is realized.
[0016] 2. The wave group position is predicted using the wave energy envelope, the wave absorption effect is enhanced in advance, and the transient reflection is significantly reduced.
[0017] 3. The relaxation zone parameters are dynamically adjusted to avoid reflection energy accumulation in long-term simulation and improve simulation stability.
[0018] 4. It is suitable for regular waves, random waves and multi-peak spectrum waves, and has wide applicability.
[0019] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.
[0020] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description
[0021] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings: Figure 1 This is a flowchart of a wave absorption method based on dynamic adjustment of the relaxation zone according to wave characteristics, as disclosed in an embodiment of the present invention. Figure 2 This is a flowchart of the reflection coefficient acquisition method disclosed in an embodiment of the present invention; Figure 3 This is a flowchart illustrating the process of adjusting the length of the relaxation region as disclosed in an embodiment of the present invention; Figure 4 This is a flowchart illustrating the process of adjusting the shape parameters of the relaxation function as disclosed in an embodiment of the present invention. Figure 5 This is a flowchart illustrating the process of adjusting the relaxation strength according to an embodiment of the present invention; Figure 6 This is a flowchart of the process for enhancing relaxation strength disclosed in an embodiment of the present invention; Figure 7 This is a flowchart of the wave-damping process disclosed in an embodiment of the present invention. Detailed Implementation
[0022] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
[0023] like Figure 1 As shown, this embodiment of the invention discloses a wave absorption method based on dynamically adjusting the relaxation zone according to wave characteristics, including: S1. Set up a numerical wave flume and arrange at least two monitoring points in the working area to obtain the wave surface time series.
[0024] Generally, to maintain robustness, the number of monitoring points should be at least two, preferably three to six. The interval between monitoring points Generally , is the dominant wavelength, the interval between two points is too small in value, and too large will affect the phase decoupling; the interval between monitoring points preferably .
[0025] The monitoring points are used to collect the wave surface time series of the free surface height, and the monitoring points are arranged upstream of the working area and downstream of the working area, so as to distinguish the incident and reflected components, and the monitoring points are linearly arranged along the wave propagation direction.
[0026] The step S1 is to obtain the basic data reflecting the wave propagation and interaction by arranging the monitoring points in the working area and collecting the wave surface time series, and to provide the input for the subsequent wave decomposition and energy analysis.
[0027] S2, the incident / reflected wave decomposition and short-time spectrum analysis are performed on the wave surface time series, and the reflection coefficients of each frequency band are obtained, wherein the wave surface time series contains multiple frequency components, and the frequency band is the entire frequency range of the wave surface time series divided into several subintervals.
[0028] The step S2 is to obtain the reflection coefficients of different frequency bands by performing incident / reflected wave decomposition and short-time spectrum analysis on the wave surface time series, so as to identify the composition and intensity of the reflected wave, and to provide the basis for dynamic control of the relaxation zone.
[0029] As Figure 2 shown, obtaining the reflection coefficients of each frequency band includes the following steps: S21, obtaining the wave surface time series.
[0030] S22, segmenting the wave surface time series to form several time periods, and performing frequency analysis in each time period to obtain the wave energy distribution under different time and frequency.
[0031] S23, establishing an incident / reflected wave decomposition model and solving the incident wave and reflected wave amplitudes.
[0032] The incident / reflected wave decomposition model is the existing technology, such as 2-point / 3-point method, time domain decomposition method, weighted least squares method, etc.
[0033] S24, calculating the reflection coefficients of each frequency band to obtain the reflection coefficients varying with time and frequency.
[0034] The reflection coefficients include the main frequency band reflection coefficient, the high frequency band reflection coefficient and the low frequency band reflection coefficient.
[0035] S3, dynamically adjusting the length of the relaxation zone, the shape parameters of the relaxation function and the relaxation strength according to the obtained reflection coefficients.
[0036] The purpose of step S3 is to dynamically adjust the length of the relaxation region, the shape parameters of the relaxation function, and the relaxation intensity according to the magnitude and distribution of the reflection coefficients of each frequency band, so as to make the wave-damping measures more targeted and improve the energy dissipation efficiency.
[0037] In step S3, the control adjustment strategy is divided into two closed loops with two time scales: a fast loop, which is used to make immediate adjustments to short-term reflection surges, with small and frequent adjustments, mainly changing the length of the relaxation region and the shape parameters of the relaxation function; and a slow loop, which is used for periodic, global optimization, evaluating and adjusting the overall length of the relaxation region and the frequency band priority, and balancing the wave attenuation effect with the computation / domain occupancy cost.
[0038] The process of dynamically adjusting the relaxation region includes: increasing the length of the relaxation region when low-frequency wave reflection is enhanced, and increasing the relaxation intensity or changing the shape of the relaxation function when high-frequency wave reflection is enhanced.
[0039] like Figure 3 As shown, the process of adjusting the length of the relaxation zone includes: Step S31a: Monitor the low-frequency band reflection coefficient and determine whether it exceeds a preset threshold within multiple consecutive time windows; Step S32a: When low-frequency reflection continues to increase and the current relaxation region length does not reach 2 to 3 times the dominant wavelength, length adjustment is triggered; Step S33a: Extend the length of the end point of the relaxation area and connect the newly added area with the original area using a smooth transition method; Step S34a: Continue monitoring the reflection. When the low-frequency reflection decreases, maintain the current length. If the reflection is still high, extend it again until the upper limit is reached. Step S35a: When the low-frequency reflection remains below the preset threshold for a long period of time, the length of the relaxation region is gradually shortened to reduce computational overhead.
[0040] like Figure 4 As shown, the process of adjusting the shape parameters of the relaxation function includes: Step S31b: Monitor the reflection distribution of each frequency band and determine whether the high-frequency reflection is enhanced or the low-frequency reflection is enhanced; Step S32b: When the high-frequency reflection coefficient continues to exceed the threshold, the relaxation function shape parameter adjustment is triggered to increase the absorption weight at the front end of the relaxation zone, so that the wave is quickly dissipated when it enters the initial stage of the relaxation zone. Step S33b: When the low-frequency reflection coefficient continues to exceed the threshold, the shape parameter adjustment of the relaxation function is triggered, and the relaxation function is modified to a smooth distribution so that the absorption process gradually unfolds throughout the relaxation region. Step S34b, monitor the adjusted result, keep the relaxation function shape parameter unchanged when the reflection drops below the threshold value; if the reflection still fluctuates greatly, continue to fine-tune the relaxation function shape parameter until the reflection is stable.
[0041] As shown in Figure 5 , the adjustment process of the relaxation strength includes: Step S31c, monitor the reflection coefficients of the main frequency band and the high frequency band, and judge whether they are significantly increased in a short time; Step S32c, trigger the relaxation strength adjustment when the high frequency or main frequency reflection exceeds the preset threshold value, or when it is identified that a wave group with energy exceeding the threshold value is about to arrive; Step S33c, gradually increase the relaxation strength coefficient in the predicted wave group propagation area or the area with enhanced reflection, and enhance the wave energy dissipation; Step S34c, preferentially increase the relaxation strength at the surface position to more effectively consume high-frequency wave energy; Step S35c, gradually reduce the relaxation strength when the wave group passes or the reflection returns to the normal range, and return to the normal level.
[0042] S4, obtain the wave energy envelope, predict the wave group propagation position, and increase the relaxation strength in advance at the position.
[0043] The function of step S4 is to identify the wave group with high energy through the energy envelope and predict its propagation position and arrival time, and to increase the relaxation strength in advance before the wave group enters the relaxation zone, so as to realize the active wave dissipation control of the strong wave group, as shown in Figure 6 . S41, obtain the wave surface time sequence collected in step S1, and obtain the wave energy envelope through signal processing method, which is used to reflect the trend of wave group energy change with time; S42, identify the local peak value in the energy envelope, and judge whether it exceeds a certain proportion of the background average value, for example, more than 30% is determined as a wave group with high energy; S43, perform short-time spectral analysis on the identified wave group with high energy to obtain its dominant period or frequency, and estimate the group velocity combined with the water depth condition; S44, calculate the time for the wave group to propagate from the monitoring point to the relaxation zone according to the group velocity, predict its arrival position and time, and set the enhancement time window in advance; S45, gradually increase the relaxation strength in the predicted propagation position range, and the enhancement amplitude is proportional to the wave group energy, and gradually return to the normal level after the wave group passes.
[0044] S5, in numerical calculation, the flow field variables are mixed and corrected by the dynamic relaxation function, so as to realize the gradual dissipation of wave energy, as shown in Figure 7 . S51, in the numerical calculation process, the current flow field solution and the target solution are obtained, and a mixed coefficient varying with position and time is generated in combination with a dynamic relaxation function; S52, the main variables such as velocity field and volume fraction are weighted and mixed to gradually transit to the target solution, so that energy dissipation is realized in the relaxation zone; S53, the relaxation coefficient is smoothed and limited in amplitude in the mixing process to avoid instability caused by numerical mutation, while maintaining mass and volume conservation; S54, the relaxation coefficient is gradually reduced after the wave group passes, so that the flow field returns to the normal state, and a closed-loop control is formed by monitoring the reflection coefficient and volume error in real time.
[0045] The role of step S5 is to mix and correct the flow field variables through the dynamic relaxation function, so that the wave gradually dissipates in the relaxation zone, thereby avoiding the reflection wave into the working area and interfering with the numerical simulation accuracy.
[0046] A computer device, comprising: a memory and a processor, which are mutually connected in communication, the memory stores computer instructions, and the processor executes the computer instructions to perform the above method.
[0047] A computer readable storage medium, the computer readable storage medium stores computer instructions, the computer instructions are used to make a computer execute the above method.
[0048] It should be understood that the specific order or hierarchy of steps in the disclosed processes is an example of an exemplary method. Based on design preferences, it is understood that the specific order or hierarchy of steps in the processes can be rearranged, and further, it is understood that some of the steps can be performed concurrently, that the amount of overlap of the steps can vary from case to case, and that some of the steps can be performed sequentially. The accompanying method claims present elements of the various steps in a sample order, and are not meant to be limited to the specific order or hierarchy presented.
[0049] In the above detailed description, various features are combined in a single embodiment for simplicity. This disclosure should not be interpreted as reflecting an intention that the claimed subject matter requires more features than are explicitly stated in each claim. Rather, as the appended claims reflect, inventive subject matter lies in fewer than all features of a single disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim acting as a separate embodiment of the application.
[0050] Those of skill would further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the embodiments disclosed herein can be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the overall system. Skilled artisans can implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present disclosure.
[0051] The steps of a method or algorithm described in connection with the embodiments disclosed herein can be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module can reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such that the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium can be integral to the processor. The processor and the storage medium can reside in an ASIC. The ASIC can reside in a user terminal. In the alternative, the processor and the storage medium can reside as discrete components in a user terminal.
[0052] For a software implementation, the techniques described herein can be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. The software codes can be stored in memory units and executed by processors. The memory unit can be implemented within the processor or external to the processor, in which case it can be communicatively coupled to the processor via various means as is known in the art.
[0053] The above description includes one or more examples of the embodiments. Of course, not all possible combinations of components or method steps described above can be claimed as embodiments. One of ordinary skill in the art can recognize that modifications and variations of the described embodiments are possible and are within the scope of the present disclosure. It is therefore intended that the embodiments described herein be considered in all respects as only illustrative and not restrictive. Specifically, the description of the embodiments should be considered to be illustrative and not exhaustive, and should be considered in the light of the claims. Further, the use of the term "comprise" in the specification is to be construed in the same way as the term "comprise" as explained in the description of the claims. Also, the use of the term "or" in the claims is to be construed as "non-exclusive or" as explained in the description of the claims.
Claims
1. A wave absorption method based on dynamically adjusting the relaxation zone according to wave characteristics, characterized in that, include: S1, Set up a numerical wave flume, arrange at least two monitoring points in the working area, and obtain the wave surface time series; S2, perform incident / reflected wave decomposition and short-time spectrum analysis on the wavefront time series to obtain the reflection coefficients of each frequency band; S3, based on the obtained reflection coefficient, dynamically adjust the length of the relaxation region, the shape parameter of the relaxation function, and the relaxation intensity; S4, Obtain the wave energy envelope, predict the wave group propagation location, and enhance the relaxation intensity in advance at the said location; S5, in numerical calculations, uses a dynamic relaxation function to mix and correct the flow field variables, thereby achieving the gradual dissipation of wave energy.
2. The method as described in claim 1, characterized in that, In step S1, monitoring points are set up upstream of the working area and downstream of the working area, and the monitoring points are arranged linearly along the wave propagation direction.
3. The method as described in claim 1, characterized in that, Step S2, obtaining the reflection coefficients for each frequency band includes the following steps: S21, Obtain the wavefront time series; S22, the wave surface time series is segmented into several time periods, and frequency analysis is performed in each time period to obtain the wave energy distribution at different times and frequencies; S23, Establish the incident / reflected wave decomposition model and solve for the amplitudes of the incident and reflected waves; S24, calculate the reflection coefficient for each frequency band to obtain the reflection coefficient that varies with time and frequency.
4. The method as described in claim 1, characterized in that, In step S3, the control adjustment strategy is divided into two closed loops with two time scales: a fast loop, which is used to make immediate adjustments to short-term reflection surges, with small and frequent adjustments, mainly changing the length of the relaxation region and the shape parameters of the relaxation function; and a slow loop, which is used for periodic, global optimization, evaluating and adjusting the overall length of the relaxation region and the frequency band priority, and balancing the wave attenuation effect with the computation / domain occupancy cost.
5. The method as described in claim 4, characterized in that, The process of dynamically adjusting the relaxation region includes: increasing the length of the relaxation region when low-frequency wave reflection is enhanced, and increasing the relaxation intensity or changing the shape of the relaxation function when high-frequency wave reflection is enhanced.
6. The method as described in claim 4, characterized in that, The process of adjusting the length of the relaxation zone includes: Step S31a: Monitor the low-frequency band reflection coefficient and determine whether it exceeds a preset threshold within multiple consecutive time windows; Step S32a: When low-frequency reflection continues to increase and the current relaxation region length does not reach 2 to 3 times the dominant wavelength, length adjustment is triggered; Step S33a: Extend the length of the end point of the relaxation area and connect the newly added area with the original area using a smooth transition method; Step S34a: Continue monitoring the reflection. When the low-frequency reflection decreases, maintain the current length. If the reflection is still high, extend it again until the upper limit is reached. Step S35a: When the low-frequency reflection remains below the preset threshold for a long period of time, the length of the relaxation region is gradually shortened to reduce computational overhead.
7. The method as described in claim 4, characterized in that, The process of adjusting the shape parameters of the relaxation function includes: Step S31b: Monitor the reflection distribution of each frequency band and determine whether the high-frequency reflection is enhanced or the low-frequency reflection is enhanced; Step S32b: When the high-frequency reflection coefficient continues to exceed the threshold, the relaxation function shape parameter adjustment is triggered to increase the absorption weight at the front end of the relaxation zone, so that the wave is quickly dissipated when it enters the initial stage of the relaxation zone. Step S33b: When the low-frequency reflection coefficient continues to exceed the threshold, the shape parameter adjustment of the relaxation function is triggered, and the relaxation function is modified to a smooth distribution so that the absorption process gradually unfolds throughout the relaxation region. Step S34b: Monitor the adjusted result. When the reflection drops below the threshold, keep the shape parameter of the relaxation function unchanged. If the reflection still fluctuates greatly, continue to fine-tune the shape parameter of the relaxation function until the reflection stabilizes.
8. The method as described in claim 4, characterized in that, The process of adjusting the relaxation strength includes: Step S31c: Monitor the reflection coefficients of the main frequency band and the high frequency band to determine whether they increase significantly in a short period of time; Step S32c: When the high-frequency or dominant frequency reflection exceeds the preset threshold, or when a wave group with energy exceeding the threshold is detected to be about to arrive, the relaxation intensity adjustment is triggered. Step S33c: In the predicted wave group propagation area or the area of enhanced reflection, gradually increase the relaxation intensity coefficient to enhance wave energy dissipation. Step S34c prioritizes increasing the relaxation intensity at the surface level to more effectively dissipate high-frequency wave energy; In step S35c, when the wave group passes through or reflects back to the normal range, the relaxation intensity is gradually reduced to return to the normal level.
9. A computer device, characterized in that, include: A memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, the processor executing the computer instructions to perform the method of any one of claims 1 to 8.
10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions for causing the computer to perform the method of any one of claims 1 to 8.