A method and system for evaluating laser irradiation damage of a detector during movement

The method employs a solid heat transfer model and neural network to assess moving laser damage in HgCdTe detectors, addressing the limitations of static laser evaluation methods, providing efficient and cost-effective damage assessment and protective strategy optimization.

CN119920384BActive Publication Date: 2025-07-15PLA PEOPLES LIBERATION ARMY OF CHINA STRATEGIC SUPPORT FORCE AEROSPACE ENG UNIV
View PDF 1 Cites 0 Cited by

Patent Information

Application Number
CN202510387839.3
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2025-03-31
Publication Date
2025-07-15
Estimated Expiration
2045-03-31

AI Technical Summary

Technical Problem

The prior art is difficult to accurately evaluate the damage of HgCdTe detectors under mobile laser radiation, and the traditional experimental methods are costly and time-consuming, and most experiments do not consider mobile laser light source radiation.

Method used

By constructing the solid heat transfer transient model of the detector, setting boundary conditions and meshing, the solver parameters are determined using the pre-trained linear neural network model, combining the domain ordinary differential and differential algebraic equation modules, the thermal damage situation of the detector is determined, and the damage area is calculated.

Benefits of technology

It realizes efficient, accurate and low-cost evaluation of the detector's mobile laser radiation damage, optimizes the protection strategy, improves the radiation resistance of the detector, and saves experimental costs.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119920384B_ABST
    Figure CN119920384B_ABST
Patent Text Reader

Abstract

The present invention belongs to the field of optoelectronic technology and relates to a method and system for evaluating the damage of a detector caused by moving laser irradiation. The evaluation method includes: determining the parameters of the moving laser and the parameters of the detector's structural material properties; constructing a transient solid heat transfer model of the detector; setting the physical field boundary conditions of the transient solid heat transfer model and performing mesh division; setting the solver parameters, specifically: inputting the moving laser pulse width and the moving laser frequency into a pre-trained linear neural network model and using the expression of the model to calculate the solver parameters; based on the solver parameters, solving the transient solid heat transfer model in the solver, and determining the thermal damage condition of the detector by calling the domain ordinary differential and differential algebraic equation module. If it is determined that the detector has been damaged, calculate the damage area and volume of the detector to obtain the spatial distribution of the damaged area of the detector.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to the field of optoelectronic technologies, and in particular, to a method and system for evaluating the damage of a detector caused by moving laser irradiation. Background Art

[0002] In the field of optoelectronic technologies, HgCdTe detectors are extremely vulnerable to laser interference damage due to their high precision, high sensitivity, high resolution, etc. Accurately determining the damage area of HgCdTe detectors under laser irradiation is crucial for evaluating the laser damage resistance of detectors, optimizing protection strategies, and improving the overall system performance.

[0003] Traditional damage determination methods mostly rely on physical inspections after experiments, such as microscopic observation, scanning electron microscope (SEM) analysis, etc. Although these methods are intuitive, they have limitations such as high cost and long time consumption. In addition, HgCdTe detectors themselves are costly, and it is difficult and costly to conduct damage experiments using laser sources. Moreover, most experiments focus on static laser interference damage and do not consider the irradiation situation of a moving laser source.

[0004] In summary, generally, the characteristics of the damage of detectors caused by moving laser irradiation cannot be obtained through a large number of experiments. Summary of the Invention

[0005] The present invention provides a method and system for evaluating the damage of a detector caused by moving laser irradiation, which evaluate the damage of the detector caused by moving laser irradiation, and obtain the temperature field change of the detector under the moving laser by constructing a model and solving it, so as to evaluate the damage situation of the detector.

[0006] To achieve the above object, the present invention adopts the following technical solutions:

[0007] In a first aspect, the present invention provides a method for evaluating the damage of a detector caused by moving laser irradiation, and the evaluation method includes the following steps:

[0008] Determine the moving laser parameters and the detector structure material property parameters;

[0009] Construct a transient solid heat transfer model of the detector;

[0010] Set the boundary conditions of the physical field of the transient solid heat transfer model and perform mesh division;

