Home Study on the impulse mechanism of optical films formed by laser plasma shock waves
Article Open Access

Study on the impulse mechanism of optical films formed by laser plasma shock waves

  • Guixia Wang EMAIL logo and Junhong Su
Published/Copyright: March 20, 2023

Abstract

In a high-power laser system, when the surface pressure of the optical film caused by laser plasma shock wave is greater than the adhesion per unit area of the film layer, the film will produce mechanical damage, and in serious cases, the whole system may not work. Therefore, studying the formation mechanism of optical film surface pressure and impulse caused by laser plasma shock wave is the key to ensure the normal operation of the high-power laser system. In this article, by studying the relaxation process of laser plasma shock wave on the surface pressure of optical film, and using the time accumulation effect of various pressures on the surface of the optical film, the calculation model of impulse on the optical film’s surface formed by laser plasma shock waves was established, and the variation rules of the impulse I st and impulse coefficient j on the unit area of single-layer Al2O3 and HfO2 optical films with different parameters were obtained. When the incident laser wavelength λ was 1,064 nm, the energy E was 0.1 J , the pulse width t p was 10 ns , the focal length of the focusing lens f was 350 mm , the distance between the film surface and the focal plane of the focusing lens z 0 was 5 mm , and the film radius R was 5 mm , the calculation and simulation results show that the impulse I st of the two films was 1 0 4 N s order of magnitude, the impulse coefficient j of the two films was 1 0 5 N s/J , the Al2O3 film with small atomic number will obtain larger I st and j , I st and j of the two films increase with the increase of E and f , and I st and j of the two films decrease with the increase of z 0 and t p . In the total impulse transfer time ( t 0 ), I st and j both increase with the increase of R .

1 Introduction

Generally speaking, in the initial stage of laser–optical film interaction, heat action heats the film and makes it locally heated, which leads to the increase of free electrons or gasification of the film. When the concentration of free electrons and ions reaches a certain amount, optical breakdown (indicating that the optical film has undergone irreversible changes, that is, failure) occurs, and then, plasma is formed [1,2,3,4,5]. The absorption of subsequent laser energy by plasma causes its own temperature to rise sharply, and its volume expands outward after reaching a high temperature and overheating state, that is, the plasma shock wave is formed [6,7,8], and the initial propagation speed of plasma shock wave reaches the order of magnitude of 104 m/s. As a result, mechanical action is generated on the optical film surface, and the film surface obtains impulse during the interaction between the plasma shock wave and the optical film surface. The magnitude of the impulse is not only related to the action laser parameters, such as energy, pulse width, and spot size, but also closely related to the properties of the action film materials and the action environment. The calculation of the impulse obtained by the film with different incident laser parameters can not only understand the mechanism of the plasma shock wave on the film impulse transfer, and the relationship between the efficiency and the laser parameters, but also have a certain reference value for the technology of using laser as the propulsion power source.

At present, much research has been done on the impact of laser-induced plasma shock wave on the target surface. In 1966, Gregy and Thomas of the University of California studied the impulse of a ruby giant pulse laser on a surface in a vacuum, and experimentally measured the impulse obtained by targets of different materials [9]. In 1972, Pirri et al. used high-power CO2 to act on carbon targets and tungsten targets, measured the impulse obtained by the targets, and tried to explain these impulses from the theory of explosive waves [10]. In 1976, Hettche et al. measured the impulse of pulsed CO2 laser on aluminum target by the target pendulum technology, and measured the change process of target surface pressure during pulse by the optical interference technology, and found out the dependence of target surface pressure and applied laser power density [11]. In 1980, Duzy et al. studied the impulse coupling coefficient of 353 nm XeF laser irradiation on aluminum target at 0.1 Torr ambient pressure [12]. Bass et al. measured the impulse coupling coefficient of the long-pulse Nd glass laser on the aluminum target and found that the impulse coupling coefficient is related to the convergence/divergence characteristics of the laser when penetrating the material, and it is inversely proportional to the migration quality of the target. At the same time, the authors also gave the relationship between the impulse coupling coefficient and the laser power density, and the relationship between the impulse coupling coefficient and the energy density [13]. Fu and Wang from the Institute of Mechanics, Chinese Academy of Sciences used the impulse matrix method to measure the impulse caused by TEA CO2-intense laser acting on the target, and also carried out the pressure measurement of shock wave induced by intense pulsed laser in the air [14,15,16]. In 2020, Liang et al. summarized the measurement methods of impulse coupling coefficient, mainly including impact pendulum method, laser interference method, laser-combined impact pendulum method, horizontal guide rail measurement method, pressure sensor method, etc., and introduced the main numerical analysis methods and influencing factors of impulse coupling coefficient [17]. In 2020, Yu et al. studied the impulse generated by ablating aluminum targets using the micro-impulse measurement method with a torsion pendulum, and the plasma plume characteristics were researched using fast photography and optical emission spectroscopy [18]. In 2022, Xu et al. used impulse measurement, spectral diagnostics, temporal evolution images, and target ablation to investigate the dynamic behaviors and parameters of Nd:YAG nanosecond laser-induced aluminum plasma at different pressures and laser fluence [19]. These works provide a lot of experimental data and results for the study of laser–target interaction, and lay a solid experimental foundation for the theoretical model research. However, these works mainly focus on the experimental measurement of target impulse, while the research on the formation mechanism of target surface impulse is relatively few; the optical film is different from the general target, and its interaction with laser is more complex. Therefore, in this article, the single-layer Al2O3 and HfO2 films are taken as examples (the single-layer Al2O3 and HfO2 film substrates are both quartz substrates, and their optical thickness is λ / 4  nm, in which the wavelength is λ = 1 , 064 nm ); this article intends to analyze the generation of optical film surface impulse from the aspects of surface pressure generation and plasma formation during the interaction between laser and optical film surface, establish the calculation model of optical film surface impulse, and calculate the optical film surface impulse under different parameter conditions. The influence of various parameters on the surface impulse of optical films is obtained, which reveals the formation mechanism of laser plasma shock wave on the surface impulse of optical films.

2 Theoretical calculation model of optical film surface impulse

The impulse on the optical film surface formed by the laser plasma shock wave is the time cumulative effect of the shock wave on various pressures. According to the generation and formation mechanism of various pressures in the process of laser and optical film, the formation of impulse and its calculation model can be discussed. The impulse calculation of optical film surface has the following two cases.

2.1 When the breakdown threshold of the optical film surface is not reached

When the applied laser power density does not reach the breakdown threshold of the optical film surface, the impulse on the unit area of the optical film surface formed by the laser has only two parts: one is the impulse caused by the light pressure itself, whereas the second is the recoil impulse caused by the material sprayed on the optical film surface.

  1. The optical pressure

    When high-power pulsed laser light acts on the surface of the optical film, considering the optical pressure p p , there are [20]

    (1) p p = P / c ,

    where P is the incident laser power density and c is the speed of light. If the incident laser is a Q-switched laser, the pulse width is generally of the order of ns . If P = 10 9 W / cm 2 , then p p = 3.33 × 10 4 N / m 2 = 33 kPa . When the laser-focusing spot diameter is 0.8 mm , the pulse width is 10 ns , the energy is 0.1 J , and the wavelength is 1 , 064 nm , the impulse caused on the optical film is the product of light pressure, the laser-focusing spot, and time, namely [20]:

    (2) P p × π × ( 0.8 ÷ 2 × 10 3 ) × 10 × 10 9 = 1.67 × 10 10 N s .

  2. Pressure produced by the ejected material

