Next Article in Journal
Second-Order PDE Constrained Controlled Optimization Problems with Application in Mechanics
Next Article in Special Issue
A Nonhomogeneous Boundary Value Problem for Steady State Navier-Stokes Equations in a Multiply-Connected Cusp Domain
Previous Article in Journal
On the Convergence of Stochastic Process Convergence Proofs
Previous Article in Special Issue
Boundary Value Problems for ψ-Hilfer Type Sequential Fractional Differential Equations and Inclusions with Integral Multi-Point Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonstationary Radiative–Conductive Heat Transfer Problem in a Semitransparent Body with Absolutely Black Inclusions

Department of Mathematical and Computer Modelling, National Research University “Moscow Power Engineering Institute”, 111250 Krasnokazarmennay St. 14, 111250 Moscow, Russia
Submission received: 11 May 2021 / Revised: 11 June 2021 / Accepted: 22 June 2021 / Published: 23 June 2021
(This article belongs to the Special Issue Nonlinear Boundary Value Problems and Their Applications)

Abstract

:
The paper is devoted to a nonstationary initial–boundary value problem governing complex heat exchange in a convex semitransparent body containing several absolutely black inclusions. The existence and uniqueness of a weak solution to this problem are proven herein. In addition, the stability of solutions with respect to data, a comparison theorem and the results of improving the properties of solutions with an increase in the summability of the data were established. All results are global in terms of time and data.

1. Introduction

Complex heat transfer problems, in which it is necessary to simultaneously take into account the transfer of energy by thermal radiation and thermal conductivity, arise in various fields of science and industry. The discussion on the properties of complex heat transfer problems and the methods for solving them constitutes an extensive physical literature (see, for example, [1,2,3,4]).
Mathematical problems of radiative–conductive heat transfer are nonstandard, interesting and rather complicated. Heat radiation is nonlinearly dependent on temperature, and integro-differential equations or nonlocal boundary conditions are used to describe radiation heat transfer. Various nonlinear nonlocal boundary and initial–boundary value problems arise in this field.
The first mathematical results in this direction were obtained by A.N. Tikhonov [5,6] in the late 1930s. The construction of the mathematical theory of radiative–conductive heat transfer problems was continued for roughly forty years [7,8,9,10,11,12,13,14]. In the early 1990s, many mathematicians were paying attention to such problems. As a result, over the past 30 years, a large number of papers have been devoted to the solvability of complex heat transfer problems (cf. [15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63]). Naturally, the above list is not exhaustive.
To date, the solvability of various statements of complex (radiative–conductive) heat exchange problems in systems consisting either only of radiation-opaque bodies or only of radiation-semitransparent bodies has been studied in sufficient detail. At the same time, the problems of radiative–conductive heat exchange in systems consisting of both radiation-opaque and of radiation-semitransparent bodies remain to date unexplored. This specific area of study, to the best of the author’s knowledge, has only been the subject of the following articles: [64,65,66].
In this paper, the existence and uniqueness of a weak solution to a nonstationary boundary value problem governing radiative–conductive heat transfer in a semitransparent body containing several absolutely black inclusions were proven. All results are global in terms of time and data. The unknown functions u and I physically represent the absolute temperature and radiation intensity. The problem was considered in a gray approximation. The technique used was developed in [41,53,57]. In the stationary version, this problem was studied in [65].
The paper is organized as follows. In Section 2, the physical sense of the problem is explained. Section 3 is devoted to notations. In Section 4, the boundary value problem for the radiative transfer equation is considered. Section 5 contains the formulation of the main results of the paper. In addition, in this section, the problem is reduced to the equivalent initial–boundary value problem for a nonlinear operator-differential equation (Problem P ) with only one unknown function u. In Section 6, a number of important auxiliary assertions are proven. In Section 7, an auxiliary problem P [ n ] is introduced and its solvabilty is proven. In Section 8, a priori estimates for weak solutions to Problems P and P [ n ] are established. Section 9 establishes the stability of weak solutions to Problem P with respect to the data. A comparison theorem, which, in particular, implies the uniqueness of a weak solution, is also proven. Section 10 contains the proof of the existence of a weak solution to the problem P . Finally, Section 11 establishes the validity of the main results of the article.

2. Physical Statement the Problem

Let G ^ be a bounded convex domain in R 3 and { G b , j } j = 1 m be a system of strictly internal subdomains of the domain G ^ . Assume that G ¯ b , i G ¯ b , j = for all i j and G s = G ^ \ m j = 1 G ¯ b , j is a domain. We put G b = m j = 1 G b , j and G = G s G b .
We assume that each of the domains G b , j , 1 j m is an absolutely black body and the domain G s is occupied by a semitransparent optically homogeneous material with a constant absorption coefficient ϰ > 0 and a scattering coefficient s 0 .
The unknown functions u ( x , t ) and I ( ω , x , t ) physically represent the absolute temperature at point x G at moment t ( 0 , T ) and the intensity of the radiation propagating at point x G s in direction ω Ω = { ω R 3 | | ω | = 1 } , respectively. The function u is defined on the set Q T = G × ( 0 , T ) . Its restrictions to the set Q s , T = G s × ( 0 , T ) and to the set Q b , T = G b × ( 0 , T ) will be denoted by u s and u b , respectively. The function I is defined on the set D s × ( 0 , T ) , where D s = Ω × G s .
To describe the nonstationary process of radiative–conductive heat transfer, a system consisting of two heat equations and radiative transfer equation is used:
c p u s t div ( λ ( x , u s ) u s ) + 4 ϰ h ( u s ) = ϰ Ω I d ω + f , ( x , t ) Q s , T ,
c p u b t div ( λ ( x , u b ) u b ) = f , ( x , t ) Q b , T ,
ω · I + ( ϰ + s ) I = s S ( I ) + ϰ π h ( u s ) , ( ω , x , t ) D s × ( 0 , T ) .
Here, c p is the heat capacity coefficient, λ ( x , u ) is the thermal conductivity coefficient, and f is the density of heat sources. The function h ( u ) = σ 0 | u | 3 u for u > 0 corresponds to the hemispherical radiation density of an absolutely black body according to the Stefan–Boltzmann law, where σ 0 > 0 is the Stefan–Boltzmann constant.
Equation (1) describes the heat transfer process in the gray semitransparent medium G s . The terms 4 ϰ h ( u s ) and ϰ Ω I d ω in it correspond to the densities of the energies emitted and absorbed in G s , respectively. Equation (2) describes the heat transfer process in opaque inclusions G s . Equation (3) describes the transfer of radiation in a radiating, absorbing and scattering medium G s . The term ω · I = i = 1 3 ω i x i I in (3) denotes the derivative of I along the direction ω . We denote by S the scattering operator:
S ( I ) ( ω , x , t ) = 1 4 π Ω θ ( ω · ω ) I ( ω , x , t ) d ω , ( ω , x , t ) D s × ( 0 , T )
with the scattering indicatrix possessing the following properties:
θ L 1 ( 1 , 1 ) , θ 0 , 1 2 1 1 θ ( μ ) d μ = 1 .
We regard R 3 as the Euclidean space of elements x = ( x 1 , x 2 , x 3 ) , y = ( y 1 , y 2 , y 3 ) equipped with the inner product x · y = i = 1 3 x i y i . Assume that the domain G s is Lipschitz. Thus, the domains G b , j , 1 j m are also Lipschitz. We denote by d ω and d σ ( x ) the measures induced by Lebesgue measure in R 3 on Ω and G s , respectively.
We also assume that the boundary G s is piecewise smooth in the following sense. There exists a closed subset G G s such that meas ( G ; d σ ) = 0 ; moreover, for each point x G s = G s \ G , there exists a neighborhood of it, in which the boundary G s is continuously differentiable.
Note that G s = G ^ G b , where G ^ G b = . We put G ^ = G ^ \ G , G b = G b \ G . It is clear that the outward normal n to the boundary G s is defined and continuous on G s and the outward normal to the boundary of the set G b coincides with n ( x ) for x G b .
We introduce the sets:
Γ = Ω × G s , Γ = { ( ω , x ) Γ | ω · n ( x ) < 0 } , Γ + = { ( ω , x ) Γ | ω · n ( x ) > 0 } , Γ ^ = { ( ω , x ) Γ | x G ^ } , Γ b = { ( ω , x ) Γ | x G b } , Ω ( x ) = { ω Ω | ω · n ( x ) < 0 } , Ω + ( x ) = { ω Ω | ω · n ( x ) > 0 } , x G b .
Denote by I | Γ ± the values (traces) of the function I on Γ ± , where I | Γ and I | Γ + are interpreted as the values of the intensity of radiations entering into G s and coming out of G s .
Endow the system (1)–(3) with the boundary conditions:
λ ( x , u s ) u s n = 0 , ( x , t ) Q ^ T ,
λ ( x , u s ) u s n + γ ( t r u s t r u b ) = 0 , ( x , t ) Q b , T ,
λ ( x , u b ) u b n + h ( t r u b ) + γ ( t r u b t r u s ) = M + ( I | Γ + ) , ( x , t ) Q b , T ,
I | Γ = 1 π h ( t r u b ) , ( ω , x , t ) Γ b × ( 0 , T ) ,
I | Γ = J * , ( ω , x , t ) Γ ^ × ( 0 , T )
and the initial condition:
u | t = 0 = u 0 , x G .
Here, Q ^ T = G ^ × ( 0 , T ) and Q b , T = G b × ( 0 , T ) are the lateral surfaces of cylinders Q ^ T = G ^ × ( 0 , T ) and Q b , T = G b × ( 0 , T ) . By t r u s and t r u b , we denote the values (traces) of u s and u b on Q b , T .
It is assumed that the body G ^ is surrounded in a vacuum. Therefore, on the boundary of the body, the boundary condition (2) is set, which means the absence of heat flux. On the boundary G b , separating the semitransparent material G s and absolutely black inclusions G b , we set two boundary conditions. They account for incoming and outgoing energy flows using a heat transfer mechanism. In addition, it is taken into account that absolutely black inclusions emit energy and absorb the incident radiation on them. Here, γ is the heat transfer coefficient and:
M + ( I | Γ + ) ( x , t ) = Ω + ( x ) I | Γ + ( ω , x , t ) ω · n ( x ) d ω
represents the flux of radiation coming out of G s and absorbed at G b . The condition (7) means that on the boundary G b , the intensity of radiation entering into G s is equal to the intensity of radiation leaving the set G b . In (8), J * denotes the intensity of external radiation incident on G ^ .

3. Function Spaces

Throughout the paper, 1 p and p denotes the conjugate exponent of p, i.e., 1 p and 1 / p + 1 / p = 1 .
Let u be a real number or a real-valued function. We put u [ N , M ] = max { N , min { u , M } } , u [ M ] = u [ M , M ] , where N 0 < M . We also put [ u ] + = max { u , 0 } and [ u ] = max { u , 0 } .
Let S be a set where the measure d μ is given. We denote by L p ( S ; d μ ) the Lebesgue space of functions f defined on Z that are measurable with respect to the measure d μ and have the finite norm:
f L p ( S ; d μ ) = S | f ( s ) | p d μ ( s ) 1 / p , 1 p < , ess sup s S | f ( s ) | , p = .

3.1. Spaces of Functions on G, G s and G b

We set:
L p ( G ) = L p ( G ; d x ) , L p ( G s ) = L p ( G s , d x ) , L p ( G b ) = L p ( G b ; d σ ( x ) ) .
We introduce the space L p ( G s , G b ) = { ( F , J ) L p ( G s ) × L p ( G b ) } equipped with the norm:
( F , J ) L p ( G s , G b ) = 4 π ϰ F L p ( G s ) p + π J L p ( G b ) p 1 / p , 1 p < , max { F L ( G s ) , J L ( G b ) } , p = .
Let functions f, g defined on G or G s are such that f g L 1 ( G ) or f g L 1 ( G s ) . In these cases, we use the notations:
( f , g ) G = G f ( x ) g ( x ) d x , ( f , g ) G s = G s f ( x ) g ( x ) d x .
Let functions f, g defined on G b are such that f g L 1 ( G b ) . In this case, we use the notation:
( f , g ) G b = G b f ( x ) g ( x ) d σ ( x ) .
Let u be a function defined on G. We denote by u s , u b and u b , j , the restrictions of u to G s , G b and G b , j , 1 j m , respectively.
By W 1 , 2 ( G ) , we understand the space:
W 1 , 2 ( G ) = { u L 2 ( G ) | u s W 1 , 2 ( G s ) , u b , j W 1 , 2 ( G b , j ) , 1 j m }
(where W 1 , 2 ( G s ) and W 1 , 2 ( G b , j ) are the classical Sobolev spaces) equipped with the norm:
u W 1 , 2 ( G ) = u s W 1 , 2 ( G s ) 2 + j = 1 m u b , j W 1 , 2 ( G b , j ) 2 1 / 2
If u W 1 , 2 ( G ) , then by t r u s and t r u b we denote the traces of the restrictions u s and u b on G b .
We remind the important multiplicative inequalities:
u L p ( G ) C 1 , p ( G ) u W 1 , 2 ( G ) 3 / 2 3 / p u L 2 ( G ) 3 / p 1 / 2 u W 1 , 2 ( G ) ,
t r u s L q ( G b ) C 2 , q ( G s ) u s W 1 , 2 ( G s ) 3 / 2 2 / q u s L 2 ( G s ) 2 / q 1 / 2 u s W 1 , 2 ( G s ) ,
t r u b L q ( G b ) C 2 , q ( G b ) u b W 1 , 2 ( G b ) 3 / 2 2 / q u b L 2 ( G b ) 2 / q 1 / 2 u b W 1 , 2 ( G b ) ,
which hold for all p [ 2 , 6 ] , q [ 2 , 4 ] .

3.2. Spaces of Functions on Q T and Q b , T

