Apparatus arrangement method and device for laser inertial confinement fusion

Abstract

The present disclosure relates to a device arrangement method and device for laser inertial confinement fusion, the method comprising obtaining an arrangement task for a target device, the target device comprising a laser launch optical assembly and / or a diagnostic device, and determining description information of a hohlraum; performing pre-processing on a structure of the hohlraum based on the description information to obtain a hohlraum model represented by a plurality of SPH particles; determining initial velocities of the SPH particles according to the description information and laser energy required for implementing a laser inertial confinement fusion process; performing SPH simulation on an explosion process of the hohlraum model based on the initial velocities of the SPH particles to obtain an explosion result of the hohlraum model; and determining an arrangement position of the target device relative to a laser inertial confinement fusion device based on the explosion result. The present disclosure can provide specific guidance for the arrangement position of the diagnostic device and the laser launch optical assembly, and reduce or avoid damage to the diagnostic device and the laser launch optical assembly caused by impact of flying debris.

Description

Technical Field

[0001] This disclosure relates to the field of laser inertial confinement fusion technology, and in particular to a method and apparatus for arranging equipment for laser inertial confinement fusion. Background Technology

[0002] Nuclear fusion is the cleanest and most efficient ultimate energy form for humankind. Currently, the main methods for achieving nuclear fusion include magnetically confined tokamak devices and inertial confinement laser fusion devices (or "laser inertial confinement fusion devices"). In laser inertial confinement fusion devices, when a laser is injected into its black cavity to drive the target pellet's fusion, it generates radiation ablation pressure and laser ablation pressure on the cavity. This causes the black cavity material to fail and fracture, resulting in extremely high-speed fragments that scatter in all directions. These fragments can impact the laser inertial confinement fusion diagnostic device (referred to as the "diagnostic device") and the laser-emitting optical components, leading to unpredictable damage.

[0003] Currently, there are two main protection schemes. The first is to adopt a strong protection scheme for the device requiring protection, but this scheme is costly and the facilities used may affect the layout and observation of the diagnostic equipment and affect the laser emission. The second is to make targeted arrangements based on the evaluation results of the preliminary test; however, this scheme also has the problem of high cost. Summary of the Invention

[0004] In view of this, this disclosure proposes a device arrangement method and apparatus for laser inertial confinement fusion, which can provide specific guidance for the arrangement of diagnostic devices and laser emitting optical components, reduce or avoid damage to diagnostic devices and laser emitting optical components from flying debris, and significantly reduce the cost of fusion implementation and diagnostic schemes.

[0005] According to an aspect of the present disclosure, a device arrangement method for laser inertial confinement fusion is provided, including: obtaining an arrangement task for a target device, the target device including a laser emission optical assembly for emitting a laser and controlling characteristics of the laser and / or a diagnostic device for detecting and analyzing a laser inertial confinement fusion process of a laser inertial confinement fusion device, and determining description information of a hohlraum in the laser inertial confinement fusion device; performing pre-processing on a structure of the hohlraum based on the description information of the hohlraum to obtain a hohlraum model represented by a plurality of smoothed particle hydrodynamics (SPH) particles; determining initial velocities of the SPH particles according to the description information of the hohlraum and laser energy required for implementing the laser inertial confinement fusion process; performing SPH simulation on an explosion process of the hohlraum model based on the initial velocities of the SPH particles to obtain an explosion result of the hohlraum model, the explosion result including particle velocities and particle positions of the SPH particles at target time points in a future time period; and determining an arrangement position of the target device relative to the laser inertial confinement fusion device based on the particle velocities and the particle positions of the SPH particles at the target time points in the future time period.

[0006] In a possible implementation, the description information of the hohlraum includes structure information and material information; and the pre-processing on the structure of the hohlraum based on the description information of the hohlraum to obtain the hohlraum model represented by the plurality of SPH particles includes: determining a computer-aided design file of the hohlraum according to the structure information; performing meshing on the hohlraum model in the computer-aided design file to obtain a hohlraum model with polyhedral meshes; and converting the hohlraum model with the polyhedral meshes into the hohlraum model represented by the plurality of SPH particles that matches the material information based on the material information.

[0007] In a possible implementation, the determination of the initial velocities of the SPH particles according to the description information of the hohlraum and the laser energy required for implementing the laser inertial confinement fusion process includes: determining laser ablation pressure and radiation ablation pressure on an inner wall material of the hohlraum according to the description information of the hohlraum and the laser energy required for implementing the laser inertial confinement fusion process; and determining the initial velocities of the SPH particles according to the description information of the hohlraum and the laser ablation pressure and the radiation ablation pressure on the inner wall material of the hohlraum.

[0008] In a possible implementation, the initial velocity of each SPH particle is determined according to the description information of the hohlraum, the laser ablation pressure and the radiation ablation pressure on the inner wall of the hohlraum, and includes: determining the laser focal region and the non-laser focal region of the hohlraum model according to the description information of the hohlraum and the laser incidence angle required for implementing the laser inertial confinement fusion process; determining the laser ablation impact velocity corresponding to the laser ablation pressure, and determining the radiation ablation impact velocity corresponding to the radiation ablation pressure; in a case where it is determined that the SPH particle belongs to the laser focal region, determining the initial velocity of the SPH particle according to the laser ablation impact velocity and the radiation ablation impact velocity; and in a case where it is determined that the SPH particle belongs to the non-laser focal region, determining the initial velocity of the SPH particle according to the radiation impact velocity.

[0009] In a possible implementation, the explosion process of the hohlraum model is simulated by SPH based on the initial velocities of the SPH particles, to obtain the explosion result of the hohlraum model, and includes: in a preset SPH simulation space, performing global search and local search of each SPH particle at each first time point based on the initial velocity and the initial position of each SPH particle, to obtain a target search result of each interaction pair of each SPH particle at each first time point, wherein each first time point is sequentially an initial time point and each target time point, and the initial position of each SPH particle is the position of the corresponding SPH particle in the SPH simulation space at the initial time point, the SPH simulation space includes a plurality of subspaces, the global search of each SPH particle is used to determine at least one target subspace matched with the SPH particle from the plurality of subspaces, the local search of each SPH particle is used to search for a sub-result in each target subspace matched with the SPH particle, each sub-result is used to indicate the particle composition and the particle parameter of each interaction pair of the SPH particle matched with the target subspace in the corresponding target subspace, and the target search result of each SPH particle includes the sub-results searched in all target subspaces matched with the SPH particle; at each first time point, the particle velocity and the particle position of each SPH particle at the next first time point are determined based on the particle composition and the particle parameter of all interaction pairs of each SPH particle at the current first time point.

[0010] In a possible implementation, in a preset SPH simulation space, global search and local search of each SPH particle at each first time point are performed based on an initial velocity and an initial position of each SPH particle, to obtain a target search result of an interaction pair of each SPH particle at each first time point, including: dividing the SPH simulation space to obtain a plurality of layers of first spaces distributed along a Z coordinate axis, each layer of first spaces including a plurality of second spaces distributed along a Y coordinate axis, and each second space including a plurality of subspaces distributed along an X coordinate axis, wherein each subspace in the SPH simulation space is uniformly distributed along the X coordinate axis, the Y coordinate axis, and the Z coordinate axis; at each first time point, performing global search on each SPH particle in each subspace in a preset subspace order to obtain a target subspace corresponding to each SPH particle; and performing the local search in the target subspace corresponding to each SPH particle at each first time point to obtain a target search result of an interaction pair corresponding to each SPH particle at each first time point.

[0011] In a possible implementation, at each of the first time points, the global search is sequentially performed on each of the SPH particles in each of the subspaces according to a preset order of the subspaces to obtain the target subspaces corresponding to the SPH particles, including: at the current first time point, starting from an initial first space where the current SPH particle is located, a first space after the initial first space is sequentially determined along a positive direction of a Z coordinate axis to determine whether the first space satisfies a first selection condition, the target first spaces are determined by combining the initial first space and each of the first spaces that satisfy the first selection condition, the first selection condition includes that a sum of a Z coordinate value of the current SPH particle at the current first time point and a smoothing length value is greater than a target Z value of the first space, and the target Z value represents a maximum boundary value of the first space in the direction of the Z coordinate axis; at the current first time point, starting from an initial second space matched by the current SPH particle in each of the target first spaces, a second space after the initial second space in the target first space is sequentially determined along a positive direction of a Y coordinate axis to determine whether the second space satisfies a second selection condition, the target second spaces are determined by combining the initial second space and each of the second spaces that satisfy the second selection condition, the second selection condition includes that a sum of a Y coordinate value of the current SPH particle at the current first time point and a smoothing length value is greater than a target Y value of the second space, and the target Y value represents a maximum boundary value of the second space in the direction of the Y coordinate axis; at the current first time point, starting from an initial subspace matched by the current SPH particle in each of the target second spaces, a subspace after the initial subspace in the target second space is sequentially determined along a positive direction of an X coordinate axis to determine whether the subspace satisfies a third selection condition, the target subspaces are determined by combining the initial subspaces and each of the subspaces that satisfy the third selection condition, and the third selection condition includes that a sum of an X coordinate value of the current SPH particle at the current first time point and a smoothing length value is greater than a target X value of the subspace, and the target X value represents a maximum boundary value of the subspace in the direction of the X coordinate axis.

[0012] In a possible implementation, at each of the first time points, the local search is performed in the target subspace corresponding to each of the SPH particles to obtain a target search result of each of the interaction pairs corresponding to each of the SPH particles at each of the first time points, including: at the current first time point, a distance between the current SPH particle and a candidate particle is determined according to particle coordinates of the current SPH particle and particle coordinates of the candidate particle, if the distance between the current SPH particle and the candidate particle is less than or equal to a sum of a smoothing length value of the current SPH particle at the current first time point and a smoothing length value of the candidate particle at the current first time point, an interaction pair including the candidate particle and the current SPH particle is determined.

[0013] In a possible implementation, the arrangement position of the target device relative to the laser inertial confinement fusion device is determined based on the particle velocity and the particle position of each of the SPH particles at each of the target time points in a future time period, including: determining the probability of occurrence of debris at each of the target positions around the target laser inertial confinement fusion device according to the particle velocity and the particle position of all the SPH particles at each of the target time points, wherein the debris is formed according to a plurality of SPH particles simultaneously present in the same region; and taking, as the arrangement position of the laser emission optical assembly relative to the laser inertial confinement fusion device, a target position that meets a preset optical assembly arrangement condition and has a probability of occurrence of debris less than a preset threshold, and / or taking, as the arrangement position of the diagnostic device relative to the laser inertial confinement fusion device, a target position that meets a preset diagnostic device arrangement condition and has a probability of occurrence of debris less than a preset threshold.

