Gravitational wave signal acquisition method and device, equipment and medium

Abstract

The invention relates to the technical field of celestial physics, and discloses a gravitational wave signal acquisition method, device and equipment and a medium, and the method comprises the steps: obtaining a first orbit parameter, a first quality parameter and a first structure parameter of double white dwarf stars; determining a target radius based on the first orbit parameter and the first quality parameter; comparing a target parameter in the first structure parameters with a target radius, and determining a second quality parameter, a second structure parameter and a second track parameter; based on the second mass parameter and the second orbit parameter, obtaining a double-white-dwarf-star gravitational wave signal or a double-white-dwarf-star gravitational wave signal and a suction disc gravitational wave signal; and judging whether a termination condition is met or not, and if the termination condition is not met, repeatedly executing from determining the target radius based on the first orbit parameter and the first quality parameter until the termination condition is met, so as to obtain a gravitational wave signal set. The gravitational wave signals can be accurately obtained by updating the orbit parameters, the quality parameters and the structure parameters.

Description

Technical Field

[0001] The present disclosure relates to the field of astrophysical technologies, and specifically relates to a method, apparatus, device, and medium for obtaining gravitational wave signals. Background Art

[0002] Space gravitational wave detectors focus on the gravitational wave signal frequency band of about 0.1 - 1000 mHz. With the development of science and technology, gravitational wave detectors have entered a new stage in data processing. Matched filtering is an important data processing method for gravitational wave detectors, and binary white dwarf gravitational wave signals are important input data for matched filtering. Therefore, efficiently calculating a large number of resolvable binary white dwarf gravitational wave signals and establishing a signal template library are the current research focuses.

[0003] Currently, the applicable frequency band of the data processing pipeline applied by the ground-based Laser Interferometer Gravitational-Wave Observatory is about 1 - 1000 Hz. For space gravitational wave detectors in the frequency band of about 0.1 - 1000 mHz, the binary white dwarf theoretical templates adopted by related technologies mainly use analytical calculations, and cannot systematically consider the influence of the orbital motion of binary white dwarfs and the structural evolution of white dwarfs on binary white dwarf gravitational wave signals when angular momentum is conserved, resulting in low accuracy of the obtained gravitational wave signals. Summary of the Invention

[0004] In view of this, the present disclosure provides a method, apparatus, device, and medium for obtaining gravitational wave signals to solve the problem of low accuracy of gravitational wave signals.

[0005] In a first aspect, the present disclosure provides a method for obtaining gravitational wave signals, the method comprising:

[0006] Obtaining first orbital parameters, first mass parameters, and first structural parameters of a binary white dwarf;

[0007] Determining a target radius based on the first orbital parameters and the first mass parameters;

[0008] Comparing target parameters in the first structural parameters with the target radius to determine second mass parameters, second structural parameters, and second orbital parameters;

[0009] Obtaining a binary white dwarf gravitational wave signal based on the second mass parameters and the second orbital parameters, or obtaining a binary white dwarf gravitational wave signal and an accretion disk gravitational wave signal based on the second mass parameters and the second orbital parameters, and storing the binary white dwarf gravitational wave signal or the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal in a gravitational wave signal set;

[0010] Determine whether the termination condition is reached. If the termination condition is not reached, start repeating from determining the target radius based on the first orbital parameter and the first mass parameter until the termination condition is reached, and obtain the gravitational wave signal set to be output. Among them, the termination condition is jointly determined by the first condition and the second condition. The first condition is jointly determined by the current mass parameter, the current structure parameter, and the current orbital parameter. The second condition is used to indicate the comparison result between the current time step and the preset time step. The gravitational wave signal set contains the double white dwarf gravitational wave signal and the accretion disk gravitational wave signal determined each time.

[0011] In the embodiments of the present disclosure, by obtaining the first orbital parameter, the first mass parameter, and the first structure parameter of the double white dwarf; determining the target radius based on the first orbital parameter and the first mass parameter; comparing the target parameter in the first structure parameter with the target radius to determine the second mass parameter, the second structure parameter, and the second orbital parameter; based on the second mass parameter and the second orbital parameter, obtaining the double white dwarf gravitational wave signal or the double white dwarf gravitational wave signal and the accretion disk gravitational wave signal; determining whether the termination condition is reached. If the termination condition is not reached, start repeating from determining the target radius based on the first orbital parameter and the first mass parameter until the termination condition is reached, and obtain the gravitational wave signal set to be output. Since the embodiments of the present disclosure can accurately obtain the gravitational wave signal by updating the orbital parameter, the mass parameter, and the structure parameter.

[0012] In an alternative embodiment, comparing the target parameter in the first structure parameter with the target radius to determine the second mass parameter, the second structure parameter, and the second orbital parameter includes:

[0013] When the target parameter is less than the target radius, determine the first mass parameter as the second mass parameter, determine the first structure parameter as the second structure parameter, and determine the first orbital parameter as the second orbital parameter;

[0014] When the target parameter is greater than or equal to the target radius, update the first mass parameter, the first structure parameter, and the first orbital parameter to obtain the second mass parameter, the second structure parameter, and the second orbital parameter.

[0015] In the embodiments of the present disclosure, by correspondingly processing the first orbital parameter, the first mass parameter, and the first structure parameter according to the size relationship between the target parameter and the target radius, the correction of the orbital parameter, the mass parameter, and the structure parameter can be realized, and the physical state changes of the double white dwarf at different evolution stages can be accurately tracked.

[0016] In an alternative embodiment, updating the first mass parameter, the first structure parameter, and the first orbital parameter to obtain the second mass parameter, the second structure parameter, and the second orbital parameter includes:

[0017] Based on the first orbital parameter, the first mass parameter, the first structural parameter, and the target radius, an accretion material proportion coefficient is obtained;

[0018] Based on the accretion material proportion coefficient, the first mass parameter and the first orbital parameter are updated to obtain a second mass parameter and a second orbital parameter;

[0019] The second mass parameter is matched with a preset mapping relationship to update the first structural parameter to obtain a second structural parameter, where the preset mapping relationship is used to represent the mapping relationship between the mass parameter and the structural parameter.

[0020] In the embodiment of the present disclosure, by updating the first mass parameter and the first orbital parameter, the correction of the mass parameter and the orbital parameter can be achieved. By using the preset mapping relationship and the second mass parameter to update the first structural parameter, the correction of the first structural parameter can be achieved, improving the efficiency of obtaining the gravitational wave signal.

[0021] In an alternative embodiment, based on the second mass parameter and the second orbital parameter, a binary white dwarf gravitational wave signal and an accretion disk gravitational wave signal are obtained, including:

[0022] In the case where the accretion material proportion coefficient is greater than a preset value, based on the second mass parameter and the second orbital parameter, an accretion disk parameter is obtained, where the preset value is used to determine whether an accretion disk is formed by the binary white dwarf;

