Next Article in Journal
Two Preconditioners for Time-Harmonic Eddy-Current Optimal Control Problems
Previous Article in Journal
On Accurate Discrete-Time Dynamic Models of an Induction Machine
Previous Article in Special Issue
An Adaptive Inertia and Damping Control Strategy Based on Enhanced Virtual Synchronous Generator Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Model of the Control Problem of the Thermal Effect of a Laser Beam on a Two-Layer Biomaterial

by
Vanya Barseghyan
1 and
Svetlana Solodusha
2,*
1
Institute of Mechanics, The National Academy of Sciences of the Republic of Armenia, Yerevan State University, Yerevan 0019, Armenia
2
Melentiev Energy Systems Institute, Siberian Branch of the Russian Academy of Sciences, 664033 Irkutsk, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(3), 374; https://doi.org/10.3390/math12030374
Submission received: 6 November 2023 / Revised: 15 January 2024 / Accepted: 22 January 2024 / Published: 24 January 2024
(This article belongs to the Special Issue Numerical Simulation and Control in Energy Systems, 2nd Edition)

Abstract

:
We consider a two-layer biological object consisting of layers with different thermophysical characteristics and subjected to laser radiation. Using the method of separation of variables and methods of control theory for finite-dimensional systems, we developed a constructive approach to constructing a control function for the thermal effect of a laser beam on a two-layer biomaterial. Under the controlled thermal influence of a laser beam, the distribution of the temperature state in a two-layer biomaterial transitions from the initial state to the final one during a given time period.

1. Introduction

The study of problems in multilayer physical systems that are under the influence of concentrated or distributed sources requires considering appropriate and adequate mathematical models. At the same time, both mathematical models and research methods must be adequate.
The work in [1] presents a review of the literature on biomedical applications of lasers. It is impossible to imagine modern medicine without lasers [2,3]. The field of applying laser radiation goes far beyond the classical concepts of a laser, in particular in energy and electronics [4,5,6,7]. One of the many areas of biomedical application of lasers is their use as a tool to influence biological objects. Lasers provide the ability to precisely deliver large amounts of energy to limited areas of a material to achieve the desired response. Laser therapy can have a positive effect on the regeneration of periodontal tissues and improve the postoperative period due to its anti-inflammatory properties. Promising technical solutions and deeper knowledge about the interaction of lasers with tissues should allow the safe and effective use of laser exposure.
With the advent of new areas for applying laser radiation for the processing of biological materials, it has become necessary to develop methods for its influence and criteria for the parameters of laser emitters. Therefore, various mathematical models are being developed to solve various problems of laser influence and evaluate the results [1], in particular, the problem of choosing the modes of thermal impact of a laser beam on a biological environment. Following the authors of [1], we note that the mechanism of the effect of a laser beam on a biological environment has not yet been sufficiently studied. So, it is necessary to carry out further diversified studies on the search for laser radiation modes to develop the possibilities of laser influence and increase the effectiveness of the impact on the biological environment.
Laser therapy is increasingly recognized as a cancer treatment method due to the localized delivery of energy to tumor tissue [8]. Thermal effects on tissue, as well as damage to adjacent tissue, are still a potential problem. Therefore, mathematical modeling of laser-tissue interaction is a necessary part of clinical treatment planning. In [8], the temperature distribution during laser-induced thermotherapy in the treatment of cancer is studied using the example of a multilayer skin with an embedded tumor model. The work in [9] presents a theoretical model that simulates the thermal effects of laser radiation incident on biological tissue. The process of thermal diffusion (as well as scattering and absorption of the laser beam) in tissue is assessed using a numerical method. The work in [10] provides a theoretical analysis of thermal damage in biological tissues caused by laser irradiation. The obtained distribution of absorbed laser energy is included in the heat transfer equation in biological tissue to solve the temperature response. The influence of laser power, exposure time and beam size, as well as tissue absorption and scattering coefficients on the process of thermal damage, is considered.
This paper considers a multilayer biological material that is exposed to laser radiation. Such an object is a system with distributed parameters [11,12,13,14,15,16,17,18]. The mathematical model of the process of laser beam action on a multilayer biological material is described using partial differential heat conduction equations with boundary conditions for the beginning and end of laser heating, boundary conditions for the interaction of the outer layer of biological material and the environment and conjugation conditions between layers. The mathematical models of these objects are characterized as heterogeneous composite systems with distributed parameters; therefore, it is advisable to use methods for investigating the control problems of composite systems (of variable structure), which are addressed, in particular, in [19,20,21,22]. The scientific area of modeling the processes of laser beam impact on a multilayer biomaterial and the study of control problems for such models have not yet been sufficiently investigated in scientific publications.
This paper considers, as a multilayer system, an object consisting of two biological layers that are inhomogeneous in their thermophysical characteristics and subjected to laser radiation. It is assumed that the process of thermal action of a laser beam on a two-layer biomaterial is controlled as follows: by changing the temperature intensity of the laser beam at the upper (left) boundary of the two-layer biomaterial, we influence the thermal state in the two-layer biomaterial. This study aims to develop an analytical approach to constructing a control function for the thermal effect of a laser beam on a two-layer biomaterial under the influence of which the distribution of the temperature state transitions from a given initial state to a given final state at a given time interval. The work is related to the research carried out in [21].

2. Mathematical Model of a Two-Layer Biomaterial and Problem Statement