[0014] According to another aspect of the present disclosure, there is provided a device arrangement apparatus for laser inertial confinement fusion, including: an acquisition module configured to acquire an arrangement task for a target device, the target device including a laser emission optical assembly for emitting laser and controlling characteristics of the laser and / or a diagnostic device for detecting and analyzing a laser inertial confinement fusion process of a laser inertial confinement fusion device, and determine description information of a hohlraum in the laser inertial confinement fusion device; a first simulation module configured to perform pre-processing on a structure of the hohlraum based on the description information of the hohlraum, to obtain a hohlraum model represented by a plurality of smoothed particle hydrodynamics (SPH) particles; a first determination module configured to determine initial velocities of the SPH particles according to the description information of the hohlraum and laser energy required for implementing the laser inertial confinement fusion process; a second simulation module configured to perform SPH simulation on an explosion process of the hohlraum model based on the initial velocities of the SPH particles, to obtain an explosion result of the hohlraum model, the explosion result including particle velocities and particle positions of the SPH particles at each of target time points in a future time period; and a second determination module configured to determine the arrangement position of the target device relative to the laser inertial confinement fusion device based on the particle velocities and the particle positions of the SPH particles at each of the target time points in the future time period.

[0015] In a possible implementation, the description information of the hohlraum includes structure information and material information; and the modeling of the structure of the hohlraum based on the description information of the hohlraum includes: determining a computer-aided design file of the hohlraum according to the structure information; performing meshing on the hohlraum model in the computer-aided design file to obtain a hohlraum model with polyhedral meshes; and converting the hohlraum model with polyhedral meshes into a hohlraum model represented by a plurality of SPH particles according to the material information.

[0016] In a possible implementation, the initial velocity of each SPH particle is determined according to the description information of the hohlraum and laser energy required for implementing the laser inertial confinement fusion process, and includes: determining laser ablation pressure and radiation ablation pressure on the inner wall material of the hohlraum according to the description information of the hohlraum and the laser energy required for implementing the laser inertial confinement fusion process; and determining the initial velocity of each SPH particle according to the description information of the hohlraum, the laser ablation pressure and the radiation ablation pressure on the inner wall material of the hohlraum.

[0017] In a possible implementation, the initial velocity of each SPH particle is determined according to the description information of the hohlraum, the laser ablation pressure and the radiation ablation pressure on the inner wall material of the hohlraum, and includes: determining a laser focal spot region and a non-laser focal spot region of the hohlraum model according to the description information of the hohlraum and a laser incident angle required for implementing the laser inertial confinement fusion process; determining a laser ablation impact velocity corresponding to the laser ablation pressure, and determining a radiation ablation impact velocity corresponding to the radiation ablation pressure; determining the initial velocity of a SPH particle according to the laser ablation impact velocity and the radiation ablation impact velocity in a case where the SPH particle belongs to the laser focal spot region; and determining the initial velocity of a SPH particle according to the radiation impact velocity in a case where the SPH particle belongs to the non-laser focal spot region.

[0018] In a possible implementation, the SPH simulation is performed on the explosion process of the black cavity model based on initial velocities of the SPH particles to obtain an explosion result of the black cavity model, including: performing global search and local search of each SPH particle at each first time point in a preset SPH simulation space based on initial velocities and initial positions of the SPH particles, to obtain target search results of interaction pairs of each SPH particle at each first time point, wherein each first time point is sequentially an initial time point and each target time point, and the initial position of each SPH particle is a position of the corresponding SPH particle in the SPH simulation space at the initial time point, the SPH simulation space includes a plurality of subspaces, the global search of each SPH particle is used to determine at least one target subspace matched with the SPH particle from the plurality of subspaces, and the local search of each SPH particle is used to search for a sub-result in each target subspace matched with the SPH particle, each sub-result is used to indicate particle composition and particle parameters of each interaction pair related to the SPH particle matched with the target subspace in the corresponding target subspace, and the target search result of each SPH particle includes the sub-results searched in all target subspaces matched with the SPH particle; at each first time point, particle velocities and particle positions of each SPH particle at a next first time point are determined based on particle composition and particle parameters of all interaction pairs of each SPH particle at the current first time point.

[0019] In a possible implementation, the global search and the local search of each SPH particle at each first time point in a preset SPH simulation space based on initial velocities and initial positions of the SPH particles are performed to obtain target search results of interaction pairs of each SPH particle at each first time point, including: dividing the SPH simulation space to obtain a plurality of layers of first spaces distributed along a Z coordinate axis, each layer of first spaces including a plurality of second spaces distributed along a Y coordinate axis, and each second space including a plurality of subspaces distributed along an X coordinate axis, and each subspace in the SPH simulation space being uniformly distributed along the X coordinate axis, the Y coordinate axis and the Z coordinate axis; at each first time point, performing the global search on each SPH particle in each subspace in a preset subspace order to obtain a target subspace corresponding to each SPH particle; and performing the local search in the target subspace corresponding to each SPH particle at each first time point to obtain the target search result of the interaction pair corresponding to each SPH particle at each first time point.

[0020] In a possible implementation, at each of the first time points, the global search is sequentially performed on each of the SPH particles in each of the subspaces according to a preset order of the subspaces to obtain the target subspaces corresponding to the SPH particles, including: at the current first time point, starting from an initial first space where the current SPH particle is located, a positive direction along a Z coordinate axis is sequentially determined whether each first space after the initial first space satisfies a first selected condition, the initial first space and each first space satisfying the first selected condition are determined as target first spaces, the first selected condition includes that a sum of a Z coordinate value of the current SPH particle at the current first time point and a smoothing length value is greater than a target Z value of the first space, and the target Z value represents a maximum boundary value of the first space in a direction of the Z coordinate axis; at the current first time point, starting from an initial second space matched by the current SPH particle in each target first space, a positive direction along a Y coordinate axis is sequentially determined whether each second space after the initial second space satisfies a second selected condition, the initial second space and each second space satisfying the second selected condition are determined as target second spaces, the second selected condition includes that a sum of a Y coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target Y value of the second space, and the target Y value represents a maximum boundary value of the second space in a direction of the Y coordinate axis; at the current first time point, starting from an initial subspace matched by the current SPH particle in each target second space, a positive direction along an X coordinate axis is sequentially determined whether each subspace after the initial subspace satisfies a third selected condition, the initial subspace and each subspace satisfying the third selected condition are determined as target subspaces, the third selected condition includes that a sum of an X coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target X value of the subspace, and the target X value represents a maximum boundary value of the subspace in a direction of the X coordinate axis.

[0021] In a possible implementation, at each of the first time points, the local search is performed in the target subspace corresponding to each of the SPH particles to obtain a target search result of each interaction pair corresponding to each of the SPH particles at each of the first time points, including: at the current first time point, a distance between the current SPH particle and a candidate particle is determined according to particle coordinates of the current SPH particle and particle coordinates of the candidate particle, if the distance between the current SPH particle and the candidate particle is less than or equal to a sum of a smoothing length value of the current SPH particle at the current first time point and a smoothing length value of the candidate particle at the current first time point, an interaction pair including the candidate particle and the current SPH particle is determined; and the candidate particle is any one of the SPH particles in the current target subspace.

[0022] In a possible implementation, the arrangement position of the target device relative to the laser inertial confinement fusion device is determined based on particle velocities and particle positions of the SPH particles at target time points in a future time period, including: determining probabilities of occurrence of debris at target positions around the target laser inertial confinement fusion device according to particle velocities and particle positions of all SPH particles at the target time points, wherein the debris is formed according to multiple SPH particles simultaneously present in the same region; taking, as the arrangement position of the laser emission optical assembly relative to the laser inertial confinement fusion device, a target position that satisfies a preset optical assembly arrangement condition and has a probability of occurrence of debris less than a preset threshold, and / or taking, as the arrangement position of the diagnostic device relative to the laser inertial confinement fusion device, a target position that satisfies a preset diagnostic device arrangement condition and has a probability of occurrence of debris less than a preset threshold.

[0023] According to another aspect of the present disclosure, a diagnostic method for a laser inertial confinement fusion device is provided, including: arranging a diagnostic device for detecting analysis of a laser inertial confinement fusion process of a laser inertial confinement fusion device at an arrangement position of the diagnostic device relative to the laser inertial confinement fusion device determined by the above device arrangement method; and diagnosing the laser inertial confinement fusion process by using the diagnostic device to obtain a diagnostic result, the diagnostic result including a detected value of a physical parameter involved in a physical phenomenon occurring in the laser inertial confinement fusion process.

[0024] According to another aspect of the present disclosure, an electronic device is provided, including: a processor; a memory for storing processor-executable instructions; wherein the processor is configured to implement the above method when executing the instructions stored in the memory.

[0025] According to another aspect of the present disclosure, a non-volatile computer-readable storage medium having computer program instructions stored thereon is provided, wherein the computer program instructions are executed by a processor to implement the above method.

[0026] According to another aspect of the present disclosure, a computer program product is provided, including computer-readable code or a non-volatile computer-readable storage medium carrying computer-readable code, when the computer-readable code is run in a processor of an electronic device, the processor in the electronic device executes the above method.

[0027] The device arrangement method and device for laser inertial confinement fusion provided by the present disclosure, by obtaining an arrangement task for a target device, the target device including a laser emission optical assembly for emitting laser and controlling characteristics of the laser and / or a diagnostic device for detecting and analyzing a laser inertial confinement fusion process of a laser inertial confinement fusion device, and determining description information of a hohlraum in the laser inertial confinement fusion device, modeling pre-processing the structure of the hohlraum based on the description information of the hohlraum to obtain a hohlraum model represented by a plurality of smoothed particle hydrodynamics (SPH) particles, determining an initial velocity of each SPH particle according to the description information of the hohlraum and laser energy required to implement the laser inertial confinement fusion process, performing SPH simulation on an explosion process of the hohlraum model based on the initial velocity of each SPH particle to obtain an explosion result of the hohlraum model, the explosion result including particle velocities and particle positions of each SPH particle at each target time point in a future time period, and determining an arrangement position of the target device relative to the laser inertial confinement fusion device based on the particle velocities and particle positions of each SPH particle at each target time point in the future time period, thus being able to provide specific guidance for the arrangement position of the diagnostic device and / or the laser emission optical assembly, reduce or avoid damage to the diagnostic device and / or the laser emission optical assembly caused by impact of scattered debris, and significantly reduce the cost of fusion implementation and diagnostic scheme.