[0011] Set the solver parameters; the solver parameters include the initial step size of the solver, the storage time step during the laser irradiation period and the laser non-irradiation period, and the maximum step size constraints during the laser irradiation period and the laser non-irradiation period; the solver parameters are determined by using a pre-trained linear neural network model. Specifically, the moving laser pulse width and the moving laser frequency are input into the pre-trained linear neural network model, and the solver parameters are obtained by calculation using the expression of the pre-trained linear neural network model.

[0012] Based on the solver parameters, solve the solid heat transfer transient model in the solver, and determine the thermal damage of the detector by calling the domain ordinary differential and differential algebraic equation module. If it is determined that the detector has been damaged, calculate the damaged area and volume of the detector to obtain the spatial distribution of the damaged area of the detector.

[0013] As a possible implementation, the pre-trained linear neural network model includes an input layer and an output layer; among them, the parameters of the input layer include the moving laser pulse width node and the moving laser frequency node; the parameters of the output layer include the storage time step node during the laser irradiation period and the laser non-irradiation period of the solver, the initial step size node, and the maximum step size constraint nodes during the laser irradiation period and the laser non-irradiation period.

[0014] As a possible implementation, the expression of the pre-trained linear neural network model is:

[0015] ;

[0016] Among them, is the weight matrix, ; is the bias vector, ; is the storage time step during the laser irradiation period of the solver; is the storage time step during the laser non-irradiation period of the solver; is the initial step size of the solver; is the maximum step size constraint during the laser irradiation period of the solver; is the maximum step size constraint during the laser non-irradiation period of the solver; is the moving laser pulse width; is the moving laser frequency.

[0017] As a possible implementation, the parameters of the input layer and the output layer are both normalized, and the normalization range is [0, 1]; the normalization expression is:

[0018] ;

[0019] Among them, is the result after normalizing the corresponding parameter values in the training dataset; is the corresponding parameter value in the training dataset; and are the minimum and maximum values of the corresponding parameter in the training dataset, respectively.

[0020] As a possible implementation, 、 and the parameters need to satisfy physical constraint conditions: ; in the case where the prediction result does not satisfy the above constraint conditions, correction processing should be performed, and the correction method is:

[0021] ; .

[0022] As a possible implementation, for mesh generation, free triangular mesh dissection is used for the upper surface of the silicon substrate and the lower surface of the indium pillar array, and the segment element size of the mesh is predefined to be finer; swept mesh dissection is performed on the remaining parts of the transient model.

[0023] As a possible implementation, the moving laser pulse width is 100 ps - 500 ns, the repetition frequency is 1K - 10K, and the moving speed is 0 - 500 mm / s.

[0024] As a possible implementation, set to use the backward difference formula to solve the time step in the solver; the initial step size of the solver, the storage time step of the solver during the laser irradiation period, the storage time step of the solver during the non - laser - irradiated period, the maximum step size constraint of the solver during the laser irradiation period, and the maximum step size constraint of the solver during the non - laser - irradiated period are all set according to the results of a pre - trained linear neural network model.

[0025] As a possible implementation, the method for determining and evaluating the thermal damage of the detector is to call the domain ordinary differential and differential - algebraic equation module to obtain the highest temperature at a certain position in the solid heat transfer transient model in the past time; if the highest temperature is greater than the melting point of the detector material, the detector is damaged; then, by calculating the area and volume of the damage, the spatial distribution of the damaged area of the detector is obtained.

[0026] In a second aspect, the present invention provides a detector moving laser irradiation damage assessment system, specifically including:

[0027] An evaluation parameter configuration unit for configuring evaluation parameters, where the evaluation parameters include moving laser parameters and detector structure material property parameters;

[0028] Detector solid heat transfer transient model construction unit, which constructs the solid heat transfer transient model of the detector and determines the boundary conditions of the heat transfer physical field of the solid heat transfer transient model;

[0029] Meshing division unit, which performs meshing division on the solid heat transfer transient model based on the heat transfer physical field;

[0030] Solver parameter configuration unit, which is used to configure solver parameters. The solver parameters include the initial step size of the solver, the laser irradiation time period, the storage time step of the laser non-irradiation time period, and the maximum step size constraint of the laser irradiation time period and the laser non-irradiation time period. The solver parameters are determined by using a pre-trained linear neural network model. Specifically, the moving laser pulse width and the moving laser frequency are input into the pre-trained linear neural network model, and the solver parameters are obtained by calculating using the expression of the pre-trained linear neural network model;