When the high-power laser is applied to the surface of the optical film, the surface temperature of the film is heated by the laser beam, and its surface temperature rises rapidly. When the vaporization temperature is reached, a Knudsen layer with high-temperature gas and liquid phases will be formed on the surface of the film, and the thickness of the layer is about 10−7 m. Because the layer continues to absorb laser pulse energy, on the one hand, the temperature continues to rise, and, on the other hand, the substances in the layer quickly leave the film surface in the form of vapor, that is, spray substances are formed [20]. The velocity of ejected material particles is Maxwell distribution, which can be characterized by a characteristic velocity. The average velocity of vapor particles can characterize the macroscopic behavior of the ejected material, namely [20]:

(3) v ¯ = 8 R g T π M mol ,

where R g is the molar gas constant, R g = 8.31432 J / ( mol K ) , T is the temperature of steam, and M mol is the molar mass of ejected material.

Suppose that when a pulsed laser with energy E acts on the optical film surface, the ejected mass of the film material is Δ m , and its mass mobility v m is defined as the mass of ejected material that can be taken away by the laser beam with unit energy. In fact, when the ejected material is ejected outward, according to the impulse theorem, the optical film surface also has a recoil impulse, so the recoil impulse I I per unit area is [20]:

(4) I I = Δ m v ¯ = v m E 8 R g T π M mol .

If the applied laser is a Q-switched laser, the mass mobility v m is about 1–10  μg / J ( μg / J indicates how much microgram ( μg ) mass can be migrated by 1 J laser energy on the surface of the film), taking the single-layer HfO2 film as an example, if E = 0.1 J , t p = 10 ns , M mol = 210.49 g / mol , T = 3 , 031 K , then I I is about 7.56 × 10 8 N s .

Obviously, when the applied laser power density does not reach the breakdown threshold of the optical film surface, the impulse caused by the light pressure is about two orders of magnitude smaller than the recoil caused by the ejected material. At this time, the impulse of the laser on the unit area of the film surface mainly depends on the recoil caused by the ejected material.

2.2 When the breakdown threshold of the optical film surface is reached

When the applied laser power density is greater than the breakdown threshold of the film surface, the impulse formed by the plasma shock wave on the unit area of the film surface consists of three parts: one is the impulse caused by the light pressure itself; second is the recoil impulse caused by ejected materials; the third is the impulse generated by the expansion and diffusion of laser plasma shock wave.

It can be seen from the previous analysis that the impulse caused by light pressure is about 1.67 × 10 10 N s , and this impulse is far less than the impulse caused by ejected material and plasma shock wave, so it can be ignored. The mass of splashed material is also very small, and the impulse I I introduced is about 7.56 × 10 8 N s . This order of magnitude impulse is not dominant in the interaction between laser and optical film surface, and can also be ignored. Therefore, when the applied laser power density is greater than the breakdown threshold of the film surface, the impulse on the optical film surface is mainly generated by the expansion and diffusion of the plasma shock wave, which will be mainly discussed below.

When the applied laser power density is greater than the breakdown threshold of the optical film surface, the optical film surface is broken down and a plasma shock wave is generated. The plasma shock wave generates great pressure on the film, making the optical film obtain impulse. If the time distribution of the pressure generated after the plasma shock wave acts on the film surface is p s ( t ) , and I s is used to represent the impulse obtained on the unit area of the film, then:

(5) I s = 0 t 0 p s ( t ) d t ,

where t 0 is the total action time of laser plasma shock wave. Therefore, the film impulse mainly depends on the pressure p s ( t ) of the expansion and diffusion of the plasma shock wave on the film surface, which is calculated below.

2.2.1 Pressure produced by the expansion and diffusion of plasma shock wave

During the pulse laser irradiation, the plasma on the thin film surface ignites, expands rapidly outward, and rapidly develops into laser supported detonation wave (LSDW). This absorption wave continues to absorb the residual pulsed laser energy and expands away from the film surface, and produces a pressure p s 10 on the film surface. However, the pulse width of the laser pulse is limited, and the disturbance behind the LSDW wave cannot catch up with LSDW during the pulse laser action. LSDW preferentially absorbs the laser pulse energy and propagates at the supersonic speed, so the intensity of LSDW would not weaken. LSDW mainly expands in the reverse direction of laser, which can be approximately regarded as one-dimensional flow field. The motion of LSDW would change after the end of the laser pulse. First, the intensity of LSDW would decrease with the propagation process. Second, the propagation of LSDW has no laser energy absorption in the opposite direction of the laser, so there is no priority. The lateral sparsity becomes important, and the motion of LSDW should be regarded as a two-dimensional flow field.

If the laser pulse width is t p and the characteristic time of LSDW two-dimensional motion is defined as t 2D , which is equal to the time required for LSDW wavefront to expand to the distance of focused spot diameter ( D s ), then [20]:

(6) D s = 2 λ f π d 1 + π z 0 d 2 λ f 2 2 1 2 .

(7) v L = 0.958 ρ 0 1 3 P 1 3 .

(8) t 2D = D s / v L .

In Eq. (6), λ is the wavelength of the incident laser beam, d is the diameter of the incident laser beam, f is the focal length of the focusing lens, and z 0 is the distance between the surface of the film and the focal plane of the focusing lens. In Eq. (7), ρ 0 is the air density, ρ 0 = 1.295 kg/m 3 , and P is the power density of the incident laser; in Eq. (8), v L is the propagation speed of the LSDW.

Use P 0 to represent the initial power density of the incident laser. When the incident laser energy is E , there are:

(9) P 0 = 4 E π D s 2 t p .

P 1 and P 2 are used to represent the radiation loss in the plasma region and the radiation loss caused by the reverse bremsstrahlung absorption of the laser passing through LSDW, respectively. Then, the actual incident laser power density P = P 0 P 1 P 2 . If the radiation loss P 1 in the plasma region is ignored, then P = P 0 P 2 . During laser irradiation, the laser acts on the surface of the film to produce a plasma with high temperature and high density. This plasma expands rapidly outward. During the expansion, the plasma continues to absorb the incident laser, so that the laser energy reaching the film will be greatly reduced, preventing the energy coupling between the laser and the film. This effect is called laser plasma shielding. Under the influence of the shielding effect of laser plasma, the laser energy reaching the film surface will gradually decrease, the plasma temperature and density will decrease, the shielding effect of laser plasma will become weak, the plasma absorbs less and less incident laser energy, the laser energy reaching the film will slowly increase, so the plasma temperature and density will rapidly increase, and the plasma shielding effect will be generated. This process is called reverse bremsstrahlung absorption, and the reverse bremsstrahlung absorption coefficient is expressed as [20]:

(10) K s = 9.4675 × 10 42 λ 2 N e 2 T 3 2 Z 2 ,

where N e represents the electron number density of the film surface material, T represents the plasma temperature, and z represents the atomic number. The actual incident laser power density after considering the reverse bremsstrahlung absorption is P [20], then

(11) P = P 0 P 2 = P 0 10 0.06 k s .

After the end of the laser pulse, considering the influence of the sparse wave on the LSDW, it is assumed that after the end of the laser pulse, a left-traveling simple wave is reflected on the LSDW, and the time of the first left-traveling simple wave reaching the surface of the film is t z , which can be expressed as [20]

(12) t z = 2 t p 2 γ b γ b 1 γ b + 1 2 ( γ b 1 ) ,

where γ b is the adiabatic index of plasma, and the value in the plasma region is 1.2 [20].

The model used for LSDW motion to study the relaxation process of film surface pressure depends on the three time parameters t p , t z , and t 2 D . There are three cases:

  1. When t p t z t 2 D and 0 t p / t 2D t p / t z , the axial rarefaction wave reached the film surface before the radial rarefaction wave reached the center of the film surface. When t z t t 2D , the plane attenuation model can be used at that time. When t t 2D , the cylindrical attenuation model was used. When t < t p , one-dimensional model was used.

  2. When t p t 2D t z and t p / t z t p / t 2D 1 , the radial rarefaction wave reached the film surface before the axial rarefaction wave reached the center of the film surface. When t p t t 2D , the plane attenuation model can be used. When t t 2D , the cylindrical attenuation model was used. When t < t p , one-dimensional model was used.

  3. When t p t 2D , the radial rarefaction wave reached the film surface before the axial rarefaction wave reached the center of the film surface, and the cylindrical attenuation model was used when t t p .

Assuming λ = 1 , 064 nm laser-focusing spot diameter D s = 0.8 mm , incident laser energy E = 0.1 J , incident laser pulse width t p = 10 ns , and the laser is incident on the surface of a single layer of hafnium oxide film, its t z equals 2.7 × 10 2 s and t 2D equals 5 × 10 8 s . Obviously, t p < t 2D < t z satisfied the second case. Therefore,

① When t = 0 , the surface pressure of the film becomes the initial atmospheric pressure p 0 = 1.013 × 10 5 N / m 2 .

② When 0 < t < t p , using one-dimensional model, the pressure p s10 of LSDW on the film surface is [20] as follows:

(13) p s10 = ρ 0 v L 2 γ b + 1 γ b + 1 2 γ b 2 γ b γ b 1 .

③ When t p t t 2D , the plane attenuation model can be used, and the pressure of LSDW on the film surface p s2 ( t ) is [20] as follows:

(14) p s 2 ( t ) = p s 1 t p t 2 3 .

When t = t 2D ,

(15) p s20 = p s10 t p t 2D 2 3 .

④ When t t 2D , t 0 is introduced; when t = t 0 , the surface pressure of LSDW becomes the initial atmospheric pressure p 0 , then the plasma flow field has diffused into two dimensions when t 2D t t 0 , and the cylindrical decay model can be used. The pressure p s3 ( t ) of LSDW on the surface of LSDW is as follows [20]:

(16) p s3 ( t ) = p s20 t 2 D t 6 5 ,

By p 0 = p s20 t 2D t 0 6 5 ,

(17) t 0 = t 2D p s20 p 0 5 6 .

From Eqs. (13)–(17), the expansion and diffusion of plasma shock wave on the surface pressure p s ( t ) of the film is as follows:

(18) p s ( t ) = p 0 , t = 0 ρ 0 v L 2 γ b + 1 γ b + 1 2 γ b 2 γ b γ b 1 , 0 < t < t p ρ 0 v L 2 γ b + 1 γ b + 1 2 γ b 2 γ b γ b 1 t p t 2 3 , t p t t 2D p s20 t 2D t 6 5 , t 2D t < t 0 p 0 , t t 0 .

2.2.2 Impulse per unit area of optical film surface formed by plasma shock wave

Since the area of the optical film is limited, the effective area of the optical film surface for impulse transmission by laser can only be the area of the optical film, and its impulse transmission will be limited by the size of the film area. Assuming that the film area is S and the radius is R , the total impulse transmission time is the total action time t 0 of the laser plasma shock wave. From the previous analysis, it can be seen that t p < t 2D < t z . Therefore, the impulse I st obtained by the film changes from Eq. (5) to the following equation:

(19) I st = 0 t 0 S p s ( t ) d S d t .

The following are discussed according to different time intervals:

  1. In the range 0 t t 2D , plasma and shock wave develop in one dimension, and their radial impulse transfer area is approximately the laser spot area π ( D s / 2 ) 2 . During this period, the impulse obtained by the film is as follows:

    (20) I st1 = π ( D s / 2 ) 2 0 t 2D p s ( t ) d t = π ( D s / 2 ) 2 3 p s10 t p p s10 p s20 1 2 2 p s10 t p .

  2. When t t 2D , considering the two-dimensional diffusion effect of plasma and shock wave, the lateral explosion wave is attenuated in the form of the cylindrical explosion wave. According to the cylindrical explosion wave theory, if the time required for the explosion wave to propagate to the film boundary R is t d , then:

(21) R = ( D s / 2 ) t d t 2D .

Therefore,

(22) t d = t 2D 2 R D s 2 .

Thus, the impulse obtained by the optical film during t 2D t t d period is as follows:

(23) I st2 = π ( D s / 2 ) 2 t 2D t d p s20 t 2D t 6 5 t t 2D d t = 5 4 π ( D s / 2 ) 2 p s20 t 2D t d t 2D 4 5 1 .

3) In the range t d t t 0 , the impulse transfer area is limited by the size of the film area, and will not increase any more, but will always be the film area S = π R 2 . During this period, the impulse obtained by the film is as follows:

(24) I st3 = π R 2 t d t 0 p s20 t 2D t 6 / 5 d t = 5 π R 2 p s20 t 2D t 2D t d 1 5 t 2D t 0 1 5 .

Therefore, the total impulse transmitted to the film during the whole action time is as follows:

(25) I st = I st1 + I st2 + I st3 = π ( D s / 2 ) 2 3 p s10 t p p s10 p s20 1 / 2 2 p s10 t p + 5 4 π ( D s / 2 ) 2 p s20 t 2D t d t 2D 4 / 5 1 + 5 π R 2 p s20 t 2D t 2D t d 1 / 5 t 2D t 0 1 / 5 .

2.2.3 Impulse coupling coefficient per unit area of optical film surface formed by plasma shock wave

Impulse energy coupling coefficient can reflect the ability of the film surface to obtain impulse, which is defined as the magnitude of the impulse obtained by the film under the action of unit laser beam energy [20], namely:

(26) j = I s t E .

The unit of j is N s / J .