[0028] Other features and aspects of the present disclosure will become apparent from the following detailed description of exemplary embodiments with reference to the drawings. BRIEF DESCRIPTION OF DRAWINGS

[0029] The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate exemplary embodiments, features, and aspects of the present disclosure and serve to explain the principles of the present disclosure.

[0030] Figure 1 A flowchart of the device arrangement method for laser inertial confinement fusion provided by the embodiment of the present disclosure is shown.

[0031] Figure 2 A schematic diagram of the hohlraum of the laser inertial confinement fusion device injected by the laser beam provided by the embodiment of the present disclosure is shown.

[0032] Figure 3 A schematic diagram of the determination process of the impact velocity in the arrangement method provided by the embodiment of the present disclosure is shown.

[0033] Figures 4a to 4d A schematic diagram of the SPH simulation process provided by the embodiment of the present disclosure is shown.

[0034] Figure 5 A distribution variation nephogram of the scattered debris at the top of the column cavity of the hohlraum provided by the embodiment of the present disclosure is shown.

[0035] Figure 6A schematic diagram showing an arrangement of a diagnostic device provided by an embodiment of the present disclosure.

[0036] Figure 7 A block diagram showing an equipment arrangement device for laser inertial confinement fusion provided by an embodiment of the present disclosure. DETAILED DESCRIPTION

[0037] Various exemplary embodiments, features, and aspects of the present disclosure will be explained in detail below with reference to the accompanying drawings. The same reference numbers in different drawings denote the same or similar elements / function. Although various aspects of the embodiments are illustrated in the drawings, the drawings are not necessarily drawn to scale unless specifically indicated.

[0038] As used herein, the terms "comprise", "comprising", "have", "having", "include", "including", "contain", "containing", or variants thereof are open-ended, and include one or more stated features, integers, elements, steps, components, or functions but do not preclude the presence or addition of one or more other features, integers, elements, steps, components, functions, or groups thereof.

[0039] When an element is referred to as being "connected", "coupled", "responsive", or variants thereof to another element, it can be directly connected, coupled, or responsive to the other element, or intervening elements can be present.

[0040] Although the terms first, second, third, etc. can be used herein to describe various elements / operations, these elements / operations should not be limited by these terms. These terms are only used to distinguish one element / operation from another element / operation. Thus, a first element / operation in some embodiments could be termed a second element / operation in other embodiments without departing from the teachings of the present inventive concept.

[0041] The word "exemplary" is used herein to mean "serving as an example, instance, or illustration." Any implementation described herein as "exemplary" is not necessarily to be construed as preferred or advantageous over other implementations.

[0042] In addition, for the purpose of convenience and brevity, detailed descriptions of well-known functions and structures incorporated herein can not be described in detail. It should be apparent to those skilled in the art that the present disclosure can be practiced without such specific details.

[0043] In order to facilitate the understanding of the technical solutions provided by the embodiments of the present disclosure, the technical environment for implementing the technical solutions will be described first.

[0044] A laser inertial confinement fusion device is a device that uses laser energy to compress and heat a target pellet containing deuterium-tritium fuel to a high energy density state and confine it for a certain time to achieve fusion energy release. The laser inertial confinement fusion device generally includes a hohlraum, a target pellet, and the like. The target pellet generally includes deuterium-tritium fuel and a shell, and the design of the target pellet needs to ensure that it can effectively absorb energy and form high-temperature and high-pressure fusion conditions under high-energy and high-power laser irradiation. The hohlraum is a place where the laser interacts with the target pellet. Laser launch optical assemblies can be used to launch lasers and control the characteristics of the lasers for generating high-energy and high-power laser beams to provide sufficient laser energy to compress and heat the target pellet. The laser inertial confinement fusion device and the laser launch optical assembly can be independently arranged, or the laser inertial confinement fusion device can include the laser launch optical assembly, and the embodiments of the present disclosure do not limit this.

[0045] For a laser inertial confinement fusion device, the loading form of laser inertial confinement fusion is divided into direct drive and indirect drive. In the indirect drive laser inertial confinement fusion device, the target pellet is generally built-in at the center of the hohlraum. When multiple high-energy laser beams are injected into the hohlraum, the laser energy is absorbed by the high-Z material of the hohlraum wall around the target pellet and partially converted into X-ray energy, and is confined in the hohlraum. Then the X-ray is absorbed by the target pellet of the fuel to generate a radiation ablation pressure to drive implosion. High-Z material refers to material with a high atomic number, which generally includes some heavy elements such as uranium, gold, and the like.

[0046] When the laser is injected into the hohlraum to be converted into X-ray for driving the target pellet fusion, the X-ray will generate tens of megabar (Mbar) level radiation ablation pressure on the cavity, and the laser focal spot area will also superimpose the pressure of the laser direct ablation, so the pressure will be stronger than the non-focal spot area (the focal spot area is about 20 Mbar-30 Mbar, and the non-focal spot area is 10 Mbar-20 Mbar). The laser radiation ablation pressure will directly cause the failure and fracture of the cavity wall metal material, and then form high-speed debris flying in all directions. The flying debris will impact the laser launch optical assembly, the laser inertial confinement fusion diagnostic device, and then cause unpredictable damage.

[0047] The laser inertial confinement fusion diagnostic device is a device for monitoring and analyzing physical phenomena in the fusion process and parameters involved in the physical phenomena. The main function is to obtain key information in the fusion reaction, such as the temperature, density, implosion speed, fusion products, etc. of the plasma, to evaluate the success of the experiment and guide subsequent research. The laser inertial confinement fusion diagnostic device and the laser inertial confinement fusion device can be independently arranged, or the laser inertial confinement fusion device includes the laser inertial confinement fusion diagnostic device. The laser inertial confinement fusion diagnostic device usually includes an optical diagnostic module, an X-ray diagnostic module, a data processing and analysis module, etc. The optical diagnostic module is mainly used for measuring the optical signal in the process of laser and plasma interaction. The X-ray diagnostic module can be a multi-channel Kirkpatrick-Baez imaging system, which can realize high-resolution X-ray space-time imaging, can image the plasma in the laser inertial confinement fusion from different angles, and obtain its information in time, space and energy spectrum, providing key experimental data for numerical simulation, and effectively promoting the understanding of key physical problems of inertial confinement fusion. The data processing and analysis module is used for processing and analyzing the data collected by the diagnostic device.

[0048] The expensive laser inertial confinement fusion diagnostic device and laser emission optical assembly are required to be protected enough during the experiment, or the diagnostic device and the laser emission optical assembly are arranged in an area with a small probability of impact of flying debris. Although the current mainstream strong protection scheme can protect the diagnostic device and the laser emission optical assembly from damage by flying debris, the strong protection scheme has obvious disadvantages, such as a sharp increase in cost, the facilities used in the strong protection scheme may affect the arrangement and observation of the diagnostic device, and affect the emission and characteristic control of the laser emission optical assembly on the laser. In addition, the strong protection scheme is a relatively blind scheme, and the arrangement of the equipment mainly relies on experience. Another mainstream scheme is to first evaluate the distribution of flying debris through several preliminary tests, and then arrange the diagnostic device and the laser emission optical assembly according to the evaluation results of the preliminary tests. The disadvantage of this scheme is that the cost of samples, manpower, etc. required for the preliminary tests is additionally increased, and the preliminary tests consume a long period of time.

[0049] In order to solve the above technical problems, the embodiment of the present disclosure provides a device arrangement method for laser inertial confinement fusion. First, the probability distribution of the black cavity debris flying position is analyzed by using a modeling simulation method, and then the target device is arranged in a place with a low probability of impact of flying debris. In this way, the arrangement position of the diagnostic device and the laser emission optical assembly can be specifically guided, the blindness in the existing scheme is avoided, and the diagnostic device and the laser emission optical assembly are reduced or avoided from being damaged by the impact of flying debris, and the cost of fusion implementation and diagnostic scheme is significantly reduced.

[0050] The present application will be described in detail below with reference to the accompanying drawings. Figures 1 to 6The device arrangement method for laser inertial confinement fusion provided by the embodiments of the present disclosure is schematically described. As shown in Figure 1 The device arrangement method can include the following steps S101 to S105.

[0051] In step S101, an arrangement task for a target device is acquired, the target device including a laser emission optical assembly for emitting laser and controlling characteristics of the laser and / or a diagnostic device for detecting and analyzing a laser inertial confinement fusion process of a laser inertial confinement fusion device, and description information of a hohlraum in the laser inertial confinement fusion device is determined.

[0052] The laser inertial confinement fusion device, especially the hohlraum, is scattered by the laser inertial confinement fusion process. By performing the arrangement task for the target device, one or more ideal arrangement positions can be determined. Arranging the target device at the arrangement position can reduce or avoid damage of the target device caused by the impact of scattered fragments. The target device can include at least one of the laser emission optical assembly and the diagnostic device, and the specific target device can be flexibly set according to actual needs, which is not limited in the embodiments of the present disclosure.

[0053] The description information of the hohlraum can include structure information and material information. The structure information can be a cavity structure, a wall thickness, a cavity size, etc., and the material information can be a cavity material, a cavity material density, etc.

[0054] The device arrangement method uses a Smoothed Particle Hydrodynamics (SPH) modeling method to simulate and analyze the scattering trajectory of the hohlraum fragments. In the SPH simulation analysis, pre-processing and post-processing of modeling can be included. The pre-processing of modeling can involve operations such as constructing a geometric model of the hohlraum, initializing attributes such as mass of SPH particles, setting boundary conditions of the hohlraum, etc. The post-processing of modeling can not only involve operations such as analyzing attributes such as running trajectory, stress distribution, energy change of the SPH particles, but also visualizing simulation results, etc.

[0055] In step S102, pre-processing of the structure of the hohlraum based on the description information of the hohlraum is performed to obtain a hohlraum model represented by a plurality of SPH particles.

[0056] The step S102 can include: determining a computer aided design file of the hohlraum according to the structure information; meshing the hohlraum model in the computer aided design file to obtain a hohlraum model with polyhedral meshes; and converting the hohlraum model with polyhedral meshes into a hohlraum model represented by a plurality of SPH particles matched with the material information based on the material information. In this way, the modeling simulation method can reduce the risk of debris impacting the target device, reduce the number of trial tests, reduce unnecessary tests and labor costs, shorten the test cycle, and the like.