[0031] Detector damage situation determination unit, which solves the solid heat transfer transient model in the solver by using the solver parameters, and determines the thermal damage situation of the detector by calling the domain ordinary differential and differential algebraic equation module. If it is determined that the detector has been damaged, the damage area and volume of the detector are calculated to obtain the spatial distribution of the damaged area of the detector.

[0032] Compared with the prior art, the beneficial effects produced by the present invention are as follows:

[0033] 1. The present invention proposes a method for evaluating the damage of a detector by moving laser irradiation. Compared with the traditional experiment of irradiating the detector at a fixed point with a static laser, the present invention can effectively solve the problem of evaluating the damage of the detector by moving laser irradiation.

[0034] 2. The present invention proposes a method for evaluating the damage of a detector by moving laser irradiation. This evaluation method is easy to implement, has high accuracy, short time consumption and low cost.

[0035] 3. The present invention proposes a method for evaluating the damage of a detector by moving laser irradiation. This method helps to optimize the laser protection strategy of the detector to improve the anti-irradiation performance of the detector, and at the same time can also save the experimental cost. Description of the Drawings

[0036] The drawings described herein are used to provide a further understanding of the present invention and constitute a part of the present invention. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0037] Figure 1 It is a flowchart of the method for evaluating the damage of a detector by moving laser irradiation;

[0038] Figure 2Schematic diagram of the detector structure material for the embodiments of the present invention;

[0039] Figure 3 Schematic diagram of the laser pulse switching function for the embodiments of the present invention;

[0040] Figure 4 Schematic diagram of the mesh division for the embodiments of the present invention;

[0041] Figure 5 Curve graph of the damage volume change of the HgCdTe layer in the embodiments of the present invention;

[0042] Figure 6 Curve graph of the damage area change of the lower surface of the HgCdTe layer in the embodiments of the present invention;

[0043] Figure 7 Temperature field distribution of the HgCdTe detector at three typical irradiation times (t1 = 1.2 us, t2 = 1201.2 us, t3 = 2334.5 us) provided in the embodiments of the present invention;

[0044] Figure 8 Damage conditions of the lower surface of the HgCdTe layer at three typical irradiation times (t1 = 1.2 us, t2 = 1201.2 us, t3 = 2334.5 us) provided in the embodiments of the present invention. Detailed implementation manners

[0045] For the convenience of clearly describing the technical solutions of the embodiments of the present invention, in the embodiments of the present invention, terms such as "first" and "second" are used to distinguish identical or similar items with basically the same functions and roles. For example, the first threshold and the second threshold are only used to distinguish different thresholds, and do not limit their sequence. Those skilled in the art can understand that terms such as "first" and "second" do not limit the quantity and execution order, and "first", "second", etc. do not necessarily mean different.

[0046] It should be noted that in the present invention, words such as "exemplary" or "for example" are used to represent examples, illustrations or explanations. Any embodiment or design solution described as "exemplary" or "for example" in the present invention should not be construed as being more preferred or having more advantages than other embodiments or design solutions. Exactly, using words such as "exemplary" or "for example" aims to present relevant concepts in a specific manner.

[0047] In the present invention, "at least one" means one or more, and "a plurality" means two or more. "And / or" describes the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist simultaneously, and B exists alone, where A and B can be singular or plural. The character " / " generally indicates that the associated objects before and after are in an "or" relationship. The following at least one (item) or its similar expression refers to any combination of these items, including any combination of single item (s) or plural item (s). For example, at least one (item) of a, b, or c can represent: a, b, c, the combination of a and b, the combination of a and c, the combination of b and c, or the combination of a, b, and c, where a, b, and c can be single or multiple.

[0048] An embodiment of the present invention aims to provide a method and system for evaluating the moving laser irradiation damage of a detector, which can efficiently and accurately evaluate the moving laser irradiation damage of the detector by constructing a model, optimizing and solving parameters. The specific implementation is as follows:

[0049] In a first aspect, an embodiment of the present invention provides a method for evaluating the moving laser irradiation damage of a detector. Specifically, refer to Figure 1 , and this evaluation method includes the following steps:

[0050] Determine the moving laser parameters and the detector structure and material property parameters;

[0051] As a possible implementation, the moving laser pulse width is 100 ps - 500 ns, the repetition frequency is 1K - 10K, and the moving speed is 0 - 500 mm / s.

[0052] As a possible implementation, the moving laser parameters include the laser pulse width, the single pulse energy of the laser, the position of the laser irradiation center point, the laser spot size, the laser repetition frequency, and the laser moving speed.

[0053] Taking a high-repetition-rate short-pulse moving laser source with a pulse width of 200 ns, a repetition frequency of 10 KHz, and a moving speed of 226.27 mm / s irradiating a HgCdTe detector as an example for illustration.

[0054] Exemplarily, the detailed parameters of the moving laser adopted in this embodiment are shown in Table 1 below:

[0055] Table 1 Moving Laser Parameter Table

[0056]

[0057] This embodiment adopts a HgCdTe detector, and this detector has a multi-layer structure, such as Figure 2As shown in the figure, the multi-layer structure of the detector from bottom to top is a CdZnTe layer, a HgCdTe layer, an indium pillar array, and a silicon substrate in sequence; among them, the thickness of the CdZnTe layer is 50um; the thickness of the HgCdTe layer is 10um; the radius of the indium pillar is 10um and the height is 10um. The indium pillars are arranged equidistantly on the upper surface of the silicon substrate in an 11×11 array, and the spacing is 15um; the thickness of the silicon substrate is 50um.

[0058] As a possible implementation, the detector material property parameters include density, thermal conductivity, and heat capacity;

[0059] Exemplarily, the detector material property parameters are shown in Table 2 below

[0060] Table 2 Detector Material Property Parameter Table

[0061]

[0062] Construct a transient solid heat transfer model of the detector, set the boundary conditions of the physical field of the transient model, and perform mesh division;

[0063] As a possible implementation, use COMSOL Multiphysics software to construct a transient solid heat transfer model of the detector. Specifically, add solid heat transfer, domain ordinary differential and differential algebraic equation modules in COMSOL Multiphysics software, and set the boundary conditions of the physical field.

[0064] Exemplarily, the expression of the solid heat transfer module is:

[0065]

[0066] Among them, T is the instantaneous temperature; is the material density; is the material heat capacity; is the material thermal conductivity.

[0067] Exemplarily, the expression of the thermal radiation of the detector surface to the environment is:

[0068] ;

[0069] Among them, is the surface emissivity, is the Stefan constant, is the environmental temperature. Since the operating temperature of the HgCdTe detector generally remains around 77K, the environmental temperature is set to 77K. Exemplarily, the expression of the convective heat transfer flux corresponding to the convective heat transfer between the detector surface and the environment is:

[0070] ;

[0071] Among them, h is the convective heat transfer coefficient; is the ambient temperature; T is the instantaneous temperature.

[0072] Exemplarily, when a moving laser irradiates a HgCdTe detector, the laser beam directly penetrates the CdZnTe layer, and most of the energy is deposited in the HgCdTe layer. Therefore, the laser heat source on the HgCdTe layer is set as a volume heat source, and the corresponding laser energy deposition expression is:

[0073] ;

[0074] ;

[0075] ;

[0076] ;

[0077] ;

[0078] Among them, is the laser heat source, is the absorption coefficient of HgCdTe to the laser, is the reflectivity of HgCdTe to the laser, is the position of any point on the laser energy deposition surface, ( ) is the spot center position of the laser energy deposition surface, is the radius of the laser spot, is the time, is the laser in directional velocity component, is the laser in directional velocity component, ( ) is the initial position of the laser spot on the laser energy deposition surface, is the laser beam in directional function, is the average output power of a single laser pulse, is the energy of a single laser pulse, is the laser pulse width, is the laser pulse switch time function, and the specific form of this function in the embodiment is as Figure 3 shown, is the laser repetition frequency, , is a Gaussian pulse function with a peak of 1 and a standard deviation equal to the pulse width.

[0079] As a possible implementation, such as Figure 4As shown, for mesh generation, free triangular mesh division is performed on the upper surface of the silicon substrate and the lower surface of the indium pillar array, and the segment element size of the mesh is predefined to be finer; swept mesh division is performed on the remaining part of the solid heat transfer transient model.