[0023] Based on the accretion disk parameter, the second mass parameter, and the second orbital parameter, an accretion disk gravitational wave signal is obtained;

[0024] Based on the second mass parameter and the second orbital parameter, a binary white dwarf gravitational wave signal is obtained.

[0025] In the embodiment of the present disclosure, by obtaining the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal based on the second mass parameter and the second orbital parameter in the case where the accretion material proportion coefficient is greater than the preset value, the gravitational wave signal can be accurately obtained, improving the accuracy and reliability of the gravitational wave signal.

[0026] In an alternative embodiment, the method further includes:

[0027] In the case where the accretion material proportion coefficient is less than or equal to the preset value, based on the second mass parameter and the second orbital parameter, a binary white dwarf gravitational wave signal is obtained.

[0028] In the embodiments of the present disclosure, when the accretion material proportion coefficient is less than or equal to a preset value, a gravitational wave signal of a double white dwarf is obtained based on a second mass parameter and a second orbital parameter, so that the gravitational wave signal can be accurately obtained, and the accuracy and reliability of the gravitational wave signal are improved.

[0029] In an alternative embodiment, obtaining the gravitational wave signal of the double white dwarf based on the second mass parameter and the second orbital parameter includes:

[0030] Based on the second mass parameter and the second orbital parameter, obtaining the mass quadrupole of the double white dwarf;

[0031] Based on the mass quadrupole of the double white dwarf, obtaining the gravitational wave amplitude and the gravitational wave frequency of the double white dwarf;

[0032] Based on the gravitational wave amplitude and the gravitational wave frequency of the double white dwarf, obtaining the gravitational wave power of the double white dwarf;

[0033] Performing frequency-domain analysis on the gravitational wave amplitude, the gravitational wave frequency, and the gravitational wave power of the double white dwarf to obtain the gravitational wave signal of the double white dwarf.

[0034] In the embodiments of the present disclosure, by obtaining the gravitational wave amplitude, the gravitational wave frequency, and the gravitational wave power of the double white dwarf based on the second mass parameter and the second orbital parameter, and performing frequency-domain analysis to obtain the gravitational wave signal of the double white dwarf, the gravitational wave signal of the double white dwarf can be accurately obtained, and the accuracy and reliability of the gravitational wave signal of the double white dwarf are improved.

[0035] In an alternative embodiment, after determining whether the termination condition is reached, the method further includes:

[0036] When the termination condition is reached, outputting a set of gravitational wave signals.

[0037] In the embodiments of the present disclosure, by outputting a set of gravitational wave signals when the termination condition is reached, the gravitational wave signal of the double white dwarf can be accurately obtained, and a complete set of gravitational wave signals can be obtained.

[0038] In a second aspect, the present disclosure provides an apparatus for obtaining a gravitational wave signal, the apparatus including:

[0039] An acquisition module, configured to acquire a first orbital parameter, a first mass parameter, and a first structure parameter of a double white dwarf;

[0040] A first determination module, configured to determine a target radius based on the first orbital parameter and the first mass parameter;

[0041] A second determination module, configured to compare a target parameter in the first structural parameter with a target radius to determine a second mass parameter, a second structural parameter, and a second orbital parameter;

[0042] A first obtaining module, configured to obtain a binary white dwarf gravitational wave signal based on the second mass parameter and the second orbital parameter, or obtain a binary white dwarf gravitational wave signal and an accretion disk gravitational wave signal based on the second mass parameter and the second orbital parameter, and store the binary white dwarf gravitational wave signal or the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal into a gravitational wave signal set;

[0043] A second obtaining module, configured to determine whether a termination condition is reached. If the termination condition is not reached, start repeating from determining the target radius based on the first orbital parameter and the first mass parameter until the termination condition is reached, to obtain a gravitational wave signal set to be output, where the termination condition is jointly determined by a first condition and a second condition, the first condition is jointly determined by a current mass parameter, a current structural parameter, and a current orbital parameter, and the second condition is used to indicate a comparison result between the current time step and a preset time step. The gravitational wave signal set contains the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal determined each time.

[0044] In a third aspect, the present disclosure provides a computer device, including: a memory and a processor, which are communicatively connected to each other. The memory stores computer instructions, and the processor executes the computer instructions to execute the method for obtaining a gravitational wave signal according to the first aspect or any corresponding implementation manner thereof.

[0045] In a fourth aspect, the present disclosure provides a computer-readable storage medium, on which computer instructions are stored, and the computer instructions are used to cause a computer to execute the method for obtaining a gravitational wave signal according to the first aspect or any corresponding implementation manner thereof.

[0046] In a fifth aspect, the present disclosure provides a computer program product, including computer instructions, and the computer instructions are used to cause a computer to execute the method for obtaining a gravitational wave signal according to the first aspect or any corresponding implementation manner thereof. Description of the Drawings

[0047] In order to more clearly illustrate the specific embodiments of the present disclosure or the technical solutions in the prior art, the following will briefly introduce the drawings required for use in the description of the specific embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of the present disclosure. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0048] Figure 1It is a schematic flowchart of a method for obtaining gravitational wave signals according to an embodiment of the present disclosure;

[0049] Figure 2 It is a schematic flowchart of another method for obtaining gravitational wave signals according to an embodiment of the present disclosure;

[0050] Figure 3 It is a schematic diagram showing the evolution of the gravitational wave amplitude of a double white dwarf over time according to an embodiment of the present disclosure;

[0051] Figure 4 It is a schematic diagram showing the evolution of the gravitational wave power spectrum of a double white dwarf over time according to an embodiment of the present disclosure;

[0052] Figure 5 It is a schematic diagram showing the evolution of the gravitational wave amplitude of another double white dwarf over time according to an embodiment of the present disclosure;

[0053] Figure 6 It is a schematic diagram showing the evolution of the gravitational wave power spectrum of another double white dwarf over time according to an embodiment of the present disclosure;

[0054] Figure 7 It is a schematic diagram showing the evolution of the gravitational wave amplitude of another double white dwarf over time according to an embodiment of the present disclosure;

[0055] Figure 8 It is a schematic diagram showing the evolution of the gravitational wave power spectrum of another double white dwarf over time according to an embodiment of the present disclosure;

[0056] Figure 9 It is a schematic diagram showing the evolution of the gravitational wave amplitude of yet another double white dwarf over time according to an embodiment of the present disclosure;

[0057] Figure 10 It is a schematic diagram showing the evolution of the gravitational wave power spectrum of yet another double white dwarf over time according to an embodiment of the present disclosure;

[0058] Figure 11 It is a block diagram of the structure of an apparatus for obtaining gravitational wave signals according to an embodiment of the present disclosure;

[0059] Figure 12 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present disclosure. Detailed implementation manners

[0060] To make the objectives, technical solutions, and advantages of the embodiments of the present disclosure clearer, the technical solutions in the embodiments of the present disclosure will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present disclosure. Apparently, the described embodiments are some, but not all, of the embodiments of the present disclosure. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present disclosure without creative efforts shall fall within the protection scope of the present disclosure.