We set:
L p ( Q T ) = L p ( Q T ; d x d t ) , L p ( Q s , T ) = L p ( Q s , T ; d x d t ) , L p ( Q b , T ) = L p ( Q b , T ; d σ ( x ) d t ) .
Note that L p ( Q T ) = L p ( 0 , T ; L p ( G ) ) and L p ( Q b , T ) = L p ( 0 , T ; L p ( G b ) ) for 1 p < .
We introduce the space V 2 ( Q T ) = L ( 0 , T ; L 2 ( G ) ) L 2 ( 0 , T ; W 1 , 2 ( G ) ) equipped with the norm:
u V 2 ( Q T ) = u L ( 0 , T ; L 2 ( G ) ) + u L 2 ( Q T )
and the space V 2 1 , 0 ( Q T ) = V 2 ( Q T ) C ( [ 0 , T ] ; L 2 ( G ) ) .
The inequalities (10)–(12) imply the estimates:
u L r ¯ 1 ( 0 , T ; L q ¯ 1 ( G ) ) C 1 , r ¯ 1 , q ¯ 1 ( G , T ) u V 2 ( Q T ) u V 2 ( Q T ) ,
t r u s L r ¯ 2 ( 0 , T ; L q ¯ 2 ( G b ) ) C 2 , r ¯ 2 , q ¯ 2 ( G , T ) u V 2 ( Q T ) u V 2 ( Q T ) ,
t r u b L r ¯ 2 ( 0 , T ; L q ¯ 2 ( G b ) ) C 2 , r ¯ 2 , q ¯ 2 ( G , T ) u V 2 ( Q T ) u V 2 ( Q T ) ,
which hold for all exponents r ¯ 1 , q ¯ 1 , r ¯ 2 , q ¯ 2 , such that:
r ¯ 1 [ 1 , ] , q ¯ 1 [ 1 , 6 ] , r ¯ 2 [ 1 , 4 ] , q ¯ 2 [ 1 , 4 ] , 2 r ¯ 1 + 3 q ¯ 1 3 2 , 2 r ¯ 2 + 2 q ¯ 2 3 2 .
From the estimates (13), (15), it follows that if u , | u | 1 / 2 u V 2 ( Q T ) , then u L 5 ( Q T ) , t r u b L 4 ( Q b , T ) and the following estimates hold:
u L 5 ( Q T ) C 3 ( G , T ) | u | 1 / 2 u V 2 ( Q T ) 2 / 3 ,
t r u b L 4 ( Q b , T ) C 4 ( G , T ) | u | 1 / 2 u V 2 ( Q T ) 2 / 3 .
We also draw attention to the following multiplicative inequality, which follows from (11), (12):
t r u s L 2 ( Q b , T ) + t r u b L 2 ( Q b , T ) C 5 ( G , T ) u L 2 ( 0 , T ; W 1 , 2 ( G ) ) 1 / 2 u L 2 ( Q T ) 1 / 2 .
This inequality, in particular, implies that if u V 2 ( Q T ) , the sequence { u k } k = 1 V 2 ( Q T ) is bounded in V 2 ( Q T ) and u k u in L 2 ( Q T ) as k , then t r u s k t r u s , t r u b k t r u b in L 2 ( Q b , T ) .

3.3. Spaces of Functions on D s and Γ

Remind that:
D s = Ω × G s , Γ = Ω × G s , Γ = { ( ω , x ) Γ | ω · n ( x ) < 0 } , Γ + = { ( ω , x ) Γ | ω · n ( x ) > 0 } , Γ ^ ± = { ( ω , x ) Γ ± | x G ^ } , Γ b ± = { ( ω , x ) Γ ± | x G b } .
We set: L p ( D s ) = L p ( D s ; d ω d x ) . We introduce the following measures on Γ and Γ ± :
d Γ ( ω , x ) = d ω d σ ( x ) , ( ω , x ) Γ , d ^ Γ ± ( ω , x ) = | ω · n ( x ) | d ω d σ ( x ) , ( ω , x ) Γ ± .
We set:
L ^ p ( Γ ± ) = L p ( Γ ± ; d ^ Γ ± ) , L ^ p ( Γ ^ ± ) = L p ( Γ ^ ± ; d ^ Γ ± ) , L ^ p ( Γ ^ b ± ) = L p ( Γ ^ b ± ; d ^ Γ ± ) , 1 p .
Note that L ^ ( Γ ± ) = L ( Γ ± ) , L ^ ( Γ ^ ± ) = L ( Γ ^ ± ) , L ^ ( Γ b ± ) = L ( Γ b ± ) .
By the weak derivative in direction ω of a function f L 1 ( D s ) , we understand a function z L 1 ( D s ) , denoted by z = ω · f and satisfying the integral identity:
D s f ( ω , x ) ω · φ ( x ) + z ( ω , x ) φ ( x ) ψ ( ω ) d ω d x = 0 φ C 0 ( G s ) , ψ L ( Ω ) .
We denote by W p ( D s ) the Banach space of functions f L p ( D s ) possessing the weak derivative ω · f L p ( D s ) and equipped with the norm:
f W p ( D s ) = f L p ( D s ) p + ω · f L p ( D s ) p 1 / p , 1 p < , max { f L ( D s ) , ω · f L ( D s ) } , p = .
We will denote by f | Γ and f | Γ + the traces of the function f W p ( D s ) on Γ and Γ + , respectively. It is known that f | Γ ± L l o c p ( Γ ± ) . Moreover, if f W p ( D s ) and f | Γ L ^ p ( Γ ) , then f | Γ + L ^ p ( Γ + ) .
We refer to [67,68,69] for more detailed information about the properties of functions f W p ( D s ) and their traces f | Γ ± .

4. Boundary Value Problem for Radiative Transfer Equation

For almost all t ( 0 , T ) , the unknown function I ( t ) , involved in the problem (1)–(9), is a solution to the following subproblem:
ω · I + ( ϰ + s ) I = s S ( I ) + ϰ F , ( ω , x ) D s ,
I | Γ = J b , ( ω , x ) Γ b ,
I | Γ = J * , ( ω , x ) Γ ^ .
where F = 1 π h ( t r u s ( t ) ) , J b = 1 π h ( t r u b ( t ) ) .
We formulate some results on the properties of the problem (20)–(22) which follow from [67,69].
Let ( F , J b , J * ) L 1 ( G s , G b ) × L ^ 1 ( Γ ^ ) . By a solution to the problem (20)–(22), we mean a function I W 1 ( D s ) that satisfies Equation (20) almost everywhere on D s , condition (21) almost everywhere on Γ b , and condition (22) almost everywhere on Γ ^ .
Theorem 1.
If ( F , J b , J * ) L p ( G s , G b ) × L ^ p ( Γ ^ ) , where 1 p , then a solution to the problem (20)–(22) exists and is unique. Moreover, I W p ( D s ) and for 1 p < the following estimates hold:
ϰ I L p ( D s ) p + I | Γ + L ^ p ( Γ + ) p ( F , J b ) L p ( G s , G b ) p + J * L ^ p ( Γ ^ ) p ,
ϰ ω · I L p ( D s ) p 2 p ϰ + s p ( F , J b ) L p ( G s , G b ) p + J * L ^ p ( Γ ^ ) p
and for p = , the following estimates hold:
I L ( D s ) max ( F , J b ) L ( G s , G b ) , J * L ( Γ ^ ) ,
ω · I L p ( D s ) 2 ϰ + s max ( F , J b ) L ( G s , G b ) , J * L ( Γ ^ ) .
In addition:
ϰ [ I ] ± L 1 ( D s ) + [ I ] ± | Γ + L ^ 1 ( Γ + ) ( [ F ] ± , [ J b ] ± ) L 1 ( G s , G b ) + [ J * ] ± L ^ 1 ( Γ ^ )
and as a consequence, if F 0 , J b 0 , J * 0 , then I 0 .
We denote by A the resolving operator for the problem:
ω · I s + ( ϰ + s ) I s = s S ( I s ) + ϰ F , ( ω , x ) D s , I s | Γ = J b , ( ω , x ) Γ b , I s | Γ = 0 , ( ω , x ) Γ ^ .
It with a pair ( F , J b ) L p ( G s , G b ) , 1 p associates the solution I s = A [ F , J b ] W p ( D s ) . This operator is linear and continuous.
We denote by A ^ the resolving operator for the problem:
ω · I * + ( ϰ + s ) I * = s S ( I * ) , ( ω , x ) D s , I * | Γ = 0 , ( ω , x ) Γ b , I * | Γ = J * , ( ω , x ) Γ ^ .
It with a function J * L ^ p ( Γ ^ ) , 1 p associates the solution I * = A ^ [ J * ] W p ( D s ) . This operator is also linear and continuous.
We introduce the operators A Ω : L p ( G s , G b ) L p ( G s ) , A ^ Ω : L ^ p ( Γ ^ ) L p ( G s ) and B : L p ( G s , G b ) L p ( G b ) , B ^ : L ^ p ( Γ ^ ) L p ( G b ) by the formulas:
A Ω [ F , J b ] ( x ) = 1 4 π Ω A [ F , J b ] ( ω , x ) d ω , A ^ Ω [ J * ] ( x ) = 1 4 π Ω A ^ [ J * ] ( ω , x ) d ω ; B [ F , J b ] ( x ) = 1 π M + ( A [ F , J b ] | Γ + ) ( x ) , B ^ [ J * ] ( x ) = 1 π M + ( A ^ [ J * ] | Γ + ) ( x ) ,
where:
M + ( ψ ) ( x ) = Ω + ( x ) ψ ( ω , x ) ω · n ( x ) d ω , x G b , ψ L ^ 1 ( Γ b + ) .
These operators are linear and continuous. Their continuity follows from the estimates (23), (25) and the estimate:
M + ( ψ ) L p ( G b ) π 1 / p ψ L ^ p ( Γ b + ) ψ L ^ p ( Γ b + ) .
We introduce the characteristic functions 1 D s , 1 G s , 1 G b and 1 ^ of sets D s , G s , G b and Γ ^ . Note that I = 1 D s is the solution to the problem (20)–(22) with ( F , J b , J ^ ) = ( 1 G s , 1 G b , 1 ^ ) . Consequently, A [ 1 G s , 1 G b ] + A ^ [ 1 ^ ] = 1 D s and:
A Ω [ 1 G s , 1 G j ] + α ^ = 1 G s , B [ 1 G s , 1 G b ] + β ^ = 1 G b ,
where α ^ = A ^ Ω [ 1 ^ ] for x G s and β ^ = B ^ [ 1 ^ ] for x G b . Since A ^ [ 1 ^ ] 0 then α ^ 0 , β ^ 0 .
It follows from (25) that:
A Ω L ( G s , G b ) L ( G s ) 1 , B L ( G s , G b ) L ( G b ) 1 .
We also introduce the operators C : L p ( G s , G b ) L p ( G s ) , D : L p ( G s , G b ) L p ( G b ) by the formulas:
C [ F , J ] = F A Ω [ F , J ] , D [ F , J ] = J B [ F , J ] .
It follows from (29) that:
C L ( G s , G b ) L ( G s ) 2 , D L ( G s , G b ) L ( G b ) 2 .
We draw attention to the following equality proven in [65]:
4 ϰ ( C [ F , J ] , F * ) G s + ( D [ F , J ] , J * ) G b = 4 ϰ ( F , C [ F * , J * ] ) G s + ( J , D [ F * , J * ] ) G b .
It holds for all ( F , J ) L p ( G s , G b ) ( F * , J * ) L p ( G s , G b ) , 1 p .

5. Main Results: Reducing the Problem (1)–(9) to Problem P

5.1. Formulation of the Main Results