[0057] For example, for the hohlraum to be checked, a three-dimensional hohlraum model matched with the structure information of the hohlraum is drawn by a CAD drawing tool to obtain a computer aided design (CAD) file and output data in stp format. After obtaining the stp data of the hohlraum, a corresponding hexahedral mesh can be generated by a mesh software to obtain a hohlraum model in hexahedral mesh. The generated hexahedral mesh is consistent in three-dimensional size, and the size change between the meshes is as small as possible. A plurality of smooth SPH particles are generated according to the generated hexahedral mesh and the material information of the hohlraum, and the hohlraum model is represented by the SPH particles. In this way, the discrete hohlraum model represented by the smooth particle fluid dynamics particles makes the hohlraum model suitable for simulation analysis by the smooth particle fluid dynamics method.

[0058] The step S103 can include: determining the initial velocity of each SPH particle according to the description information of the hohlraum and the laser energy required to implement the laser inertial confinement fusion process.

[0059] The step S103 can include: determining the initial velocity of each SPH particle according to the description information of the hohlraum and the laser energy required to implement the laser inertial confinement fusion process.

[0060] As Figure 2As shown, the multiple laser beams are injected into the inner surface of the cavity of the hohlraum through the injection holes at both ends of the hohlraum. The energy of the laser beams is related to the laser power, waveform and duration. The determination of the laser parameters such as the laser power, waveform and duration depends on the laser inertial confinement fusion process to be achieved and the description of the hohlraum for achieving the laser inertial confinement fusion process. Simply speaking, the laser energy E can be calculated by E=Pxt, P represents the laser power, and t represents the laser duration, both of which can be determined by the description of the hohlraum and the laser inertial confinement fusion process to be achieved. After the laser energy is absorbed by the high-Z material of the inner wall of the hohlraum and partially converted into X-ray energy, the X-ray will generate a radiation ablation pressure (hereinafter referred to as "radiation pressure") of tens of megabars on the cavity of the hohlraum, and the laser focal spot area will also superimpose the pressure of laser direct ablation (i.e. laser ablation pressure, hereinafter referred to as "ablation pressure"). Therefore, in order to achieve the laser inertial confinement fusion process, the pressure on the laser focal spot area and the non-laser focal spot area of the inner wall of the hohlraum is different. In terms of the hohlraum model, the SPH particles in the laser focal spot area and the non-laser focal spot area of the hohlraum model are subjected to different pressures, which makes the initial velocities of different SPH particles different.

[0061] As shown in Figure 3 , the process of converting laser energy into ablation pressure and radiation pressure can be scaled, i.e. quantified and standardized, to ensure accurate measurement of the pressure change from laser energy to compression of fusion fuel. Generally, based on the scaling, the energy conversion efficiency (i.e. the efficiency of converting laser energy into ablation pressure and radiation pressure) can be determined, and the compression process of the target fuel affected by the laser energy can be determined according to the target design parameters (such as the material, shape and size of the target), so as to determine the ablation pressure and the radiation pressure.

[0062] The material of the inner wall of the hohlraum is generally gold. By using the rain-gony parameter of gold, a conversion relationship between pressure (i.e. laser ablation pressure, radiation ablation pressure) and impact velocity can be established, as shown in the following formula 1:

[0063] P' = p0Du, D = C0 + λu Formula 1

[0064] In formula 1, P' is the ablation pressure or the radiation pressure, u is the impact velocity, p0 is the material density of gold, p0 = 19.3 g / cm3, C0 is the propagation speed of the inner wave of the material of the inner wall of the hohlraum, C0 = 3.08 km / s, and λ is a constant, usually λ = 1.56. In this way, after the ablation pressure or the radiation pressure suffered by each SPH particle of the hohlraum model is determined, the impact velocity corresponding to the ablation pressure or the radiation pressure can be determined by the conversion relationship between pressure and impact velocity such as formula 1, so as to determine the initial velocity of the corresponding SPH particle according to the pressure condition suffered by each SPH particle.

[0065] Specifically, step S103 above, determining the initial velocity of each SPH particle based on the description information of the black cavity, the laser ablation pressure and radiation ablation pressure on the inner wall material of the black cavity, can further include: determining the laser focal spot region and non-laser focal spot region of the black cavity model based on the description information of the black cavity and the laser incident angle required to achieve the laser inertial confinement fusion process, wherein the SPH particles in the laser focal spot region are subjected to both radiation ablation pressure and laser ablation pressure, while the SPH particles in the non-laser focal spot region are only subjected to radiation ablation pressure; determining the impact velocity corresponding to the laser ablation pressure based on the laser ablation pressure on the inner wall material of the black cavity, and determining the impact velocity corresponding to the radiation ablation pressure based on the radiation ablation pressure on the inner wall material of the black cavity. The impact velocity corresponding to the ablation pressure can, for example, be established through the material's equation of state and the Hugoniot equation to establish a conversion relationship between pressure and impact velocity. Thus, given the laser ablation pressure, the impact velocity corresponding to the laser ablation pressure can be determined based on this conversion relationship. The same principle applies to radiation ablation pressure, and will not be elaborated further. The initial velocity of each SPH particle is determined according to the region type of the black cavity model in which each SPH particle is located. That is, if it is determined that the SPH particle belongs to the laser focal spot region, the initial velocity of the SPH particle is determined based on the ablation impact velocity and the radiation impact velocity; and if it is determined that the SPH particle belongs to the non-laser focal spot region, the initial velocity of the SPH particle is determined based on the radiation impact velocity. In this way, the device arrangement method determines the laser focal spot region and non-laser focal spot region based on the description information of the black cavity and the laser incident angle, which can more accurately simulate the laser inertial confinement fusion process. This precise region division makes the simulation results closer to the actual physical process. Furthermore, by processing the SPH particles in the black cavity according to the region type, the simulation process can perform more accurate calculations on the pressure received by different regions. For particles in the laser focal spot region, the effects of both laser ablation pressure and radiation ablation pressure are considered, while for particles in the non-laser focal spot region, only the effect of radiation ablation pressure is considered. This avoids unnecessary computational redundancy, improves simulation efficiency, and also improves the accuracy of simulation results.

[0066] In addition, such as Figure 3 As shown, the process of converting the laser ablation pressure and radiation ablation pressure applied to the black cavity by laser energy into impact velocity can be simplified, that is, directly finding the quantitative relationship between laser energy and impact velocity, so as to convert the laser energy into initial impact velocity to load the inner surface of the cavity.

[0067] Step S104: performing SPH simulation on the explosion process of the hohlraum model based on the initial velocities of the SPH particles to obtain an explosion result of the hohlraum model, the explosion result including particle velocities and particle positions of the SPH particles at each target time point in a future time period;

[0068] The method simulates and analyzes the explosion process of a hohlraum model generated by laser radiation hohlraum by using a smoothed particle hydrodynamics method. In general, when simulating by using the SPH method, at a current time point, global search is first performed in a preset SPH simulation space to determine each target subspace matched by each SPH particle, and then local search is performed to determine interaction pairs of each SPH particle in each target subspace, so as to determine all forces (such as pressure, viscous force, and gravity) borne by the SPH particle according to the interaction pairs of the SPH particle, and calculate particle parameters such as acceleration, velocity, and displacement of the SPH particle at a next time point according to the action force conditions.

[0069] The hohlraum model in which the SPH particles loaded with initial velocities are located is solved by using the smoothed particle hydrodynamics method, and the essence is to solve the Navier-Stokes equation in the Lagrangian framework by using the SPH method. The control equation of the SPH method can include the continuity equation shown in the following formula 2:

[0070]

[0071] In formula 2, ρ represents the density of the SPH particle, t represents time, x represents the position vector of the SPH particle, v represents the velocity vector of the SPH particle, and β represents a three-dimensional space index. The control equation of the SPH method can also include the momentum equation shown in the following formula 3:

[0072]

[0073] In formula 3, v represents the velocity vector of the SPH particle, ∈ represents a second-order stress tensor, α represents a three-dimensional space index, and the remaining parameters are described above. For brevity, they will not be described herein again. The control equation of the SPH method can also include the energy equation shown in the following formula 4:

[0074]

[0075] In formula 4, e represents the unit mass energy of the SPH particle, and the remaining parameters are described above. For brevity, they will not be described herein again. The control equation of the SPH method can also include the position vector equation shown in the following formula 5:

[0076]

[0077] In formula 5, x represents a SPH particle position vector, and the remaining parameters are as described above, and will not be repeated here for brevity. The control equation of the SPH method can further include a smoothing length equation shown in the following formula 6:

[0078]

[0079] In formula 6, h represents a smoothing length, dim represents a problem dimension, and dim=3 since it is a three-dimensional problem, and the remaining parameters are as described above, and will not be repeated here for brevity. The control equation of the SPH method can further include a material constitutive equation shown in the following formula 7:

[0080] σ αβ = -Pδ αβ +τ αβ Formula 7

[0081] In formula 7, σ represents a second-order stress tensor, δ represents a Kronecker symbol, τ represents a viscous stress tensor, P represents a pressure received by a SPH particle, and the remaining parameters are as described above, and will not be repeated here for brevity. The control equation of the SPH method can further include a Gruneisen state equation in a compression state shown in the following formula 8, and a Gruneisen state equation in an expansion state shown in the following formula 9:

[0082]

[0083] In formula 8 and formula 9, P represents a pressure received by a SPH particle, ρ0 represents an initial density of a material, C0 represents an initial sound speed of the material, C represents an intercept of a Hugoniot impact test curve, and can be considered as a sound speed (i.e., C0), μ represents an indicator of a compression or expansion state, μ>1 represents compression, μ<1 represents expansion, η represents a ratio of a density after change to an initial state, γ0 represents a Gruneisen parameter, a represents a volume correction coefficient, a is generally 0, e represents an internal energy per unit mass, and the remaining parameters are as described above, and will not be repeated here for brevity.

[0084] Step S104, which involves performing an SPH simulation of the explosion process of the black cavity model based on the initial velocities of each SPH particle to obtain the explosion result of the black cavity model, may include: in a preset SPH simulation space, performing a global and local search of each SPH particle at each first time point based on the initial velocity and initial position of each SPH particle, to obtain the target search results for the interaction pairs of each SPH particle at each first time point. Here, each first time point is sequentially the initial time point and the target time point. The initial position of each SPH particle is the position of the corresponding SPH particle in the SPH simulation space at the initial time point. The SPH simulation space includes multiple subspaces, and the global search of each SPH particle is used for... At least one target subspace matching the SPH particle is determined from multiple subspaces. The local search of each SPH particle is used to search for sub-results in each target subspace matching the SPH particle. Each sub-result is used to indicate the particle composition and particle parameters of each interaction pair related to the SPH particle matching the target subspace in the corresponding target subspace. The target search results of each SPH particle include the sub-results searched in all target subspaces matching the SPH particle. At each first time point, based on the particle composition and particle parameters of all interaction pairs of each SPH particle at the current first time point, the particle velocity and particle position of each SPH particle at the next first time point are determined.