[0061] Space gravitational wave detectors focus on the gravitational wave signal frequency band of about 0.1 - 1000 mHz. With the development of science and technology, gravitational wave detectors have entered a new stage in data processing. Matched filtering is an important data processing method for gravitational wave detectors, and the gravitational wave signal of double white dwarfs is an important input data for matched filtering. Therefore, efficiently calculating a large number of resolvable gravitational wave signals of double white dwarfs and establishing a signal template library are the current research focuses.

[0062] Currently, the applicable frequency band of the data processing pipeline applied by ground-based laser interferometer gravitational wave observatories is about 1 - 1000 Hz. For space gravitational wave detectors in the frequency band of about 0.1 - 1000 mHz, the theoretical templates of double white dwarfs used in related technologies mainly use analytical calculations and cannot systematically consider the influence of the orbital motion of double white dwarfs and the structural evolution of white dwarfs on the gravitational wave signal of double white dwarfs when considering angular momentum conservation, resulting in low accuracy of the obtained gravitational wave signal. In addition, the numerical relativity method has a long simulation calculation time and low efficiency for the merger of double white dwarfs.

[0063] To solve the above problems, according to the embodiments of the present disclosure, an embodiment of a method for obtaining a gravitational wave signal is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although the logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in a different order than here.

[0064] In this embodiment, a method for obtaining a gravitational wave signal is provided, as Figure 1 shown, Figure 1 is a schematic flowchart of a method for obtaining a gravitational wave signal according to the embodiments of the present disclosure. This process can be applied to a server and includes the following steps:

[0065] Step S101, obtain the first orbital parameters, first mass parameters, and first structural parameters of the double white dwarfs.

[0066] Optionally, in the embodiments of the present disclosure, the first orbital parameters are included in the orbital parameters, and the orbital parameters include the orbital spacing, orbital period, ellipticity, orbital phase, and orbital inclination (when the normal of the orbital plane faces the observer, the orbital inclination is zero) of the double white dwarfs, etc. The first mass parameters are included in the mass parameters, and the mass parameters include the masses of the double white dwarfs respectively. The first structural parameters are included in the structural parameters, and the structural parameters include the radii, densities, and pressures of the double white dwarfs respectively.

[0067] Specifically, as Figure 2 shown, the server can use the observed parameters of the double white dwarf identification source as the parameter values, or use the default parameter values, or set the parameter values according to the parameter range in the program.

[0068] It should be noted that the double white dwarf identification source refers to a double white dwarf celestial system that has been observed and for which observation parameters have been given. Based on the observation parameters and some assumptions, the gravitational wave signals of these identification sources can be obtained, providing a reference basis for future space gravitational wave detection. Among them, the observation parameters of the double white dwarf identification source are obtained through observations of the double white dwarf celestial system using telescopes in the optical, infrared, radio, or X-ray bands.

[0069] In addition, after the server obtains the first orbital parameter, the first mass parameter, and the first structure parameter, it can set the initial conditions for the evolution of the double white dwarf according to these parameters, that is, whether mass transfer has occurred initially. The server can also pre-obtain some physical constants, including the solar mass, solar radius, solar luminosity, gravitational constant, pi, speed of light, electron mass, and Planck constant, etc.

[0070] Step S102, determine the target radius based on the first orbital parameter and the first mass parameter.

[0071] Optionally, in the embodiments of the present disclosure, the target radius refers to the Roche lobe. In a double white dwarf system, the Roche lobe is a teardrop-shaped region that surrounds each star in the double white dwarf. When the stellar matter exceeds the range of its Roche lobe, this matter is no longer completely bound by the gravity of that star, and thus mass transfer occurs.

[0072] Specifically, the server obtains the mass ratio of the double white dwarf based on the first mass parameter, and then obtains the Roche lobe, that is, the target radius, based on the mass ratio of the double white dwarf and the orbital separation in the first orbital parameter.

[0073] Step S103, compare the target parameter in the first structure parameter with the target radius to determine the second mass parameter, the second structure parameter, and the second orbital parameter.

[0074] Optionally, in the embodiments of the present disclosure, the target parameter refers to the respective radii of the double white dwarfs in the first structure parameter. The second mass parameter is included in the mass parameter, the second structure parameter is included in the structure parameter, and the second orbital parameter is included in the orbital parameter.

[0075] Specifically, as Figure 2 shown, the server compares the target parameter in the first structure parameter with the target radius to determine whether the radii of the double white dwarfs fill the Roche lobe. When the radii of all the stars in the double white dwarf do not exceed their equivalent Roche lobes, the server determines the first mass parameter as the second mass parameter, the first structure parameter as the second structure parameter, and the first orbital parameter as the second orbital parameter;

[0076] When the radius of any one of the two white dwarfs exceeds its equivalent Roche lobe, the part of the star that exceeds the Roche lobe will break away from the star, resulting in mass transfer between the two white dwarfs. The server updates the first mass parameter, the first structure parameter, and the first orbital parameter to obtain the second mass parameter, the second structure parameter, and the second orbital parameter.

[0077] Step S104: Based on the second mass parameter and the second orbital parameter, obtain the gravitational wave signal of the double white dwarf, or based on the second mass parameter and the second orbital parameter, obtain the gravitational wave signal of the double white dwarf and the gravitational wave signal of the accretion disk, and store the gravitational wave signal of the double white dwarf or the gravitational wave signal of the double white dwarf and the gravitational wave signal of the accretion disk into the gravitational wave signal set.

[0078] Optionally, in the embodiments of the present disclosure, the gravitational wave signal of the double white dwarf includes the evolution of the gravitational wave amplitude, frequency, and power of the double white dwarf over time. The gravitational wave signal of the accretion disk includes the evolution of the gravitational wave amplitude, frequency, and power of the accretion disk over time. The gravitational wave signal set contains the gravitational wave signal of the double white dwarf and the gravitational wave signal of the accretion disk.

[0079] As Figures 3 - 10 shown, the server selects 4 stars from the double white dwarf sources obtained from astronomical observations, sets different double white dwarf masses and initial orbital periods for these four stars, and sets the same ellipticity (such as 0), distance (such as 10 kiloparsecs), magnetic field (such as 10 10 Gauss), azimuth (such as 0°), pitch (such as 60°), orbital phase (such as 0°), and angular momentum transfer coefficient (such as 1.0).