In what follows, it is assumed that the following conditions on the data are satisfied.
( A 1 ) The function λ ( x , u ) is defined on G × R , and for any u R , it is measurable with respect to x. Furthermore:
0 < λ min λ ( x , u ) λ max ( x , u ) G × R ,
where λ min and λ max are constants.
In addition, the following holder condition holds:
| λ ( x , u + v ) λ ( x , u ) | L | v | 1 / 2 ( x , u ) G × R , v [ 1 , 1 ] ,
where L is a constant.
( A 2 ) c p L ( G ) , 0 < c ̲ p c p ( x ) c ¯ p for x G , where c ̲ p and c ¯ p are constants.
( A 3 ) γ L ( G b ) , γ 0 .
( A 4 ) u 0 L 3 ( G ) , f L r ( 0 , T ; L q ( G ) ) , J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) , where:
r [ 1 , ] , q [ 9 / 7 , ] , r * [ 3 / 2 , ] , q * [ 3 / 2 , ] , 2 r + 3 q 3 , 2 r * + 2 q * 2 .
We introduce the spaces:
V ( Q T ) = { u V 2 1 , 0 ( Q T ) L 4 ( Q s , T ) | t r u b L 4 ( Q b , T ) } , V = W 1 , 2 ( G ) L ( G ) , C * [ 0 , T ] = { η C [ 0 , T ] | η ( T ) = 0 } .
By a weak solution to the problem (1)–(9), we mean a pair of functions ( u , I ) V ( Q T ) × L 1 ( 0 , T ; W 1 ( D s ) ) such that:
(1) The following identity holds:
0 T ( c p u ( t ) , v ) G d η d t ( t ) d t + 0 T a ( u ( t ) , v ) η ( t ) d t + 0 T b ( u ( t ) , I ( t ) , v ) η ( t ) d t = ( c p u 0 , v ) G · η ( 0 ) + 0 T ( f ( t ) , v ) G η ( t ) d t d t v V , η C * [ 0 , T ] ,
where:
a ( u , v ) = a 0 ( u , v ) + a 1 ( u , v ) , a 0 ( u , v ) = ( λ ( · , u ) u , v ) G = G λ ( x , u ) u · v d x , a 1 ( u , v ) = ( γ [ t r u s t r u b ] , t r v s t r v b ) G b , b ( u , I , v ) = 4 ϰ h ( u s ) ϰ Ω I d ω , v s G s + h ( t r u b ) M + ( I | Γ + ) , t r v b G b .
Here and below, by u s ( t ) , v s and u b ( t ) , v b we denote the restrictions of u ( t ) , v to G s and G b , respectively, by t r u s ( t ) , t r u b ( t ) , t r v s , t r v b we denote the traces of the functions u s ( t ) , u b ( t ) , v s , v b on G b .
(2) For almost all t ( 0 , T ) , the function I ( t ) satisfies Equation (3) almost everywhere on D s and the conditions (7), (8) almost everywhere on Γ b , Γ ^ , respectively. meaning that:
I ( t ) = I s ( t ) + I * ( t ) ,
where I s ( t ) = 1 π A [ h ( u s ( t ) ) , h ( t r u b ( t ) ) ] , I * ( t ) = A ^ [ J * ( t ) ] for almost all t ( 0 , T ) .
Remark 1.
The fulfillment of the identity (34) is equivalent to the fact that
( c p u , v ) G W 1 , 1 ( 0 , T ) for all v V ; moreover:
d d t ( c p u , v ) G + a ( u , v ) + b ( u , I , v ) = ( f , v ) G for almost all t ( 0 , T ) , ( c p u ( t ) , v ) G ( c p u 0 , v ) G as t 0 .
Remark 2.
It follows from Theorem 1 that I * L r * ( 0 , T ; W q * ( D s ) ) and:
I * L r * ( 0 , T ; W q * ( D s ) ) ϰ 1 / q * [ 1 + 2 ( ϰ + s ) ] J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) .
In addition, for almost all t ( 0 , T ) , the following estimates hold:
ϰ I * ( t ) L q * ( D s ) q * + I * ( t ) | Γ + L ^ q * ( Γ + ) q * 1 / q * J * ( t ) L ^ q * ( Γ ^ ) , 1 q * < ,
max { I * ( t ) L ( D s ) , I * ( t ) | Γ + L ( Γ + ) } J * ( t ) L ( Γ ^ ) , q * = .
In what follows, the following notations are used:
( f , J * ) r , q , r * , q * = f L r ( 0 , T ; L q ( G ) ) + J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) , ( u 0 , f , J * ) p , r , q , r * , q * = u 0 L p ( G ) + f L r ( 0 , T ; L q ( G ) ) + J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) .
The main results of this paper are the following theorems.
Theorem 2.
A weak solution to the problem (1)–(9) exists and is unique.
Theorem 3.
Let ( u 1 , I 1 ) and ( u 2 , I 2 ) be two weak solutions to the problem (1)–(9) with ( u 0 , 1 , f 1 , J * 1 ) and ( u 0 , 2 , f 2 , J * 2 ) instead of ( u 0 , f , J * ) . Then, the following estimates hold:
c p Δ u C ( [ 0 , T ] ; L 1 ( G ) ) c p Δ u 0 L 1 ( G ) + Δ f L 1 ( Q T ) + Δ J * L 1 ( 0 , T ; L ^ 1 ( Γ ^ ) ) ,
c p [ Δ u ] ± C ( [ 0 , T ] ; L 1 ( G ) ) c p [ Δ u 0 ] ± L 1 ( G ) + [ Δ f ] ± L 1 ( Q T ) + [ Δ J * ] ± L 1 ( 0 , T ; L ^ 1 ( Γ ^ ) ) ,
where Δ u = u 1 u 2 , Δ u 0 = u 0 , 1 u 0 , 2 , Δ f = f 1 f 2 , Δ J * = J * 1 J * 2 .
Theorem 4
(Comparison theorem). Let ( u 1 , I 1 ) and ( u 2 , I 2 ) be two weak solutions to the problem (1)–(9), with ( u 0 , 1 , f 1 , J * 1 ) and ( u 0 , 2 , f 2 , J * 2 ) instead of ( u 0 , f , J * ) . If u 0 , 1 u 0 , 2 , f 1 f 2 and J * 1 J * 2 , then u 1 u 2 and I 1 I 2 .
Note that the uniqueness of the solution to the problem (1)–(9) and Theorem 4 are direct consequences of Theorem 3.
Consider that u and I are interpreted as the absolute temperature and the radiation intensity. Therefore, it is important to show that u and I are nonnegative under some natural assumptions on the data. It is clear that ( u , I ) = ( 0 , 0 ) is a solution to the problem (1)–(9) with u 0 = 0 , f = 0 and J * = 0 . Thus, Theorem 4 implies the following result.
Corollary 1.
Let ( u , I ) be a weak solution to the problem (1)–(9). If u 0 0 , f 0 and J * 0 , then u 0 and I 0 .
The following three theorems show that an increase in summability exponents of f and g leads to improved properties of a weak solution.
Theorem 5.
Let ( u , I ) be a weak solution to the problem (1)–(9). If:
u 0 L p ( G ) , f L r ( 0 , T ; L q ( G ) ) , J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) ,
p [ 3 , ) , r [ 1 , ] , q 3 2 + 1 / p , , r * 2 1 + 1 / p , , q * 2 1 + 1 / p , ,
2 / r + 3 / q 2 + 3 / p , 2 / r * + 2 / q * 1 + 3 / p ,
then | u | γ 1 u V 2 ( Q T ) for all γ [ 1 , p / 2 ] , I s L r s ( 0 , T ; W q s ( D ) ) for all r s [ 1 , p / 2 ] , q s [ 1 , p / 2 ] such that 1 / r s + 1 / q s 6 / p and:
| u | γ 1 u V 2 ( Q T ) 1 / γ C ( u 0 , f , g ) p , r , q , r * , q * ,
I s L r s ( 0 , T ; W q s ( D s ) ) C ( u 0 , f , g ) p , r , q , r * , q * 4 .
In addition, u C ( [ 0 , T ] ; L s ( G ) ) for all s [ 1 , p ) .
By C (with or without indices), we denote various positive constants that may depend on G, T, c ̲ p , c ¯ p , λ min , λ max , σ 0 , ϰ , s and p , r , q , r * , q * .
Theorem 6.
Let ( u , I ) be a weak solution to the problem (1)–(9). If:
e β 0 | u 0 | L 1 ( G ) for some β 0 > 0 , f L r ( 0 , T ; L q ( G ) ) , J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) ,
r [ 1 , ] , q [ 3 / 2 , ] , r * [ 2 , ) , q * [ 2 , ) , 2 r + 3 q 2 , 2 r * + 2 q * 1 ,
then there exists a constant β ( 0 , β 0 / 2 ] such that e β | u | V 2 ( Q T ) and:
e β | u | V 2 ( Q T ) C ( e β | u 0 | L 2 ( G ) + 1 ) ,
where C depends on ( f , J * ) r , q , r * , q * , but does not depend on u 0 .
In addition, u C ( [ 0 , T ] ; L s ( G ) ) for all s [ 1 , ) and I s L r s ( 0 , T ; W q s ( D s ) ) for all r s , q s [ 1 , ) .
Theorem 7.
Let ( u , I ) be a weak solution to the problem (1)–(9). If:
u 0 L ( G ) , f L r ( 0 , T ; L q ( G ) ) , J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) ,
r ( 1 , ] , q ( 3 / 2 , ] , r * ( 2 , ] , q * ( 2 , ] , 2 r + 3 q < 2 , 2 r * + 2 q * < 1 ,
then u L ( Q T ) and I s ( t ) W ( D s ) for almost all t ( 0 , T ) . Moreover:
u L ( Q T ) C ( u 0 , f , J * ) , r , q , r * , q * ,
ess sup t ( 0 , T ) I ( t ) W ( D s ) C ( u 0 , f , g ) , r , q , r * , q * 4 .

5.2. Reducing the Problem (1)–(9) to Problem P