[0085] In one example, such as Figure 4a As shown, the main process of SPH simulation analysis includes input data, initialization data, performing global and local searches at each first time point and calculating particle parameters (including but not limited to force, particle acceleration, particle velocity, and displacement) at the next first time point based on the search results, outputting data, and finally determining whether the program should terminate by calculating the physical time. If the current time is greater than the end time t... end If the SPH simulation ends, the analysis is terminated; otherwise, the simulation continues. For example... Figure 4a As shown, calculating the forces can further include calculating internal forces (such as pressure difference force and velocity difference viscosity force), calculating external forces (such as the gravity acting on SPH particles), and calculating contact forces (such as the collision forces between particles). In this example, as... Figure 4b As shown, the main program for SPH simulation analysis can include reading in data and completing data initialization, then starting the program's main loop. During each step of the main loop, the time step, particle forces, and data output will be calculated. For example... Figure 4c As shown, the data initialization program architecture mainly performs operations such as initializing SPH data, initializing the smoothness length of SPH particles, initializing the search, and allocating corresponding memory for SPH particles, action pairs, materials, and components. The main program loop is the most important part, and its main execution architecture is as follows: Figure 4dAs shown, the main process includes operations such as time integration of physical quantities and calculation of physical quantities, wherein the process of calculation of physical quantities can mainly include searching and updating of interaction pairs, density updating, calculation of internal forces (energy calculation is in the same subprogram), smoothing length updating, etc.

[0086] In the process of SPH simulation analysis, the global search for each SPH particle at each target time point can refer to determining at least one target sub-space matching the current SPH particle from a plurality of sub-spaces, which is actually to find the neighbor particles of the current SPH particle around it that have potential interaction, so that a local search is performed based on these neighbor particles, which can refer to determining which neighbor particles can form an interaction pair with the current SPH particle from the neighbor particles with potential interaction found in the global search. Whether one of the two SPH particles falls within the influence domain of the other SPH particle can be determined according to the smoothing length of the two SPH particles, and if it falls within, the two SPH particles are considered to constitute a real interaction pair, otherwise not. In this way, the accuracy of the SPH simulation result can be ensured.

[0087] In the preset SPH simulation space, the global search and the local search of each SPH particle at each first time point are performed based on the initial velocity and the initial position of each SPH particle, and the target search result of the interaction pair of each SPH particle at each first time point is obtained, which can further include: dividing the SPH simulation space to obtain a plurality of first spaces distributed along the Z coordinate axis direction, each first space including a plurality of second spaces distributed along the Y coordinate axis direction, and each second space including a plurality of sub-spaces distributed along the X coordinate axis direction, wherein each sub-space in the SPH simulation space is uniformly distributed along the X coordinate axis direction, the Y coordinate axis direction and the Z coordinate axis direction. This step can be understood as a grid processing of the SPH simulation space, and the SPH simulation space is preset to be very large to ensure that the running of each SPH particle is located in the SPH simulation space; at each first time point, the global search of each SPH particle in each sub-space is performed in turn according to a preset sub-space order, and the target sub-space corresponding to each SPH particle is obtained; at each first time point, the local search is performed in the target sub-space corresponding to each SPH particle, and the target search result of the interaction pair corresponding to each SPH particle at each first time point is obtained.

[0088] The above-mentioned local search in the target subspace corresponding to each SPH particle at each first time point to obtain the target search result of the interaction pair corresponding to each SPH particle at each first time point can include: at the current first time point, determining the distance between the current SPH particle and the candidate particle according to the particle coordinates of the current SPH particle and the particle coordinates of the candidate particle; if the distance between the current SPH particle and the candidate particle is less than or equal to the sum of the smoothing length value of the current SPH particle at the current first time point and the smoothing length value of the candidate particle at the current first time point, determining that an interaction pair including the candidate particle and the current SPH particle is formed, the particle composition of the interaction pair is the candidate particle and the current SPH particle, and the particle parameters of the interaction pair are the forces, particle accelerations, particle velocities, and particle displacements, etc. experienced by the candidate particle and the current SPH particle respectively; if the distance between the current SPH particle and the candidate particle is greater than the sum of the smoothing length value of the current SPH particle at the current first time point and the smoothing length value of the candidate particle at the current first time point, it is considered that the candidate particle and the current SPH particle do not form an interaction pair. Wherein, the candidate particle is any one SPH particle in the current target subspace, and the calculation method of the smoothing length value of the SPH particle at the current first time point can be to first calculate the density according to the particle velocity and the particle position of the corresponding SPH particle at the current first time point, and this step can specifically use the above-mentioned formula 2 to calculate the density; and then further calculate the smoothing length value of the SPH particle at the current first time point according to the density, and this step can specifically use the above-mentioned formula 6 to calculate the smoothing length value.

[0089] The global search can be performed along the directions of the X coordinate axis, the Y coordinate axis and the Z coordinate axis corresponding to the SPH simulation space, respectively. The above global search of each SPH particle in each sub-space at each first time point in the preset sub-space order can include: at the current first time point, determining whether each first space after an initial first space in which the current SPH particle is located satisfies a first selection condition along the Z coordinate axis in a positive direction, the initial first space being the initial first space, the first selection condition including that the sum of the Z coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target Z value of the first space, the target Z value representing the maximum boundary value of the first space in the Z coordinate axis direction; at the current first time point, determining whether each second space after an initial second space in which the current SPH particle is located in each target first space satisfies a second selection condition along the Y coordinate axis in a positive direction, the initial second space being the initial second space in the target first space, the second selection condition including that the sum of the Y coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target Y value of the second space, the target Y value representing the maximum boundary value of the second space in the Y coordinate axis direction; at the current first time point, determining whether each sub-space after an initial sub-space in which the current SPH particle is located in each target second space satisfies a third selection condition along the X coordinate axis in a positive direction, the initial sub-space being the initial sub-space in the target second space, the third selection condition including that the sum of the X coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target X value of the sub-space, the target X value representing the maximum boundary value of the sub-space in the X coordinate axis direction.

[0090] Exemplarily, the space of the SPH computing system in this example (i.e. the preset SPH simulation space) is divided into 3*3*3 total 27 subspaces (or lattices), i.e. there are 3 layers of first spaces (denoted as z=1, z=2, z=3 respectively) distributed along the Z coordinate axis direction, each layer of first space includes 3 second spaces (denoted as y=1, y=2, y=3 respectively) distributed along the Y coordinate axis direction, and each second space includes 3 subspaces (denoted as x=1, x=2, x=3 respectively) distributed along the X coordinate axis direction. The 27 subspaces are preset with serial numbers, and each SPH particle in each subspace also has a serial number. At the initial time point, the initial velocity of each SPH particle in each subspace can be obtained by the above step S103, and the initial position of each SPH particle in the SPH simulation space (i.e. the coordinate of the center point of the SPH particle) can also be determined in advance.

[0091] At the current first time point, the search can be performed according to the preset subspace order. For the first subspace x=1, y=1, z=1, the SPH particle with the smallest serial number in the subspace can be taken as the current SPH particle, and the initial first space (i.e. z=1) where the current SPH particle is located is determined as the target first space. The initial second space (i.e. y=1, z=1) where the current SPH particle is located is determined as the target second space. In y=1, z=1, the initial subspace (i.e. x=1, y=1, z=1) where the current SPH particle is located is determined as the target subspace and a local search is performed in the target subspace to obtain the sub-result corresponding to the target subspace, wherein the process of the local search is described in detail above and will not be repeated here. Then the migration search is performed along the positive direction of the X coordinate axis, i.e. starting from the initial subspace where the current SPH particle is matched in y=1, z=1, each subspace after the initial subspace in the target second space is determined in turn along the positive direction of the X coordinate axis to determine whether it satisfies the third selected condition. The third selected condition includes that the sum of the X coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than the target X value of the subspace (i.e. the maximum boundary value of the subspace in the X coordinate axis direction). First, it is determined whether the subspace x=2, y=1, z=1 satisfies the third selected condition. If the subspace x=2, y=1, z=1 does not satisfy the third selected condition, the migration search along the positive direction of the X coordinate axis in y=1, z=1 is stopped. If the subspace x=2, y=1, z=1 satisfies the third selected condition, the subspace is taken as the target subspace and a local search is performed in the target subspace to obtain the sub-result corresponding to the target subspace. Then it is determined whether the subspace x=3, y=1, z=1 satisfies the third selected condition. If the subspace does not satisfy the third selected condition, the migration search in y=1, z=1 is stopped; if the subspace satisfies the third selected condition, the subspace is taken as the target subspace and a local search is performed in the target subspace to obtain the sub-result corresponding to the target subspace.

[0092] After determining the target sub-space matched by the current SPH particle in y = 1, z = 1, the migration search is performed in the positive direction along the Y coordinate axis. According to the original sub-space in which the current SPH particle is located, i.e., x = 1, y = 1, z = 1, it can be determined that if the migration search is performed in the Y direction, the first second space after the initial second space (i.e., y = 1, z = 1) of the SPH particle is y = 2, z = 1. Next, it is determined whether y = 2, z = 1 satisfies the second selected condition. The second selected condition includes that the sum of the Y coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than the target Y value (i.e., the maximum boundary value of the second space in the Y coordinate axis direction) of the second space. It should be noted that the determination of whether y = 2, z = 1 satisfies the second selected condition can be performed at the beginning or in the process of migration search, and the embodiments of the present disclosure do not limit this. If y = 2, z = 1 does not satisfy the second selected condition, the migration search in the positive direction along the Y coordinate axis is stopped. If y = 2, z = 1 satisfies the second selected condition, y = 2, z = 1 is determined as the target second space, and x = 1, y = 2, z = 1 is determined as the initial sub-space matched by the current SPH particle in y = 2, z = 1. Then, the initial sub-space x = 1, y = 2, z = 1 is determined as the target sub-space, and the local search is performed in the target sub-space to obtain the sub-result corresponding to the target sub-space. Then, in y = 2, z = 1, it is determined in sequence whether each sub-space after the initial sub-space x = 1, y = 2, z = 1 in the target second space satisfies the third selected condition along the X coordinate axis in the positive direction, and the process is the same as that of the current SPH particle in y = 1, z = 1. Here, details are not repeated.