[0080] Among them, Figure 3 and Figure 4 are respectively schematic diagrams of the evolution of the gravitational wave amplitude and power spectrum of the double white dwarf over time. The corresponding double white dwarf masses are 0.55 times the solar mass and 0.27 times the solar mass, and the orbital period is about 321.6 seconds. Figure 5 and Figure 6 are respectively schematic diagrams of the evolution of the gravitational wave amplitude and power spectrum of another double white dwarf over time. The corresponding double white dwarf masses are 1.17 times the solar mass and 0.17 times the solar mass, and the orbital period is about 285.2 seconds. Figure 7 and Figure 8 are respectively schematic diagrams of the evolution of the gravitational wave amplitude and power spectrum of another double white dwarf over time. The corresponding double white dwarf masses are 1.38 times the solar mass and 1.35 times the solar mass, and the orbital period is about 5.7 seconds. Figure 9 and Figure 10 are respectively schematic diagrams of the evolution of the gravitational wave amplitude and power spectrum of yet another double white dwarf over time. The corresponding double white dwarf masses are 0.55 times the solar mass and 0.27 times the solar mass, and the orbital period is about 321.6 seconds.

[0081] It should be noted that an accretion disk is a disk-shaped structure composed of diffuse matter around celestial bodies, and these substances move in orbital motion under the action of the gravitational field of the central body. The central body includes young stars, protostars, white dwarfs, neutron stars, and black holes, etc. In a binary white dwarf system, if matter transfer occurs between the binary white dwarfs and the proportion coefficient of the accreted matter is greater than the preset value, part of the accreted matter will form an accretion disk. For example, in a close binary white dwarf system, when the more massive star first evolves into a white dwarf and the less massive companion star evolves to the giant stage and exceeds its Roche lobe, gas will flow from the companion star to the white dwarf, thereby forming an accretion disk.

[0082] Specifically, as Figure 2 shown, when the radius of the binary white dwarfs does not fill the Roche lobe, the server obtains the gravitational wave signal of the binary white dwarfs based on the second mass parameter and the second orbital parameter, and stores the gravitational wave signal of the binary white dwarfs in the gravitational wave signal set;

[0083] When the radius of the binary white dwarfs fills the Roche lobe and no accretion disk is formed, the server obtains the gravitational wave signal of the binary white dwarfs based on the second mass parameter and the second orbital parameter, and stores the gravitational wave signal of the binary white dwarfs in the gravitational wave signal set;

[0084] When the radius of the binary white dwarfs fills the Roche lobe and an accretion disk is formed, the server obtains the gravitational wave signal of the binary white dwarfs and the gravitational wave signal of the accretion disk based on the second mass parameter and the second orbital parameter, and stores the gravitational wave signal of the binary white dwarfs and the gravitational wave signal of the accretion disk in the gravitational wave signal set.

[0085] Step S105, determine whether the termination condition is reached. If the termination condition is not reached, start repeating from determining the target radius based on the first orbital parameter and the first mass parameter until the termination condition is reached, and obtain the gravitational wave signal set to be output. Among them, the termination condition is jointly determined by the first condition and the second condition. The first condition is jointly determined by the current mass parameter, the current structure parameter, and the current orbital parameter. The second condition is used to indicate the comparison result between the current time step and the preset time step. The gravitational wave signal set contains the gravitational wave signal of the binary white dwarfs and the gravitational wave signal of the accretion disk determined each time.

[0086] Optionally, in the embodiments of the present disclosure, the termination condition includes the first condition and the second condition. Among them, the first condition refers to a phase change between the binary white dwarfs or the merger of the binary white dwarfs, and the second condition refers to that the current time step reaches the preset time step. The preset time step refers to the time step preset to end the program loop.

[0087] Among them, the phase changes include four situations: (1) Under normal circumstances, the electron degeneracy pressure inside a white dwarf is balanced with the gravitational pressure, so the white dwarf can be in an equilibrium state. If accretion causes the gravitational pressure to be too large, electrons may be forced into the atomic nucleus and combine with protons to become neutrons. When the neutron degeneracy pressure can balance the stellar gravitational pressure, the white dwarf will become a neutron star, and at this time, the double white dwarf system will become a binary system of a neutron star and a white dwarf; (2) If the star loses mass too quickly in a short period of time or undergoes a nuclear reaction inside, the degenerate state of the electrons will end, that is, the electrons will become non-degenerate. At this time, the properties of the white dwarf will change and it will no longer be the original white dwarf; (3) After a white dwarf steadily loses too much mass, it will gradually evolve into a planet; (4) The double white dwarfs collide.

[0088] The merger of double white dwarfs refers to the fusion phenomenon that occurs when two white dwarfs orbit each other. When the sum of the radii of the double white dwarfs is greater than the orbital spacing, it is determined that the double white dwarfs have merged.

[0089] The time unit step size can be selected as a value from seconds to thousands of years. The smaller the step size, the higher the accuracy. Generally, a small step size corresponds to a white dwarf with a large mass, and a large time step size corresponds to a white dwarf with a small mass. For example, when the mass of the white dwarf is set to 1.2 times the solar mass, the time unit step size can be taken as about 300 seconds; when the mass of the white dwarf is set to 0.05 times the solar mass, the time unit step size can be taken as about 10 years.

[0090] The server accumulates the time unit step size in each loop to obtain the current time step size. When the current time step size reaches the preset time step size, it is determined that the termination condition is reached.

[0091] Specifically, the server determines whether the termination condition is reached. If none of the above termination conditions are reached, it starts to repeat the execution from determining the target radius based on the first orbital parameter and the first mass parameter, obtains the corresponding gravitational wave signals of the double white dwarfs and the accretion disk gravitational wave signals, and stores the corresponding gravitational wave signals of the double white dwarfs and the accretion disk gravitational wave signals in the gravitational wave signal set until one of the above termination conditions is reached, obtaining the gravitational wave signal set to be output.

[0092] In the embodiments of the present disclosure, by obtaining the first orbital parameters, the first mass parameters, and the first structural parameters of the double white dwarf; determining the target radius based on the first orbital parameters and the first mass parameters; comparing the target parameters in the first structural parameters with the target radius to determine the second mass parameters, the second structural parameters, and the second orbital parameters; based on the second mass parameters and the second orbital parameters, obtaining the gravitational wave signal of the double white dwarf or the gravitational wave signal of the double white dwarf and the gravitational wave signal of the accretion disk; determining whether the termination condition is reached, and in the case where the termination condition is not reached, repeating the execution from determining the target radius based on the first orbital parameters and the first mass parameters until the termination condition is reached, obtaining the set of gravitational wave signals to be output. Since the embodiments of the present disclosure can accurately obtain the gravitational wave signal by updating the orbital parameters, the mass parameters, and the structural parameters.

[0093] In some alternative embodiments, comparing the target parameters in the first structural parameters with the target radius to determine the second mass parameters, the second structural parameters, and the second orbital parameters includes:

[0094] In the case where the target parameters are less than the target radius, determining the first mass parameters as the second mass parameters, determining the first structural parameters as the second structural parameters, and determining the first orbital parameters as the second orbital parameters;

[0095] In the case where the target parameters are greater than or equal to the target radius, updating the first mass parameters, the first structural parameters, and the first orbital parameters to obtain the second mass parameters, the second structural parameters, and the second orbital parameters.