[0080] Set the solver parameters; the solver parameters include the initial step size of the solver, the storage time step during the laser irradiation period and the laser non-irradiation period, and the maximum step size constraints during the laser irradiation period and the laser non-irradiation period; the solver parameters are determined by using a pre-trained linear neural network model. Specifically, the moving laser pulse width and the moving laser frequency are input into the pre-trained linear neural network model, and the solver parameters are obtained by calculating using the expression of the pre-trained linear neural network model.

[0081] As a possible implementation, the pre-trained linear neural network model includes an input layer and an output layer; among them, the parameters of the input layer include the moving laser pulse width node and the moving laser frequency node; the parameters of the output layer include the storage time step node of the solver during the laser irradiation period and the laser non-irradiation period, the initial step size node, and the maximum step size constraint nodes during the laser irradiation period and the laser non-irradiation period.

[0082] As a possible implementation, the expression of the pre-trained linear neural network model is:

[0083] ;

[0084] Among them, is the weight matrix, ; is the bias vector, ; is the storage time step of the solver during the laser irradiation period; is the storage time step of the solver during the laser non-irradiation period; is the initial step size of the solver; is the maximum step size constraint of the solver during the laser irradiation period; is the maximum step size constraint of the solver during the laser non-irradiation period; is the moving laser pulse width; is the moving laser frequency.

[0085] As a possible implementation, the parameters of the input layer and the output layer are both normalized, and the normalization range is [0, 1]; the normalization expression is:

[0086] ;

[0087] Among them, It is the result after normalizing the corresponding parameter values in the training dataset; It is the corresponding parameter value in the training dataset; and are respectively the minimum and maximum values of the corresponding parameter in the training dataset.

[0088] As a possible implementation, 、 and the parameters need to satisfy the physical constraint conditions: ; In the case where the prediction result does not meet the above constraint conditions, correction processing should be performed, and the correction method is:

[0089] ; .

[0090] Based on the solver parameters, solve the transient solid heat transfer model in the solver, and determine the thermal damage of the detector by calling the domain ordinary differential and differential algebraic equation module. If it is determined that the detector is damaged, calculate the damaged area and volume of the detector to obtain the spatial distribution of the damaged area of the detector.

[0091] During specific implementation, set the use of backward difference formula to solve the time step in the solver; the solver step size is set to be adaptive; the initial step size , set to ; During the laser irradiation period, the maximum step size is set to one-tenth of the laser pulse; during the period when the laser is not irradiated, the maximum step size is set to one-third of the pulse period.

[0092] Exemplarily, call the domain ordinary differential and differential algebraic equation module to obtain the highest temperature at a certain position in the transient solid heat transfer model in the past time step. Specifically, construct a second-order differential equation about the field variable u with respect to time t as shown in the following formula, where the field variable is used to store the highest temperature; g is the source term, and it is set to judge whether the current temperature is the highest temperature of the grid node.

[0093] ;

[0094] Exemplarily, judge whether the highest temperature exceeds the melting point of the detector material. If the highest temperature exceeds the melting point of the detector, it is determined that the detector is damaged, and then calculate the damaged volume of the HgCdTe layer of the detector and the damaged area of the lower surface of the HgCdTe layer respectively. Specifically, perform volume integration on the HgCdTe layer and surface integration on the lower surface of the HgCdTe layer in post-processing, and enter in the expression column of the integration:

[0095] if(u>993[K],1,0);

[0096] This expression sets 993 K (the melting point of HgCdTe) as the thermal damage threshold for the HgCdTe layer; this expression can also be set as the expression for the plotting group to analyze and evaluate the spatial distribution of the damaged area.

[0097] The corresponding change curves of the damaged volume of the HgCdTe layer and the damaged area on the lower surface of the HgCdTe layer are as Figure 5 and Figure 6 shown. As Figure 7 shown, in this embodiment, a three-dimensional plotting group is added to plot the temperature distribution maps of the HgCdTe detector at t1 = 1.2 us, t2 = 1201.2 us, and t3 = 2334.5 us respectively; as Figure 8 shown, in this embodiment, a two-dimensional plotting group is added, the lower surface of the HgCdTe layer is selected, and in the expression column, input: if(u>993[K],1,0), and the damaged area distribution maps at t1 = 1.2 us, t2 = 1201.2 us, and t3 = 2334.5 us are plotted respectively.