Let us consider a two-layer biological material infinite in coordinates x and y (Figure 1) with different thermophysical characteristics (thermal conductivity coefficients, mass density and heat capacity) of the layers.
According to the multilayer structure of the biomaterial [15,16,17,18], in the case when the thermal conductivity coefficients are constant, the differential equation of thermal conductivity is transformed into the system of the following differential equations of thermal conductivity:
ρ 1 c 1 T 1 ( z , t ) t = K 1 2 T 1 ( z , t ) z 2 , z [ 0 , l 1 ] , ρ 2 c 2 T 2 ( z , t ) t = K 2 2 T 2 ( z , t ) z 2 , z [ l 1 , l 1 + l 2 ] ,
where ρ j —mass density coefficient for the j-th layer of the biological material, j = 1 , 2 (kg/m3); c j —heat capacity coefficient for the j-th layer of the biological material (J/(kg·K)); T j ( z , t ) —temperature field for the j-th layer of the biological material (K); z—penetration depth of laser radiation in biological material (m); t—thermal exposure time (s); and K j —thermal conductivity factor for the j-th layer of the biological material (W/(m·K)).
Let us assume that the boundary conditions of the thermal effect on a two-layer biological material, respectively, are as follows
T 1 ( z , t ) | z = 0 = u ( t ) , T 2 ( z , t ) | z = l 1 + l 2 = P ( t ) ,
where u ( t ) —the temperature of the laser beam on the upper border of the two-layer biomaterial, which varies over time; P ( t ) —the temperature of the laser beam on the lower border of the two-layer biological material, which is considered known.
Let us introduce the layers matching conditions, which express the equalities of the continuity of the temperature fields along the time coordinate, and the conditions for the ideal thermal contact of the layers as follows:
T 1 ( z , t ) | z = l 1 0 = T 2 ( z , t ) | z = l 1 + 0 , K 1 T 1 ( z , t ) z | z = l 1 0 = K 2 T 2 ( z , t ) z | z = l 1 + 0 , t [ t 0 , t 2 ] .
It is assumed that the initial condition (at t = t 0 )
T 1 ( z , t ) | t = t 0 = T H ( z )
and final condition (at t = t 2 )
T 2 ( z , t ) | t = t 2 = T K ( z )
are set.
The thermal effect of a laser beam on a two-layer biomaterial is controlled as follows: by changing the intensity (temperature) of the laser beam on the upper border of the two-layer biomaterial, we influence the thermal state of the two-layer biomaterial. The boundary function u ( t ) is a control action (boundary control).
It is assumed that the admissible control u ( t ) belongs to L 2 ( t 0 , t 2 ) . The function T j ( z , t ) L 2 ( Ω ) , j = 1 , 2 , where we set Ω = ( z , t ) : z [ 0 , l 1 + l 2 ] , t [ t 0 , t 2 ] , and the functions T H ( z ) , T K ( z ) belong to L 2 ( 0 , l 1 + l 2 ) . It is also assumed that all functions are such that the following matching conditions are satisfied:
u ( t 0 ) = T H ( 0 ) , P ( t 0 ) = T H ( l 1 + l 2 ) , u ( t 2 ) = T K ( 0 ) , P ( t 2 ) = T K ( l 1 + l 2 ) .
The control problem of the thermal effect of a laser beam on a two-layer biomaterial can be formulated as follows: find a control law u ( t ) , t [ t 0 , t 2 ] of the thermal effect of a laser beam on a two-layer biomaterial under the influence of which the distribution of the thermal state (1) at the time t = t 0 transitions from the initial state (4) to the specified end state (5) on a given time interval [ t 0 , t 2 ] .

3. Reduction of the Problem to a Problem with Zero Boundary Conditions

We introduce the notation a j 2 = K j c j ρ j , j = 1 , 2 . To solve the problem posed, we consider it expedient to introduce a new variable [21]
ξ = z , z [ 0 , l 1 ] , a 1 a 2 z + l 1 1 a 1 a 2 , z [ l 1 , l 1 + l 2 ] .
The replacement of variable (7) leads to the expansion or contraction of the segment [ l 1 , l 1 + l 2 ] with respect to point z = l 1 . In this case, instead of segment [ l 1 , l 1 + l 2 ] , we consider segment [ l 1 , L ] , where L = l 1 + a 1 a 2 l 2 .
Note that, for convenience, after the replacement of variable (7), we keep all the above functions in their original notation. Thus, Equation (1) can be written in the form
T 1 ( ξ , t ) t = a 1 2 2 T 1 ( ξ , t ) ξ 2 , ξ [ 0 , l 1 ] , T 2 ( ξ , t ) t = a 1 2 2 T 2 ( ξ , t ) ξ 2 , ξ [ l 1 , L ] .
Denote
T ( ξ , t ) = T 1 ( ξ , t ) , ξ [ 0 , l 1 ] , T 2 ( ξ , t ) , ξ [ l 1 , L ] .
Therefore, two identical Equations (8) with the function T ( ξ , t ) , ξ [ 0 , L ] , t [ t 0 , t 2 ] introduced in (9) will be written with the equation
T ( ξ , t ) t = a 1 2 2 T ( ξ , t ) ξ 2 , ξ [ 0 , L ] , t [ t 0 , t 2 ] ,
with corresponding boundary conditions
T ( 0 , t ) = u ( t ) , T ( L , t ) = P ( t ) , t 0 t t 2 ,
with initial conditions
T ( ξ , t 0 ) = T H ( ξ ) , ξ [ 0 , L ] ,
with final conditions
T ( ξ , t 2 ) = T K ( ξ ) , ξ [ 0 , L ] ,
and with matching conditions at joining point ξ = l 1 of the areas
T 1 ( ξ , t ) | ξ = l 1 0 = T 2 ( ξ , t ) | ξ = l 1 + 0 , a 2 K 1 T 1 ( ξ , t ) ξ | ξ = l 1 0 = a 1 K 2 T 2 ( ξ , t ) ξ | ξ = l 1 + 0 ,
t [ t 0 , t 2 ] .
Given the inhomogeneity of the boundary conditions (11), the solution to Equation (10) will be in the form of the sum
T ( ξ , t ) = V ( ξ , t ) + W ( ξ , t ) ,
where V ( ξ , t ) is a function with homogeneous boundary conditions
V ( 0 , t ) = V ( L , t ) = 0
to be determined, and the function W ( ξ , t ) is the solution to (8), subject to
W ( 0 , t ) = u ( t ) , W ( L , t ) = P ( t ) ,
which has the form
W ( ξ , t ) = u ( t ) + ξ L P ( t ) u ( t ) .
To find the function V ( ξ , t ) , from Formulas (8), (15) and (18), we obtain
V ( ξ , t ) t = a 1 2 2 V ( ξ , t ) ξ 2 + F ( ξ , t ) , ξ [ 0 , L ] , t [ t 0 , t 2 ] ,
where
F ( ξ , t ) = ξ L u ˙ ( t ) P ˙ ( t ) u ˙ ( t ) .
Function V ( ξ , t ) satisfies the matching condition corresponding to (14) at joining point ξ = l 1 of the areas. Note that, according to (7), from condition (6), we have
T H ( l 1 + l 2 ) = T H ( L ) , T K ( l 1 + l 2 ) = T K ( L ) .
Using the approaches given in [21,22], from the initial (12) and final (13) conditions, following the matching conditions, we obtain that the function V ( ξ , t ) should meet the following initial:
V ( ξ , t 0 ) = T H ( ξ ) u ( t 0 ) ξ L P ( t 0 ) u ( t 0 ) ,
and final
V ( ξ , t 2 ) = T K ( ξ ) u ( t 2 ) ξ L P ( t 2 ) u ( t 2 )
conditions.
Given conditions (6) and (21), conditions (22) and (23) are written, respectively, as follows:
V ( ξ , t 0 ) = T H ( ξ ) T H ( 0 ) ξ L P ( t 0 ) T H ( 0 ) ,
V ( ξ , t 2 ) = T K ( ξ ) T K ( 0 ) ξ L P ( t 2 ) T K ( 0 ) .
Solving the original problem is reduced to solving the problem of controlling the thermal effect of a laser beam on a two-layer biomaterial, which is defined by Equation (19) with homogeneous boundary conditions (16). The resulting problem is formulated as follows: find a control law u ( t ) , t [ t 0 , t 2 ] under the influence of which the distribution of the thermal state defined by Equation (19) with boundary conditions (16) transitions from the given initial state (24) to the final state (25) at a given time interval [ t 0 , t 2 ] .