[0096] Optionally, in the embodiments of the present disclosure, as Figure 2 shown, the server compares the target parameters (i.e., the radius of the double white dwarf) with the target radius (i.e., the Roche lobe) to determine whether the radius of the double white dwarf fills the Roche lobe. In the case where the target parameters are less than the target radius (i.e., the radii of all the stars in the double white dwarf do not exceed their equivalent Roche lobes), the server determines the first mass parameters as the second mass parameters, determines the first structural parameters as the second structural parameters, and determines the first orbital parameters as the second orbital parameters;

[0097] In the case where the target parameters are greater than or equal to the target radius (i.e., the radius of any one of the stars in the double white dwarf exceeds its equivalent Roche lobe), the server updates the first mass parameters, the first structural parameters, and the first orbital parameters to obtain the second mass parameters, the second structural parameters, and the second orbital parameters.

[0098] In the embodiments of the present disclosure, by correspondingly processing the orbital parameters, mass parameters, and structural parameters according to the magnitude relationship between the target parameters and the target radius, it is possible to correct the orbital parameters, mass parameters, and structural parameters, and accurately track the changes in the physical states of double white dwarfs at different evolutionary stages.

[0099] In some alternative embodiments, updating the first mass parameter, the first structural parameter, and the first orbital parameter to obtain the second mass parameter, the second structural parameter, and the second orbital parameter includes:

[0100] Based on the first orbital parameter, the first mass parameter, the first structural parameter, and the target radius, obtain the proportion coefficient of the accreted matter;

[0101] Based on the proportion coefficient of the accreted matter, update the first mass parameter and the first orbital parameter to obtain the second mass parameter and the second orbital parameter;

[0102] Match the second mass parameter with the preset mapping relationship to update the first structural parameter to obtain the second structural parameter, where the preset mapping relationship is used to represent the mapping relationship between the mass parameter and the structural parameter.

[0103] Optionally, in the embodiments of the present disclosure, the proportion coefficient of the accreted matter refers to the proportion of the mass of the accreted matter of the primary star to the mass of the overflowing matter of the companion star, and the value range is from 0 to 1. Among them, the companion star is the star that provides matter and overflows matter outward when the radius exceeds the Roche lobe radius; the primary star is the accreted star that receives part of the matter overflowing from the companion star. The preset mapping relationship refers to the mapping relationship between the mass parameter and the structural parameter, including but not limited to forms such as numerical tables, scatter plots, line charts, or surface charts, function equations, and databases.

[0104] It should be noted that the result after the white dwarf fills the Roche lobe is that matter overflows from the companion star, part of it is transferred to the primary star, and part of it is transferred to the region outside the binary star. When the white dwarf fills the Roche lobe, the server can calculate what proportion of the matter is transferred to the primary star and what proportion of the matter is transferred outside the binary star by comparing the binary star luminosity and the Eddington luminosity: If the binary star luminosity is less than the Eddington luminosity, all the detached matter will be accreted by the primary star; if the binary star luminosity is greater than the Eddington luminosity, calculate how much matter is accreted according to the law of conservation of energy, and the remaining matter will enter other regions outside the binary star.

[0105] Specifically, as Figure 2As shown, the server first obtains the mass of the overflowed matter of the companion star in the binary white dwarf system according to the first orbital parameter, the first mass parameter, the first structural parameter and the target radius, and calculates the mass transfer rate in the binary white dwarf system, i.e., the binary star luminosity, by the overflowed matter mass of the companion star and the time unit step. Then, the server compares the binary star luminosity with the Eddington luminosity to obtain the ratio of the accreted matter mass of the primary star to the overflowed matter mass of the companion star, i.e., the accretion matter ratio coefficient. Next, the server calculates the mass changes of the companion star and the primary star respectively according to the overflowed matter mass of the companion star and the accretion matter ratio coefficient to obtain the second mass parameter, and updates the first orbital parameter to obtain the second orbital parameter.

[0106] In addition, the server can pre-establish a mapping relationship between mass parameters and structural parameters, that is, a preset mapping relationship. The server first obtains the general relationship between the internal pressure and density of a white dwarf based on the laws of quantum mechanics and relativity derived by Chandrasekhar, that is, the equation of state. Then, the server obtains the fluid statics equation based on the equation of state, and writes the fluid statics equation as the differential of pressure with respect to radius equal to the differential of gravitational potential with respect to radius multiplied by the negative value of density. Next, the server combines the fluid statics equation and the Poisson equation. The combined equation is that the second-order derivative of the gravitational potential with respect to radius is equal to 4 times pi multiplied by the universal gravitational constant multiplied by the density. By introducing a function related to the momentum of the electron, the gravitational potential is eliminated in the combined equation, and the degenerate second-order differential equation for radius is obtained. Finally, the server uses the difference method to solve the second-order differential equation, dividing the star into multiple layers along the radius. The physical properties of each layer are approximately the same. The solution starts from the center of the star (layer 0 and layer 1), and the second layer is solved according to the relationship of the second-order differential equation. Then, it iterates layer by layer until the solution is reached on the surface of the star. At this time, the mass of the star basically no longer changes with the increase of radius, and the preset mapping relationship is obtained.

[0107] like Figure 2 As shown, after obtaining the second quality parameter and the preset mapping relationship, the server obtains the structural parameter matching the second quality parameter in the preset mapping relationship, updates the first structural parameter, and obtains the second structural parameter.

[0108] In the embodiment of the present disclosure, by updating the first mass parameter and the first orbital parameter, the mass parameter and the orbital parameter can be corrected. By updating the first structural parameter using a preset mapping relationship and the second mass parameter, the first structural parameter can be corrected, thereby improving the efficiency of acquiring gravitational wave signals.

[0109] In some optional implementations, obtaining a binary white dwarf gravitational wave signal and an accretion disk gravitational wave signal based on the second mass parameter and the second orbital parameter includes:

[0110] When the accretion mass ratio coefficient is greater than a preset value, based on the second mass parameter and the second orbital parameter, obtain the accretion disk parameters, where the preset value is used to determine whether an accretion disk is formed in the double white dwarf;

[0111] Based on the accretion disk parameters, the second mass parameter, and the second orbital parameter, obtain the accretion disk gravitational wave signal;

[0112] Based on the second mass parameter and the second orbital parameter, obtain the double white dwarf gravitational wave signal.

[0113] Optionally, in the embodiments of the present disclosure, the preset value is used to determine whether an accretion disk is formed in the double white dwarf, and the value can be 0. An accretion disk is a disk-shaped structure composed of diffuse matter that rotates around a central body. If the accretion mass ratio coefficient is greater than the preset value, part of the accreted matter will form an accretion disk, and the other part will collide with the surface of the primary star. It is defaulted that the mass of this part of the matter is very small and can be ignored. The accretion disk parameters include the accretion disk mass parameter, the accretion disk radius parameter, etc.