Let ( u , I ) be a weak solution to the problem (1)–(9) and I * = A [ J * ] . We put:
f * = ϰ Ω I * d ω = 4 π ϰ A ^ Ω [ J * ] for ( x , t ) Q s , T , g * = M + ( I * | Γ + ) = π B ^ [ J * ] for ( x , t ) Q b , T .
It follows from (37), (38) that f * L r * ( 0 , T ; L q * ( G s ) ) , g * L r * ( 0 , T ; L q * ( G b ) ) and
f * L r * ( 0 , T ; L q * ( G s ) ) + g * L r * ( 0 , T ; L q * ( G b ) ) ( ( 4 π ϰ ) 1 / q * + π 1 / q * ) J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) .
It follows from (35) that:
4 ϰ h ( u s ) ϰ Ω I d ω = 4 ϰ h ( u s ) 4 ϰ A Ω [ h ( u s ) , h ( t r u b ) ] f * = 4 ϰ C [ h ( u s ) , h ( t r u b ) ] f * , h ( t r u b ) M + ( I | Γ + ) = h ( t r u b ) 1 π M + ( A [ h ( u s ) , h ( t r u b ) ] | Γ b + ) g * = h ( t r u b ) B [ h ( u s ) , h ( t r u b ) ] g * = D [ h ( u s ) , h ( t r u b ) ] g * .
Using these formulas, we exclude the function I from the problem (1)–(9) and arrive at the problem:
c p u s t div ( λ ( x , u s ) u s ) + 4 ϰ C [ h ( u s ) , h ( t r u b ) ] = f + f * , ( x , t ) Q s , T , c p u b t div ( λ ( x , u b ) u b ) = f , ( x , t ) Q b , T , λ ( x , u s ) u s n = 0 , ( x , t ) Q ^ T , λ ( x , u s ) u s n + γ ( t r u s t r u b ) = 0 , ( x , t ) Q b , T , λ ( x , u b ) u b n + D [ h ( u s ) , h ( t r u b ) ] + γ ( t r u b t r u s ) = g * , ( x , t ) Q b , T , u | t = 0 = u 0 , x G ,
in which only one function u is unknown. This problem will be called Problem P .
Remind that V ( Q T ) = { u V 2 1 , 0 ( Q T ) L 4 ( Q s , T ) | t r u b L 4 ( Q b , T ) } . Therefore, it follows from u V ( Q T ) and the boundedness of the operators C : L 1 ( G s , G b ) L 1 ( G s ) , D : L 1 ( G s , G b ) L 1 ( G b ) that C [ h ( u s ) , h ( t r u b ) ] L 1 ( Q s , T ) , D [ h ( u s ) , h ( t r u b ) ] L 1 ( Q b , T ) .
By a weak solution to Problem P , we mean a function u V ( Q T ) satisfying the identity:
0 T ( c p u ( t ) , v ) G d η d t ( t ) d t + 0 T A ( u ( t ) , v ) η ( t ) d t = ( c p u 0 , v ) G · η ( 0 ) + 0 T F ( t ) , v η ( t ) d t v V , η C * [ 0 , T ] ,
where:
A ( u , v ) = a ( u , v ) + b ( u , v ) , b ( u , v ) = ( 4 ϰ C [ h ( u s ) ) , h ( t r u b ) ] , v s ) G s + ( D [ h ( u s ) , h ( t r u b ) ] , t r v b ) G b , F ( t ) , v = ( f ( t ) , v ) G + ( f * ( t ) , v s ) G s + ( g * ( t ) , t r v b ) G b .
Remark 3.
Due to the equality (31) instead of the formula (55), it is possible to use an alternative formula:
b ( u , v ) = ( h ( u s ) , 4 ϰ C [ v s , t r v b ] ) G s + ( h ( t r u b ) , D [ v s , t r v b ] ) G b .
Remark 4.
It is easy to see that if ( u , I ) V ( Q T ) × L 1 ( 0 , T ; W 1 ( D s ) is a weak solution to the problem (1)–(9), then u is a weak solution to Problem P . On the other hand, if u V ( Q T ) is a weak solution to Problem P , then defining I by the formula (35) for almost all t ( 0 , T ) , we obtain the pair ( u , I ) V ( Q T ) × L 1 ( 0 , T ; W 1 ( D s ) ) that is a weak solution to the problem (1)–(9). The fact that I L 1 ( 0 , T ; W 1 ( D s ) ) follows from the continuity of the operator A : L 1 ( G s , G b ) W 1 ( D s ) and properties h ( u s ) L 1 ( Q s , T ) , h ( t r u b ) L 1 ( Q b , T ) that the function u V ( Q T ) possesses.

6. Auxiliaries

6.1. Forms D ^ ( U i , U I i , V i , V I i ) , D ( U i , U I i , V i , V I i ) and Some of Their Properties

We set:
d ^ ( u I , u I I , v I , v I I ) = ( 4 ϰ C [ u I , u I I ] , v I ) G s + ( D [ u I , u I I ] , v I I ) G b , d ( u I , u I I , v I , v I I ) = ( 4 ϰ A Ω [ 1 G s , 1 G b ] u I , v I ) G s ( 4 ϰ A Ω [ u I , u I I ] , v I ) G s + B [ 1 G s , 1 G b ] u I I , v I I G b B [ u I , u I I ] , v I I G b ,
where ( u I , u I I , v I , v I I ) L 1 ( G s , G b ) × L ( G s , G b ) or ( u I , u I I , v I , v I I ) L 2 ( G s , G b ) × L 2 ( G s , G b ) .
It follows from (28) that:
d ^ ( u I , u I I , v I , v I I ) = d ( u I , u I I , v I , v I I ) + 4 ϰ α ^ u I , v I G s + ( β ^ u I I , v I I ) G b .
The following three statements are proven in [65].
Lemma 1.
Assume that ( u I , u I I ) L 1 ( G s , G b ) , ( v I , v I I ) L ( G s , G b ) , { ( u I N , u I I N ) } N = 1   L 1 ( G s , G b ) , { ( v I N , v I I N ) } N = 1 L ( G s , G b ) . Assume also that ( u I N , u I I N ) ( u I , u I I ) in L 1 ( G s , G b ) , v I N v I almost everywere on G s , v I I N v I I almost everywhere on G b as N and sup N 1 ( v I N , v I I N ) L ( G s , G b ) < . Then:
lim N d ^ ( u I N , u I I N , v I N , v I I N ) = d ^ ( u I , u I I , v I , v I I ) .
Lemma 2.
Assume that ( u I , u I I ) L 2 ( G s , G b ) , { ( u I N , u I I N ) } N = 1 L 2 ( G s , G b ) and ( u I N , u I I N ) ( u I , u I I ) in L 2 ( G s , G b ) as N . Then:
lim N d ^ ( u I N , u I I N , u I N , u I I N ) = d ^ ( u I , u I I , u I , u I I ) .
Let { E s i } i = 1 N be a system of measurable pairwise disjoint subsets of G s , such that G s = i = 1 N E s i and let { S b } = 1 N be a system of measurable pairwise disjoint subsets of G b , such that G b = = 1 N S b . We denote by 1 s i and 1 b the characteristic functions of sets E s i and S b , respectively. We set:
α k i = 4 ϰ A Ω [ 1 s k , 0 ] , 1 s i G s , β k = B [ 0 , 1 s ] , 1 s k G b , δ i = 4 ϰ A Ω [ 0 , 1 b ] , 1 s i G s .
Note that α k i 0 , β k 0 , δ i 0 .
Lemma 3.
Let:
u I N ( x ) = i = 1 N u s i 1 s i ( x ) , u I I N ( x ) = = 1 N u b 1 b ( x )
be simple functions defined on G s , G b , respectively, and let:
v I N ( x ) = i = 1 N v s i 1 s i ( x ) , v I I N ( x ) = = 1 N v b 1 b ( x )
be other simple functions defined on G s , G b , respectively. Then:
d ( u I N , u I I N , v I N , v I I N ) = 1 2 i = 1 N k = 1 N α k i [ u s i u s k ] [ v s i v s k ] + 1 2 = 1 N k = 1 N β k [ u b u b k ] [ v b v b k ] + i = 1 N = 1 N δ i [ u s i u b ] [ v s i v b ] .

6.2. Forms b ( u , v ) , b [ n ] ( u , v ) and Some of Their Properties

Consider that:
b ( u , v ) = ( 4 ϰ C [ h ( u s ) , h ( t r u b ) ] , v s ) G s + ( D [ h ( u s ) , h ( t r u b ) ] , t r v b ) G b , ( u , v ) W 1 , 2 ( G ) × V .
We also set:
b [ n ] ( u , v ) = ( 4 ϰ C [ h [ n ] ( u s ) , h [ n ] ( t r u b ) ] , v s ) G s + ( D [ h [ n ] ( u s ) , h [ n ] ( t r u b ) ] , t r v b ) G b , ( u , v ) W 1 , 2 ( G ) × W 1 , 2 ( G ) ,
where h [ n ] = min { max { h ( u ) , n } , n } .
Note that:
b ( u , v ) = d ^ ( h ( u s ) , h ( t r u b ) , v s , t r v b ) = d ( h ( u s ) , h ( t r u b ) , v s , t r v b ) + 4 ϰ α ^ h ( u s ) , v s G s + ( β ^ h ( t r u b ) , t r v s ) G b , b [ n ] ( u , v ) = d ^ ( h [ n ] ( u s ) , h [ n ] ( t r u b ) , v s , t r v b ) = d ( h [ n ] ( u s ) , h [ n ] ( t r u b ) , v s , t r v b ) + 4 ϰ α ^ h [ n ] ( u s ) , v s G s + ( β ^ h [ n ] ( t r u b ) , t r v s ) G b .
Lemma 4.
For all n 1 , the following inequality holds:
b [ n ] ( u , u ) 0 u W 1 , 2 ( G ) .
Proof of Lemma 4.
We construct sequences of simple functions { u I N } N = 1 and { u I I N } N = 1 of the forms (58) such that:
| u I N ( x ) | | u s ( x ) | and lim N u I N ( x ) = u s ( x ) for almost all x G s , | u I I N ( x ) | | t r u b ( x ) | and lim N u I I N ( x ) = t r u b ( x ) for almost all x G b .
It follows from Lemma 3 and the monotonicity of the function h that:
d ^ ( h [ n ] ( u I N ) , h [ n ] ( u I I N ) , u I N , u I I N ) = 1 2 i = 1 N k = 1 N α k i [ h [ n ] ( u s i N ) h [ n ] ( u s k N ) ] [ u s i N u s k N ] + 1 2 = 1 N k = 1 N β k [ h [ n ] ( u b N ) h [ n ] ( u b k N ) ] [ u b N u b k N ] + i = 1 N = 1 N δ i [ h [ n ] ( u s i N ) h [ n ] ( u b N ) ] [ u s i N u b N ] + ( 4 ϰ α ^ h [ n ] ( u I N ) , u I N ) G s + ( β ^ h [ n ] ( u I I N ) , u I I N ) G b 0 .
It is clear that ( h [ n ] ( u I N ) , h [ n ] ( u I I N ) ) ( h [ n ] ( u s ) , h [ n ] ( t r u b ) ) in L 2 ( G s , G b ) and ( u I N , u I I N ) ( u s , t r u b ) in L 2 ( G s , G b ) . Therefore, by Lemma 2:
lim N d ^ ( h [ n ] ( u I N ) , h [ n ] ( u I I N ) , u I N , u I I N ) d ^ ( h [ n ] ( u s ) , h [ n ] ( t r u b ) , u s , t r u b ) = b [ n ] ( u , u ) .
Passing in the inequality (61) to the limit as N , we arrive at the inequality (60).  ☐
Lemma 5.
Assume that w C ( R ) , w be a non-decreasing function such that w ( 0 ) = 0 . Then, for all n 1 , M > 0 , the following inequalities hold:
b ( u , w ( u [ M ] ) ) 0 u W 1 , 2 ( G ) ,
b [ n ] ( u , w ( u [ M ] ) ) 0 u W 1 , 2 ( G ) .
Proof of Lemma 5.
Let { u I N } N = 1 and { u I I N } N = 1 be the same sequences of simple functions as in the proof of the previous lemma.
It follows from Lemma 3 and the monotonicity of the function h that:
d ^ h ( u I N ) , h ( u I I N ) , w ( ( u I N ) [ M ] ) , w ( ( u I I N ) [ M ] ) = 1 2 i = 1 N k = 1 N α k i h ( u s i N ) h ( u s k N ) w ( ( u s i N ) [ M ] ) w ( ( u s k N ) [ M ] ) + 1 2 = 1 N k = 1 N β k h ( u b N ) h ( u b k N ) w ( ( u b N ) [ M ] ) w ( ( u b k N ) [ M ] ) + i = 1 N = 1 N δ i [ h ( u s i N ) h ( u b N ) ] [ w ( ( u s i N ) [ M ] ) w ( ( u b N ) [ M ] ) ] + ( 4 ϰ α ^ h ( u I N ) , w ( ( u I N ) [ M ] ) ) G s + ( β ^ h ( u I I N ) , w ( ( u I I N ) [ M ] ) ) G b 0 .
It is clear that ( h ( u I N ) , h ( u I I N ) ) ( h ( u s ) , h ( t r u b ) ) in L 1 ( G s , G b ) as N ; w ( ( u I N ) [ M ] ) w ( u s [ M ] ) almost everywhere on G s , w ( ( u I I N ) [ M ] ) w ( t r u b [ M ] ) almost everywhere on G b as N ; in addition, w ( ( u I N ) [ M ] ) L ( G s ) C M , w ( ( u I I N ) [ M ] ) L ( G b ) C M , where C M = max { w ( M ) , | w ( M ) | } . Therefore, by Lemma 1:
lim N d ^ h ( u I N ) , h ( u I I N ) , w ( ( u I N ) [ M ] ) , w ( ( u I I N ) [ M ] ) = d ^ h ( u s ) , h ( t r u b ) , w ( u s [ M ] ) , w ( t r u b [ M ] ) = b u , w ( u [ M ] ) .
Passing in (64) to the limit as N , we arrive at the inequality (62). The proof of the inequality (63) is quite the same.  ☐

7. An Auxiliary Problem P [ N ] and Its Solvability

Consider an auxiliary Problem P [ n ] , which differs from Problem P only in that in its formulation, the function h ( u ) is replaced by h [ n ] = min { max { h ( u ) , n } , n } , where 1 n is a natural parameter.
By a weak solution to Problem P [ n ] , we mean a function u V 2 1 , 0 ( Q T ) satisfying the identity:
0 T ( c p u ( t ) , v ) G d d t η ( t ) d t + 0 T A [ n ] ( u ( t ) , v ) η ( t ) d t = ( c p u 0 , v ) G · η ( 0 ) + 0 T F ( t ) , v η ( t ) d t v W 1 , 2 ( G ) , η C * [ 0 , T ] ,
where A [ n ] ( u , v ) = a ( u , v ) + b [ n ] ( u , v ) and b [ n ] ( u , v ) is given by the formula (59).
Theorem 8.
Assume that u 0 L 2 ( G ) , f L r ( 0 , T ; L q ( G ) ) , J * L r * ( 0 , T ; L q * ( Γ ^ ) ) , where r ( 1 , ] , q [ 6 / 5 , ] , r * [ 4 / 3 , ] , q * [ 4 / 3 , ] , 2 / r + 3 / q 7 / 2 , 2 / r * + 2 / q * 5 / 2 . Then, a weak solution to Problem P [ n ] exists.
Proof of Theorem 8.
Let { e } = 1 be a basis in W 1 , 2 ( G ) that is orthonormal in L 2 ( G ) with weight c p .
We set V k = span { e 1 , , e k } , k 1 and will seek an approximate solution to Problem P [ n ] in the form u ( k ) ( t ) = = 1 k d ( k ) ( t ) e , determining the coefficients d ( k ) from the Galerkin method:
c p d d t u ( k ) ( t ) , v G + A [ n ] ( u ( k ) ( t ) , v ) = F ( t ) , v , t ( 0 , T ) v V k , u ( k ) ( 0 ) = u 0 , k = = 1 k ( c p u 0 , e ) G e .
Note that u 0 , k u 0 in L 2 ( G ) as k , moreover c p 1 / 2 u 0 , k L 2 ( G ) c p 1 / 2 u 0 L 2 ( G ) .
The Caratheodory theorem implies the existence of a time-local solution u ( k ) . It is defined on the whole interval ( 0 , T ) by virtue of the global to time a priori estimate:
u ( k ) V 2 ( Q T ) C 1 ( u 0 , f , J * ) 2 , r , q , r * , q * .
To obtain this estimate, we substitute v = u ( k ) ( t ) in (66) and use the inequalities:
λ min u ( k ) L 2 ( G ) 2 a ( u ( k ) , u ( k ) ) , 0 b [ n ] ( u ( k ) , u ( k ) ) ,
which follow from the condition (32) and Lemma 4, and arrive at the inequality:
1 2 d d t c p 1 / 2 u ( k ) ( t ) L 2 ( G ) 2 + λ min u ( k ) ( t ) L 2 ( G ) 2 F ( t ) , u ( k ) ( t ) , t ( 0 , T ) .
Integrating it, we deduce the inequality:
1 2 c p 1 / 2 u ( k ) ( t ) L 2 ( G ) 2 + λ min u ( k ) L 2 ( Q t ) 2 1 2 c p 1 / 2 u 0 , k L 2 ( G ) 2 + 0 t F ( t ) , u ( k ) ( t ) d t 1 2 c p 1 / 2 u 0 L 2 ( G ) 2 + f L r ( 0 , T ; L q ( G ) ) u ( k ) L r ( 0 , T ; L q ( G ) ) + f * L r * ( 0 , T ; L q * ( G ) u s ( k ) L r * ( 0 , T ; L q * ( G s ) ) + g * L r * ( 0 , T ; L q * ( G b ) ) t r u b ( k ) L r * ( 0 , T ; L q * ( G b ) ) , t ( 0 , T )
Applying the inequality (13) with r , q and r * , q * instead of r ¯ 1 , q ¯ 1 , the inequality (14) with r * , q * instead of r ¯ 2 , q ¯ 2 and using the inequality (53), we arrive at the estimate (67).
Let us derive one more estimate. Since h [ n ] ( u s ( k ) ) , h [ n ] ( t r u b ( k ) ) L ( G s , G b ) n , then it follows from (30) that:
C [ h [ n ] ( u s ( k ) ) , h [ n ] ( t r u b ( k ) ) ] L ( G s ) 2 n , D [ h [ n ] ( u s ( k ) ) , h [ n ] ( t r u b ( k ) ) ] L ( G b ) 2 n .
Consequently:
| b [ n ] ( u ( k ) , v ) | 2 n 4 ϰ v s L 1 ( G s ) + t r v b L 1 ( G b ) .
Integrating (66) over ( t , t + τ ) , where 0 < τ < T , and taking into account the estimate (67), we have:
( c p Δ ( τ ) u ( k ) ( t ) , v ) G = t t + τ A [ n ] ( u ( k ) ( t ) , v ) + F ( t ) , v d t τ 1 / 2 λ max u ( k ) L 2 ( Q T ) v L 2 ( G ) + τ 1 / 2 γ L ( G b ) t r u b ( k ) t r u s ( k ) L 2 ( Q b , T ) t r v b t r v s L 2 ( G b ) + 2 τ n 4 ϰ v s L 1 ( G s ) + t r v b L 1 ( G b ) + τ 1 1 / r f L r ( 0 , T ; L q ( G ) ) v L q ( G ) + τ 1 1 / r * f * L r * ( 0 , T ; L q * ( G ) ) v s L q * ( G ) + τ 1 1 / r * g * L r * ( 0 , T ; L q * ( G b ) ) t r v b L q * ( G b ) ,
where Δ ( τ ) u ( k ) ( t ) = u ( k ) ( t + τ ) u ( k ) ( t ) .
Taking v = Δ ( τ ) u ( k ) ( t ) , integrating the resulting inequality over t from 0 to T τ and using the inequalities (53), (67), we obtain:
c ̲ p Δ ( τ ) u ( k ) L 2 ( Q T τ ) 2 2 τ 1 / 2 T 1 / 2 λ max u ( k ) L 2 ( Q T ) 2 + γ L ( G b ) t r u b ( k ) t r u s ( k ) L 2 ( Q b , T ) 2 + 4 τ n 4 ϰ u s ( k ) L 1 ( Q s , T ) + t r u b ( k ) L 1 ( Q b , T ) + 2 τ 1 1 / r f L r ( 0 , T ; L q ( G ) ) u ( k ) L 1 ( 0 , T ; L q ( G ) ) + 2 τ 1 1 / r * f * L r * ( 0 , T ; L q * ( G s ) ) u s ( k ) L 1 ( 0 , T ; L q * ( G s ) ) + 2 τ 1 1 / r * | g L r * ( 0 , T ; L q * ( G b ) ) t r u b ( k ) L 1 ( 0 , T ; L q * ( G b ) ) ] τ ν C 2 ( u 0 , f , J * ) 2 , r , g , r * , q * + n 2 ,
where ν = min { 1 / 2 , 1 1 / r , 1 1 / r * } . Thus:
Δ ( τ ) u ( k ) L 2 ( Q T τ ) C 3 τ ν / 2 ( u 0 , f , J * ) 2 , r , g , r * , q * + n .
It follows from (67) that there exist a function u V 2 ( Q T ) and the subsequence { u ( k ) } = 1 such that u ( k ) u weakly in L 2 ( 0 , T ; W 1 , 2 ( G ) ) and weakly stars in L ( 0 , T ; L 2 ( G ) ) as k .
By virtue of the Riesz precompactness criterion for L 2 ( Q T ) , the estimates (67) and (68) allow us to select a subsequence, such that u ( k ) u strongly in L 2 ( Q T ) and almost everywhere on Q T .
It is clear that λ ( · , u ( k ) ) u ( k ) λ ( · , u ) u weakly in L 2 ( Q T ) and therefore
a 0 ( u ( k ) , v ) a 0 ( u , v ) weakly in L 1 ( 0 , T ) for all v W 1 , 2 ( G ) .
From the estimate (19) applied to u ( k ) u , it follows that t r u s ( k ) t r u s and t r u b ( k ) t r u b in L 2 ( Q b , T ) as . So a 1 ( u ( k ) , v ) a 1 ( u , v ) in L 1 ( 0 , T ) .
It is also easy to see that h [ n ] ( u s ( k ) ) h [ n ] ( u s ) in L 1 ( Q s , T ) and h [ n ] ( t r u b ( k ) ) h [ n ] ( t r u b ) in L 1 ( Q b , T ) . Using the formula (59), we have:
b [ n ] ( u ( k ) , v ) = 4 ϰ h [ n ] ( u s ( k ) ) , C [ v s , t r v b ] G s + h [ n ] ( t r u b ( k ) ) , D [ v s , t r v b ] G b 4 ϰ ( h [ n ] ( u s ) , C [ v s , t r v b ] ) G s + ( h [ n ] ( t r u b ) , D [ v s , t r v b ] ) G b = b [ n ] ( u , v )
in L 1 ( 0 , T ) for all v W 1 , 2 ( G ) .
Multiplying (66) on η ( t ) , where η C * [ 0 , T ] , and integrating the result over t from 0 to T, we have:
0 T ( c p u ( k ) ( t ) , v ) G d d t η ( t ) d t + 0 T A [ n ] ( u ( k ) ( t ) , v ) η ( t ) d t = ( c p u 0 , k , v ) G · η ( 0 ) + 0 T F ( t ) , v η ( t ) d t .
Passing to the limit at k = k , we establish the validity of the identity (65) for an arbitrary function v k = 1 V k . Since the set k = 1 V k is dense everywhere in W 1 , 2 ( G ) , then the identity (65) holds for all v W 1 , 2 ( G ) .
Since the function u satisfies this identity, it follows (see, for example [41], Lemma 4.1) that u C ( [ 0 , T ] ; L 2 ( G ) ) . Thus, u V 2 1 , 0 ( Q T ) .  ☐

8. Estimates for Weak Solutions to Problems P and P [ N ]

We need the following statement, following from [41], in Lemma 4.4.
Lemma 6.
Assume that a function u V 2 1 , 0 ( Q T ) satisfies the identity:
0 T ( c p u ( t ) , v ) G d d t η ( t ) d t = 0 T F ^ ( t ) , v η ( t ) d t v V , η C 0 [ 0 , T ] ,
where:
F ^ ( t ) , v = ( F ( t ) , v ) G + ( f ^ ( t ) , v ) G + ( g s ( t ) , t r v s ) G b + ( g b ( t ) , t r v b ) G b , F L 2 ( 0 , T ; ( L 2 ( G ) ) 3 ) , f ^ L 1 ( Q T ) , g s , g b L 1 ( Q b , T ) .
Assume also that w C 1 ( R ) , w 0 , w ( 0 ) = 0 and W ( M ) ( u ) = 0 u w ( s [ M ] ) d s , where M > 0 .
Then:
c p W ( M ) ( u ( t ) ) L 1 ( G ) = c p W ( M ) ( u ( 0 ) ) L 1 ( G ) + 0 t F ^ ( t ) , w ( u [ M ] ( t ) ) d t , t [ 0 , T ] ,
c p W ( M ) ( [ u ] + ( t ) ) L 1 ( G ) = c p W ( M ) ( [ u ] + ( 0 ) ) L 1 ( G ) + 0 t F ^ ( t ) , w ( u [ 0 , M ] ( t ) ) d t , t [ 0 , T ] .
Lemma 7.
Let u be a weak solution to Problem P or to Problem P [ n ] .
Suppose that U C 1 ( R ) , U 0 , w ( u ) = 0 u ( U ( s ) ) 2 d s , W ( u ) = 0 u w ( s ) d s . Then:
c p W ( u [ M ] ( t ) ) L 1 ( G ) + λ min U ( u [ M ] ) L 2 ( Q t ) 2 c p W ( u 0 ) L 1 ( G ) + 0 t F ( t ) , w ( u [ M ] ( t ) ) d t , t [ 0 , T ] .
Proof of Lemma 7.
A weak solution to Problem P satisfies the identity (69) with:
F ˜ ( t ) , v = A ( u ( t ) , v ) + F ( t ) , v .
Using Lemma 6, we arrive at the equality:
c p W ( M ) ( u ( t ) ) L 1 ( G ) + 0 t A u ( t ) , w ( u [ M ] ( t ) ) d t = c p W ( M ) ( u 0 ) L 1 ( G ) + 0 t F ( t ) , w ( u [ M ] ( t ) ) d t .
Note that:
a 0 ( u , w ( u [ M ] ) ) = ( λ ( · , u ) u , w ( u [ M ] ) ) G = ( λ ( · , u ) u , [ U ( u [ M ] ) ] 2 u [ M ] ) G = ( λ ( · , u ) U ( u [ M ] ) , U ( u [ M ] ) ) G λ min U ( u [ M ] ) L 2 ( G ) 2 , a 1 ( u , w ( u [ M ] ) ) = ( γ [ t r u b ( t ) t r u s ( t ) ] , w ( t r u b [ M ] ) w ( t r u s [ M ] ) ) G b 0 .
Using these inequalities, the inequality (62) and the estimates:
W ( u [ M ] ) W ( M ) ( u ) W ( u ) ,
we arrive from (73) at (72).
The inequality (72) for the weak solution to Problem P [ n ] is established in the same way. The only difference is that the inequality (63) is used instead of the inequality (62).  ☐
Theorem 9.
Let u be a weak solution to Problem P or to Problem P [ n ] . If the assumptions (41)–(43) are satisfied, then | u | γ 1 u V 2 ( Q T ) for all γ [ 1 , p / 2 ] ; moreover, the estimate (44) holds. In addition, u C ( [ 0 , T ] ; L s ( G ) ) for all s [ 1 , p ) .
Proof of Theorem 9.
Let γ [ 1 , p / 2 ] . We set:
U γ ( u ) = | u | γ 1 u , w γ ( u ) = 0 u ( U γ ( s ) ) 2 d s = γ 2 2 γ 1 | u | 2 γ 2 u , W γ ( u ) = 0 u w γ ( s ) d s = γ 2 ( 2 γ 1 ) | u | 2 γ .
Since:
1 4 [ U γ ( u ) ] 2 W γ ( u ) 1 2 [ U γ ( u ) ] 2 , | w γ ( u ) | γ | U γ ( u ) | 2 1 / γ ,
the inequality (72) implies the inequality:
1 4 U γ ( u [ M ] ( t ) ) L 2 ( G ) 2 + λ min U γ ( u [ M ] ) L 2 ( Q t ) 2 1 2 U γ ( u 0 ) L 2 ( G ) 2 + γ f L r ( 0 , T ; L q ( G ) ) U γ ( u [ M ] ) L r ¯ ( 0 , T ; L q ¯ ( G ) ) 2 1 / γ + γ f * L r * ( 0 , T ; L q * ( G s ) ) U γ ( u s [ M ] ) L r ¯ * ( 0 , T ; L q ¯ * ( G s ) ) 2 1 / γ + γ g * L r * ( 0 , T ; L q * ( G b ) ) U γ ( t r u b [ M ] ) L r ¯ * ( 0 , T ; L q ¯ * ( G b ) ) 2 1 / γ ,
where r ¯ = ( 2 1 / γ ) r , q ¯ = ( 2 1 / γ ) q , r ¯ * = ( 2 1 / γ ) r * , q ¯ * = ( 2 1 / γ ) q * .
As a consequence, we have:
U γ ( u [ M ] ) V 2 ( Q T ) 2 C 1 [ u 0 L 2 γ ( G ) 2 γ + γ f L r ( 0 , T ; L q ( G ) ) U γ ( u [ M ] ) L r ¯ ( 0 , T ; L q ¯ ( G ) ) 2 1 / γ + γ J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) U γ ( u s [ M ] ) L r ¯ * ( 0 , T ; L q ¯ * ( G s ) ) 2 1 / γ + U γ ( t r u b [ M ] ) L r ¯ * ( 0 , T ; L q ¯ * ( G b ) ) 2 1 / γ ] .
Using (13) and (14), we derive from (74) the estimate:
| u [ M ] | γ 1 u [ M ] V 2 ( Q T ) 1 / γ C ( u 0 , f , J * ) p , r , q , r * , q *
with a constant C that does not depend on M. Since u [ M ] u in L 2 ( 0 , T ; W 1 , 2 ( G ) ) and u [ M ] ( t ) u ( t ) in L 2 ( G ) for all t [ 0 , T ] as M , the estimate (75) implies that | u | γ 1 u V 2 ( Q T ) and the estimate (44) holds.
Since | u | p / 2 1 u V 2 ( Q T ) then u L ( 0 , T ; L p ( G ) ) . Taking into account that u C ( [ 0 , T ] ; L 2 ( G ) ) , we come to the conclusion that u C ( [ 0 , T ] ; L s ( G ) ) for all s [ 1 , p ) .  ☐
Theorem 10.
Let u be a weak solution to Problem P or to Problem P [ n ] . If the assumptions (49), (50) are satisfied, then u L ( Q T ) and the estimate (51) holds.
Proof of Theorem 10.
We put A = ( u 0 , f , J * ) , r , q , r * , q * . If A = 0 , then it follows from (44) that u = 0 and the estimate (51) holds.
Let A > 0 . We divide both sides of the inequality (74) by A 2 γ . Taking into account that γ 1 and (74) holds for all M > 0 , we obtain:
U γ ( u ¯ [ M ] ) V 2 ( Q T ) 2 C 1 γ [ meas G + U γ ( u ¯ [ M ] ) L r ¯ ( 0 , T ; L q ¯ ( G ) ) 2 1 / γ + U γ ( u ¯ [ M ] ) L r ¯ * ( 0 , T ; L q ¯ * ( G s ) ) 2 1 / γ + U γ ( t r u ¯ b [ M ] ) L r ¯ * ( 0 , T ; L q ¯ * ( G b ) ) 2 1 / γ ] C 2 γ [ U γ ( u ¯ [ M ] ) L 2 r ( 0 , T ; L 2 q ( G ) ) 2 + U γ ( u ¯ [ M ] ) L 2 r * ( 0 , T ; L 2 q * ( G ) ) 2 + U γ ( t r u ¯ b [ M ] ) L 2 r * ( 0 , T ; L 2 q * ( G b ) ) 2 + 1 ] ,
where u ¯ = u / A .
By condition (50), we can assume that 2 / r + 3 / q 2 3 δ , 2 / r * + 2 / q * 1 3 δ with some δ ( 0 , 1 ) .
Setting γ = γ k = ( 1 + δ ) k , k 0 , we obtain the inequality:
U γ k ( u ¯ [ M ] ) V 2 ( Q T ) 2 C 2 ( 1 + δ ) k [ U γ k 1 ( u ¯ [ M ] ) L r ˜ ( 0 , T ; L q ˜ ( G ) ) 2 ( 1 + δ ) + U γ k 1 ( u ¯ [ M ] ) L r ˜ * ( 0 , T ; L q ˜ * ( G ) ) 2 ( 1 + δ ) + U γ k 1 ( t r u ¯ b [ M ] ) L r ˜ * ( 0 , T ; L q ˜ * ( G ^ ) ) 2 ( 1 + δ ) + 1 ] , k 1 ,
where r ˜ = 2 ( 1 + δ ) r , q ˜ = 2 ( 1 + δ ) q , r ˜ * = 2 ( 1 + δ ) r * , q ˜ * = 2 ( 1 + δ ) q * . It is easy to check that r ˜ , q ˜ and r ˜ * , q ˜ * satisfy (16) in the role of r ¯ 1 , q ¯ 1 and r ˜ * , q ˜ * satisfy (16) in the role of r ¯ 2 , q ¯ 2 . Using (13) and (15), we arrive at the inequality:
U γ k ( u ¯ [ M ] ) V 2 ( Q T ) 2 + 1 C 3 ( 1 + δ ) k U γ k 1 ( u ¯ [ M ] ) V 2 ( Q T ) 2 ( 1 + δ ) + 1 , k 1 ,
which implies the inequality:
d k C 3 1 / γ k ( 1 + δ ) k / γ k d k 1 , k 1 ,
where d k = | u ¯ [ M ] | γ k V 2 ( Q T ) 2 + 1 1 / γ k .
Iterating these inequalities, we find:
d k C 4 d 0 = C 4 ( u ¯ [ M ] V 2 ( Q T ) 2 + 1 ) C 4 ( u ¯ V 2 ( Q T ) 2 + 1 ) , k 1 .
Thus:
u ¯ [ M ] L 2 γ k ( Q T ) 2 T 1 / γ k | u ¯ [ M ] ) | γ k 1 u ¯ [ M ] V 2 ( Q T ) 2 / γ k T 1 / γ k y k T 1 / γ k C 4 u ¯ V 2 ( Q T ) 2 + 1 , k 1 .
The limit passage as k leads to the estimate:
u ¯ [ M ] L ( Q T ) 2 C 4 u ¯ V 2 ( Q T ) 2 + 1 M 1 ,
which implies that u L ( Q T ) and:
u L ( Q T ) 2 C 6 u V 2 ( Q T ) 2 + A 2 .
Taking into account the estimate (44) with γ = 1 , we obtain the estimate (51).  ☐
Theorem 11.
Let u be a weak solution to Problem P . If the assumptions (46), (47) are satisfied, then there exists a constant β ( 0 , β 0 / 2 ) such that e β | u | V 2 ( Q T ) , the estimate (48) holds and u C ( [ 0 , T ] ; L s ( G ) for all s [ 1 , ) .
Proof of Theorem 11.
Let:
U β ( u ) = [ e β | u | 1 ] sgn u , w β ( u ) = 0 u ( U β ( s ) ) 2 d s = β 2 e 2 β | u | 1 sgn u , W β ( u ) = 0 u w β ( s ) d s = 1 4 e 2 β | u | 2 β u ,
where β > 0 . Since:
1 4 ( U β ( u ) ) 2 W β ( u ) 1 2 ( U β ( u ) ) 2 , w β ( u ) β U β ( u ) 2 + 1 ,
inequality (72) implies the inequality
U β ( u [ M ] ( t ) ) L 2 ( G ) 2 + U β ( u [ M ] ) L 2 ( Q t ) 2 C 1 U β ( u [ M ] ( 0 ) ) L 2 ( G ) 2 + C 1 β [ f L r ( 0 , T ; L q ( G ) ) U β ( u [ M ] ) L r ¯ 1 ( 0 , T ; L q ¯ 1 ( G ) ) 2 + 1 + f * L r * ( 0 , T ; L q * ( G ) ) U β ( u [ M ] ) L r ¯ * ( 0 , T ; L q ¯ * ( G ) ) 2 + 1 + g * L r * ( 0 , T ; L q * ( G ^ ) ) U β ( t r u b [ M ] ) L r ¯ * ( 0 , T ; L q ¯ * ( G ^ ) ) + 1 ] , t ( 0 , T ] ,
where r ¯ 1 = 2 r , q ¯ 1 = 2 q , r ¯ * = 2 r * , q ¯ * = 2 q * .
Using (13) and (14), we derive the estimate:
U β ( u [ M ] ) V 2 ( Q T ) 2 2 C 1 U β ( u 0 ) L 2 ( G ) 2 + β C f , J * U β ( u [ M ] ) V 2 ( Q T ) 2 + 1 ,
where C f , J * = C 2 f L r ( 0 , T ; L q ( G ) ) + J * L r * ( 0 , T ; L ^ q * ( Γ ^ ) ) .
Taking β = min { β 0 , C f , J * 1 } / 2 , we obtain the estimate:
e β | u [ M ] | V 2 ( Q T ) C e β | u 0 | L 1 ( G ) + 1 .
Since u [ M ] u in L 2 ( 0 , T ; W 1 , 2 ( G ) ) and u [ M ] ( t ) u ( t ) in L 2 ( G ) for all t [ 0 , T ] as M , the estimate (76) implies that e β | u | V 2 ( Q T ) and the estimate (48) holds.
Note that the conditions of Theorem 9 are satisfied for all p [ 3 , ) . Therefore, u C ( [ 0 , T ] ; L s ( G ) ) for all s [ 1 , ) .  ☐

9. Stability and Uniqueness of Weak Solutions to Problem P : Comparison Theorem

The proof given in this section uses some ideas of the method [70] proposed for proving comparison theorems for quasilinear elliptic equations. Special modifications of this method for some nonstationary radiative–conductive heat transfer problems were used in [41,53,57].
The following theorem concerns the stability of weak solutions to Problem P with respect to data.
Theorem 12.
Let u 1 and u 2 be two weak solutions to Problem P with ( u 0 , 1 , f 1 , J * 1 ) and ( u 0 , 2 , f 2 , J * 2 ) instead of ( u 0 , f , J * ) . Then, the estimates (39), (40) hold.
Proof of Theorem 12.
We put Δ u = u 1 u 2 , Δ u 0 = u 0 , 1 u 0 , 2 Δ f = f 1 f 2 , Δ J * = J * 1 J * 2 .
Let 0 < t T , 0 < δ < 1 , δ is a parameter. We introduce the sets:
Q t + = { ( x , t ) Q t | Δ u ( x , t ) > 0 } , Q t = { ( x , t ) Q t | Δ u ( x , t ) 0 } , Q t δ = { ( x , t ) Q t | Δ u ( x , t ) δ } , Q t ( 0 , δ ) = { ( x , t ) Q t | 0 < Δ u ( x , t ) < δ } , Q s , t + = Q s , t Q t + , Q s , t = Q s , t Q t , Q s , t δ = Q s , t Q t δ , Q s , t ( 0 , δ ) = Q s , t Q t ( 0 , δ ) .
We introduce the function v δ = δ 1 ( Δ u ) [ 0 , δ ] = min { δ 1 [ Δ u ] + , 1 } . It is clear that 0 v δ 1 ; moreover v δ ( x , t ) = 0 for ( x , t ) Q T , v δ ( x , t ) = 1 for ( x , t ) Q T δ and lim δ 0 v δ ( x , t ) = 1 for ( x , t ) Q T + .
Subtracting from each other the identities (34) corresponding to the definitions of the weak solutions u 1 and u 2 leads to the identity:
0 T ( c p Δ u ( t ) , v ) G d d t η ( t ) d t + 0 T [ a ( u 1 ( t ) , v ) a ( u 2 ( t ) , v ) ] η ( t ) d t + 0 T [ b ( u 1 ( t ) , v ) b ( u 2 ( t ) , v ) ] η ( t ) d t = ( c p Δ u 0 , v ) G · η ( 0 ) + 0 T ( Δ f ( t ) , v ) G η ( t ) d t + 0 T ( Δ f * ( t ) , v s ) G s η ( t ) d t + 0 T ( Δ g * ( t ) , t r v b ) G b η ( t ) d t v V , η C * [ 0 , T ] ,
where Δ f * = f * 1 f * 2 = 4 ϰ A ^ Ω [ Δ J * ] , Δ g * = g * 1 g * 2 = π B ^ [ Δ J * ] .
Using Lemma 6 with w ( u ) = δ 1 u , M = δ , W ( δ ) ( u ) = δ 1 0 u s [ δ ] d s and taking into account that:
W ( δ ) ( [ Δ u 0 ] + ) [ Δ u 0 ] + , 0 w ( ( Δ u ) [ 0 , δ ] ) = v δ 1 , [ Δ f * ] + L 1 ( Q t ) + [ Δ g * ] + L 1 ( Q b , t ) [ Δ J * ] L 1 ( 0 , T ; L ^ 1 ( Γ ^ ) )
(the inequality (78) follows from (27)), we have:
c p W ( δ ) ( [ Δ u ] + ( t ) ) L 1 ( G ) + ( λ ( · , u 1 ) u 1 λ ( · , u 2 ) u 2 , v δ ) L 2 ( Q t ) + 0 t [ a 1 ( u 1 ( t ) , v δ ( t ) ) a 1 ( u 2 ( t ) , v δ ( t ) ) ] d t + 0 t [ b ( u 1 ( t ) , v δ ( t ) ) b ( u 2 ( t ) , v δ ( t ) ) ] d t = c p W ( δ ) ( [ Δ u 0 ] + ) L 1 ( G ) + 0 t ( Δ f ( t ) , v δ ( t ) ) G + ( Δ f * ( t ) , v s δ ( t ) ) G s + ( Δ g * ( t ) , t r v b δ ( t ) ) G b d t c p [ Δ u 0 ] + L 1 ( G ) + [ Δ f ] + L 1 ( Q t ) + [ Δ J * ] + L 1 ( 0 , t ; L ^ 1 ( Γ ^ ) ) t ( 0 , T ] .
Using the fact that v δ = δ 1 ( u 1 u 2 ) almost everywhere on Q T ( 0 , δ ) and v δ = 0 almost everywhere on Q T Q T δ and taking into account assumptions (32), (33), we find that:
( λ ( · , u 1 ) u 1 λ ( · , u 2 ) u 2 , v δ ) L 2 ( Q t ) = δ ( λ ( · , u 1 ) v δ , v δ ) L 2 ( Q t ) + ( [ λ ( · , u 1 ) λ ( · , u 2 ) ] u 2 , v δ ) L 2 ( Q t ( 0 , δ ) ) δ λ min v δ L 2 ( Q t ) 2 L δ 1 / 2 u 2 L 2 ( Q t ( 0 , δ ) ) v δ L 2 ( Q t ) L 2 4 λ min u 2 L 2 ( Q t ( 0 , δ ) ) 2 .
We note also that:
a 1 ( u 1 , v δ ) a 1 ( u 2 , v δ ) = ( γ [ t r Δ u s t r Δ u b ] , t r v s δ t r v b δ ) G b 0 .
We set Δ h ( u s ) = h ( u s 1 ) h ( u s 2 ) for ( x , t ) Q s , T , Δ h ( t r u b ) = h ( t r u b 1 ) h ( t r u b 2 ) for ( x , t ) Q b , T and introduce the sets:
+ Q b , T = { ( x , t ) Q b , T | t r Δ u b ( x , t ) > 0 } , Q b , T = { ( x , t ) Q b , T | t r Δ u b ( x , t ) 0 } , δ Q b , T = { ( x , t ) Q b , T | t r Δ u b ( x , t ) δ } , ( 0 , δ ) Q b , T = { ( x , t ) Q b , T | 0 < t r Δ u b ( x , t ) < δ } .
From the formula (56), it follows that:
b ( u 1 , v δ ) b ( u 2 , v δ ) = ( 4 ϰ Δ h ( u s ) , C [ v s δ , t r v b δ ] ) G s + ( Δ h ( t r u b ) , D [ v s δ , t r v b δ ] ) G b .
Noticing that:
C [ v s δ , t r v b δ ] = v s δ A Ω [ v s δ , t r v b δ ] = A Ω [ v s δ , t r v b δ ] 0 , ( x , t ) Q s , T , C [ v s δ , t r v b δ ] = v s δ A Ω [ v s δ , t r v b δ ] v s δ 1 , ( x , t ) Q s , T + , D [ v s δ , t r v b δ ] = t r v b δ B [ v s δ , t r v b δ ] = B [ v s δ , t r v b δ ] 0 , ( x , t ) Q b , T , D [ v s δ , t r v b δ ] = t r v b δ B [ v s δ , t r v b δ ] t r v b δ 1 , ( x , t ) + Q b , T ,
we find that:
0 t [ b ( u 1 ( t ) , v δ ( t ) ) b ( u 2 ( t ) , v δ ( t ) ) ] d t = ( 4 ϰ Δ h ( u s ) , C [ v s δ , t r v b δ ] ) Q s , t + ( 4 ϰ ( Δ h ( u s ) , C [ v s δ , t r v b δ ] ) Q s , t + + ( Δ h ( t r u b ) , D [ v s δ , t r v b δ ] ) Q b , t + ( Δ h ( t r u b ) , D [ v s δ , t r v b δ ] ) + Q b , t ( 4 ϰ Δ h ( u s ) , v s δ 1 ) Q s , t + + ( Δ h ( t r u b ) , t r v b δ 1 ) + Q b , t = ( 4 ϰ Δ h ( u s ) , v s δ 1 ) Q s , t ( 0 , δ ) + ( Δ h ( t r u b ) , t r v b δ 1 ) ( 0 , δ ) Q b , t 4 ϰ Δ h ( u s ) L 1 ( Q s , t ( 0 , δ ) ) Δ h ( t r u b ) L 1 ( ( 0 , δ ) Q b , t ) .
It follows from the inequality (79) and the estimates (80)–(82) that:
c p W ( δ ) ( [ Δ u ] + ( t ) ) L 1 ( G ) L 2 4 λ min u 2 L 2 ( Q T ( 0 , δ ) ) 2 + 4 ϰ Δ h ( u s ) L 1 ( Q s , T ( 0 , δ ) ) + Δ h ( t r u b ) L 1 ( ( 0 , δ ) Q b , T ) + c p [ Δ u 0 ] + L 1 ( G ) + ( [ Δ f ] + , [ Δ J * ] + ) 1 , 1 , 1 , 1 t ( 0 , T ] .
We pass to the limit as δ 0 in this inequality. Since:
W ( δ ) ( [ Δ u ] + ) [ Δ u ] + δ ,
then:
c p W ( δ ) ( [ Δ u ] + ( t ) ) L 1 ( G ) c p [ Δ u ] + ( t ) L 1 ( G ) .
The first three terms on the right hand side of (83) tend to zero as δ 0 , since:
meas ( Q T ( 0 , δ ) ; d x d t ) 0 , meas ( Q s , T ( 0 , δ ) ; d x d t ) 0 , meas ( ( 0 , δ ) Q b , T ; d σ ( x ) d t ) 0 .
Thus, (83) implies the inequality:
c p [ Δ u ] + ( t ) L 1 ( G ) c p [ Δ u 0 ] + L 1 ( G ) + ( [ Δ f ] + , [ Δ J * ] + ) 1 , 1 , 1 , 1 t ( 0 , T ] .
The following inequality can be established in the same way:
c p [ Δ u ] ( t ) L 1 ( G ) c p [ Δ u 0 ] L 1 ( G ) + ( [ Δ f ] , [ Δ J * ] ) 1 , 1 , 1 , 1 t ( 0 , T ] .
Adding (84) and (85), we obtain the inequality:
c p Δ u ( t ) L 1 ( G ) c p Δ u 0 L 1 ( G ) + ( Δ f , Δ J * ) 1 , 1 , 1 , 1 t ( 0 , T ] .
The inequalities (84), (86) imply the estimates (39), (40).  ☐
Corollary 2
(Comparison theorem). If u 0 , 1 u 0 , 2 , f 1 f 2 and J * 1 J * 2 , then u 1 u 2 .
Corollary 3
(Uniqueness theorem). If a weak solution to Problem P exists, then it is unique.

10. Solvability of Problem P

Theorem 13.
A weak solution to Problem P exists and is unique.
Proof of Theorem 13.
Firstly, we suppose that assumptions (49), (50) hold. By Theorems 8 and 10, for all n > 0 , there exists a function u V 2 1 , 0 ( Q T ) L ( Q T ) V ( Q T ) , which is a weak solution to Problem P [ n ] and satisfies the estimate:
u L ( Q T ) M = C ( u 0 , f , J * ) , r , q , r * , q * ,
where M does not depend on n. By this estimate, h [ n ] ( u ) = h ( u ) for n > h ( M ) . Therefore, a weak solution to Problem P [ n ] with n > h ( M ) is simultaneously a weak solution to Problem P .
Now, we prove the existence of a solution without additional assumptions (49), (50). Let N,M be natural numbers. Since ( u 0 ) [ N , M ] L ( G ) , f [ N , M ] L ( Q T ) L ( 0 , T ; L 2 ( G ) ) , J * [ N , M ] L ( Γ ^ × ( 0 , T ) ) L ( 0 , T ; L ^ 3 ( Γ ^ ) ) then, by the first part of the proof, Problem P with ( u 0 ) [ N , M ] , f [ N , M ] and J * [ N , M ] in the role of u 0 , f, and J * has a weak solution u ( N , M ) such that:
0 T ( c p u ( N , M ) ( t ) , v ) G d d t η ( t ) d t + 0 T A ( u ( N , M ) ( t ) , v ) η ( t ) d t = ( c p ( u 0 ) [ N , M ] , v ) G · η ( 0 ) + 0 T ( f [ N , M ] ( t ) , v ) G η ( t ) d t + 0 T ( f * [ N , M ] , v s ) G s + ( g * ( N , M ) ( t ) , t r v b ) G b η ( t ) d t v V , η C * [ 0 , T ] .
Here, f * ( N , M ) = 4 π ϰ A ^ Ω [ J * [ N , M ] ] , g * ( N , M ) = π B ^ [ J * [ N , M ] ] .
Note that:
( u 0 ) [ N , M ] , f [ N , M ] , J * [ N , M ] ) p , r , q , r * , q * C 0 = ( u 0 , f , J * ) p , r , q , r * , q * .
So, by Theorem 9, the following uniform parameters N and M estimates hold:
u ( N , M ) V 2 ( Q T ) C ( u 0 ) [ N , M ] , f [ N , M ] , J * [ N , M ] ) p , r , q , r * , q * C 1 = C C 0 ,
| u [ N , M ] | 1 / 2 u [ N , M ] V 2 ( Q T ) 2 / 3 C ( u 0 ) [ N , M ] , f [ N , M ] , J * [ N , M ] ) p , r , q , r * , q * C 1 .
The estimate (89) and the inequalities (17), (18) imply the estimates
u ( N , M ) L ( 0 , T ; L 3 ( G ) ) C 1 ,
u ( N , M ) L 5 ( Q T ) C 2 ,
t r u b ( N , M ) L 4 ( Q b , T ) C 3 .
We fix N. Since:
( u 0 ) [ N , M ] ( u 0 ) [ N , M + 1 ] , f [ N , M ] f [ N , M + 1 ] , J * [ N , M ] J * [ N , M + 1 ] ,
then by Corollary 2, the sequence { u ( N , M ) } M = 1 is non-decreasing with respect to M. Therefore, from the estimate (91), by virtue of Levi’s monotone convergence theorem, there exists a function u ( N ) L 5 ( Q T ) such that u ( N , M ) u ( N ) in L 5 ( Q T ) and almost everywhere on Q T as M . From (88), it follows that u ( N ) V 2 ( Q T ) , u ( N , M ) u ( N ) weakly in L 2 ( 0 , T ; W 1 , 2 ( G ) ) and weakly stars in L ( 0 , T ; L 3 ( G ) ) as M . From the multiplicative inequality (19), it follows that t r u s ( N , M ) t r u s ( N ) , t r u b ( N , M ) t r u b ( N ) in L 2 ( Q b , T ) .
Since the sequence t r u b ( N , M ) M = 1 does not decrease with respect to M, then it follows from the estimate (92) and Levi’s monotone convergence theorem that t r u b ( N ) L 4 ( Q b , T ) and t r u b ( N , M ) t r u b ( N ) in L 4 ( Q b , T ) . So h ( t r u b ( N , M ) ) h ( t r u b ( N ) ) in L 1 ( Q b , T ) .
Passage to the limit as M in (88), (90)–(92) leads to the inequalities:
u ( N ) V 2 ( Q T ) C 1 , u ( N ) L ( 0 , T ; L 3 ( G ) ) C 1 ,
u ( N ) L 5 ( Q T ) C 2 , t r u b ( N ) L 4 ( Q b , T ) C 2 .
Let v V . Since u ( N , M ) u ( N ) weakly in L 2 ( Q T ) , λ ( · , u ( N , M ) ) λ ( · , u ( N ) ) almost everywhere on Q T , then a 0 ( u ( N , M ) , v ) a 0 ( u ( N ) , v ) weakly in L 1 ( 0 , T ) as M . It is also clear that a 1 ( u ( N , M ) , v ) a 1 ( u ( N ) , v ) in L 1 ( 0 , T ) . Thus, a ( u ( N , M ) , v ) a ( u ( N ) , v ) weakly in L 1 ( 0 , T ) as M .
Since h ( u ( N , M ) ) h ( u ( N ) ) in L 1 ( Q T ) and h ( t r u b ( N , M ) ) h ( t r u b ( N ) ) in L 1 ( Q b , T ) , then:
( h ( u ( N , M ) ) , C [ v s , t r v b ] ) G s ( h ( u ( N ) ) , C [ v s , t r v b ] ) G s in L 1 ( 0 , T ) , ( h ( t r u b ( N , M ) ) , D [ v s , t r v b ] ) G b ( h ( t r u b ( N ) ) , D [ v s , t r v b ] ) G b in L 1 ( 0 , T ) .
Thus, b ( u ( N , M ) , v ) b ( u ( N ) , v ) in L 1 ( 0 , T ) .
Passing to the limit as M in the identity (87), we arrive at the identity:
0 T ( c p u ( N ) ( t ) , v ) G d d t η ( t ) d t + 0 T A ( u ( N ) ( t ) , v ) η ( t ) d t
= ( c p ( u 0 ) [ N , ] , v ) G · η ( 0 ) + 0 T ( f [ N , ] ( t ) , v ) G η ( t ) d t + 0 T ( f * ( N , ) ( t ) , v ) G s η ( t ) d t
+ ( g * ( N , ) ( t ) , v ) G b η ( t ) d t v V , η C * [ 0 , T ] ,
where ( u 0 ) [ N , ] = max { u 0 , N } , f [ N , ] = max { f , N } , J * [ N , ] = max { J * , N } , f * ( N , ) = 4 π ϰ A ^ Ω [ J * [ N , ] ] , g * ( N , ) = π B ^ [ J * [ N , ] ] .
Theorem 12 implies the estimate:
c ̲ p u ( N , M 1 ) u ( N , M 2 ) C ( [ 0 , T ] ; L 1 ( G ) ) c ¯ p ( u 0 ) [ N , M 1 ] ( u 0 ) [ N , M 2 ] L 1 ( G ) + ( f [ N , M 1 ] f [ N , M 2 ] , J * [ N , M 1 ] J * [ N , M 2 ] ) 1 , 1 , 1 , 1 M 1 1 , M 2 1 .
This estimate means that { u ( N , M ) } M = 1 is a Cauchy sequence in C ( [ 0 , T ] ; L 1 ( G ) ) . Hence, u ( N ) C ( [ 0 , T ] ; L 1 ( G ) ) . Taking into account that u ( N ) L ( 0 , T ; L 3 ( G ) ) , we have u ( N ) C ( [ 0 , T ] ; L 2 ( G ) ) .
Thus, the function u ( N ) is a weak solution to Problem P corresponding to the data ( u 0 ) [ N , ] , f [ N , ] and J * [ N , ] in the role of u 0 , f and J * .
Since:
( u 0 ) [ ( N + 1 ) , ] ( u 0 ) [ N , ] , f [ ( N + 1 ) , ] f [ N , ] , J * [ ( N + 1 ) , ] g [ N , ] ,
then by virtue of Corollary 2, the sequence { u ( N ) } N = 1 is non-increasing. Therefore, from the estimates (93), (94) it follows that there exists a function u V 2 ( Q T ) L 5 ( Q T ) such that u ( N ) u weakly in L 2 ( 0 , T ; W 1 , 2 ( G ) ) , weakly stars in L ( 0 , T ; L 3 ( G ) ) , strongly in L 5 ( Q T ) and almost everywhere on Q T as N . In addition, t r u s ( N ) t r u s in L 2 ( Q b , T ) and t r u b ( N ) t r u b in L 4 ( Q b , T ) . As a consequence, h ( u ( N ) ) h ( u ) in L 1 ( Q T ) , h ( t r u b ( N ) ) h ( u b ) in L 1 ( Q ^ b , T ) .
Therefore, a ( u ( N ) , v ) a ( u , v ) weakly in L 1 ( 0 , T ) and b ( u ( N ) , v ) b ( u , v ) in L 1 ( 0 , T ) for all v V . Passage to the limit as N in the identity (96) gives the identity (34).
Theorem 12 implies the estimate:
c ̲ p u ( N 1 ) u ( N 2 ) C ( [ 0 , T ] ; L 1 ( G ) ) c ¯ p ( u 0 ) [ N 1 , ] ( u 0 ) [ N 2 , ] L 1 ( G ) + ( f [ N 1 , ] f [ N 2 , ] , J * [ N 1 , ] J * [ N 2 , ] ) 1 , 1 , 1 , 1 N 1 1 , N 2 1 .
This inequality means that { u ( N ) } N = 1 is a Cauchy sequence in C ( [ 0 , T ] ; L 1 ( G ) ) . Hence, u C ( [ 0 , T ] ; L 1 ( G ) ) . Taking into account that u L ( 0 , T ; L 3 ( G ) ) , we have u C ( [ 0 , T ] ; L 2 ( G ) ) .
We proved the existence of a weak solution to Problem P . Its uniqueness follows from Corollary 3.  ☐