When the incident laser wavelength λ = 1 , 064 nm , the energy E = 0.1 J , the pulse width t p = 10 ns , the diameter of incident laser beam d = 28 mm , the focal length of the focusing lens f = 350 mm , the distance between the film surface and the focal plane of the focusing lens z 0 = 5 mm , and the film radius R = 5 mm , taking single-layer Al2O3 and HfO2 thin films as examples, the impulse I st and impulse coupling coefficient j of optical thin films per unit area are calculated according to Eqs. (25) and (26). The calculation results are shown in Table 1.

Table 1

Calculation results

I st (N s) j (N s/J)
Al2O3 thin film 2.60 × 10−4 2.60 × 10−3
HfO2 thin film 2.27 × 10−4 2.27 × 10−3

It can be seen from Table 1 that under the same laser action parameters, a single layer of Al2O3 thin film with a smaller atomic number Z has a larger impulse and impulse coupling coefficient than HfO2 thin film, indicating that the single layer of Al2O3 thin film has a stronger ability to obtain impulse. In the literature [21], an experimental study on the impulse formation of laser target was carried out. It was found that copper target can obtain more impulse than aluminum target under the same conditions, and the impulse coupling coefficient of copper target and aluminum target is about 10−4 N s/J, which is basically consistent with the conclusion in this article. In this article, the impulse coupling coefficient obtained by single-layer aluminum oxide and hafnium oxide film is about 10−3 N s/J. The difference in data is mainly caused by the difference in laser parameters and materials.

Different materials have different impulses. The main reason for this phenomenon is that under the action of the same laser action parameters, the effects on plasma and shock waves generated by different films are different. In the process of reverse bremsstrahlung absorption during plasma formation, the reverse bremsstrahlung absorption coefficient is related to the atomic number Z of the film material because the atomic number Z of Al2O3 thin film is smaller than that of HfO2 thin film. As a result, the plasma absorbs more laser energy, expands more violently, and obtains greater impulse.

3 Analysis of factors affecting the impulse of optical film surface formed by shock wave

The variation of I s t and j would be studied when the incident laser energy ( E ), the distance between the film surface and the focal plane of the focusing lens ( z 0 ), the radius of film (R), the focal length of the focusing lens (f), and the incident laser pulse width ( t p ) change, respectively.

3.1 Incident laser energy

When the incident laser energy E is changed and other parameters are unchanged ( z 0 = 5 mm , f = 350 mm , t p = 10 ns , R = 5 mm ), the curves of I st under different E actions are obtained, as shown in Figure 1. As shown in Figure 1, I st increases with increasing E . And the impulse of the single-layer Al2O3 film is larger than that of HfO2 thin film.