[0114] Specifically, as Figure 2 shown, the server compares the accretion mass ratio coefficient with the preset value (such as 0) to determine whether an accretion disk is formed. When the accretion mass ratio coefficient is greater than the preset value (i.e., an accretion disk is formed), the server first obtains the accretion disk parameters based on the second mass parameter and the second orbital parameter, then obtains the accretion disk gravitational wave signal based on the accretion disk parameters, the second mass parameter, and the second orbital parameter, and finally obtains the double white dwarf gravitational wave signal based on the second mass parameter and the second orbital parameter.

[0115] In the embodiments of the present disclosure, by obtaining the double white dwarf gravitational wave signal and the accretion disk gravitational wave signal based on the second mass parameter and the second orbital parameter when the accretion mass ratio coefficient is greater than the preset value, the gravitational wave signal can be accurately obtained, improving the accuracy and reliability of the gravitational wave signal.

[0116] In some alternative embodiments, the method further includes:

[0117] When the accretion mass ratio coefficient is less than or equal to the preset value, based on the second mass parameter and the second orbital parameter, obtain the double white dwarf gravitational wave signal.

[0118] Optionally, in the embodiments of the present disclosure, the server compares the accretion mass ratio coefficient with the preset value (such as 0) to determine whether an accretion disk is formed. When the accretion mass ratio coefficient is less than or equal to the preset value (i.e., no accretion disk is formed), the server obtains the double white dwarf gravitational wave signal based on the second mass parameter and the second orbital parameter.

[0119] In the embodiments of the present disclosure, by obtaining the gravitational wave signal of a binary white dwarf based on the second mass parameter and the second orbital parameter when the accretion material proportion coefficient is less than or equal to a preset value, the gravitational wave signal can be accurately obtained, improving the accuracy and reliability of the gravitational wave signal.

[0120] In some alternative embodiments, obtaining the gravitational wave signal of a binary white dwarf based on the second mass parameter and the second orbital parameter includes:

[0121] Based on the second mass parameter and the second orbital parameter, obtaining the mass quadrupole moment of the binary white dwarf;

[0122] Based on the mass quadrupole moment of the binary white dwarf, obtaining the gravitational wave amplitude and the gravitational wave frequency of the binary white dwarf;

[0123] Based on the gravitational wave amplitude and the gravitational wave frequency of the binary white dwarf, obtaining the gravitational wave power of the binary white dwarf;

[0124] Performing frequency domain analysis on the gravitational wave amplitude, the gravitational wave frequency, and the gravitational wave power of the binary white dwarf to obtain the gravitational wave signal of the binary white dwarf.

[0125] Optionally, in the embodiments of the present disclosure, the mass quadrupole moment of the binary white dwarf is a physical quantity that describes the degree of deviation of the mass distribution of the binary white dwarf system from spherical symmetry. According to the theory of general relativity, the mass distribution and orbital motion state of the binary white dwarf will affect the mass quadrupole moment of the binary white dwarf. The greater the mass of the binary white dwarf and the more complex the change in the orbital spacing, the more significant the change characteristics of the mass quadrupole moment of the binary white dwarf. The server substitutes the second mass parameter and the second orbital parameter into the relevant formula to deduce the mass quadrupole moment of the binary white dwarf.

[0126] During the mutual revolution of the binary white dwarfs, the change in the mass quadrupole moment of the binary white dwarf is the root cause of the generation of gravitational waves. The speed of change determines the frequency of the gravitational waves, and the magnitude of the change corresponds to the amplitude of the gravitational waves. The server substitutes the mass quadrupole moment of the binary white dwarf into the relevant formula to deduce the gravitational wave amplitude and the gravitational wave frequency of the binary white dwarf.

[0127] The gravitational wave power represents the energy size radiated by the gravitational waves per unit time. The gravitational wave amplitude and frequency are important bases for calculating the gravitational wave power. Generally speaking, the greater the gravitational wave amplitude and the higher the frequency, the greater the gravitational wave power, and the stronger the gravitational wave energy radiated by the binary white dwarf system. The server substitutes the gravitational wave amplitude and the gravitational wave frequency of the binary white dwarf into the relevant formula to deduce the gravitational wave power of the binary white dwarf.

[0128] After that, the server performs a frequency-domain analysis on the gravitational wave amplitude, frequency, and power of the double white dwarf, using methods such as discrete Fourier transform to convert the time-domain signal to the frequency domain, obtaining the spectral distribution of the time-domain gravitational wave intensity in the frequency domain. Then, the time-domain gravitational wave intensity is segmented according to the time series, and the spectral distribution of each segment is calculated, thereby obtaining the evolution of the spectral distribution over time, that is, the power spectrum evolution, and obtaining the evolution of the gravitational wave amplitude, frequency, and power of the double white dwarf over time, that is, the gravitational wave signal of the double white dwarf.

[0129] In the embodiments of the present disclosure, by obtaining the gravitational wave amplitude, frequency, and power of the double white dwarf based on the second mass parameter and the second orbital parameter, and performing a frequency-domain analysis to obtain the gravitational wave signal of the double white dwarf, the gravitational wave signal of the double white dwarf can be accurately obtained, improving the accuracy and reliability of the gravitational wave signal of the double white dwarf.

[0130] In some alternative embodiments, after determining whether the termination condition is reached, the method further includes:

[0131] In the case where the termination condition is reached, output the gravitational wave signal set.

[0132] Optionally, in the embodiments of the present disclosure, the gravitational wave signal set includes the gravitational wave signal of the double white dwarf and the gravitational wave signal of the accretion disk. When the server determines that any of the termination conditions in the above embodiments is reached, it outputs the gravitational wave signal set and ends the loop.

[0133] In the embodiments of the present disclosure, by outputting the gravitational wave signal set in the case where the termination condition is reached, the gravitational wave signal of the double white dwarf can be accurately obtained, and a complete gravitational wave signal set can be obtained.

[0134] In this embodiment, a device for obtaining a gravitational wave signal is further provided. This device is used to implement the above embodiments and preferred embodiments, and those that have been described will not be repeated. As used below, the term "module" can be a combination of software and / or hardware that can achieve a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation in hardware, or a combination of software and hardware is also possible and contemplated.

[0135] This embodiment provides a device for obtaining a gravitational wave signal, as Figure 11 shown, including:

[0136] An acquisition module 1101, configured to acquire the first orbital parameter, the first mass parameter, and the first structural parameter of the double white dwarf;

[0137] A first determination module 1102, configured to determine the target radius based on the first orbital parameter and the first mass parameter;

[0138] A second determination module 1103, configured to compare a target parameter in the first structure parameter with a target radius, and determine a second mass parameter, a second structure parameter, and a second orbital parameter;

[0139] A first obtaining module 1104, configured to obtain a binary white dwarf gravitational wave signal based on the second mass parameter and the second orbital parameter, or obtain a binary white dwarf gravitational wave signal and an accretion disk gravitational wave signal based on the second mass parameter and the second orbital parameter, and store the binary white dwarf gravitational wave signal or the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal into a gravitational wave signal set;