4. Solving the Problem

The solution to Equation (19), subject to the boundary conditions (16) and the consistency condition, is sought in the form
V ( ξ , t ) = k = 1 V k ( t ) sin π k ξ L , V k ( t ) = 2 L 0 L V ( ξ , t ) sin π k ξ L d ξ .
By presenting the functions F ( ξ , t ) , V ( ξ , t 0 ) , V ( ξ , t 2 ) in the form of Fourier series in the basis sin π k ξ L ( k = 1 , 2 , ), and substituting their values together with V ( ξ , t ) in Equations (19) and (20) and in conditions (21) and (22), we obtain that the Fourier coefficients V k ( t ) satisfy a countable number of systems of ordinary differential equations
V ˙ k ( t ) + λ k V k ( t ) = F k ( t ) , λ k = a 1 π k L 2 , k = 1 , 2 , ,
F k ( t ) = 2 π k ( 1 ) k P ˙ ( t ) u ˙ ( t ) ,
V k ( t 0 ) = T k ( H ) 2 π k T H ( 0 ) ( 1 ) k P ( t 0 ) ,
V k ( t 2 ) = T k ( K ) 2 π k T K ( 0 ) ( 1 ) k P ( t 2 ) .
Here, the Fourier coefficients of functions F ( ξ , t ) , V ( ξ , t 0 ) , V ( ξ , t 2 ) , T H ( ξ ) and T K ( ξ ) are denoted by F k ( t ) , V k ( t 0 ) , V k ( t 2 ) , T k ( H ) and T k ( K ) , respectively.
The general solution to Equation (27) with initial condition (29) has the form [12]
V k ( t ) = V k ( t 0 ) e λ k ( t t 0 ) + t 0 t F k ( τ ) e λ k ( t τ ) d τ .
Now, taking into account the final conditions (30), we obtain that the functions F k ( τ ) , τ [ t 0 , t 2 ] for each k = 1 , 2 , must meet the following relation:
t 0 t 2 F k ( τ ) e λ k ( t 2 τ ) d τ = V k ( t 2 ) V k ( t 0 ) e λ k ( t 2 t 0 ) .
Using the approaches given in [18,19,20], we obtain that the control function u ( t ) for each k = 1 , 2 , must satisfy the integral relation:
t 0 t 2 u τ e λ k τ d τ = C k ,
C k = 1 λ k π k 2 V k ( t 2 ) e λ k t 2 V k ( t 0 ) e λ k t 0 + T K ( 0 ) e λ k t 2 T H ( 0 ) e λ k t 0 ( 1 ) k t 0 t 2 P ˙ ( τ ) e λ k τ d τ .
In practice, the first few n ( k = 1 , 2 , , n ) relations (33) are usually chosen and the problem of control synthesis is solved using methods of control theory for finite-dimensional systems [23,24]. We will solve the problem following this approach. Therefore, for the first n relations, from (33), we will have
t 0 t 2 H n ( τ ) u n ( τ ) d τ = η n ,
H n ( τ ) = e λ 1 τ e λ 2 τ e λ n τ T , η n = C 1 C 2 C n T .
Here and below, “n” in the lower index of the letter means “for the first n” modes.
Relation (34) implies the validity of the following statement about complete controllability [23,24].
Proposition 1.
The first n modes of the dynamic process described by equality (27) with Equations (28)–(30) are completely controllable if and only if for any vector η n one can find control u n ( t ) , t [ t 0 , t 2 ] , which meets condition (34).
The control action u n ( t ) satisfying the integral relation (34) can be represented in the form [22,23,24]
u n ( t ) = H n T ( t ) Q n 1 η n + f n ( t ) ,
where H n T ( t ) is transposed matrix, f n ( t ) is vector function for which
t 0 t 2 H n ( t ) f n ( t ) d t = 0 , Q n = t 0 t 2 H n ( t ) H n T ( t ) d t .
Here,
Q n = t 0 t 2 e 2 λ 1 τ e ( λ 1 + λ 2 ) τ e ( λ 1 + λ n ) τ e ( λ 1 + λ 2 ) τ e 2 λ 2 τ e ( λ 2 + λ n ) τ e ( λ 1 + λ n ) τ e ( λ 2 + λ n ) τ e 2 λ n τ d τ
is a symmetric matrix for which det Q n 0 . It follows from Formula (35) that there are many control functions that solve the boundary control problem.
Having function expressions u n ( t ) , t t 0 , t 2 , from (28) and (31), we obtain an explicit expression for the function V k ( t ) in the form
V k ( t ) = = V k ( t 0 ) e λ k ( t t 0 ) + 2 ( 1 ) k π k t 0 t P ˙ ( τ ) e λ k ( t τ ) d τ 2 π k t 0 t H ˙ n T ( τ ) e λ k ( t τ ) d τ Q n 1 η n , k = 1 , 2 , , n ,
where
H ˙ n T ( τ ) = λ 1 e λ 1 τ λ 2 e λ 2 τ λ n e λ n τ .
Note, in order not to use more cumbersome formulas, it is assumed that the vector function f n ( t ) = 0 . However, in the case f n ( t ) 0 , from Formulas (35), (28) and (31), we obtain a formula similar to (37), which differs only in one term associated with the function f n ( t ) . We believe that the case f n ( t ) 0 differs from the one considered only in a technical sense, and the obtained formulas have more cumbersome notations. Therefore, all further conclusions are given under the assumption f n ( t ) = 0 .
From (26), we obtain an explicit expression for the function V n ( ξ , t ) , t [ t 0 , t 2 ] .
Further, with the help of (15) and (18), the oscillation function Q n ( ξ , t ) , 0 ξ L for the first n harmonics will be written as
T n ( ξ , t ) = k = 1 n V k ( t ) sin π k ξ L + u n ( t ) + ξ L P ( t ) u n ( t ) , ξ [ 0 , L ] , t [ t 0 , t 2 ] .
Given the notation (9) for the functions T ( ξ , t ) , at ξ [ 0 , L ] , t [ t 0 , t 2 ] , we will have
T 1 n ( ξ , t ) = k = 1 n V k ( t 0 ) e λ k ( t t 0 ) + 2 ( 1 ) k π k t 0 t P ˙ ( τ ) e λ k t τ d τ 2 π k t 0 t H ˙ n T ( τ ) e λ k ( t τ ) d τ Q n 1 η n sin π k ξ L + u n ( t ) + ξ L P ( t ) u n ( t ) , ξ 0 , l 1 ,
T 2 n ( ξ , t ) = k = 1 n V k ( t 1 ) e λ k t t 1 + 2 ( 1 ) k π k t 1 t P ˙ ( τ ) e λ k t τ d τ 2 π k t 1 t H ˙ n T ( τ ) e λ k t τ d τ Q n 1 η n sin π k ξ L + u n ( t ) + ξ L P ( t ) u n ( t ) , ξ [ l 1 , L ] ,
where
V k ( t 1 ) = V k ( t 0 ) e λ k ( t 1 t 0 ) + 2 ( 1 ) k π k t 0 t 1 P ˙ ( τ ) e λ k ( t 1 τ ) d τ 2 π k t 0 t 1 H ˙ n T ( τ ) e λ k ( t 1 τ ) d τ Q n 1 η n .
Taking into account notation (7), the functions T 1 n ( z , t ) at z [ 0 , l 1 ] , t [ t 0 , t 1 ] and T 2 n ( z , t ) at z [ l 1 , l 1 + l 2 ] , t [ t 1 , t 2 ] are presented in the following form:
T 1 n ( z , t ) = k = 1 n V k ( t 0 ) e λ k ( t t 0 ) + 2 ( 1 ) k π k t 0 t P ˙ ( τ ) e λ k ( t τ ) d τ 2 π k t 0 t H ˙ n T ( τ ) e λ k ( t τ ) d τ Q n 1 η n sin π k z L + u n ( t ) + z L P ( t ) u n ( t ) , z [ 0 , l 1 ] ,
T 2 n ( z , t ) = k = 1 n V k ( t 1 ) e λ k t t 1 + 2 ( 1 ) k π k t 1 t P ˙ ( τ ) e λ k ( t τ ) d τ 2 π k t 1 t H ˙ n T ( τ ) e λ k ( t τ ) d τ Q n 1 η n sin π k L a 1 a 2 z + l 1 1 a 1 a 2 + + u n ( t ) + 1 L P ( t ) u n ( t ) a 1 a 2 z + l 1 1 a 1 a 2 , z l 1 , l 1 + l 2 .
Note that if we assume that under boundary conditions (2), the known function of the temperature field P t is constant, then Formulas (39) and (40) take a simpler form.