Figure 1 
                  Curves of 
                        
                           
                           
                              
                                 
                                    I
                                 
                                 
                                    st
                                 
                              
                           
                           {I}_{\text{st}}
                        
                      
                     versus 
                     
                        
                           
                           
                              E
                           
                           E
                        
                      (
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                              =
                              5
                              
                              mm
                           
                           {z}_{0}=5\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              f
                              =
                              350
                              
                              mm
                           
                           f=350\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                              =
                              10
                              
                              ns
                           
                           {t}_{\text{p}}=10\hspace{.25em}\text{ns}
                        
                     , 
                        
                           
                           
                              R
                              =
                              5
                              
                              mm
                           
                           R=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 1

Curves of I st versus E ( z 0 = 5 mm , f = 350 mm , t p = 10 ns , R = 5 mm ).

The curve of impulse coupling coefficient j versus E is shown in Figure 2. With the increase of E , j decreases, and the impulse coupling coefficient of the single-layer Al2O3 thin film is larger than that of HfO2 thin film, indicating that the ability of the single-layer Al2O3 thin film to obtain impulse is stronger than that of HfO2 thin film.

Figure 2 
                  Curves of 
                        
                           
                           
                              j
                           
                           j
                        
                      
                     versus E (
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                              =
                              5
                              
                              mm
                           
                           {z}_{0}=5\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              f
                              =
                              350
                              
                              mm
                           
                           f=350\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                              =
                              10
                              
                              ns
                           
                           {t}_{\text{p}}=10\hspace{.25em}\text{ns}
                        
                     , 
                        
                           
                           
                              R
                              =
                              5
                              
                              mm
                           
                           R=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 2

Curves of j versus E ( z 0 = 5 mm , f = 350 mm , t p = 10 ns , R = 5 mm ).

In 1972, Pirri et al. carried out experiments on impulse coupling coefficient of laser-supported detonation wave acting on a flat plate. The experimental results show that the impulse coupling coefficient decreases with the increase of laser incident energy [22]. In 2011, Chen et al. also studied the change rule of the impulse coupling coefficient of the plate under the action of laser-supported detonation wave with the laser incident energy, and obtained a conclusion through experiments: the impulse obtained by the plate increases with the increase of the laser incident energy, while the impulse coupling coefficient decreases with the increase of the laser incident energy [23]. The results in these works are basically consistent with the analysis results in this article, but the results given in the works are the change curves of data fitting obtained from experimental measurements, while this study is the theoretical curve obtained from modeling.

3.2 Distance between the film surface and the focal plane of the focusing lens z 0

By changing the distance between the film surface and the focal plane of the focusing lens z 0 , other parameters are unchanged ( E = 0.1 J , f = 350 mm , t p = 10 ns , R = 5 mm ), and the curves of I st versus z 0 are obtained, and I st decreases as z 0 increases, as shown in Figure 3. The impulse coupling coefficient j varies with z 0 , and j decreases with the increase of z 0 , as shown in Figure 4.

Figure 3 
                  Curves of 
                        
                           
                           
                              
                                 
                                    I
                                 
                                 
                                    st
                                 
                              
                           
                           {I}_{\text{st}}
                        
                      
                     versus 
                     
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                           
                           {z}_{0}
                        
                      (
                        
                           
                           
                              E
                              =
                              0.1
                              
                              J
                           
                           E=0.1\hspace{.25em}\text{J}
                        
                     , 
                        
                           
                           
                              f
                              =
                              350
                              
                              mm
                           
                           f=350\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                              =
                              10
                              
                              ns
                           
                           {t}_{\text{p}}=10\hspace{.25em}\text{ns}
                        
                     , 
                        
                           
                           
                              R
                              =
                              5
                              
                              mm
                           
                           R=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 3

Curves of I st versus z 0 ( E = 0.1 J , f = 350 mm , t p = 10 ns , R = 5 mm ).

Figure 4 
                  Curves of 
                        
                           
                           
                              j
                           
                           j
                        
                      
                     versus 
                     
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                           
                           {z}_{0}
                        
                      (
                        
                           
                           
                              E
                              =
                              0.1
                              
                              J
                           
                           E=0.1\hspace{.25em}\text{J}
                        
                     , 
                        
                           
                           
                              f
                              =
                              350
                              
                              mm
                           
                           f=350\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                              =
                              10
                              
                              ns
                           
                           {t}_{\text{p}}=10\hspace{.25em}\text{ns}
                        
                     , 
                        
                           
                           
                              R
                              =
                              5
                              
                              mm
                           
                           R=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 4

Curves of j versus z 0 ( E = 0.1 J , f = 350 mm , t p = 10 ns , R = 5 mm ).

In 2007, Zhang et al. [24] conducted an experimental study on the influence of defocusing amount z 0 on the mechanical effects of plasma shock waves. The experimental results show that the impulse coupling coefficient decreases with the increase of z 0, which is consistent with the theoretical analysis results in this study. In 2011, Chen et al. [23] obtained the curve of impulse and impulse coupling coefficient changing with defocusing amount z 0, and found that impulse and impulse coupling coefficient decrease with increasing z 0 , which is consistent with the conclusion of this study.

This is because as the distance z 0 between the film surface and the focal plane of the focusing lens decreases, the radius of the focusing spot on the film surface will decrease, so the laser power density on the film surface will increase, and the impulse transferred to the film will be greater.

3.3 Radius of film R

By changing the radius of film R , other parameters are unchanged ( E = 0.1 J , f = 350 mm , t p = 10 ns , R = 5 mm ), and the curves of I st versus R are obtained, and I st increases as R increases, as shown in Figure 5. The impulse coupling coefficient j varies with R , and j increases with the increase of R , as shown in Figure 6.

Figure 5 
                  Curves of 
                        
                           
                           
                              
                                 
                                    I
                                 
                                 
                                    st
                                 
                              
                           
                           {I}_{\text{st}}
                        
                      
                     versus 
                     
                        
                           
                           
                              R
                           
                           R
                        
                      (
                        
                           
                           
                              E
                              =
                              0.1
                              
                              J
                           
                           E=0.1\hspace{.25em}\text{J}
                        
                     , 
                        
                           
                           
                              f
                              =
                              350
                              
                              mm
                           
                           f=350\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                              =
                              10
                              
                              ns
                           
                           {t}_{\text{p}}=10\hspace{.25em}\text{ns}
                        
                     , 
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                              =
                              5
                              
                              mm
                           
                           {z}_{0}=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 5

Curves of I st versus R ( E = 0.1 J , f = 350 mm , t p = 10 ns , z 0 = 5 mm ).

Figure 6 
                  Curves of 
                        
                           
                           
                              j
                           
                           j
                        
                      
                     versus 
                     
                        
                           
                           
                              R
                           
                           R
                        
                      (
                        
                           
                           
                              E
                              =
                              0.1
                              
                              J
                           
                           E=0.1\hspace{.25em}\text{J}
                        
                     , 
                        
                           
                           
                              f
                              =
                              350
                              
                              mm
                           
                           f=350\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                              =
                              10
                              
                              ns
                           
                           {t}_{\text{p}}=10\hspace{.25em}\text{ns}
                        
                     , 
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                              =
                              5
                              
                              mm
                           
                           {z}_{0}=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 6

Curves of j versus R ( E = 0.1 J , f = 350 mm , t p = 10 ns , z 0 = 5 mm ).

In 2011, Chen et al. [23] obtained the curve of impulse and impulse coupling coefficient changing with target area when the critical area of the laser target is 77.33, 113.91, and 128.85 cm2, respectively. The results showed that impulse and impulse coupling coefficient increase with the increase of target area, and the larger the critical area of the target, the larger the impulse and impulse coupling coefficient. Because the target area is proportional to R, this is consistent with the conclusion in this study.

The size of film radius R has a significant impact on impulse transfer. When the film radius increases, the film area increases, and the impulse transferred to the film also increases, because the impulse transferred to the film is the time accumulation effect of the force, and the action area of the force increases with the expansion of the plasma. Within a certain range, when the film surface area increases, the impulse transferred to the film by the laser also increases, but not unlimited. It can be seen from the previous analysis that the impulse has a total transmission time t 0 . When the transmission time is t 0 , the corresponding film area is the limit area of impulse transmission. Only the film surface within this limit area has impulse transmission, and there is no impulse transmission beyond this limit area. In addition, when the applied laser energy is small, there is no significant difference in the impulse obtained by the films with different film areas for the same energy. This is because at this time, the limit area of impulse transmission decreases with the weakening of the film surface plasma, so that all the impulses caused by laser plasma can be transmitted to the film.

3.4 Focal length of focusing lens f

By changing the focal length of focusing lens f , other parameters are unchanged ( E = 0.1 J , z 0 = 5 mm , t p = 10 ns , R = 5 mm ), and the curves of I st versus f are obtained, and I st increases as f increases, as shown in Figure 7. The impulse coupling coefficient j varies with f , and j increases with the increase of f , as shown in Figure 8.

Figure 7 
                  Curves of 
                        
                           
                           
                              
                                 
                                    I
                                 
                                 
                                    st
                                 
                              
                           
                           {I}_{\text{st}}
                        
                      
                     versus 
                     
                        
                           
                           
                              f
                           
                           f
                        
                      (
                        
                           
                           
                              E
                              =
                              0.1
                              
                              J
                           
                           E=0.1\hspace{.25em}\text{J}
                        
                     , 
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                              =
                              5
                              
                              mm
                           
                           {z}_{0}=5\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                              =
                              10
                              
                              ns
                           
                           {t}_{\text{p}}=10\hspace{.25em}\text{ns}
                        
                     , 
                        
                           
                           
                              R
                              =
                              5
                              
                              mm
                           
                           R=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 7

Curves of I st versus f ( E = 0.1 J , z 0 = 5 mm , t p = 10 ns , R = 5 mm ).

Figure 8 
                  Curves of 
                        
                           
                           
                              j
                           
                           j
                        
                      
                     versus 
                     
                        
                           
                           
                              f
                           
                           f
                        
                      (
                        
                           
                           
                              E
                              =
                              0.1
                              
                              J
                           
                           E=0.1\hspace{.25em}\text{J}
                        
                     , 
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                              =
                              5
                              
                              mm
                           
                           {z}_{0}=5\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                              =
                              10
                              
                              ns
                           
                           {t}_{\text{p}}=10\hspace{.25em}\text{ns}
                        
                     , 
                        
                           
                           
                              R
                              =
                              5
                              
                              mm
                           
                           R=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 8

Curves of j versus f ( E = 0.1 J , z 0 = 5 mm , t p = 10 ns , R = 5 mm ).

Because the focal length f of the focusing lens directly affects the size of the laser focus spot on the film surface, that is, it directly affects the laser power density in the focal spot area, thus affecting the characteristics of plasma generated on the film surface. With the increase of f, the obtained impulse and the impulse coupling coefficient increase.

3.5 Incident laser pulse width t p

By changing the incident laser pulse width t p , other parameters are unchanged ( E = 0.1 J , z 0 = 5 mm , f = 350 mm , R = 5 mm ), and the curves of I st versus t p are obtained, and I st decreases as t p increases, as shown in Figure 9. The impulse coupling coefficient j varies with t p , and j decreases with the increase of t p , as shown in Figure 10.

Figure 9 
                  Curves of 
                        
                           
                           
                              
                                 
                                    I
                                 
                                 
                                    st
                                 
                              
                           
                           {I}_{\text{st}}
                        
                      
                     versus 
                     
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                           
                           {\text{t}}_{p}
                        
                      (
                        
                           
                           
                              E
                              =
                              0.1
                              
                              J
                           
                           E=0.1\hspace{.25em}\text{J}
                        
                     , 
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                              =
                              5
                              
                              mm
                           
                           {z}_{0}=5\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              f
                              =
                              350
                              
                              mm
                           
                           f=350\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              R
                              =
                              5
                              
                              mm
                           
                           R=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 9

Curves of I st versus t p ( E = 0.1 J , z 0 = 5 mm , f = 350 mm , R = 5 mm ).

Figure 10 
                  Curves of 
                        
                           
                           
                              j
                           
                           j
                        
                      
                     versus 
                     
                        
                           
                           
                              
                                 
                                    t
                                 
                                 
                                    p
                                 
                              
                           
                           {\text{t}}_{p}
                        
                      (
                        
                           
                           
                              E
                              =
                              0.1
                              
                              J
                           
                           E=0.1\hspace{.25em}\text{J}
                        
                     , 
                        
                           
                           
                              
                                 
                                    z
                                 
                                 
                                    0
                                 
                              
                              =
                              5
                              
                              mm
                           
                           {z}_{0}=5\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              f
                              =
                              350
                              
                              mm
                           
                           f=350\hspace{.25em}\text{mm}
                        
                     , 
                        
                           
                           
                              R
                              =
                              5
                              
                              mm
                           
                           R=5\hspace{.25em}\text{mm}
                        
                     ).
Figure 10

Curves of j versus t p ( E = 0.1 J , z 0 = 5 mm , f = 350 mm , R = 5 mm ).

This is because the larger the t p is, the smaller the laser power acting on the unit area of the film surface is. Therefore, the smaller the impulse I st and impulse coupling coefficient j formed on the film surface are.

4 Conclusion

  1. In this study, the impulse calculation model of plasma shock wave acting on the film surface is established. The impulse I st and impulse coefficient j of the single-layer Al2O3 and HfO2 films on the unit area are calculated, respectively, under the specified parameters. The results show that the impulse obtained by Al2O3 films is larger, which indicates that I st and j are related to the atomic number Z of the film material. The smaller Z is, the larger I st and j are.

  2. The factors affecting the size of I st and j are analyzed. The results show that I st and j of the two films increase with the increase of E and f, and I st and j decrease with the increase of z 0 and t p . In the total impulse transmission time t 0 , I st and j increase with the increase of R , and when the action time exceeds t 0 , I st and j are both zero.

The above research shows that the pressure and impulse of the plasma shock wave on the optical film are the factors that cannot be ignored in the process of laser damage to most of the films. Considering that there are also thermal effects, field effects and plasma flash splashing away vaporized film material and other destructive effects in the process of laser damage to the film, it can be asserted that the mechanical effect of plasma shock wave on multilayer optical film, combined with the stress existing in the film, can strengthen the destructive effect of other laser effects on the film system, leading to the collapse of a large area of the film.

In addition, since the laser-supported detonation wave drives the light ship, its propulsion performance is mainly reflected by thrust, impulse, and impulse coupling coefficient. Therefore, the research in this study can also provide theoretical basis for the impulse coupling in laser propulsion technology and the selection of light ship, and has reference significance for the selection of laser parameters, focusing system parameters, and target parameters in laser propulsion research. Since the impulse coupling coefficient is an important parameter to calculate the laser clearance of space debris, this work can also provide a way for the research of laser clearance of space debris.

  1. Funding information: This work was supported in part by the National Natural Science Foundation of China (NSFC) (No. 61378050 and No. 62205263).

  2. Author contributions: All authors have accepted the responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: All data generated or analyzed during this study are included in this published article.

References

[1] Lai QY, Feng GY, Yan J. Damage threshold of substrates for nanoparticles removal using a laser-induced plasma shockwave. Appl Surf Sci. 2021;539:148282.10.1016/j.apsusc.2020.148282Search in Google Scholar

[2] Zhu CQ, Dyomin V, Yudin N. Laser induced damage threshold of nonlinear GaSe and GaSe: In crystals upon exposure to pulsed radiation at a wavelength of 2.1 mum. APPL SCI-BASEL. 2021;11(3):1208.10.3390/app11031208Search in Google Scholar

[3] Xie LY, Zhang JL, Zhang ZY. Rectangular multilayer dielectric gratings with broadband high diffraction efficiency and enhanced laser damage resistance. OPT EXPRESS. 2021;29(2):2669–78.10.1364/OE.415847Search in Google Scholar PubMed

[4] Lian X, Yao WD, Liu WL. KNa2ZrF7:A mixed-metal fluoride exhibits phase matchable second-harmonic-generation effect and high laser induced damage threshold. Inorg Chem. 2021;60(1):19–23.10.1021/acs.inorgchem.0c03198Search in Google Scholar PubMed

[5] Shan C, Zhao YA, Zhang XH, Hu GH, Wang YL, Peng XC, et al. Study on laser damage threshold of optical element surface based on Gaussian pulsed laser spatial resolution. Chin J Lasers. 2018;45(1):0104002.10.3788/CJL201845.0104002Search in Google Scholar

[6] Ling XL, Liu SH, Liu XF. Enhancement of laser-induced damage threshold of optical coatings by ion-beam etching in vacuum environment. Optik. 2020;200:163429.10.1016/j.ijleo.2019.163429Search in Google Scholar

[7] Kumar S, Shankar A, Kishore N, Mukherjee C, Kamparath R, Thakur S. Laser-induced damage threshold study on TiO2/SiO2 multilayer reflective coatings. Indian J Phys. 2020;94:105–15.10.1007/s12648-019-01445-3Search in Google Scholar

[8] Li YY, Zhang WY, Liu Z, Li MX, Fu XH. Cumulative effect of thin film laser damage under S-on-1 measurement mode. Laser Technol. 2018;42(1):39–42.Search in Google Scholar

[9] Gregg DW, Thomas SJ. Momentum transfer and plasma formation above a surface with a high-power CO2 laser. J Appl Phys. 1966;37(7):2787–9.10.1063/1.1782123Search in Google Scholar

[10] Piri AN, Schlier R, Northam D. Momentum transfer and plasma formation above a surface with a high-power CO2 laser. Appl Phys Lett. 1972;21(3):79–81.Search in Google Scholar

[11] Hettche LR, Tucker TR, Schriempf JT, Stegman RL, Metz SA. Mechanical response and thermal coupling of metallic targets to high-intensity 1.06μm laser radiation. J Appl Phys. 1976;47(4):1415–21.10.1063/1.322802Search in Google Scholar

[12] Duzy C, Woodroffe JA, Hsia JC, Ballantyne A. Interaction of a pulsed XeF laser with an aluminum surface. Appl Phys Lett. 1980;37(6):542–4.10.1063/1.91979Search in Google Scholar

[13] Bass M, Nassar MA, Swimm RT. Impulse coupling to aluminum resulting from Nd: glass laser irradiation induced material removal. J Appl Phys. 1987;61(3):1137–44.10.1063/1.338158Search in Google Scholar

[14] Fu YS, Wang CK. Investigation on the laser-generated stress waves and its application. Laser J. 1984;5(1):1-6,11.Search in Google Scholar

[15] Fu YS, Wang CK. The pressure attenuation of the air spot explode by the laser focusing. Laser J. 1983;4(3):174–5.Search in Google Scholar

[16] Fu YS, Wang CK. Measurement of shock wave pressure induced by intense pulse laser in air - an application of KL-I Sensor. Mech Eng. 1980;2:61–3.Search in Google Scholar

[17] Liang XB, Zhang GP, Fang YW. Research progress of impulse coupling coefficient of pulsed laser to aluminum target fragments. Laser & Infrared. 2020;50(9):1027–34.Search in Google Scholar

[18] Yu CH, Zhou WJ, Chang H, Chen YF. Experimental research on impulse coupling characteristics and plasma plume dynamics of a nanosecond pulsed laser irradiated aluminum target. IEEE Access. 2020;8:205272–81.10.1109/ACCESS.2020.3037302Search in Google Scholar

[19] Xu YF, Yang L, Zhou DJ, Liu BK, Li QW, Shi WB, et al. Experimental study on the dynamics and parameters of nanosecond laser-induced aluminum plasma. J Phys D Appl Phys. 2022;32(55):325201.10.1088/1361-6463/ac6a27Search in Google Scholar

[20] Lu J, Ni XW, He AZ. Physics of interaction between laser and materials. Beijing: China Machine Press; 1996.Search in Google Scholar

[21] Lu J, Ni XW, He AZ. Momentum formation on a target surface with a laser. J Nanjing Univ Sci Technol. 1994;3:61–4.Search in Google Scholar

[22] Pirri AN, Schilier R, Northam. D. Momentum transfer and plasma formation above a surface with a high power CO2 laser. Appl Phys Lett. 1972;21(3):79–81.10.1063/1.1654297Search in Google Scholar

[23] Chen L, Lu JY, Wu JY, Feng CG. Laser supported detonation wave. BeiJing: National Defence Industry Press; 2011.Search in Google Scholar

[24] Zhang YZ, Wang GA, Shen ZH, Ni XW, Lu J. Influence of intensity on mechanical effect of laser plasma shock wave. Laser Technol. 2007;(6):659–62.Search in Google Scholar

Received: 2022-10-19
Revised: 2023-02-02
Accepted: 2023-02-23
Published Online: 2023-03-20

© 2023 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Dynamic properties of the attachment oscillator arising in the nanophysics
  3. Parametric simulation of stagnation point flow of motile microorganism hybrid nanofluid across a circular cylinder with sinusoidal radius
  4. Fractal-fractional advection–diffusion–reaction equations by Ritz approximation approach
  5. Behaviour and onset of low-dimensional chaos with a periodically varying loss in single-mode homogeneously broadened laser
  6. Ammonia gas-sensing behavior of uniform nanostructured PPy film prepared by simple-straightforward in situ chemical vapor oxidation
  7. Analysis of the working mechanism and detection sensitivity of a flash detector
  8. Flat and bent branes with inner structure in two-field mimetic gravity
  9. Heat transfer analysis of the MHD stagnation-point flow of third-grade fluid over a porous sheet with thermal radiation effect: An algorithmic approach
  10. Weighted survival functional entropy and its properties
  11. Bioconvection effect in the Carreau nanofluid with Cattaneo–Christov heat flux using stagnation point flow in the entropy generation: Micromachines level study
  12. Study on the impulse mechanism of optical films formed by laser plasma shock waves
  13. Analysis of sweeping jet and film composite cooling using the decoupled model
  14. Research on the influence of trapezoidal magnetization of bonded magnetic ring on cogging torque
  15. Tripartite entanglement and entanglement transfer in a hybrid cavity magnomechanical system
  16. Compounded Bell-G class of statistical models with applications to COVID-19 and actuarial data
  17. Degradation of Vibrio cholerae from drinking water by the underwater capillary discharge
  18. Multiple Lie symmetry solutions for effects of viscous on magnetohydrodynamic flow and heat transfer in non-Newtonian thin film
  19. Thermal characterization of heat source (sink) on hybridized (Cu–Ag/EG) nanofluid flow via solid stretchable sheet
  20. Optimizing condition monitoring of ball bearings: An integrated approach using decision tree and extreme learning machine for effective decision-making
  21. Study on the inter-porosity transfer rate and producing degree of matrix in fractured-porous gas reservoirs
  22. Interstellar radiation as a Maxwell field: Improved numerical scheme and application to the spectral energy density
  23. Numerical study of hybridized Williamson nanofluid flow with TC4 and Nichrome over an extending surface
  24. Controlling the physical field using the shape function technique
  25. Significance of heat and mass transport in peristaltic flow of Jeffrey material subject to chemical reaction and radiation phenomenon through a tapered channel
  26. Complex dynamics of a sub-quadratic Lorenz-like system
  27. Stability control in a helicoidal spin–orbit-coupled open Bose–Bose mixture
  28. Research on WPD and DBSCAN-L-ISOMAP for circuit fault feature extraction
  29. Simulation for formation process of atomic orbitals by the finite difference time domain method based on the eight-element Dirac equation
  30. A modified power-law model: Properties, estimation, and applications
  31. Bayesian and non-Bayesian estimation of dynamic cumulative residual Tsallis entropy for moment exponential distribution under progressive censored type II
  32. Computational analysis and biomechanical study of Oldroyd-B fluid with homogeneous and heterogeneous reactions through a vertical non-uniform channel
  33. Predictability of machine learning framework in cross-section data
  34. Chaotic characteristics and mixing performance of pseudoplastic fluids in a stirred tank
  35. Isomorphic shut form valuation for quantum field theory and biological population models
  36. Vibration sensitivity minimization of an ultra-stable optical reference cavity based on orthogonal experimental design
  37. Effect of dysprosium on the radiation-shielding features of SiO2–PbO–B2O3 glasses
  38. Asymptotic formulations of anti-plane problems in pre-stressed compressible elastic laminates
  39. A study on soliton, lump solutions to a generalized (3+1)-dimensional Hirota--Satsuma--Ito equation
  40. Tangential electrostatic field at metal surfaces
  41. Bioconvective gyrotactic microorganisms in third-grade nanofluid flow over a Riga surface with stratification: An approach to entropy minimization
  42. Infrared spectroscopy for ageing assessment of insulating oils via dielectric loss factor and interfacial tension
  43. Influence of cationic surfactants on the growth of gypsum crystals
  44. Study on instability mechanism of KCl/PHPA drilling waste fluid
  45. Analytical solutions of the extended Kadomtsev–Petviashvili equation in nonlinear media
  46. A novel compact highly sensitive non-invasive microwave antenna sensor for blood glucose monitoring
  47. Inspection of Couette and pressure-driven Poiseuille entropy-optimized dissipated flow in a suction/injection horizontal channel: Analytical solutions
  48. Conserved vectors and solutions of the two-dimensional potential KP equation
  49. The reciprocal linear effect, a new optical effect of the Sagnac type
  50. Optimal interatomic potentials using modified method of least squares: Optimal form of interatomic potentials
  51. The soliton solutions for stochastic Calogero–Bogoyavlenskii Schiff equation in plasma physics/fluid mechanics
  52. Research on absolute ranging technology of resampling phase comparison method based on FMCW
  53. Analysis of Cu and Zn contents in aluminum alloys by femtosecond laser-ablation spark-induced breakdown spectroscopy
  54. Nonsequential double ionization channels control of CO2 molecules with counter-rotating two-color circularly polarized laser field by laser wavelength
  55. Fractional-order modeling: Analysis of foam drainage and Fisher's equations
  56. Thermo-solutal Marangoni convective Darcy-Forchheimer bio-hybrid nanofluid flow over a permeable disk with activation energy: Analysis of interfacial nanolayer thickness
  57. Investigation on topology-optimized compressor piston by metal additive manufacturing technique: Analytical and numeric computational modeling using finite element analysis in ANSYS
  58. Breast cancer segmentation using a hybrid AttendSeg architecture combined with a gravitational clustering optimization algorithm using mathematical modelling
  59. On the localized and periodic solutions to the time-fractional Klein-Gordan equations: Optimal additive function method and new iterative method
  60. 3D thin-film nanofluid flow with heat transfer on an inclined disc by using HWCM
  61. Numerical study of static pressure on the sonochemistry characteristics of the gas bubble under acoustic excitation
  62. Optimal auxiliary function method for analyzing nonlinear system of coupled Schrödinger–KdV equation with Caputo operator
  63. Analysis of magnetized micropolar fluid subjected to generalized heat-mass transfer theories
  64. Does the Mott problem extend to Geiger counters?
  65. Stability analysis, phase plane analysis, and isolated soliton solution to the LGH equation in mathematical physics
  66. Effects of Joule heating and reaction mechanisms on couple stress fluid flow with peristalsis in the presence of a porous material through an inclined channel
  67. Bayesian and E-Bayesian estimation based on constant-stress partially accelerated life testing for inverted Topp–Leone distribution
  68. Dynamical and physical characteristics of soliton solutions to the (2+1)-dimensional Konopelchenko–Dubrovsky system
  69. Study of fractional variable order COVID-19 environmental transformation model
  70. Sisko nanofluid flow through exponential stretching sheet with swimming of motile gyrotactic microorganisms: An application to nanoengineering
  71. Influence of the regularization scheme in the QCD phase diagram in the PNJL model
  72. Fixed-point theory and numerical analysis of an epidemic model with fractional calculus: Exploring dynamical behavior
  73. Computational analysis of reconstructing current and sag of three-phase overhead line based on the TMR sensor array
  74. Investigation of tripled sine-Gordon equation: Localized modes in multi-stacked long Josephson junctions
  75. High-sensitivity on-chip temperature sensor based on cascaded microring resonators
  76. Pathological study on uncertain numbers and proposed solutions for discrete fuzzy fractional order calculus
  77. Bifurcation, chaotic behavior, and traveling wave solution of stochastic coupled Konno–Oono equation with multiplicative noise in the Stratonovich sense
  78. Thermal radiation and heat generation on three-dimensional Casson fluid motion via porous stretching surface with variable thermal conductivity
  79. Numerical simulation and analysis of Airy's-type equation
  80. A homotopy perturbation method with Elzaki transformation for solving the fractional Biswas–Milovic model
  81. Heat transfer performance of magnetohydrodynamic multiphase nanofluid flow of Cu–Al2O3/H2O over a stretching cylinder
  82. ΛCDM and the principle of equivalence
  83. Axisymmetric stagnation-point flow of non-Newtonian nanomaterial and heat transport over a lubricated surface: Hybrid homotopy analysis method simulations
  84. HAM simulation for bioconvective magnetohydrodynamic flow of Walters-B fluid containing nanoparticles and microorganisms past a stretching sheet with velocity slip and convective conditions
  85. Coupled heat and mass transfer mathematical study for lubricated non-Newtonian nanomaterial conveying oblique stagnation point flow: A comparison of viscous and viscoelastic nanofluid model
  86. Power Topp–Leone exponential negative family of distributions with numerical illustrations to engineering and biological data
  87. Extracting solitary solutions of the nonlinear Kaup–Kupershmidt (KK) equation by analytical method
  88. A case study on the environmental and economic impact of photovoltaic systems in wastewater treatment plants
  89. Application of IoT network for marine wildlife surveillance
  90. Non-similar modeling and numerical simulations of microploar hybrid nanofluid adjacent to isothermal sphere
  91. Joint optimization of two-dimensional warranty period and maintenance strategy considering availability and cost constraints
  92. Numerical investigation of the flow characteristics involving dissipation and slip effects in a convectively nanofluid within a porous medium
  93. Spectral uncertainty analysis of grassland and its camouflage materials based on land-based hyperspectral images
  94. Application of low-altitude wind shear recognition algorithm and laser wind radar in aviation meteorological services
  95. Investigation of different structures of screw extruders on the flow in direct ink writing SiC slurry based on LBM
  96. Harmonic current suppression method of virtual DC motor based on fuzzy sliding mode
  97. Micropolar flow and heat transfer within a permeable channel using the successive linearization method
  98. Different lump k-soliton solutions to (2+1)-dimensional KdV system using Hirota binary Bell polynomials
  99. Investigation of nanomaterials in flow of non-Newtonian liquid toward a stretchable surface
  100. Weak beat frequency extraction method for photon Doppler signal with low signal-to-noise ratio
  101. Electrokinetic energy conversion of nanofluids in porous microtubes with Green’s function
  102. Examining the role of activation energy and convective boundary conditions in nanofluid behavior of Couette-Poiseuille flow
  103. Review Article
  104. Effects of stretching on phase transformation of PVDF and its copolymers: A review
  105. Special Issue on Transport phenomena and thermal analysis in micro/nano-scale structure surfaces - Part IV
  106. Prediction and monitoring model for farmland environmental system using soil sensor and neural network algorithm
  107. Special Issue on Advanced Topics on the Modelling and Assessment of Complicated Physical Phenomena - Part III
  108. Some standard and nonstandard finite difference schemes for a reaction–diffusion–chemotaxis model
  109. Special Issue on Advanced Energy Materials - Part II
  110. Rapid productivity prediction method for frac hits affected wells based on gas reservoir numerical simulation and probability method
  111. Special Issue on Novel Numerical and Analytical Techniques for Fractional Nonlinear Schrodinger Type - Part III
  112. Adomian decomposition method for solution of fourteenth order boundary value problems
  113. New soliton solutions of modified (3+1)-D Wazwaz–Benjamin–Bona–Mahony and (2+1)-D cubic Klein–Gordon equations using first integral method
  114. On traveling wave solutions to Manakov model with variable coefficients
  115. Rational approximation for solving Fredholm integro-differential equations by new algorithm
  116. Special Issue on Predicting pattern alterations in nature - Part I
  117. Modeling the monkeypox infection using the Mittag–Leffler kernel
  118. Spectral analysis of variable-order multi-terms fractional differential equations
  119. Special Issue on Nanomaterial utilization and structural optimization - Part I
  120. Heat treatment and tensile test of 3D-printed parts manufactured at different build orientations
Downloaded on 10.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2022-0237/html
Scroll to top button