[0140] A second obtaining module 1105, configured to determine whether a termination condition is reached. If the termination condition is not reached, start repeating from determining the target radius based on the first orbital parameter and the first mass parameter until the termination condition is reached, and obtain a gravitational wave signal set to be output, where the termination condition is jointly determined by a first condition and a second condition, the first condition is jointly determined by a current mass parameter, a current structure parameter, and a current orbital parameter, and the second condition is used to indicate a comparison result between the current time step and a preset time step, and the gravitational wave signal set contains the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal determined each time.

[0141] In an embodiment of the present disclosure, by obtaining a first orbital parameter, a first mass parameter, and a first structure parameter of a binary white dwarf; determining a target radius based on the first orbital parameter and the first mass parameter; comparing the target parameter in the first structure parameter with the target radius to determine a second mass parameter, a second structure parameter, and a second orbital parameter; obtaining a binary white dwarf gravitational wave signal or a binary white dwarf gravitational wave signal and an accretion disk gravitational wave signal based on the second mass parameter and the second orbital parameter; determining whether a termination condition is reached, and if the termination condition is not reached, starting to repeat from determining the target radius based on the first orbital parameter and the first mass parameter until the termination condition is reached, a gravitational wave signal set to be output is obtained. Since the embodiment of the present disclosure updates the orbital parameter, the mass parameter, and the structure parameter, the gravitational wave signal can be accurately obtained.

[0142] In some optional implementation manners, the second determination module 1103 includes:

[0143] A determination sub-module, configured to, when the target parameter is less than the target radius, determine the first mass parameter as the second mass parameter, determine the first structure parameter as the second structure parameter, and determine the first orbital parameter as the second orbital parameter;

[0144] The first obtaining sub-module is configured to update the first mass parameter, the first structure parameter, and the first orbital parameter to obtain a second mass parameter, a second structure parameter, and a second orbital parameter when the target parameter is greater than or equal to the target radius.

[0145] In some alternative embodiments, the first obtaining sub-module includes:

[0146] The first obtaining unit is configured to obtain an accretion matter proportion coefficient based on the first orbital parameter, the first mass parameter, the first structure parameter, and the target radius;

[0147] The second obtaining unit is configured to update the first mass parameter and the first orbital parameter based on the accretion matter proportion coefficient to obtain a second mass parameter and a second orbital parameter;

[0148] The third obtaining unit is configured to match the second mass parameter with a preset mapping relationship to update the first structure parameter to obtain a second structure parameter, where the preset mapping relationship is used to represent the mapping relationship between the mass parameter and the structure parameter.

[0149] In some alternative embodiments, the first obtaining module 1104 includes:

[0150] The obtaining sub-module is configured to obtain an accretion disk parameter based on the second mass parameter and the second orbital parameter when the accretion matter proportion coefficient is greater than a preset value, where the preset value is used to determine whether an accretion disk is formed in a double white dwarf;

[0151] The second obtaining sub-module is configured to obtain an accretion disk gravitational wave signal based on the accretion disk parameter, the second mass parameter, and the second orbital parameter;

[0152] The third obtaining sub-module is configured to obtain a double white dwarf gravitational wave signal based on the second mass parameter and the second orbital parameter.

[0153] In some alternative embodiments, the apparatus further includes:

[0154] The third obtaining module is configured to obtain a double white dwarf gravitational wave signal based on the second mass parameter and the second orbital parameter when the accretion matter proportion coefficient is less than or equal to the preset value.

[0155] In some alternative embodiments, the first obtaining module 1104 includes:

[0156] The fourth obtaining sub-module is configured to obtain a double white dwarf mass quadrupole based on the second mass parameter and the second orbital parameter;

[0157] The fifth obtaining sub-module is configured to obtain a double white dwarf gravitational wave amplitude and a double white dwarf gravitational wave frequency based on the double white dwarf mass quadrupole;

[0158] A sixth obtaining sub-module, configured to obtain the gravitational wave power of the double white dwarf based on the gravitational wave amplitude of the double white dwarf and the gravitational wave frequency of the double white dwarf;

[0159] A seventh obtaining sub-module, configured to perform frequency domain analysis on the gravitational wave amplitude of the double white dwarf, the gravitational wave frequency of the double white dwarf, and the gravitational wave power of the double white dwarf, so as to obtain the gravitational wave signal of the double white dwarf.

[0160] In some alternative embodiments, the apparatus further includes:

[0161] An output module, configured to output a set of gravitational wave signals.

[0162] The further function descriptions of the foregoing modules and units are the same as those in the corresponding foregoing embodiments, and will not be elaborated herein.

[0163] The apparatus for obtaining a gravitational wave signal in this embodiment is presented in the form of functional units. Here, the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and a memory that execute one or more software or fixed programs, and / or other devices that can provide the foregoing functions.

[0164] This embodiment of the present disclosure further provides a computer device having the foregoing Figure 11 shown apparatus for obtaining a gravitational wave signal.

[0165] Please refer to Figure 12 , Figure 12 which is a schematic structural diagram of a computer device provided by an alternative embodiment of the present disclosure. As Figure 12 shown, the computer device includes: one or more processors 10, a memory 20, and an interface for connecting each component, including a high-speed interface and a low-speed interface. Each component communicates with each other using different buses and can be installed on a common motherboard or installed in other ways as needed. The processor can process instructions executed within the computer device, including instructions stored in the memory or on the memory to display graphical information of the GUI on an external input / output device (such as a display device coupled to the interface). In some alternative embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (for example, as a server array, a set of blade servers, or a multi-processor system). Figure 12 In

[0166] The processor 10 may be a central processing unit, a network processor, or a combination thereof. Among them, the processor 10 may further include a hardware chip. The above-mentioned hardware chip may be an application-specific integrated circuit, a programmable logic device, or a combination thereof. The above-mentioned programmable logic device may be a complex programmable logic device, a field programmable gate array, a generic array logic, or any combination thereof.

[0167] Among them, the memory 20 stores instructions executable by at least one processor 10, so that the at least one processor 10 executes the method shown in the above embodiments.

[0168] The memory 20 may include a program storage area and a data storage area. Among them, the program storage area may store an operating system and application programs required for at least one function; the data storage area may store data created according to the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory, and may further include a non-transitory memory, such as at least one magnetic disk storage device, a flash memory device, or other non-transitory solid-state storage devices. In some alternative embodiments, the memory 20 may optionally include a memory remotely disposed relative to the processor 10, and these remote memories may be connected to the computer device through a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an enterprise intranet, a local area network, a mobile communication network, and combinations thereof.

[0169] The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk, or a solid-state drive; the memory 20 may further include a combination of the above types of memories.

[0170] The computer device further includes a communication interface 30 for the computer device to communicate with other devices or a communication network.