[0093] After determining the target sub-space matched by the current SPH particle in y=2, z=1, the migration search is continued along the positive direction of the Y coordinate axis, and the search process in y=3, z=1 is similar to that in y=2, z=1. For brevity, the search process in y=3, z=1 is not described herein. Thus, the migration search of the current SPH particle in z=1 is completed. Next, the migration search is continued along the positive direction of the Z coordinate axis. Starting from z=1 where the current SPH particle is located, it is determined whether each first space after z=1 along the positive direction of the Z coordinate axis satisfies a first selection condition. The first selection condition includes that the sum of the Z coordinate value and the smoothing length value of the current SPH particle at the current first time point is greater than the target Z value (i.e., the maximum boundary value of the first space in the direction of the Z coordinate axis) of the first space. Each first space satisfying the first selection condition is determined as a target first space. Specifically, it is first determined whether z=2 satisfies the first selection condition. If z=2 does not satisfy the first selection condition, the migration search along the positive direction of the Z coordinate axis is stopped. If z=2 satisfies the first selection condition, z=2 is determined as the target first space, y=1, z=2 is determined as the initial second space matched by the current SPH particle in z=2, and x=1, y=1, z=2 is determined as the initial sub-space of the current SPH particle in y=1, z=2. Then, the initial sub-space x=1, y=1, z=2 is determined as the target sub-space, and a local search is performed in the target sub-space to obtain a sub-result corresponding to the target sub-space. Next, in y=1, z=2, it is determined whether each sub-space after the initial sub-space x=1, y=1, z=2 in the target second space satisfies a third selection condition along the positive direction of the X coordinate axis. The process is similar to that in y=1, z=1, and is not described herein.

[0094] After determining the target subspaces matched by the current SPH particle in y = 1, z = 2, the migration search is continued along the positive direction of the Y coordinate axis, and the search processes in y = 2, z = 2 and y = 3, z = 2 are similar to the search process in y = 1, z = 2. For brevity, the search processes in y = 2, z = 2 and y = 3, z = 2 are not described herein. Thus, the migration search of the current SPH particle in z = 2 is completed. Next, the migration search is continued along the positive direction of the Z coordinate axis, and the search process in the first space z = 3 is similar to the search process in z = 2. For brevity, the search process in z = 3 is not described herein. In this way, the global search and the local search of the SPH particle with the smallest serial number (denoted as particle A) in the first subspacex = 1, y = 1, z = 1 are completed, and thus the corresponding sub-results in all target subspaces matched by the particle A can be searched, and all the sub-results are the target search results of the interaction pair of the particle A at the current first time point. The search processes of the remaining SPH particles in the first subspacex = 1, y = 1, z = 1 are similar to the search process of the particle A, and thus the search processes of the remaining SPH particles in the first subspacex = 1, y = 1, z = 1 are not described herein. The search processes in the remaining subspaces are similar to the search process in the first subspacex = 1, y = 1, z = 1, and thus the search processes in the remaining subspaces are not described herein.

[0095] In fact, the search (including the global search and the local search) for each SPH particle at each first time point can be performed simultaneously or in a certain order (i.e., the order in the above example along the subspaces and the serial numbers of the SPH particles in the subspaces), and the embodiments of the present disclosure do not limit the order.

[0096] The target X value, target Y value, and target Z value are described. Assume that the length values of each subspace in the X coordinate axis direction, Y coordinate axis direction, and Z coordinate axis direction are all 1. For the subspace of x=3, y=1, and z=1, the maximum boundary value of the subspace in the X coordinate axis direction is x*1=3*1=3. If it is to be determined whether the subspace satisfies the third selected condition, it is to be determined whether the sum of the X coordinate value in the center point coordinate of the current SPH particle at the current first time point and the smoothing length value of the SPH particle at the current first time point is greater than the maximum boundary value 3. For the second space of y=2 and z=1, the maximum boundary value of the second space in the Y coordinate axis direction is y*1=2*1=2. If it is to be determined whether the second space satisfies the second selected condition, it is to be determined whether the sum of the Y coordinate value in the center point coordinate of the current SPH particle at the current first time point and the smoothing length value of the SPH particle at the current first time point is greater than the maximum boundary value 2. For the first space of z=3, the maximum boundary value of the first space in the Z coordinate axis direction is z*1=3*1=3. If it is to be determined whether the first space satisfies the first selected condition, it is to be determined whether the sum of the Z coordinate value in the center point coordinate of the current SPH particle at the current first time point and the smoothing length value of the SPH particle at the current first time point is greater than the maximum boundary value 3. The target X value of the remaining subspace, the target Y value of the second space, and the target Z value of the first space are the same as the above example, and are not described again.

[0097] In summary, the migration search process for each SPH particle can include migration out and migration in processes in the X coordinate axis positive direction (hereinafter referred to as the X axis direction), Y coordinate axis positive direction (hereinafter referred to as the Y axis direction), and Z coordinate axis positive direction (hereinafter referred to as the Z axis direction) in order. After the SPH simulation space is divided into multiple grids, multiple layers of grids can be formed in the Z axis direction, each layer can form multiple rows of grids in the Y axis direction, and each row can form multiple grids in the X axis direction. The migration strategy can perform the following three steps: first, the X axis direction migration can include migrating the SPH particle to the next grid in the X axis direction to continue detecting interaction pairs if it is determined that the X axis coordinate of the SPH particle after adding the smoothing length value exceeds the boundary of the grid in the X axis direction. Then, after each row of grids in the X axis direction detects interaction pairs, the Y axis direction migration can include migrating the SPH particle to the next row in the Y axis direction to continue detecting interaction pairs if it is determined that the Y axis coordinate of the SPH particle after adding the smoothing length value exceeds the boundary of the grid in the Y axis direction. Finally, after each layer of grids detects interaction pairs, the Z axis direction migration can include migrating the SPH particle to the next layer in the Z axis direction to continue detecting interaction pairs if it is determined that the Z axis coordinate of the SPH particle after adding the smoothing length value exceeds the boundary of the grid in the Z axis direction.

[0098] Thus, for each SPH particle, the search result of the interaction pair of the SPH particle at the current first time point can be obtained through the global search and the local search of the SPH particle at the current first time point, and the search result can be used to indicate the particle composition and the particle parameters of each interaction pair, the particle composition represents two different SPH particles forming the interaction pair, and the particle parameters include but are not limited to the force, the particle acceleration, the particle velocity and the particle displacement of the SPH particle, so that the force analysis result of the SPH particle can be determined according to the particle composition and the particle parameters of all the interaction pairs of the SPH particle at the current first time point, and the particle velocity and the particle position of the SPH particle at the next first time point can be calculated according to the force analysis result.

[0099] The migration search algorithm complexity of the SPH simulation analysis provided by the device arrangement method is O(n), and the algorithm complexity of the traditional search method is O(NLogN), since O(n) < O(nlogn), the SPH simulation analysis provided by the device arrangement method can realize linear search efficiency O(N), which is higher than the efficiency of the traditional search method.

[0100] In step S105, the arrangement position of the target device relative to the laser inertial confinement fusion device is determined based on the particle velocity and the particle position of each SPH particle at each target time point in a future time period.

[0101] Step S105 can include: determining the probability of the occurrence of debris at each target position around the target laser inertial confinement fusion device according to the particle velocity and the particle position of all SPH particles at each target time point, wherein the debris is formed according to a plurality of SPH particles simultaneously appearing in the same region; taking the target position satisfying the preset optical assembly arrangement condition and having a debris occurrence probability less than a preset threshold as the arrangement position of the laser emission optical assembly relative to the laser inertial confinement fusion device, and / or taking the target position satisfying the preset diagnostic device arrangement condition and having a debris occurrence probability less than a preset threshold as the arrangement position of the diagnostic device relative to the laser inertial confinement fusion device. The optical assembly arrangement condition can refer to a first position region required to be satisfied for realizing the laser inertial confinement fusion process, and the diagnostic device arrangement condition can refer to a second position region required to be satisfied for realizing the detection and analysis of the laser inertial confinement fusion process. The first position region and the second position region can be set according to actual needs, and the embodiments of the present disclosure are not limited in this regard.

[0102] In this equipment arrangement method, the target angle can also be determined based on the velocity direction of the free expansion of SPH particles in the fragments and the motion trajectory of the SPH particles. This allows the diagnostic device to be arranged towards the laser inertial confinement fusion device at the target angle, enabling the diagnostic device to better monitor the physical phenomena and parameters involved in the fusion process, thereby more accurately analyzing the laser inertial confinement fusion process.

[0103] SPH simulation analysis can simulate the temporal changes and velocity of fragment dispersion within the black cavity. For example... Figure 5 As shown, based on the distribution and variation cloud map of the scattered debris, the probability of debris occurrence in each area surrounding the black cavity can be statistically analyzed. This allows for the identification of areas where scattered debris is more likely to occur, as well as areas where debris with greater arrival time and mass is more likely to appear. Furthermore, the optimal layout scheme for the target equipment can be determined based on the debris distribution results. For example... Figure 6 As shown, the detection windows of the diagnostic device can be placed in locations with a low probability of debris occurrence. Specifically, the statistical analysis process may include first determining the velocity of each debris at each target time point to determine the trajectory of each debris at each target time point, and then determining the spatial distribution of debris at all locations at one or more target time points of interest based on the trajectory of each debris at each target time point, thereby selecting the ideal target location. In this way, the device placement method can obtain the difference in the probability of each debris occurrence through the results of a single SPH simulation. If there is a high distribution of SPH particles, the probability of debris occurrence is high, and vice versa. By comparing the differences in the number of debris occurrences in different regions, the optimal placement scheme of the target device relative to the laser inertial confinement fusion device can be determined.

[0104] The device arrangement method for laser inertial confinement fusion proposed in this disclosure has many advantages. It first simulates debris dispersion and statistically analyzes debris characteristics using smooth particle hydrodynamics, and then arranges the target equipment based on the analysis results. For example, the arrangement scheme is more targeted, and it can determine areas with low debris dispersion risk through simulation, resulting in high safety for the arrangement of the target equipment. It saves the high costs of equipment, samples, and labor consumed by strong protection schemes and preliminary experiments. It avoids the time consumed by preliminary experiments and can accelerate the overall implementation efficiency.