11. Justification of the Main Results

Proof of Theorem 2.
Note that (see Remark 4) the pair ( u , I ) V ( Q T ) × L 1 ( 0 , T ; W 1 ( D ) ) is a weak solution to the problem (1)–(9) if and only if u V ( Q T ) is a weak solution to problem P and I is expressed by the formula:
I = 1 π A [ h ( u s ) , h ( t r u b ) ] + A ^ [ J * ] .
Therefore, the existence and uniqueness of a weak solution to the problem (1)–(9) follow directly from Theorem 13.  ☐
Proof of Theorem 4.
Assume that the conditions of Theorem 4 be satisfied.
By Corallary 2 we have u 1 u 2 . Thus, h ( u s 1 ) h ( u s 2 ) , h ( t r u b 1 ) h ( t r u b 2 ) and
I 1 = 1 π A [ h ( u s 1 ) , h ( t r u b 1 ) ] + A ^ [ J * 1 ] 1 π A [ h ( u s 2 ) , h ( t r u b 2 ) ] + A ^ [ J * 2 ] = I 2 .
 ☐
Proof of Theorem 5.
Assume that the conditions of Theorem 5 are satisfied.
By Theorem 9, | u | γ 1 u V 2 ( Q T ) for all γ [ 1 , p / 2 ] , the estimate (44) is valid and u C ( [ 0 , T ] ; L s ( G ) ) for all s [ 1 , p ) .
From | u | p / 2 1 u V 2 ( Q T ) , it follows (see (13), (14)) that h ( u s ) L r s ( 0 , T ; L q s ( G s ) ) , h ( t r u b ) L r s ( 0 , T ; L q s ( G b ) ) for all r s [ 1 , p / 2 ] , q s [ 1 , p / 2 ] such that 1 / r s + 1 / q s 6 / p . In addition, the following estimate holds:
h ( u s ) L r s ( 0 , T ; L q s ( G s ) ) + h ( t r u b ) L r s ( 0 , T ; L q s ( G b ) ) C 1 | u | p / 2 1 u V 2 ( Q T ) 8 / p .
From this estimate, the boundedness of the operator A : L q s ( G s , G b ) W q s ( D s ) and the estimate (44) with γ = p / 2 , it follows that I s L r s ( 0 , T ; W q s ( D s ) ) and:
I s L r s ( 0 , T ; W q s ( D s ) ) C 2 | u | p / 2 1 u V 2 ( Q T ) 8 / p C ( u 0 , f , J * ) p , r , q , r * , q * 4 .
 ☐