5. Problem Solving for n = 2

Let us illustrate the above for n = 2 . In this case, assuming that f 2 ( t ) = 0 from (35), it follows that
u 2 ( t ) = H 2 T ( t ) Q 2 1 η 2 ,
where
H 2 ( τ ) = e λ 1 τ e λ 2 τ T , η n = C 1 C 2 T ,
Q 2 = q 11 q 12 q 21 q 22 = t 0 t 2 e 2 λ 1 τ e λ 1 + λ 2 τ e λ 1 + λ 2 τ e 2 λ 2 τ d τ , Q 2 1 = q ^ 11 q ^ 12 q ^ 21 q ^ 22 .
Note that det Q 2 = q 11 q 22 q 12 q 21 0 , where
q j j = 1 2 λ j e 2 λ j t 2 e 2 λ j t 0 , j = 1 , 2 ; q 12 = q 21 = 1 λ 1 + λ 2 e λ 1 + λ 2 t 2 e λ 1 + λ 2 t 0 .
The elements of the inverse matrix Q 2 1 have the following explicit form:
q ^ 11 = 1 det Q 1 2 λ 2 e 2 λ 2 t 2 e 2 λ 2 t 0 , q ^ 22 = 1 det Q 1 2 λ 1 e 2 λ 1 t 2 e 2 λ 1 t 0 ,
q ^ 12 = q ^ 21 = 1 det Q · 1 λ 1 + λ 2 e λ 1 + λ 2 t 2 e λ 1 + λ 2 t 0 .
Thus, we obtain
u 2 ( t ) = ( C 1 q ^ 11 + C 2 q ^ 12 ) e λ 1 t + ( C 1 q ^ 21 + C 2 q ^ 22 ) e λ 2 t ,
C 1 = 1 λ 1 π 2 V 1 ( t 2 ) e λ 1 t 2 V 1 ( t 0 ) e λ 1 t 0 + T K ( 0 ) e λ 1 t 2 T H ( 0 ) e λ 1 t 0 + t 0 t 2 P ˙ τ e λ 1 τ d τ ,
C 2 = 1 λ 2 π V 2 ( t 2 ) e λ 2 t 2 V 2 ( t 0 ) e λ 2 t 0 + T K ( 0 ) e λ 2 t 2 T H ( 0 ) e λ 2 t 0 t 0 t 2 P ˙ τ e λ 2 τ d τ .
T 2 ( ξ , t ) = V 1 ( t ) sin π ξ L + V 2 ( t ) sin 2 π ξ L + u 2 ( t ) + ξ L P ( t ) u 2 ( t ) , ξ [ 0 , L ] , t [ t 0 , t 2 ] .
V k ( t ) = V k ( t 0 ) e λ k t t 0 + t 0 t F k ( τ ) e λ k t τ d τ , k = 1 , 2 ,
F k ( t ) = 2 π k ( 1 ) k P ˙ ( t ) λ 1 ( C 1 q ^ 11 + C 2 q ^ 12 ) e λ 1 t λ 2 ( C 1 q ^ 21 + C 2 q ^ 22 ) e λ 2 t .
Given (39) and (40), function T 2 ( z , t ) has the following form:
T 2 ( z , t ) = T 12 ( z , t ) , z [ 0 , l 1 ] , T 22 ( z , t ) , z [ l 1 , l 1 + l 2 ] ,
where functions T 12 ( z , t ) at z [ 0 , l 1 ] and T 22 ( z , t ) at z [ l 1 , l 1 + l 2 ] are presented as follows:
T 12 ( z , t ) = V 1 ( t ) sin π z L + V 2 ( t ) sin 2 π z L + u 2 ( t ) + z L P ( t ) u 2 ( t ) , z [ 0 , l 1 ] ,
T 22 ( z , t ) = V 1 ( t ) sin π L a 1 a 2 z + l 1 1 a 1 a 2 + V 2 ( t ) sin 2 π L a 1 a 2 z + l 1 1 a 1 a 2 +
+ u 2 ( t ) + 1 L P ( t ) u 2 ( t ) a 1 a 2 z + l 1 1 a 1 a 2 , z [ l 1 , l 1 + l 2 ] .
Thus, using the proposed approach, we constructed explicit expressions for the control function of the thermal process for n = 2 , which solves the problem posed, and an explicit expression for the corresponding temperature distribution functions in the two-layer biomaterial.