[0098] From Figure 5 and Figure 6 it can be seen that the damaged volume and the damaged area increase in a "step-like" manner, which well demonstrates the damage process of the detector by the pulsed laser with a certain repetition frequency during the movement; from Figure 7 the spatio-temporal distribution of the temperature field of the detector during the movement of the laser spot can be obtained, and the thermal influence of the pulsed laser on the detector during the movement at different times can be obtained; from Figure 8 it can be seen that the damaged surface area during the movement of the pulsed laser can provide a basis for the spatio-temporal distribution evolution characteristics of the damaged area size of the detector.

[0099] In this invention, COMSOL Multiphysics is used to establish a transient model of solid heat transfer of the detector to analyze the damage of the moving laser to the detector. Combining physical processes such as heat conduction, convective heat transfer with the environment, and thermal radiation during the laser irradiation process, the distribution of the temperature field at different irradiation times is obtained. Through the domain ordinary differential and differential algebraic equation modules in the software, the thermal damage area can be judged and calculated. The evaluation method for the damage of the moving laser irradiation to the detector helps to optimize the laser protection strategy of the detector to improve the anti-irradiation performance of the detector, and at the same time can also save experimental costs.

[0100] Although the present invention has been described in connection with the various embodiments, however, in the process of implementing the claimed invention, those skilled in the art can understand and achieve other variations of the disclosed embodiments by viewing the drawings, the disclosure, and the description of the drawings. In the specification, the word "comprising" does not exclude other components or steps, and "a" or "an" does not exclude a plurality. A single processor or other unit can implement several functions recited in the specification. Certain measures are recited in mutually different embodiments, but this does not mean that these measures cannot be combined to produce good results.

[0101] Although the present invention has been described in connection with specific features and their embodiments, it is obvious that various modifications and combinations can be made without departing from the spirit and scope of the present invention. Accordingly, this specification and the drawings are merely exemplary descriptions of the present invention and are considered to have covered any and all modifications, variations, combinations, or equivalents within the scope of the present invention. Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the present invention and its equivalent technologies, the present invention is also intended to include these changes and modifications.

Claims

1. A method for evaluating the laser irradiation damage of a detector during movement, characterized in that The evaluation method includes the following steps: Determine the mobile laser parameters and the detector structure material property parameters; Construct a transient solid heat transfer model of the detector; Set the boundary conditions of the physical field of the transient solid heat transfer model and perform mesh generation; Set the solver parameters; the solver parameters include the initial step size of the solver, the storage time step during the laser irradiation period and the laser non-irradiation period, and the maximum step size constraints during the laser irradiation period and the laser non-irradiation period; the solver parameters are determined by using a pre-trained linear neural network model. Specifically, the mobile laser pulse width and the mobile laser frequency are input into the pre-trained linear neural network model, and the solver parameters are obtained by calculation using the expression of the pre-trained linear neural network model; the pre-trained linear neural network model includes an input layer and an output layer; among them, the parameters of the input layer include the mobile laser pulse width node and the mobile laser frequency node; the parameters of the output layer include the storage time step node of the solver during the laser irradiation period and the laser non-irradiation period, the initial step size node, and the maximum step size constraint nodes during the laser irradiation period and the laser non-irradiation period; the expression of the pre-trained linear neural network model is: ; Among them, is the weight matrix, ; is the bias vector, ; is the storage time step of the solver during the laser irradiation period; is the storage time step of the solver during the period when the laser is not irradiated; is the initial step size of the solver; is the maximum step size constraint of the solver during the laser irradiation period; is the maximum step size constraint of the solver during the period when the laser is not irradiated; is the moving laser pulse width; is the moving laser frequency; Based on the solver parameters, solve the transient solid heat transfer model in the solver, and determine the thermal damage condition of the detector by calling the domain ordinary differential and differential algebraic equation module. If it is determined that the detector is damaged, calculate the damaged area and volume of the detector to obtain the spatial distribution of the damaged area of the detector.