[0105] The embodiment of the present disclosure further provides an apparatus arrangement device for laser inertial confinement fusion, comprising: an acquisition module configured to acquire a layout task for a target apparatus, the target apparatus comprising a laser emission optical assembly for emitting laser and controlling characteristics of the laser and / or a diagnostic device for detecting and analyzing a laser inertial confinement fusion process of a laser inertial confinement fusion device, and determine description information of a hohlraum in the laser inertial confinement fusion device; a first simulation module configured to perform pre-processing on a structure of the hohlraum based on the description information of the hohlraum, to obtain a hohlraum model represented by a plurality of smoothed particle hydrodynamics (SPH) particles; a first determination module configured to determine initial velocities of the SPH particles according to the description information of the hohlraum and laser energy required for implementing the laser inertial confinement fusion process; a second simulation module configured to perform SPH simulation on an explosion process of the hohlraum model based on the initial velocities of the SPH particles, to obtain explosion results of the hohlraum model, the explosion results comprising particle velocities and particle positions of the SPH particles at each target time point in a future time period; and a second determination module configured to determine a layout position of the target apparatus relative to the laser inertial confinement fusion device based on the particle velocities and the particle positions of the SPH particles at each target time point in the future time period.

[0106] In a possible implementation, the description information of the hohlraum comprises structure information and material information; and the pre-processing on the structure of the hohlraum based on the description information of the hohlraum to obtain the hohlraum model represented by the plurality of SPH particles comprises: determining a computer-aided design file of the hohlraum according to the structure information; performing meshing on the hohlraum model in the computer-aided design file to obtain a hohlraum model with polyhedral meshes; and converting the hohlraum model with polyhedral meshes into the hohlraum model represented by the plurality of SPH particles that matches the material information based on the material information.

[0107] In a possible implementation, the determination of the initial velocities of the SPH particles according to the description information of the hohlraum and the laser energy required for implementing the laser inertial confinement fusion process comprises: determining laser ablation pressure and radiation ablation pressure on an inner wall material of the hohlraum according to the description information of the hohlraum and the laser energy required for implementing the laser inertial confinement fusion process; and determining the initial velocities of the SPH particles according to the description information of the hohlraum and the laser ablation pressure and the radiation ablation pressure on the inner wall material of the hohlraum.

[0108] In a possible implementation, the initial velocity of each SPH particle is determined according to the description information of the hohlraum, the laser ablation pressure and the radiation ablation pressure on the inner wall of the hohlraum, and includes: determining the laser focal region and the non-laser focal region of the hohlraum model according to the description information of the hohlraum and the laser incidence angle required for implementing the laser inertial confinement fusion process; determining the laser ablation impact velocity corresponding to the laser ablation pressure, and determining the radiation ablation impact velocity corresponding to the radiation ablation pressure; in a case where it is determined that the SPH particle belongs to the laser focal region, determining the initial velocity of the SPH particle according to the laser ablation impact velocity and the radiation ablation impact velocity; and in a case where it is determined that the SPH particle belongs to the non-laser focal region, determining the initial velocity of the SPH particle according to the radiation impact velocity.

[0109] In a possible implementation, the explosion process of the hohlraum model is simulated by SPH based on the initial velocities of the SPH particles, to obtain the explosion result of the hohlraum model, and includes: performing global search and local search of each SPH particle at each first time point in a preset SPH simulation space based on the initial velocity and the initial position of each SPH particle, to obtain a target search result of each interaction pair of each SPH particle at each first time point, wherein each first time point is sequentially an initial time point and each target time point, and the initial position of each SPH particle is the position of the corresponding SPH particle in the SPH simulation space at the initial time point, the SPH simulation space includes a plurality of subspaces, the global search of each SPH particle is used to determine at least one target subspace matched with the SPH particle from the plurality of subspaces, the local search of each SPH particle is used to search for a sub-result in each target subspace matched with the SPH particle, each sub-result is used to indicate the particle composition and the particle parameter of each interaction pair of the SPH particle matched with the target subspace in the corresponding target subspace, and the target search result of each SPH particle includes the sub-results searched in all target subspaces matched with the SPH particle; at each first time point, the particle velocity and the particle position of each SPH particle at the next first time point are determined based on the particle composition and the particle parameter of all interaction pairs of each SPH particle at the current first time point.

[0110] In a possible implementation, in a preset SPH simulation space, global search and local search of each SPH particle at each first time point are performed based on an initial velocity and an initial position of each SPH particle, to obtain a target search result of an interaction pair of each SPH particle at each first time point, including: dividing the SPH simulation space to obtain a plurality of first spaces distributed along a Z coordinate axis, each first space including a plurality of second spaces distributed along a Y coordinate axis, and each second space including a plurality of subspaces distributed along an X coordinate axis, wherein each subspace in the SPH simulation space is uniformly distributed along the X coordinate axis, the Y coordinate axis, and the Z coordinate axis; at each first time point, performing global search on each SPH particle in each subspace in a preset subspace order to obtain a target subspace corresponding to each SPH particle; and performing the local search in the target subspace corresponding to each SPH particle at each first time point to obtain a target search result of an interaction pair corresponding to each SPH particle at each first time point.

[0111] In a possible implementation, at each of the first time points, the global search is sequentially performed on each of the SPH particles in each of the subspaces according to a preset order of the subspaces to obtain the target subspaces corresponding to the SPH particles, including: at the current first time point, starting from an initial first space where the current SPH particle is located, a positive direction along a Z coordinate axis is sequentially determined whether each first space after the initial first space satisfies a first selected condition, the initial first space and each first space satisfying the first selected condition are determined as target first spaces, the first selected condition includes that a sum of a Z coordinate value of the current SPH particle at the current first time point and a smoothing length value is greater than a target Z value of the first space, and the target Z value represents a maximum boundary value of the first space in a direction of the Z coordinate axis; at the current first time point, starting from an initial second space matched by the current SPH particle in each target first space, a positive direction along a Y coordinate axis is sequentially determined whether each second space after the initial second space satisfies a second selected condition, the initial second space and each second space satisfying the second selected condition are determined as target second spaces, the second selected condition includes that a sum of a Y coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target Y value of the second space, and the target Y value represents a maximum boundary value of the second space in a direction of the Y coordinate axis; at the current first time point, starting from an initial subspace matched by the current SPH particle in each target second space, a positive direction along an X coordinate axis is sequentially determined whether each subspace after the initial subspace satisfies a third selected condition, the initial subspace and each subspace satisfying the third selected condition are determined as target subspaces, the third selected condition includes that a sum of an X coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target X value of the subspace, and the target X value represents a maximum boundary value of the subspace in a direction of the X coordinate axis.

[0112] In a possible implementation, at each of the first time points, the local search is performed in the target subspace corresponding to each of the SPH particles to obtain a target search result of each interaction pair corresponding to each of the SPH particles at each of the first time points, including: at the current first time point, a distance between the current SPH particle and a candidate particle is determined according to particle coordinates of the current SPH particle and particle coordinates of the candidate particle, if the distance between the current SPH particle and the candidate particle is less than or equal to a sum of a smoothing length value of the current SPH particle at the current first time point and a smoothing length value of the candidate particle at the current first time point, an interaction pair including the candidate particle and the current SPH particle is determined; and the candidate particle is any one of the SPH particles in the current target subspace.

[0113] In a possible implementation, the arrangement position of the target device relative to the laser inertial confinement fusion device is determined based on particle velocities and particle positions of the SPH particles at target time points in a future time period, and includes: determining probabilities of occurrence of debris at target positions around the target laser inertial confinement fusion device according to particle velocities and particle positions of all the SPH particles at the target time points, wherein the debris is formed according to a plurality of SPH particles simultaneously present in the same region; taking, as the arrangement position of the laser emission optical assembly relative to the laser inertial confinement fusion device, a target position that meets a preset optical assembly arrangement condition and has a probability of occurrence of debris less than a preset threshold, and / or taking, as the arrangement position of the diagnostic device relative to the laser inertial confinement fusion device, a target position that meets a preset diagnostic device arrangement condition and has a probability of occurrence of debris less than a preset threshold.

[0114] In some embodiments, the apparatus provided by the embodiments of the present disclosure has functions or includes modules that can be used to perform the methods described in the above method embodiments, and specific implementations can refer to the descriptions of the above arrangement method embodiments. For brevity, details are not described here.

[0115] The embodiments of the present disclosure also provide a diagnostic method for a laser inertial confinement fusion device, including: arranging a diagnostic device for detecting and analyzing a laser inertial confinement fusion process of a laser inertial confinement fusion device at an arrangement position of the diagnostic device relative to the laser inertial confinement fusion device determined by the above device arrangement method; and using the diagnostic device to diagnose the laser inertial confinement fusion process to obtain a diagnostic result, the diagnostic result including a detected value of a physical parameter involved in a physical phenomenon occurring in the laser inertial confinement fusion process. In some embodiments, specific implementations of the diagnostic method provided by the embodiments of the present disclosure can refer to the descriptions of the above arrangement method embodiments. For brevity, details are not described here.

[0116] The embodiments of the present disclosure also provide a computer-readable storage medium having computer program instructions stored thereon, the computer program instructions being executed by a processor to implement the above method. The computer-readable storage medium can be a volatile or non-volatile computer-readable storage medium.

[0117] The embodiments of the present disclosure also provide an electronic device, including: a processor; a memory for storing processor-executable instructions; wherein the processor is configured to implement the above method when executing the instructions stored in the memory.

[0118] The embodiment of the present disclosure further provides a computer program product, comprising computer readable code or a nonvolatile computer readable storage medium carrying computer readable code, when the computer readable code is executed in a processor of an electronic device, the processor in the electronic device performs the above method.

[0119] Figure 7 A block diagram of an apparatus for arranging device of laser inertial confinement fusion is shown. For example, the device 1900 can be provided as a server or terminal device. Referring to Figure 7 , the device 1900 includes a processing component 1922, which further includes one or more processors, and a memory resource represented by a memory 1932, for storing instructions executable by the processing component 1922, such as an application program. The application program stored in the memory 1932 can include one or more than one module each corresponding to a set of instructions. In addition, the processing component 1922 is configured to execute the instructions to perform the above method.

[0120] The device 1900 can also include a power supply component 1926 configured to perform power management of the device 1900, a wired or wireless network interface 1950 configured to connect the device 1900 to a network, and an input output interface 1958 (I / O interface). The device 1900 can operate based on an operating system stored in the memory 1932, such as Windows Server TM , MacOS X TM , Unix TM , Linux TM , FreeBSD TM or the like.

[0121] In an exemplary embodiment, a non-volatile computer readable storage medium, such as a memory 1932 including computer program instructions, is also provided, which can be executed by the processing component 1922 of the device 1900 to complete the above method.

[0122] The present disclosure can be a system, a method, and / or a computer program product. The computer program product can include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present disclosure.