6. Computational Experiment

Let, at t = t 0 = 0 and t = t 2 , the corresponding initial condition
T H ( z ) = z 2 , 0 z l 1 , 1 2 4 l 1 5 + 1 5 z 2 l 1 , l 1 z l 1 + l 2 ,
and final condition
T K ( z ) = l 1 3 + z 2 , 0 z l 1 , l 1 2 + z l 1 2 , l 1 z l 1 + l 2 ,
be specified.
Assume that the values l 1 = 0.02 (m), l 2 = 0.03 (m) are given. Then, we have T H ( 0 ) = 0 , T K ( 0 ) = 6.6667 × 10 3 , T H ( L ) = 0.0205 , T K ( L ) = 0.0109 . Let us choose P ( t ) in the form
P ( t ) = 0.0205 α t .
As can be seen from the above formulas, in particular (39) and (40), they are applicable for both linear and nonlinear forms of the function P ( t ) . For simplicity of numerical calculations, we take the function P ( t ) in the form (41).
We carry out calculations in two cases (for different parameter values α 1 and α 2 ).

6.1. Case 1

For a 1 = 0.001 (m2/s), a 2 = 0.02 (m2/s), we have L = l 1 + a 1 a 2 l 2 = 0.0215 (m). The values a 1 and a 2 were chosen for mathematical reasons as test ones. After performing variable substitution for the functions
T H ( ξ ) = ξ 2 , 0 ξ l 1 , 2 l 1 5 + 1 10 l 1 a 2 a 1 2 ξ l 1 1 a 1 a 2 2 , l 1 ξ L ,
and
T K ( ξ ) = l 1 3 + ξ 6 , 0 ξ l 1 , l 1 2 + a 2 a 1 ξ l 1 1 a 1 a 2 l 1 2 , l 1 ξ L ,
Fourier coefficients (at n = 2 ) are equal T 1 ( H ) = 3.4409 × 10 3 , T 2 ( H ) = 1.7488 × 10 3 , T 1 ( K ) = 5.3852 × 10 3 , T 2 ( K ) = 5.7132 × 10 4 . Then
λ 1 = a 1 π L 2 = 0.0214 , λ 2 = 2 a 2 π L 2 = 0.0854 ,
The control function is
u 2 ( t ) = β 1 e 0.0214 t + β 2 e 0.0854 t ,
and
F 1 ( t ) = γ 11 + γ 21 e 0.0214 t γ 31 e 0.0854 t ,
F 2 ( t ) = γ 12 + γ 22 e 0.0214 t γ 32 e 0.0854 t ,
V 1 ( t ) = χ 11 χ 21 e 0.0214 t + χ 31 e 0.0214 t + χ 41 e 0.0854 t ,
V 2 ( t ) = χ 12 + χ 22 e 0.0854 t + χ 32 e 0.0214 t χ 42 e 0.0854 t .
For a comparative analysis of the results obtained, further calculations are carried out for different time periods of laser exposure.
Calculations were performed for the values t = 30 (c), t = 60 (c), and t = 90 (c). The values of the main parameters obtained as a result of calculations are given in Appendix A.
The behavior of the obtained temperature distribution function T 22 ( t , L ) and the given function P ( t ) for the values of time instants t 2 is presented in Figure 2.
To estimate the deviation modulus of the functions T 22 ( t , L ) and P ( t ) , we denote
ε t 2 = max 0 t t 2 T 22 ( t , L ) P ( t ) .
A comparative analysis of the results obtained is given in the Table 1.
Thus, the results of the numerical experiment show (see Table 1) that the value of the error ε ( t 2 ) decreases when the duration ( t 2 ) of laser exposure increases.

6.2. Case 2