[0171] The embodiments of the present disclosure also provide a computer-readable storage medium. The method according to the embodiments of the present disclosure may be implemented in hardware, firmware, or may be implemented as computer code recorded on a storage medium, or may be implemented as computer code originally stored in a remote storage medium or a non-temporary machine-readable storage medium and downloaded through a network and to be stored in a local storage medium, so that the method described herein may be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium may be a magnetic disk, an optical disc, a read-only memory, a random access memory, a flash memory, a hard disk, or a solid-state drive, etc.; further, the storage medium may further include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component capable of storing or receiving software or computer code, and when the software or computer code is accessed and executed by the computer, the processor, or the hardware, the method shown in the above embodiments is implemented.

[0172] A part of the present disclosure can be applied as a computer program product, such as computer program instructions, which, when executed by a computer, can call or provide the methods and / or technical solutions according to the present disclosure through the operations of the computer. Those skilled in the art should understand that the forms of existence of computer program instructions in a computer-readable medium include, but are not limited to, source files, executable files, installation package files, etc. Correspondingly, the ways for a computer to execute computer program instructions include, but are not limited to: the computer directly executes the instructions, or the computer compiles the instructions and then executes the corresponding compiled program, or the computer reads and executes the instructions, or the computer reads and installs the instructions and then executes the corresponding installed program. Herein, the computer-readable medium can be any available computer-readable storage medium or communication medium accessible to the computer.

[0173] Although the embodiments of the present disclosure have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the present disclosure, and such modifications and variations all fall within the scope defined by the appended claims.

Claims

1. A method for acquiring a gravitational wave signal, characterized in that: The method comprises: Obtaining the first orbital parameter, the first mass parameter and the first structural parameter of the binary white dwarf; determining a target radius based on the first orbit parameter and the first mass parameter; Compare the target parameter in the first structural parameter with the target radius to determine a second mass parameter, a second structural parameter and a second orbital parameter; Based on the second mass parameter and the second orbital parameter, a binary white dwarf gravitational wave signal is obtained, or based on the second mass parameter and the second orbital parameter, the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal are obtained, and the binary white dwarf gravitational wave signal or the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal are stored in a gravitational wave signal set; Determine whether a termination condition is met. If the termination condition is not met, repeat the process starting from determining the target radius based on the first orbital parameter and the first mass parameter until the termination condition is met, and obtain the gravitational wave signal set to be output, wherein the termination condition is jointly determined by a first condition and a second condition, the first condition is jointly determined by a current mass parameter, a current structural parameter, and a current orbital parameter, the second condition is used to indicate a comparison result between a current time step and a preset time step, and the gravitational wave signal set includes the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal determined each time.

2. The method according to claim 1, characterized in that The step of comparing the target parameter in the first structural parameter with the target radius to determine the second mass parameter, the second structural parameter and the second orbital parameter comprises: When the target parameter is smaller than the target radius, the first mass parameter is determined as the second mass parameter, the first structure parameter is determined as the second structure parameter, and the first orbit parameter is determined as the second orbit parameter; When the target parameter is greater than or equal to the target radius, the first mass parameter, the first structure parameter, and the first orbit parameter are updated to obtain the second mass parameter, the second structure parameter, and the second orbit parameter.

3. The method according to claim 2, characterized in that The updating of the first quality parameter, the first structure parameter and the first orbit parameter to obtain the second quality parameter, the second structure parameter and the second orbit parameter includes: Obtaining an accretion matter proportion coefficient based on the first orbital parameter, the first mass parameter, the first structural parameter, and the target radius; The first mass parameter and the first orbital parameter are updated based on the accretion matter proportion coefficient to obtain the second mass parameter and the second orbital parameter; The second quality parameter is matched with a preset mapping relationship, and the first structural parameter is updated to obtain the second structural parameter, wherein the preset mapping relationship is used to characterize the mapping relationship between the quality parameter and the structural parameter.

4. The method according to claim 3, characterized in that: The obtaining of the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal based on the second mass parameter and the second orbital parameter comprises: When the accretion matter proportion coefficient is greater than a preset value, acquiring an accretion disk parameter based on the second mass parameter and the second orbit parameter, wherein the preset value is used to determine whether the binary white dwarf forms an accretion disk; Based on the accretion disk parameter, the second mass parameter and the second orbit parameter, obtaining the accretion disk gravitational wave signal; Based on the second mass parameter and the second orbital parameter, the binary white dwarf gravitational wave signal is obtained.

5. The method according to claim 4, characterized in that The method further comprises: When the accretion matter proportion coefficient is less than or equal to the preset value, the binary white dwarf gravitational wave signal is obtained based on the second mass parameter and the second orbital parameter.

6. The method according to claim 1, 4 or 5, characterized in that: The step of obtaining a binary white dwarf gravitational wave signal based on the second mass parameter and the second orbital parameter comprises: Based on the second mass parameter and the second orbital parameter, obtaining the double white dwarf mass quadrupole moment; Based on the double white dwarf mass quadrupole moment, the double white dwarf gravitational wave amplitude and the double white dwarf gravitational wave frequency are obtained; Based on the binary white dwarf gravitational wave amplitude and the binary white dwarf gravitational wave frequency, obtaining the binary white dwarf gravitational wave power; The double white dwarf gravitational wave amplitude, the double white dwarf gravitational wave frequency and the double white dwarf gravitational wave power are analyzed in the frequency domain to obtain the double white dwarf gravitational wave signal.

7. The method according to claim 1, characterized in that After determining whether the termination condition is met, the method further includes: When the termination condition is met, the gravitational wave signal set is output.

8. A device for acquiring gravitational wave signals, characterized in that: The device comprises: An acquisition module, used to acquire a first orbital parameter, a first mass parameter, and a first structural parameter of the binary white dwarf; A first determination module, configured to determine a target radius based on the first orbit parameter and the first mass parameter; A second determination module, configured to compare the target parameter in the first structural parameter with the target radius to determine a second mass parameter, a second structural parameter and a second orbital parameter; A first obtaining module is used to obtain a binary white dwarf gravitational wave signal based on the second mass parameter and the second orbital parameter, or to obtain the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal based on the second mass parameter and the second orbital parameter, and store the binary white dwarf gravitational wave signal or the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal into a gravitational wave signal set; The second obtaining module is used to determine whether the termination condition is met. If the termination condition is not met, the module repeats the process from determining the target radius based on the first orbital parameter and the first mass parameter until the termination condition is met, thereby obtaining the gravitational wave signal set to be output, wherein the termination condition is jointly determined by the first condition and the second condition, the first condition is jointly determined by the current mass parameter, the current structural parameter and the current orbital parameter, the second condition is used to indicate the comparison result between the current time step and the preset time step, and the gravitational wave signal set includes the binary white dwarf gravitational wave signal and the accretion disk gravitational wave signal determined each time.

9. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the method for acquiring gravitational wave signals according to any one of claims 1 to 7 by executing the computer instructions.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the method for acquiring gravitational wave signals according to any one of claims 1 to 7.

Images