[0123] Computer readable storage media can be tangible storage media which can retain and store instructions for use by an instruction execution device. Computer readable storage media can be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of computer readable storage media include the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or raised structures in a groove having instructions recorded thereon, and any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.

[0124] Computer readable program instructions described herein can be downloaded to respective computing / processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and / or a wireless network. The network can comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and / or edge servers. A network adapter card or network interface in each computing / processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing / processing device.

[0125] Computer readable program instructions for carrying out operations of the present disclosure can be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++ or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The computer readable program instructions can execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer can be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection can be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate array (FPGA), or programmable logic array (PLA) can execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present disclosure.

[0126] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide processes for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0127] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide processes for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0128] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable data processing apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer, other programmable data processing apparatus, or other device implement the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0129] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable data processing apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer, other programmable data processing apparatus, or other device implement the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0130] Embodiments of the present disclosure have been described above, and the description is intended to be illustrative, and not restrictive, of the disclosed embodiments. Many modifications and variations of the described embodiments are possible, and all such modifications and variations are intended to be within the scope of the described embodiments. The description used herein is intended to best explain the principles of the various embodiments, the practical application, and the best mode of using the present disclosure, and to enable others skilled in the art to understand the disclosure, various embodiments, and / or various implementations detailed herein. Any terminology used herein should not be considered limiting of the disclosure, various embodiments, and / or various implementations detailed herein.

Claims

1. A method of arranging equipment for laser inertial confinement fusion, characterized by, The method comprises the following steps: obtaining a layout task for a target device, the target device comprising a laser emission optical assembly for emitting laser and controlling characteristics of the laser and / or a diagnostic device for detecting a laser inertial confinement fusion process of a laser inertial confinement fusion device, and determining description information of a hohlraum in the laser inertial confinement fusion device; modeling the structure of the hohlraum based on the description information of the hohlraum to obtain a hohlraum model represented by a plurality of smoothed particle hydrodynamics (SPH) particles; determining initial velocities of the SPH particles according to the description information of the hohlraum and laser energy required for implementing the laser inertial confinement fusion process; performing SPH simulation on an explosion process of the hohlraum model based on the initial velocities of the SPH particles to obtain explosion results of the hohlraum model, the explosion results comprising particle velocities and particle positions of the SPH particles at target time points in a future time period; determining a layout position of the target device relative to the laser inertial confinement fusion device based on the particle velocities and the particle positions of the SPH particles at the target time points in the future time period.

2. The method of claim 1, wherein, The description information of the hohlraum comprises structure information and material information; wherein the modeling the structure of the hohlraum based on the description information of the hohlraum to obtain a hohlraum model represented by a plurality of smoothed particle hydrodynamics (SPH) particles comprises: determining a computer-aided design file of the hohlraum according to the structure information; performing mesh division on the hohlraum model in the computer-aided design file to obtain a hohlraum model with polyhedral meshes; converting the hohlraum model with polyhedral meshes into a hohlraum model represented by a plurality of SPH particles that matches the material information based on the material information.

3. The method of claim 1, wherein, The determining initial velocities of the SPH particles according to the description information of the hohlraum and laser energy required for implementing the laser inertial confinement fusion process comprises: determining laser ablation pressure and radiation ablation pressure on an inner wall material of the hohlraum according to the description information of the hohlraum and the laser energy required for implementing the laser inertial confinement fusion process; determining the initial velocities of the SPH particles according to the description information of the hohlraum, the laser ablation pressure and the radiation ablation pressure on the inner wall material of the hohlraum.

4. The method of claim 3, wherein, The determining the initial velocities of the SPH particles according to the description information of the hohlraum, the laser ablation pressure and the radiation ablation pressure on the inner wall material of the hohlraum comprises: determining a laser focal spot region and a non-laser focal spot region of the hohlraum model according to the description information of the hohlraum and a laser incidence angle required for implementing the laser inertial confinement fusion process; determining a laser ablation impact velocity corresponding to the laser ablation pressure and a radiation ablation impact velocity corresponding to the radiation ablation pressure; in a case where it is determined that a SPH particle belongs to the laser focal spot region, determining the initial velocity of the SPH particle according to the laser ablation impact velocity and the radiation ablation impact velocity; and In a case where it is determined that the SPH particle belongs to the non-laser focal spot region, an initial velocity of the SPH particle is determined according to the radiation ablation impact velocity.

5. The method according to any one of claims 1 to 4, characterized in that, An explosion process of the hohlraum model is simulated based on the initial velocities of the SPH particles to obtain an explosion result of the hohlraum model, including: In the preset SPH simulation space, global search and local search of each SPH particle at each first time point are performed based on the initial velocities and initial positions of the SPH particles to obtain target search results of interaction pairs of each SPH particle at each first time point, wherein each first time point is sequentially an initial time point and each target time point, and the initial position of each SPH particle is a position of the corresponding SPH particle in the SPH simulation space at the initial time point, the SPH simulation space includes a plurality of subspaces, the global search of each SPH particle is used to determine at least one target subspace matched with the SPH particle from the plurality of subspaces, the local search of each SPH particle is used to search for a sub-result in each target subspace matched with the SPH particle, each sub-result is used to indicate particle composition and particle parameters of each interaction pair of the SPH particle matched with the target subspace in the corresponding target subspace, and the target search result of each SPH particle includes the sub-results searched in all target subspaces matched with the SPH particle; At each first time point, particle velocities and particle positions of each SPH particle at a next first time point are determined based on particle composition and particle parameters of all interaction pairs of each SPH particle at the current first time point.

6. The method of claim 5, wherein, In the preset SPH simulation space, global search and local search of each SPH particle at each first time point are performed based on the initial velocities and initial positions of the SPH particles to obtain target search results of interaction pairs of each SPH particle at each first time point, including: The SPH simulation space is divided to obtain a plurality of layers of first spaces distributed along a Z coordinate axis direction, each layer of first spaces includes a plurality of second spaces distributed along a Y coordinate axis direction, and each second space includes a plurality of subspaces distributed along an X coordinate axis direction, and each subspace in the SPH simulation space is uniformly distributed in the X coordinate axis direction, the Y coordinate axis direction and the Z coordinate axis direction; At each first time point, the global search of each SPH particle in each subspace is sequentially performed according to a preset subspace order to obtain a target subspace corresponding to each SPH particle; At each first time point, the local search is performed in the target subspace corresponding to each SPH particle to obtain a target search result of interaction pairs corresponding to each SPH particle at each first time point.

7. The method of claim 6, wherein, At each first time point, the global search of each SPH particle in each subspace is sequentially performed according to a preset subspace order to obtain a target subspace corresponding to each SPH particle, including: At the current first time point, whether each first space after an initial first space where the current SPH particle is located satisfies a first selected condition is determined along the positive direction of the Z coordinate axis in sequence from the initial first space, and each first space satisfying the first selected condition and the initial first space are determined as target first spaces, the first selected condition including that the sum of the Z coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target Z value of the first space, the target Z value representing the maximum boundary value of the first space in the direction of the Z coordinate axis; At the current first time point, whether each second space after an initial second space matched by the current SPH particle in each target first space satisfies a second selected condition is determined along the positive direction of the Y coordinate axis in sequence from the initial second space, and each second space satisfying the second selected condition and the initial second space are determined as target second spaces, the second selected condition including that the sum of the Y coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target Y value of the second space, the target Y value representing the maximum boundary value of the second space in the direction of the Y coordinate axis; At the current first time point, whether each sub-space after an initial sub-space matched by the current SPH particle in each target second space satisfies a third selected condition is determined along the positive direction of the X coordinate axis in sequence from the initial sub-space, and each sub-space satisfying the third selected condition and each initial sub-space are determined as target sub-spaces, the third selected condition including that the sum of the X coordinate value of the current SPH particle at the current first time point and the smoothing length value is greater than a target X value of the sub-space, the target X value representing the maximum boundary value of the sub-space in the direction of the X coordinate axis.

8. The method according to claim 6 or 7, characterized in that, At each first time point, the local search is performed in the target sub-space corresponding to each SPH particle to obtain a target search result of an interaction pair corresponding to each SPH particle at each first time point, including that at the current first time point, The distance between the current SPH particle and the candidate particle is determined according to the particle coordinates of the current SPH particle and the particle coordinates of the candidate particle, and if the distance between the current SPH particle and the candidate particle is less than or equal to the sum of the smoothing length value of the current SPH particle at the current first time point and the smoothing length value of the candidate particle at the current first time point, an interaction pair including the candidate particle and the current SPH particle is determined. The candidate particle is any one of the SPH particles in the current target sub-space.

9. The method according to any one of claims 1 to 4, characterized in that, The arrangement position of the target device relative to the laser inertial confinement fusion device is determined based on the particle speed and the particle position of each SPH particle at each target time point in a future time period, including that Determine a probability of debris appearing at each target position around the inertial confinement fusion device according to particle velocities and particle positions of all SPH particles at each target time point, wherein the debris is formed according to multiple SPH particles appearing in the same region at the same time; Take a target position satisfying a preset optical component arrangement condition and having a probability of debris appearing less than a preset threshold as an arrangement position of the laser emission optical component relative to the inertial confinement fusion device, and / or take a target position satisfying a preset diagnostic device arrangement condition and having a probability of debris appearing less than a preset threshold as an arrangement position of the diagnostic device relative to the inertial confinement fusion device.

10. An apparatus arrangement for laser inertial confinement fusion, characterized by Comprise: An acquisition module configured to acquire an arrangement task for a target device, and determine description information of a hohlraum in the inertial confinement fusion device, wherein the target device comprises a laser emission optical component for emitting laser and controlling characteristics of the laser, and / or a diagnostic device for detecting and analyzing a laser inertial confinement fusion process of the inertial confinement fusion device; A first simulation module configured to perform pre-processing on a structure of the hohlraum based on the description information of the hohlraum, to obtain a hohlraum model represented by multiple smoothed particle hydrodynamics (SPH) particles; A first determination module configured to determine initial velocities of the SPH particles according to the description information of the hohlraum and laser energy required for implementing the laser inertial confinement fusion process; A second simulation module configured to perform SPH simulation on an explosion process of the hohlraum model based on the initial velocities of the SPH particles, to obtain explosion results of the hohlraum model, wherein the explosion results comprise particle velocities and particle positions of the SPH particles at each target time point in a future time period; A second determination module configured to determine arrangement positions of the target device relative to the inertial confinement fusion device based on the particle velocities and the particle positions of the SPH particles at each target time point in the future time period.

Images