Leaving the same values for l 1 , l 2 , we chose a 1 = 4.96 × 10 4 (m2/s), a 2 = 4.40 × 10 4 (m2/s). For the calculations, the following values of the initial parameters were specified: c 1 = 2.4 × 10 3 , c 2 = 1.9 × 10 3 , ρ 1 = 2.2 × 10 3 , ρ 2 = 1.9 × 10 3 , K 1 = 1.3 , K 2 = 0.7 , l 1 = 0.4 , l 2 = 0.7 . These values were chosen to be close to the values from [25]. Then, L = 0.0538 (m). The Fourier coefficients for the functions T H ( ξ ) and T K ( ξ ) in this case are
T 1 ( H ) = 7 . 1190 × 10 3 , T 2 ( H ) = 2 . 7900 × 10 3 , T 1 ( K ) = 6 . 2227 × 10 3 , T 2 ( K ) = 4 . 7874 × 10 4 .
Then, we obtain
λ 1 = a 1 π L 2 = 8.3939 × 10 4 , λ 2 = 2 a 2 π L 2 = 3.3576 × 10 3 .
The control function is
u 2 ( t ) = β 1 e 8.3939 × 10 4 t + β 2 e 3.3576 × 10 3 t ,
and
F 1 ( t ) = γ 11 + γ 12 e 3.3576 × 10 3 t γ 13 e 8.3939 × 10 4 t ,
F 2 ( t ) = γ 12 + γ 22 e 3.3576 × 10 3 t γ 32 e 3.3576 × 10 3 t ,
V 1 ( t ) = χ 11 χ 21 e 8.3939 × 10 4 t + χ 31 e 3.3576 × 10 3 t χ 41 e 8.3939 × 10 4 t ,
V 2 ( t ) = χ 12 χ 22 e 8.3939 × 10 4 t + χ 32 e 3.3576 × 10 3 t χ 42 e 3.3576 × 10 3 t .
In this case, calculations were performed for values t = 180 (c) and t = 210 (c). The obtained values of the main parameters are given in Appendix B.
The behavior of the obtained temperature distribution function T 22 ( t , L ) and the given function P ( t ) is presented in Figure 3.
Note that in this case, all calculations were also performed for the values 30, 60, 90. We do not give the parameter values for them so as not to overload the article with numerical values. A comparative analysis for these cases was performed and is given below.
Using the previously introduced notation ε t 2 , we choose t 2 [ 30 , 210 ] and conduct a comparative analysis of the results. Some of the values ε t 2 for the selected range t 2 are presented in Table 2. We found that the value of ε t 2 decreases with increasing t 2 = 30, 60, 90, 180, 210.
In addition, calculations showed that it is advisable to further develop work related to the search for the value of t 2 that provides the minimum deviation of T 22 ( t , L ) from P ( t ) at a fixed value of t = t 2 . We do not consider this extremal problem; this requires a separate study.
Thus, using the proposed approach, at n = 2 , a comparative analysis was carried out for the constructed control function for the thermal process of laser exposure and some selected numerical values of the biomaterial parameters. A comparison of the results of a computational experiment showed that, under the found law of laser action, the value of the resulting temperature distribution function T 22 ( t , L ) at the end of a two-layer biological material is quite close to the values of the given function P ( t ) . The computational experiment also showed that with increasing duration of laser exposure, the maximum modulus of deviation of the functions T 22 ( t , L ) and P ( t ) decreases.

7. Conclusions

This paper considers the process of the effect of a laser beam on a two-layer biological material, the mathematical model of which is described using differential equations of thermal conductivity. Such a model is characterized as a composite dynamic system (of variable structure) with distributed parameters. By performing some substitution of variables, and then using the method of separation of variables and methods of the control theory for finite-dimensional systems, we developed a constructive approach to constructing a control function for the thermal effect of a laser beam on a two-layer biomaterial. Under the constructed controlled thermal effect of a laser beam, the distribution of the temperature state of a two-layer biomaterial transitions from the initial state to the final state at a given time interval. The constructiveness of the developed approach is illustrated by a specific example with a performed computational experiment and subsequent analysis of the results. The results obtained, even for the first two modes, show that under the influence of the constructed boundary controls of the thermal effect of the laser beam on a two-layer biomaterial, the distribution function of the temperature state at the border of the two-layer biological material is quite close to the given desired function. The proposed approach, using the Fourier method, can be extended to other non-one-dimensional processes.

Author Contributions

Conceptualization, V.B.; methodology, V.B.; software, S.S.; validation, V.B. and S.S.; formal analysis, V.B. and S.S.; investigation, V.B. and S.S.; resources, V.B. and S.S.; data curation, V.B. and S.S.; writing—original draft preparation, V.B. and S.S.; writing—review and editing, V.B. and S.S.; visualization, V.B. and S.S.; supervision, V.B. and S.S.; project administration, V.B. and S.S.; funding acquisition, V.B. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

The research by S.S. was carried out under State Assignment Project (no. FWEU-2021-0006) of the Fundamental Research Program of the Russian Federation 2021–2030 using the resources of the High-Temperature Circuit Multi-Access Research Center (Ministry of Science and Higher Education of the Russian Federation, project no 13.CKP.21.0038).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

For t 2 = 30 (s), we obtain C 1 = 0.5966 , C 2 = 1.8454 , q 11 = 60.8982 , q 22 = 978.0206 , q 12 = q 21 = 221.0502 .
Therefore, the matrix Q 2 1 has the form
Q 2 1 = 9.1433 × 10 2 2.0666 × 10 2 2.0666 × 10 2 5.6933 × 10 3 .
The constants in (41)–(46) have the following values:
α = 3.2 × 10 4 , β 1 = 9.2684 × 10 2 , β 2 = 2.2835 × 10 2 , γ 11 = 6.4 × 10 4 π ,
γ 21 = 1.2589 × 10 3 , γ 31 = 1.2416 × 10 3 , γ 12 = 3.2 × 10 4 π , γ 22 = 6.2991 × 10 4 ,
γ 32 = 6.2078 × 10 4 , χ 11 = 9.5413 × 10 3 , χ 21 = 3.0816 × 10 2 , χ 31 = 2.9502 × 10 2 , χ 41 = 1.1630 × 10 2 , χ 12 = 1.1927 × 10 3 , χ 22 = 5.9956 × 10 3 , χ 32 = 5.9004 × 10 3 , χ 42 = 3.6343 × 10 3 .
For t 2 = 60 (s), the values of C 1 , C 2 are C 1 = 1.0740 , C 2 = 21.5196 , the elements of matrix Q 2 are equal q 11 = 280.1624 , q 22 = 1.6534 × 10 5 , q 12 = q 21 = 5.6585 × 10 3 , and the matrix Q 2 1 has the following form:
Q 2 1 = 1.1560 × 10 2 3.9562 × 10 4 3.9562 × 10 4 1.9588 × 10 5 .
In this case, the constants from (41)–(46) have the following values: α = 1.6 × 10 5 , β 1 = 2.0929 × 10 2 , β 2 = 8.4641 × 10 4 , γ 11 = 3.2 × 10 4 π , γ 21 = 2.8448 × 10 4 ,
γ 31 = 4.6020 × 10 5 , γ 12 = 1.6 × 10 4 π , γ 22 = 1.4224 × 10 4 , γ 32 = 2.3010 × 10 5 ,
χ 11 = 4.7706 × 10 3 , χ 21 = 1.4404 × 10 2 , χ 31 = 6.6618 × 10 3 , χ 41 = 4.3107 × 10 4 , χ 12 = 5.9633 × 10 4 , χ 22 = 1.0717 × 10 3 , χ 32 = 1.3324 × 10 3 , χ 42 = 1.3471 × 10 4 .
At the value of t 2 = 90 (s), the calculations showed the following results: C 1 = 1.8953 , C 2 = 264.0481 , q 11 = 1.6962 × 10 3 , q 22 = 2.7787 × 10 7 , q 12 = q 21 = 1.3941 × 10 5 ,
Q 2 1 = 2.7015 × 10 3 1.3554 × 10 5 1.3554 × 10 5 1.0399 × 10 7 ,
α = 1.0667 × 10 4 , β 1 = 8.6987 × 10 3 , β 2 = 5.3145 × 10 5 , γ 11 = 2.12 × 10 4 π ,
γ 21 = 2.7689 × 10 3 , γ 31 = 2.7067 × 10 5 , γ 12 = 1.07 × 10 4 π , γ 22 = 5.9119 × 10 5 , γ 32 = 1.4448 × 10 6 , χ 11 = 3.1804 × 10 3 , χ 21 = 9.3250 × 10 3 , χ 31 = 2.7689 × 10 3 , χ 41 = 2.7067 × 10 5 , χ 12 = 3.9755 × 10 4 , χ 22 = 1.5252 × 10 3 , χ 32 = 5.5378 × 10 4 , χ 41 = 8.4584 × 10 6 .