2. The method for evaluating the laser irradiation damage of a detector movement according to claim 1, wherein The parameters of the input layer and the output layer are both normalized, and the normalization range is [0, 1]; the normalization expression is: ; Among them, is the result after normalizing the corresponding parameter values in the training dataset; is the corresponding parameter value in the training dataset; and are respectively the minimum and maximum values of the corresponding parameter in the training dataset.

3. The method for evaluating the laser irradiation damage of a detector during movement according to claim 1, characterized in that, , and The parameters need to satisfy the physical constraints: ; In the case where the prediction result does not meet the above constraints, correction processing should be performed, and the correction method is: 。 4. The method for evaluating the laser irradiation damage of a detector movement according to claim 1, wherein For the mesh generation, free triangular mesh division is adopted for the upper surface of the silicon substrate and the lower surface of the indium pillar array, and the segment element size of the mesh is predefined to be finer; swept mesh division is performed on the rest of the transient solid heat transfer model.

5. The method for evaluating laser irradiation damage of a detector movement according to claim 1, wherein The mobile laser pulse width is 100 ps - 500 ns, the repetition frequency is 1K - 10K, and the mobile speed is 0 - 500 mm / s.

6. The method for evaluating the laser irradiation damage of a detector movement according to claim 1, wherein Set the use of the backward difference formula to solve the time step in the solver; the initial step size of the solver, the storage time step of the solver during the laser irradiation period, the storage time step of the solver during the laser non-irradiation period, the maximum step size constraint of the solver during the laser irradiation period, and the maximum step size constraint of the solver during the laser non-irradiation period are all set according to the results of the pre-trained linear neural network model.

7. The method for evaluating the laser irradiation damage of a detector movement according to claim 1, wherein The method for determining and evaluating the thermal damage condition of the detector is to call the domain ordinary differential and differential algebraic equation module to obtain the highest temperature at a certain position in the transient solid heat transfer model in the past time; if the highest temperature is greater than the melting point of the detector material, the detector is damaged; Then, by calculating the damaged area and volume, the spatial distribution of the damaged area of the detector is obtained.

8. A detector moving laser irradiation damage assessment system, characterized in that, Including: An evaluation parameter configuration unit for configuring evaluation parameters, and the evaluation parameters include mobile laser parameters and detector structure material property parameters; A transient solid heat transfer model construction unit for the detector constructs a transient solid heat transfer model of the detector and determines the boundary conditions of the heat transfer physical field of the transient solid heat transfer model; A meshing unit meshes the transient solid heat transfer model based on the heat transfer physical field; A solver parameter configuration unit is used to configure solver parameters, where the solver parameters include the initial step size of the solver, the laser irradiation time period, the storage time step for the laser non-irradiation time period, and the maximum step size constraints for the laser irradiation time period and the laser non-irradiation time period; the solver parameters are determined by applying a pre-trained linear neural network model. Specifically, the moving laser pulse width and the moving laser frequency are input into the pre-trained linear neural network model, and the solver parameters are obtained by calculating using the expression of the pre-trained linear neural network model; the pre-trained linear neural network model includes an input layer and an output layer; among them, the parameters of the input layer include a moving laser pulse width node and a moving laser frequency node; the parameters of the output layer include a storage time step node for the solver during the laser irradiation time period and the laser non-irradiation time period, an initial step size node, and a maximum step size constraint node for the laser irradiation time period and the laser non-irradiation time period; the expression of the pre-trained linear neural network model is: ; Among them, is the weight matrix, ; is the bias vector, ; is the storage time step of the solver during the laser irradiation period; is the storage time step of the solver during the period when the laser is not irradiated; is the initial step size of the solver; is the maximum step size constraint of the solver during the laser irradiation period; is the maximum step size constraint of the solver during the period when the laser is not irradiated; is the moving laser pulse width; is the moving laser frequency; A detector damage situation determination unit applies the solver parameters to solve the transient solid heat transfer model in the solver, and determines the thermal damage situation of the detector by calling the domain ordinary differential and differential algebraic equation module. If it is determined that the detector is damaged, the damage area and volume of the detector are calculated to obtain the spatial distribution of the damaged area of the detector.

Citation Information

Patent Citations

  • Method and system for predicting damage time of laser irradiation metal material

    CN117852361A