Proof of Theorem 6.
Assume that the conditions of Theorem 6 are satisfied.
By Theorem 11, there exists a constant β ( 0 , β 0 / 2 ) such that e β | u | V 2 ( Q T ) and the estimate (48) holds.
Note that the conditions of Theorem 5 are satisfied for all p [ 3 , ) . Thus, I s L r s ( 0 , T ; W q s ( D s ) ) for all r s [ 1 , ) , q s [ 1 , ) .  ☐
Proof of Theorem 7.
Assume that the conditions of Theorem 7 are satisfied.
By Theorem 10, u L ( Q T ) and the estimate (51) holds. Consequently,
( h ( u s ( t ) ) , h ( t r u b ( t ) ) L ( G s , G b ) for almost all t ( 0 , T ) and:
ess sup t ( 0 , T ) ( h ( u s ( t ) ) , h ( t r u b ( t ) ) L ( G s , G b ) σ 0 ess sup t ( 0 , T ) u ( t ) L ( G ) 4 C 1 ( u 0 , f , J * ) , r , q , r * , q * 4 .
Therefore, I s ( t ) = 1 π A [ h ( u s ( t ) ) , h ( t r u b ( t ) ) ] W ( D s ) for almost all t ( 0 , T ) and by virtue of the estimate (98):
ess sup t ( 0 , T ) I s ( t ) W ( D s ) C ( u 0 , f , J * ) , r , q , r * , q * 4 ,
where C = C 1 1 π A L ( G s , G b ) W ( D s ) .  ☐

12. Conclusions

In this paper, the author continues to construct a mathematical theory of complex heat transfer problems.
A nonstationary initial–boundary value problem governing a radiative–conductive heat transfer in a convex semitransparent body with an absolutely black inclusions was considered. To describe the process, a system consisting of two heat equations and the integro-differential radiative transfer equation was used. This system is supplied by boundary conditions, which describe the energy exchange between semitransparent body, external media and opaque inclusions.
The unique solvability of this problem was proven. In addition, the stability of solutions with respect to the data was proven, which established a comparison theorem. Besides, results on improving the properties of solutions with an increase in the summability of the data were established. All results are global in terms of time and data.
The considered mathematical model of radiative–conductive heat transfer contains a number of simplifying assumptions. One should consider the process of heat transfer in a system of bodies, and not in one convex body. In a more complex model, it should be taken into account that the properties of the semitransparent medium and the radiation intensity depend on the radiation frequency. In addition, inclusions may not be completely black, but gray or even “colored”. The author expects to study the more complex corresponding models in the near future.

Funding

The results were obtained within the framework of the state assignments of the Russian Ministry of Education and Science (project FSWF-2020-0022).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Sparrow, E.M.; Cess, R.D. Radiation Heat Transfer; Hemisphere: New York, NY, USA, 1978. [Google Scholar]
  2. Siegel, R.; Howell, J.R. Thermal Radiation Heat Transfer; Taylor and Francis-Hemishpere: Washington, DC, USA, 2001. [Google Scholar]
  3. Ozisik, M.N. Radiative Transfer and Interactions with Conduction and Convection; John Willey & Sons: New York, NY, USA, 1973. [Google Scholar]
  4. Modest, F.M. Radiative Heat Transfer; Academic Press: New York, NY, USA, 2003. [Google Scholar]
  5. Tikhonov, A.N. Cooling of bodies under radiation following to the Stefan–Boltzmann law. Izv. AN SSSR Ser. Geogr. Geophiz. 1937, 3, 461–479. (In Russian) [Google Scholar]
  6. Tikhonov, A.N. On Volterra type functional equations and their applications to some problems of mathematical physics. Bull. MGU Sec. A Ser. Mat. Mekh. 1938, 8, 1–25. (In Russian) [Google Scholar]
  7. Amosov, A.A. A positive solution of an elliptic equation with nonlinear integral boundary condition of the radiation type. Math. Notes 1977, 22, 555–561. [Google Scholar] [CrossRef]
  8. Forste, J. Steady flow of incompressible volatile fluid under action of radiation. ZAMM 1977, 57, 265–267. [Google Scholar]
  9. Kuraev, G.N. The problem of stationary free convection under non-linear boundary conditions. USSR Comput. Math. Math. Phys. 1978, 18, 784–789. (In Russian) [Google Scholar] [CrossRef]
  10. Amosov, A.A. The solvability of a problem of radiation heat transfer. Sov. Phys. Dokl. 1979, 24, 261–262. [Google Scholar]
  11. Amosov, A.A. The limit connection between two problems of radiation heat transfer. Sov. Phys. Dokl. 1979, 24, 439–441. [Google Scholar]
  12. Amosov, A.A. Solvability of the radiative heat transfer problem in accordance with the Stefan–Boltzmann law. Mosc. Univ. Comput. Math. Cybern. 1980, 3, 18–27. [Google Scholar]
  13. Amosov, A.A. Well-posedness “in the large” of initial-and-boundary-value problems for the system of dynamical equations of a viscous radiating gas. Sov. Phys. Dokl. 1985, 30, 129–131. [Google Scholar]
  14. Kuraev, G.N. On the solvability of the problem of stationary heat convection with radiant convective heat exchange on the boundary. USSR Comput. Math. Math. Phys. 1986, 26, 184–188. [Google Scholar] [CrossRef]
  15. Perret, C.; Witomski, P. Equation de la chaleur et reflections multiples. Ann. de l’Institut Henri Poincare (C) Non Lineare Anal. 1991, 8, 677–689. [Google Scholar] [CrossRef] [Green Version]
  16. Saldanha da Gama, R.M. Existence, uniqueness and construction of the solution of the energy transfer problem in a rigid and non-convex blackbody. Z. Angew. Math. Phys. 1991, 42, 334–347. [Google Scholar] [CrossRef]
  17. Saldanha da Gama, R.M. An alternative mathematical modelling for coupled conduction/radiation energy transfer phenomenon in a system of n grey bodies surrounded by a vacuum. Int. J. Non-Linear Mech. 1995, 30, 433–447. [Google Scholar] [CrossRef]
  18. Kelley, C. Existence and uniqueness of solutions of nonlinear systems of conductive radiative heat transfer equations. Transp. Theory Stat. Phys. 1996, 25, 249–260. [Google Scholar] [CrossRef] [Green Version]
  19. Qatanani, N. Loesungsverfahren und Analysis der Integral-Gleichung Fuer das Hohlraum-Strahlungs-Problem. Ph.D. Thesis, University of Stuttgart, Stuttgart, Germany, 1996. [Google Scholar]
  20. Tiihonen, T. Stefan-Boltzmann radiation on non-convex surfaces. Math. Methods Appl. Sci. 1997, 20, 47–57. [Google Scholar] [CrossRef]
  21. Tiihonen, T. A nonlocal problem arising from heat radiation in non-convex surfaces. Eur. J. Appl. Math. 1997, 8, 403–416. [Google Scholar] [CrossRef]
  22. Laitinen, M.; Tiihonen, T. Heat transfer in conducting and radiating bodies. Appl. Math. Lett. 1997, 10, 5–8. [Google Scholar] [CrossRef] [Green Version]
  23. Gergo, L.; Stoyan, G. On a mathematical model of a radiating, viscous, heatconducting fluid: Remarks on a paper by J. Förste. ZAMM Zeitschrift für Angewandte Mathematik und Mechanik 1997, 77, 367–375. [Google Scholar] [CrossRef]
  24. Laitinen, M.T.; Tiihonen, T. Integro—Differential equation modelling heat transfer in conducting, radiating and semitransparant materials. Math. Methods Appl. Sci. 1998, 21, 375–392. [Google Scholar] [CrossRef]
  25. Metzger, M. Existence for a time-dependent heat equation with non-local radiation terms. Math. Methods Appl. Sci. 1999, 22, 1101–1119. [Google Scholar] [CrossRef]
  26. Laitinen, M.; Tiihonen, T. Conductive-radiative heat transfer in grey materials. Q. Appl. Math. 2001, 59, 737–768. [Google Scholar] [CrossRef] [Green Version]
  27. Laitinen, M.T. Asymptotic analysis of conductive-radiative heat transfer. Asymptot. Anal. 2002, 29, 323–342. [Google Scholar]
  28. Saldanha da Gama, R.M. On the conduction/radiation heat transfer problem in a body with wavelength-dependent properties. Appl. Math. Model. 2004, 28, 795–816. [Google Scholar] [CrossRef]
  29. Thompson, M.; Segatto, C.; De Vilhena, M.T. Existence theory for the solution of a stationary nonlinear conductive-radiative heat-transfer problem in three space dimensions. Transp. Theory Stat. Phys. 2004, 33, 563–576. [Google Scholar] [CrossRef]
  30. Amosov, A.A. Global solvability of nonlinear nonstationary problem with nonlocal boundary condition of radiative heat transfer type. Differ. Equ. 2005, 41, 96–109. [Google Scholar] [CrossRef]
  31. Qatanani, N. Analysis of the heat equation with non-local radiation terms in a non-convex diffuse and grey surfaces. Eur. J. Sci. Rec. 2006, 15, 245–254. [Google Scholar]
  32. Qatanani, N. Qualitative analysis of the radiative energy transfer model. Eur. J. Sci. Rec. 2007, 17, 379–391. [Google Scholar]
  33. Qatanani, N.; Barham, R.; Heeh, Q. Existence and uniqueness of the solution of the coupled conduction-radiation energy transfer on diffuse-grey surfaces. Surv. Math. Appl. 2007, 2, 43–58. [Google Scholar]
  34. Qatanani, N.M.; Heeh, Q.M. On existence and uniqueness theorem concerning time-dependent heat transfer model. Appl. Appl. Math. 2008, 2, 235–253. [Google Scholar]
  35. Thompson, M.; De Vilhena, M.T.; Bodmann, B.E.J. Existence theory for radiative flows. Transp. Theory Stat. Phys. 2008, 37, 307–326. [Google Scholar] [CrossRef]
  36. Druet, P.-E. Weak solutions to a stationary heat equation with nonlocal radiation boundary condition and right-hand side in Lp (p≥1). Math. Methods Appl. Sci. 2009, 32, 135–166. [Google Scholar] [CrossRef]
  37. Druet, P.-E. Existence for the stationary MHD equations coupled to heat transfer with nonlocal radiation effects. Czechoslov. Math. J. 2009, 59, 791–825. [Google Scholar] [CrossRef] [Green Version]
  38. Druet, P.-E. Existence of weak solution to time-dependent MHD equations coupled to heat transfer with nonlocal radiation boundary conditions. Nonlinear Anal. Real. World Appl. 2009, 10, 2914–2936. [Google Scholar] [CrossRef]
  39. Druet, P.-E. Weak solutions to a time-dependent heat equation with nonlocal radiation boundary condition and arbitrary p-summable right-hand side. Appl. Math. 2010, 55, 111–149. [Google Scholar] [CrossRef] [Green Version]
  40. Amosov, A.A. Stationary nonlinear nonlocal problem of radiative-conductive heat transfer in a system of opaque bodies with properties depending on radiation frequency. J. Math. Sci. 2010, 164, 309–344. [Google Scholar] [CrossRef]
  41. Amosov, A.A. Nonstationary nonlinear nonlocal problem of radiative-conductive heat transfer in a system of opaque bodies with properties depending on the radiation frequency. J. Math. Sci. 2010, 165, 1–41. [Google Scholar] [CrossRef]
  42. Amosov, A.A. Nonstationary radiative-conductive heat transfer problem in a periodic system of grey heat shields. J. Math. Sci. 2010, 169, 1–45. [Google Scholar] [CrossRef]
  43. Druet, P.-E.; Philip, P. Noncompactness of integral operators modeling diffuse-gray radiation in polyhedral and transient settings. Integr. Equ. Oper. Theory 2011, 69, 101–111. [Google Scholar] [CrossRef]
  44. Sauter, E.; De Azevedo, F.S.; Thompson, M. Existence theory for one-dimensional quasilinear coupled conductive-radiative flows. Appl. Math. Comput. 2014, 233, 545–556. [Google Scholar] [CrossRef]
  45. Kovtanyuk, A.E.; Chebotarev, A.Y. Steady—State problem of complex heat transfer. Comp. Math. Math. Phys. 2014, 54, 719–726. [Google Scholar] [CrossRef]
  46. Kovtanyuk, A.E.; Chebotarev, A.Y.; Botkin, N.D.; Hoffmann, K.-H. Solvability of P1 approximation of a conductive-radiative heat transfer problem. Appl. Math. Comput. 2014, 249, 247–252. [Google Scholar] [CrossRef]
  47. Kovtanyuk, A.E.; Chebotarev, A.Y. Stationary free convection problem with radiative heat exchange. Differ. Equ. 2014, 50, 1592–1599. [Google Scholar] [CrossRef]
  48. Kovtanyuk, A.E.; Chebotarev, A.Y.; Botkin, N.D.; Hoffmann, K.-H. The unique solvability of a complex 3d heat transfer problem. J. Math. Anal. Appl. 2014, 409, 808–815. [Google Scholar] [CrossRef]
  49. Grenkin, G.V.; Chebotarev, A.Y. A nonstationary problem of complex heat transfer. Comput. Math. Mathem. Phys. 2014, 54, 1737–1747. [Google Scholar] [CrossRef]
  50. Kovtanyuk, A.E.; Chebotarev, A.Y.; Botkin, N.D.; Hoffmann, K.-H. Unique solvability a steady-state complex heat transfer model. Commun. Nonlinear Sci. Numer. Simulat. 2015, 20, 776–784. [Google Scholar] [CrossRef]
  51. Amosov, A. Solvability of a nonstationary problem of radiative—Conductive heat transfer in a system of semitransparent bodies. In Integral Methods in Science and Engineering; Constanda, C., Kirsch, A., Eds.; Birkhauser: Basel, Switzerland, 2015; pp. 1–13. [Google Scholar]
  52. Grenkin, G.V.; Chebotarev, A.Y. Nonstationary problem of free convection with radiative heat transfer. Comput. Math. Mathem. Phys. 2016, 56, 278–285. [Google Scholar] [CrossRef]
  53. Amosov, A. Unique solvability of a nonstationary problem of radiative—Conductive heat exchange in a system of semitransparent bodies. Russ. J. Math. Phys. 2016, 23, 309–334. [Google Scholar] [CrossRef]
  54. Amosov, A.A. Unique solvability of stationary radiative-conductive heat transfer problem in a system of semitransparent bodies. J. Math. Sci. 2017, 224, 618–646. [Google Scholar] [CrossRef]
  55. Amosov, A.A. Stationary problem of complex heat transfer in a system of semitransparent bodies with boundary conditions of diffuse reflection and refraction of radiation. Comput. Math. Mathem. Phys. 2017, 57, 515–540. [Google Scholar] [CrossRef]
  56. Chebotarev, A.Y.; Grenkin, G.V.; Kovtanyuk, A.E. Inhomogeneous steady-state problem of complex heat transfer. ESAIM Math. Model. Numer. Anal. 2017, 51, 2511–2519. [Google Scholar] [CrossRef]
  57. Amosov, A.A. Nonstationary problem of complex heat transfer in a system of semitransparent bodies with boundary-value conditions of diffuse reflection and refraction of radiation. J. Math. Sci. 2018, 233, 777–806. [Google Scholar] [CrossRef]
  58. Ghattassi, M.; Roche, J.R.; Schmitt, D. Existence and uniqueness of a transient state for the coupled radiative—conductive heat transfer problem. Comput. Math. Appl. 2018, 75, 3918–3928. [Google Scholar] [CrossRef] [Green Version]
  59. Chebotarev, A.Y.; Grenkin, G.V.; Kovtanyuk, A.E.; Botkin, N.D.; Hoffmann, K.-H. Diffusion approximation of the radiative-conductive heat transfer model with Fresnel matching conditions. Commun. Nonlinear Sci. Numer. Simul. 2018, 57, 290–298. [Google Scholar] [CrossRef]
  60. Grenkin, G.V.; Chebotarev, A.Y. Stability of stationary dslutins of the radiative heat transfer equations. Comput. Math. Mathem. Phys. 2018, 58, 1420–1425. [Google Scholar] [CrossRef]
  61. Kolobov, A.G.; Pak, T.V.; Chebotarev, A.Y. Stationary problem of radiative heat transfer with Cauchy boundary conditions. Comput. Math. Mathem. Phys. 2019, 59, 1199–1203. [Google Scholar] [CrossRef]
  62. Grenkin, G.V.; Chebotarev, A.Y. Inverse problem for equations of complex heat transfer. Comput. Math. Mathem. Phys. 2019, 59, 1361–1371. [Google Scholar] [CrossRef]
  63. Chebotarev, A.Y.; Kovtanyuk, A.E.; Botkin, N.D. Problem of radiation heat exchange with boundary conditions of the Cauchy type. Commun. Nonlinear Sci. Numer. Simul. 2019, 75, 262–269. [Google Scholar] [CrossRef]
  64. Amosov, A. Unique solvability of a stationary radiative-conductive heat transfer problem in a system consisting of an absolutely black body and several semitransparent bodies. Math. Methods Appl. Sci. 2021. [Google Scholar] [CrossRef]
  65. Amosov, A. Unique solvability of a stationary radiative-conductive heat transfer problem in a semitransparent body with absolutely black inclusions. Z. Angew. Math. Phys. 2021, 72, 104. [Google Scholar] [CrossRef]
  66. Amosov, A.A. Unique solvability of the stationary complex heat transfer problem in a system of gray bodies with semitransparent inclusions. J. Math. Sci. 2021, 255, 353–388. [Google Scholar] [CrossRef]
  67. Amosov, A.A. Radiative transfer equation with diffuse reflection and refraction conditions in a system of bodies with piecewise smooth boundaries. J. Math. Sci. 2016, 216, 155–216. [Google Scholar] [CrossRef]
  68. Amosov, A.A. Radiative transfer equation with Fresnel reflection and refraction conditions in a system of bodies with piecewise smooth boundaries. J. Math. Sci. 2016, 219, 821–849. [Google Scholar] [CrossRef]
  69. Amosov, A.A. Boundary Value Problems for the Radiation Transfer Equation with Reflection and Refraction Conditions; Tamara Rozhkovskaya Publisher: Novosibirsk, Russia, 2017. (In Russian) [Google Scholar]
  70. Křižek, M.; Liu, L. On a comparison principle for a quasilinear elliptic boundary value problem of a nonmonotone type. Appl. Math. 1996, 24, 97–107. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Amosov, A. Nonstationary Radiative–Conductive Heat Transfer Problem in a Semitransparent Body with Absolutely Black Inclusions. Mathematics 2021, 9, 1471. https://0-doi-org.brum.beds.ac.uk/10.3390/math9131471

AMA Style

Amosov A. Nonstationary Radiative–Conductive Heat Transfer Problem in a Semitransparent Body with Absolutely Black Inclusions. Mathematics. 2021; 9(13):1471. https://0-doi-org.brum.beds.ac.uk/10.3390/math9131471

Chicago/Turabian Style

Amosov, Andrey. 2021. "Nonstationary Radiative–Conductive Heat Transfer Problem in a Semitransparent Body with Absolutely Black Inclusions" Mathematics 9, no. 13: 1471. https://0-doi-org.brum.beds.ac.uk/10.3390/math9131471

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