Appendix B

For t 2 = 180 (s), we have C 1 = 4.8224 , C 2 = 3.6465 , q 11 = 210.1558 , q 22 = 349.8378 , q 12 = q 21 = 268.9029 .
The inverse matrix Q 2 1 has the following form:
Q 2 1 = 2.8872 × 10 1 2.2193 × 10 1 2.2193 × 10 1 1.7344 × 10 1 .
The constants in (41), (47)–(51) have the following values:
α = 5.3 × 10 5 , β 1 = 5.8309 × 10 1 , β 2 = 4.3777 × 10 1 , γ 11 = 1.67 × 10 4 π ,
γ 21 = 9.3572 × 10 4 , γ 31 = 3.1159 × 10 4 , γ 12 = 5.3 × 10 5 π , γ 22 = 4.6786 × 10 4 ,
γ 32 = 1.4698 × 10 4 , χ 11 = 4.0450 × 10 2 , χ 21 = 8.7808 × 10 2 , χ 31 = 2.2295 × 10 1 , χ 41 = 1.8560 × 10 1 , χ 12 = 5.0562 × 10 3 , χ 22 = 3.7120 × 10 2 , χ 32 = 6.9673 × 10 2 , χ 42 = 2.1723 × 10 2 .
For t 2 = 210 (s), the calculations showed the following results: C 1 = 4.6219 , C 2 = 4.4147 , q 11 = 251.7793 , q 22 = 461.1499 , q 12 = q 21 = 336.9543 ,
Q 2 1 = 1.7945 × 10 1 1.3112 × 10 1 1.3112 × 10 1 9.7977 × 10 2 .
α = 4.57 × 10 5 , β 1 = 2.5055 × 10 1 , β 2 = 1.7350 × 10 1 , γ 11 = 9.14 × 10 5 π ,
γ 21 = 3.7085 × 10 4 , γ 31 = 1.3389 × 10 4 , γ 12 = 4.57 × 10 5 π , γ 22 = 1.8543 × 10 4 ,
γ 32 = 6.6943 × 10 5 , χ 11 = 3.0467 × 10 2 , χ 21 = 5.3289 × 10 2 , χ 31 = 8.8362 × 10 2 , χ 41 = 7.9752 × 10 2 , χ 12 = 4.3334 × 10 3 , χ 22 = 1.5950 × 10 2 , χ 32 = 2.7613 × 10 3 , χ 42 = 1.5555 × 10 3 .

References

  1. Shangina, O.R.; Gaynutdinova, R.D. Interaction of laser radiation with biological tissues. Pract. Med. 2019, 17, 24–27. (In Russian) [Google Scholar] [CrossRef]
  2. Elagin, V.V.; Shakhova, M.A.; Karabut, M.M.; Kuznetsova, D.S.; Bredikhin, V.I.; Prodanets, N.N.; Snopova, L.B.; Baskina, O.S.; Shakhov, A.V.; Kamensky, V.A. Evaluation of Cutting Properties of a Laser Scalpel with Heavily Absorbing Coatings of an Optical Fiber. Sovrem. Tehnol. Med. 2015, 7, 55–60. [Google Scholar] [CrossRef]
  3. Skobelkin, O.K. Lasers in Surgery; Meditsina: Moscow, Russia, 1989; 256p. (In Russian) [Google Scholar]
  4. Privalov, V.E.; Turkin, V.A.; Shemanin, V.G. Lasers and green energy. Lasers Meas. Inf. 2022, 2, 5–11. (In Russian) [Google Scholar] [CrossRef]
  5. Kumar, R.; Pérez del Pino, A.; Sahoo, S.; Kumar Singh, R.; Tan, W.K.; Kar, K.K.; Matsuda, A.; Joanni, E. Laser processing of graphene and related materials for energy storage: State of the art and future prospects. Prog. Energy Combust. Sci. 2022, 91, 100981. [Google Scholar] [CrossRef]
  6. Gupta, M.C.; Carlson, D.E. Laser processing of materials for renewable energy applications. MRS Energy Sustain. 2015, 2, 2. [Google Scholar] [CrossRef]
  7. Son, E.E.; Sheindlin, M.A. Lasers in high-temperature energy materials. Proc. Russ. Acad. Sci. Energy 2011, 5, 88–103. (In Russian) [Google Scholar]
  8. Wongchadakul, P.; Rattanadecho, P. Mathematical modeling of multilayered skin with embedded tumor through combining laser ablation and nanoparticles: Effects of laser beam area, wavelength, intensity, tumor absorption coefficient and its position. Int. J. Heat Technol. 2021, 39, 89–100. [Google Scholar] [CrossRef]
  9. Cummins, L.; Nauenberg, M. Thermal effects of laser radiation in biological tissue. Biophys. J. 1983, 42, 99–102. [Google Scholar] [CrossRef]
  10. Zhou, J.; Chen, J.K.; Zhang, Y. Theoretical analysis of thermal damage in biological tissues caused by laser irradiation. Mol. Cell. Biomech. 2007, 4, 27–39. [Google Scholar] [PubMed]
  11. Belozerov, L.G.; Kireev, V.A. Composite Shells under Force and Thermal Effects; Fizmatlit: Moscow, Russia, 2003; 388p. (In Russian) [Google Scholar]
  12. Gubin, A.I.; Malaya, Y.A. Mathematical modeling of thermal processes during laser processing of materials based on the nonlinear hyperbolic equation of heat conduction. In Tekhnichna Teplofizika i Promislova Teploenergetika; New Ideology: Dnipropetrovsk, Ukraine, 2011; Volume 3, pp. 72–85. (In Russian) [Google Scholar]
  13. Svet, E.V. The nonstationary problem of heat conduction in three-dimensional formulation for laminated plates with complex shape. Bull. NTU KhPI 2013, 63, 122–131. (In Russian) [Google Scholar]
  14. Yamaoka, N.; Sugie, J. Multilayer structures of second-order linear differential equations of Euler type and their application to nonlinear oscillations Ukr. Math. J. 2006, 58, 1935–1949. [Google Scholar]
  15. Megel, Y.E.; Levkin, D.A. Mathematical model of thermal heating of a multilayer microbiological object. East Eur. J. Adv. Technol. 2012, 3/4, 4–7. (In Russian) [Google Scholar]
  16. Shupikov, A.N.; Buzko, Y.P.; Smetankina, N.V.; Ugrimov, S.V. Non-Stationary Oscillations of Multilayer Plates and Shells and Their Optimization; Scientific Publication; Publishing House of KhNEU: Kharkov, Ukraine, 2004; 252p. (In Russian) [Google Scholar]
  17. Pyatkov, S.G. Certain inverse problems for parabolic equations. J. Math. Sci. 2008, 150, 2422–2433. [Google Scholar] [CrossRef]
  18. Kantor, B.Y.; Smetankina, N.V.; Shupikov, A.N. Analysis of non-stationary temperature fields in laminated strips and plates. Int. J. Solids Struct. 2001, 38, 8673–8684. [Google Scholar] [CrossRef]
  19. Barseghyan, V.R. Control of Stage by Stage Changing Linear Dynamic Systems. Yugosl. J. Oper. Res. 2012, 22, 31–39. [Google Scholar] [CrossRef]
  20. Barseghyan, V.R. On the controllability and observability of linear dynamic systems with variable structure. In Proceedings of the 2016 International Conference “Stability and Oscillations of Nonlinear Control Systems” (Pyatnitskiy’s Conference), STAB 2016, Moscow, Russia, 1–3 June 2016. [Google Scholar]
  21. Barseghyan, V.R. The Problem of Boundary Control of Displacement at Two Ends by the Process of Oscillation of a Rod Consisting of Two Sections of Different Density and Elasticity. Mech. Solids 2023, 58, 483–491. [Google Scholar] [CrossRef]
  22. Barseghyan, V.; Solodusha, S. On the Optimal Control Problem for Vibrations of the Rod/String Consisting of Two Non-Homogeneous Sections with the Condition at an Intermediate Time. Mathematics 2022, 10, 4444. [Google Scholar] [CrossRef]
  23. Barseghyan, V.R.; Barseghyan, T.V. On an Approach to the Problems of Control of Dynamic System with Nonseparated Multipoint Conditions. Autom. Remote Control 2015, 76, 549–559. [Google Scholar] [CrossRef]
  24. Zubov, V.I. Lectures on Control Theory; Nauka: Moscow, Russia, 1975. (In Russian) [Google Scholar]
  25. Mispahov, I.S.; Magomedov, D.A. Mathematical model of the thermoelectric system for short-term storage and transportation of biological material. Bull. Dagestan State Tech. Univ. Tech. Sci. 2012, 26, 7–15. (In Russian) [Google Scholar]
Figure 1. Block diagram of the impact of a laser beam on a two-layer biological material.
Figure 1. Block diagram of the impact of a laser beam on a two-layer biological material.
Mathematics 12 00374 g001
Figure 2. Graph of functions P ( t ) (red line) and T 22 ( t , L ) (green line): (a) Case of t 2 = 30 (s). (b) Case of t 2 = 60 (s). (c) t 2 = 90 (s).
Figure 2. Graph of functions P ( t ) (red line) and T 22 ( t , L ) (green line): (a) Case of t 2 = 30 (s). (b) Case of t 2 = 60 (s). (c) t 2 = 90 (s).
Mathematics 12 00374 g002
Figure 3. Graph of functions P ( t ) (red line) and T 22 ( t , L ) (green line): (a) Case of t 2 = 180 (s). (b) Case of t 2 = 210 (s).
Figure 3. Graph of functions P ( t ) (red line) and T 22 ( t , L ) (green line): (a) Case of t 2 = 180 (s). (b) Case of t 2 = 210 (s).
Mathematics 12 00374 g003
Table 1. Values ε t 2 .
Table 1. Values ε t 2 .
t 2 (s)306090
ε t 2 7.3684 × 10 3 4.0699 × 10 3 3.3119 × 10 3
Table 2. Values ε t 2 .
Table 2. Values ε t 2 .
t 2 (s)306090180210
ε t 2 1.0522 0.2235 8.0783 × 10 2 4.1035 × 10 3 2.0224 × 10 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Barseghyan, V.; Solodusha, S. A Model of the Control Problem of the Thermal Effect of a Laser Beam on a Two-Layer Biomaterial. Mathematics 2024, 12, 374. https://doi.org/10.3390/math12030374

AMA Style

Barseghyan V, Solodusha S. A Model of the Control Problem of the Thermal Effect of a Laser Beam on a Two-Layer Biomaterial. Mathematics. 2024; 12(3):374. https://doi.org/10.3390/math12030374

Chicago/Turabian Style

Barseghyan, Vanya, and Svetlana Solodusha. 2024. "A Model of the Control Problem of the Thermal Effect of a Laser Beam on a Two-Layer Biomaterial" Mathematics 12, no. 3: 374. https://doi.org/10.3390/math12030374

APA Style

Barseghyan, V., & Solodusha, S. (2024). A Model of the Control Problem of the Thermal Effect of a Laser Beam on a Two-Layer Biomaterial. Mathematics, 12(3), 374. https://doi.org/10.3390/math12030374

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop