## Generalized Solutions to Hyperbolic Systems with Nonlinear Conditions and Strongly Singular Data

### I. Kmit Institute for Applied Problems of Mechanics and Mathematics, Ukrainian Academy of Sciences Naukova St. 3b 79060 Lviv, Ukraine E-mail: kmit@ov.litech.net

Abstract
Using the framework of Colombeau algebras of generalized functions, we prove the existence and uniqueness results for global generalized solvability of semilinear hyperbolic systems with nonlinear nonlocal boundary conditions. We admit strong singularities in the differential equations as well as in the initial and boundary conditions. Our analysis covers the case of non-Lipshitz nonlinearities both in the differential equations and in the boundary conditions.

1 Introduction

We study existence and uniqueness of global generalized solutions to mixed problems for semilinear hyperbolic systems with nonlinear nonlocal boundary conditions.
Specifically, in the domain $\Pi =\left\{\left(x,t\right)|0  , $t>0\right\}$  we study the following problem:
 $\begin{array}{ccc}\left({\partial }_{t}+\Lambda \left(x,t\right){\partial }_{x}\right)U& =& F\left(x,t,U\right),\left(x,t\right)\in \Pi \end{array}$ (1)
 $\begin{array}{ccc}U\left(x,0\right)& =& A\left(x\right),x\in \left(0,l\right)\end{array}$ (2)
 $\begin{array}{ccc}{U}_{i}\left(0,t\right)& =& {H}_{i}\left(t,V\left(t\right)\right),k+1\le i\le n,t\in \left(0,\infty \right)\end{array}$
 $\begin{array}{ccc}{U}_{i}\left(l,t\right)& =& {H}_{i}\left(t,V\left(t\right)\right),1\le i\le k,t\in \left(0,\infty \right),\end{array}$ (3)
where $U$  , $F$  , and $A$  are real $n$  -vectors, $\Lambda =diag\left({\Lambda }_{1},...,{\Lambda }_{n}\right)$  is a diagonal matrix, ${\Lambda }_{1},...,{\Lambda }_{k}<0$  , ${\Lambda }_{k+1},...,{\Lambda }_{n}>0$  for some $1\le k\le n$  , and $V\left(t\right)=\left({U}_{1}\left(0,t\right),...,$  ${U}_{k}\left(0,t\right),{U}_{k+1}\left(l,t\right),...,{U}_{n}\left(l,t\right)\right)$  . Due to the conditions imposed on $\Lambda$  , the system ( 1 ) is non-strictly hyperbolic. Note also that the boundary of $\Pi$  is not characteristic.
We will denote $H=\left({H}_{1},...,{H}_{n}\right)$  .
Special cases of ( 1 )–( 3 ) arise in laser dynamics [8, 19, 21, 22and chemical kinetics [23.
All the data of the problem are allowed to be strongly singular, namely, they can be of any desired order of singularity. This entails nonlinear superpositions for (strongly singular) distributions in the right-hand sides of ( 1 )–( 3 ) (including compositions of the singular initial data and the singular characteristic curves). To tackle this complication, we use the framework of Colombeau algebra of generalized functions $\mathcal{G}\left(\overline{\Pi }\right)$  [1, 2, 16. We show that all superpositions appearing here are well defined in $\mathcal{G}\left(\overline{\Pi }\right)$  .
We establish a positive existence-uniqueness result in $\mathcal{G}\left(\overline{\Pi }\right)$  for the problem ( 1 )–( 3 ) with strongly singular initial data and with nonlinearities of the following type (more detailed description is given in Section 3): The functions $F$  and $H$  are either Lipshitz with Colombeau generalized numbers as Lipshitz constants or non-Lipshitz with less than quadratic growth in $U$  and $V$  .
For different aspects of the subject we refer the reader to sources [4, 5, 9, 12, 13, 16, 17, 14.
The essential assumption made on $F$  in papers [12, 16is that ${\text{grad}}_{U}F$  is globally bounded uniformly over $\left(x,t\right)$  varying in any compact set. The main complication with the non-Lipshitz nonlinearities when investigating Colombeau solutions lies in the following. The solutions in $\mathcal{G}$  are nets of smooth functions which are classical solutions to the associated problems with smooth initial data. To guarantee the existence of such solutions, in general, one has to assume that the nonlinearities have bounded gradients. In the present paper we tackle this complication combining classical and nonclassical approaches.
Papers [14, 15deal with Cauchy problems for semilinear hyperbolic systems ( 1 ) with $F$  slowly increasing at infinity. The nonlinear term is replaced by a suitable regularization ${F}_{\varepsilon }$  having a bounded gradient with respect to $U$  for every fixed $\varepsilon$  and converging to $F$  as $\varepsilon \to 0$  . The regularized system is solved in $\mathcal{G}\left({\mathbb{R}}^{2}\right)$  . Moreover, in [14the components of $\Lambda$  are allowed to be 1-tempered generalized functions. The authors replace $\Lambda$  by its regularization (the regularization procedure is similar to the one for $F$  ) which is 1-tempered generalized function of bounded growth and solve the regularized problem. In [4, 5the author investigates weak limits for semilinear hyperbolic systems and nonlinear superpositions for strongly singular distributions appearing in these systems.
He establishes an optimal link between the singularity of the initial data and the growth of the nonlinear term. Weak limits of strongly singular Cauchy problems for semilinear hyperbolic systems with bounded, sublinear, and superlinear growth are investigated in [3, 10, 18, 20.
Existence-uniqueness results within Colombeau algebras for two-dimensional hyperbolic problems with discontinuous coefficients ${\Lambda }_{i}$  are obtained in [9, 13, 16, 17. These papers impose an essential restriction on the coefficients, which allows one to avoide the negative effect of infinite propagation speed. Namely, the coefficients are assumed to be globally bounded in the Colombeau algebra $\mathcal{G}$  . At the present paper we do not assume global boundedness of the ${\Lambda }_{i}$  in ( 1 ), thereby allowing them to be strongly singular (this is also the case for [11).
A novelty of the paper is that it treats (within $\mathcal{G}\left(\overline{\Pi }\right)$  ) strongly singular initial data of the problem (including coefficients ${\Lambda }_{i}$  in ( 1 )), nonlinear boundary conditions, and non-Lipshitz nonlinearities in ( 1 ) and ( 3 ).
The plan of our exposition is as follows. In Section 2 we compile some facts about Colombeau algebra of generalized functions. In Section 3 we state and prove our main results. We prove the existence-uniqueness result within $\mathcal{G}\left(\overline{\Pi }\right)$  for Colombeau Lipshitz nonlinearities (Subsection 3.1) and extend its to the case of non-Lipshitz nonlinearities (Subsection 3.2).

2 Preliminaries

In this section we summarize the relevant material on the full version of Colombeau algebras of generalized functions.
Let $\Omega \subset {\mathbb{R}}^{n}$  be a domain in ${\mathbb{R}}^{n}$  . We denote by $\mathcal{G}\left(\Omega \right)$  and $\mathcal{G}\left(\overline{\Omega }\right)$  the full version of Colombeau algebra of generalized functions over $\Omega$  and $\overline{\Omega }$  , respectively. To define $\mathcal{G}\left(\Omega \right)$  and $\mathcal{G}\left(\overline{\Omega }\right)$  , we first introduce the mollifier spaces used to parametrize the regularizing sequences of generalized functions. Given $q\in {\mathbb{N}}_{0}$  denote
 $\begin{array}{c}{\mathcal{A}}_{q}\left(\mathbb{R}\right)=\left\{\phi \in \mathcal{D}\left(\mathbb{R}\right)|\int \phi \left(x\right)dx=1,\int {x}^{k}\phi \left(x\right)dx=0\text{for}1\le k\le q\right\},\end{array}$
 $\begin{array}{c}{\mathcal{A}}_{q}\left({\mathbb{R}}^{n}\right)=\left\{\phi \left({x}_{1},...,{x}_{n}\right)={\prod }_{i=1}^{n}{\phi }_{0}\left({x}_{i}\right)|{\phi }_{0}\in {\mathcal{A}}_{q}\left(\mathbb{R}\right)\right\}.\end{array}$
If $\phi \in {\mathcal{A}}_{0}\left({\mathbb{R}}^{n}\right)$  , let ${\phi }_{\varepsilon }\left(x\right)=\frac{1}{{\varepsilon }^{n}}\phi \left(\frac{x}{\varepsilon }\right).$  Set $\mathcal{ℰ}\left(\overline{\Omega }\right)=\left\{u:{\mathcal{A}}_{0}×\overline{\Omega }\to \mathbb{R}|u\left(\phi ,.\right)\in {\text{C}}^{\infty }\left(\overline{\Omega }\right)\forall \phi \in {\mathcal{A}}_{0}\left(\mathbb{R}\right)\right\}.$  We define the algebra of moderate elements $\mathcal{ℰ}\text{M}\left(\overline{\Omega }\right)$  to be the subalgebra of $\mathcal{ℰ}\left(\overline{\Omega }\right)$  consisting of the elements $u\in \mathcal{ℰ}\left(\overline{\Omega }\right)$  such that
 $\begin{array}{c}\forall K\subset \overline{\Omega }\text{compact},\forall \alpha \in {\mathbb{N}}_{0}^{n},\exists N\in \mathbb{N}\text{such that}\forall \phi \in {\mathcal{A}}_{N}\left({\mathbb{R}}^{n}\right)\end{array}$
 $\begin{array}{c}\exists C>0,\exists \eta >0\text{with}{sup}_{x\in K}|{\partial }^{\alpha }u\left({\phi }_{\varepsilon },x\right)|\le C{\varepsilon }^{-N},0<\varepsilon <\eta .\end{array}$
The ideal $\mathcal{N}\left(\overline{\Omega }\right)$  consists of all $u\in \mathcal{ℰ}\text{M}\left(\overline{\Omega }\right)$  such that
 $\begin{array}{c}\forall K\subset \overline{\Omega }\text{compact},\forall \alpha \in {\mathbb{N}}_{0}^{n},\exists N\in \mathbb{N}\text{such that}\forall q\ge N,\forall \phi \in {\mathcal{A}}_{q}\left({\mathbb{R}}^{n}\right)\end{array}$
 $\begin{array}{c}\exists C>0,\exists \eta >0\text{with}{sup}_{x\in K}|{\partial }^{\alpha }u\left({\phi }_{\varepsilon },x\right)|\le C{\varepsilon }^{q-N},0<\varepsilon <\eta .\end{array}$
Finally, $\mathcal{G}\left(\overline{\Omega }\right)=\mathcal{ℰ}\text{M}\left(\overline{\Omega }\right)/\mathcal{N}\left(\overline{\Omega }\right).$  This is an associative and commutative differential algebra. The algebra $\mathcal{G}\left(\Omega \right)$  on open set is constructed in the same manner, with $\Omega$  in place of $\overline{\Omega }$  . Note that $\mathcal{G}\left(\Omega \right)$  admits a canonical embedding of ${\mathcal{D}}^{\prime }\left(\Omega \right)$  . We will use the notation $U=\left[\left(u\left(\phi ,x\right){\right)}_{\phi \in {\mathcal{A}}_{0}\left({\mathbb{R}}^{n}\right)}\right]$  for the elements $U$  of $\mathcal{G}\left(\Omega \right)$  with representative $u\left(\phi ,x\right)$  .
One of the advantages of using Colombeau algebra of generalized functions $\mathcal{G}$  lies in the fact that in a variety of important cases the division by generalized functions, in particular the division by discontinuous functions and measures, is defined in $\mathcal{G}$  .
Complete description of the cases when the division is possible in the full version of Colombeau algebras is given by the following criterion of invertibility [11(the criterion of invertibility for special version of Colombeau algebras ${\mathcal{G}}_{s}\left(\Omega \right)$  is proved in [7):
Theorem 1 Let $U\in \mathcal{G}\left(\Omega \right)$  ( $U\in \mathcal{G}\left(\overline{\Omega }\right)$  ). Then the following two conditions are equivalent:
(i) $U$  is invertible in $\mathcal{G}\left(\Omega \right)$  (in $\mathcal{G}\left(\overline{\Omega }\right)$  ), i.e., there exists $V\in \mathcal{G}\left(\Omega \right)$  ( $V\in \mathcal{G}\left(\overline{\Omega }\right)$  ) such that $UV=1$  in $\mathcal{G}\left(\Omega \right)$  (in $\mathcal{G}\left(\overline{\Omega }\right)$  ).
(ii) For each representative $\left(u\left(\phi ,x\right){\right)}_{\phi \in {\mathcal{A}}_{0}\left({\mathbb{R}}^{n}\right)}$  of $U$  and each compact set $K\subset \Omega$  ( $K\subset \overline{\Omega }$  ) there exists $p\in \mathbb{N}$  such that for all $\phi \in {\mathcal{A}}_{p}\left({\mathbb{R}}^{n}\right)$  there is $\eta >0$  with ${inf}_{K}|u\left({\phi }_{\varepsilon },x\right)|\ge {\varepsilon }^{p}$  for all $0<\varepsilon <\eta$  .

3 Existence-uniqueness results in the Colombeau algebra of generalized functions

3.1 Colombeau Lipshitz nonlinearities

We here develop some results of [11and [12to the case of nonlinear nonlocal boundary conditions and Colombeau Lipshitz nonlinearities in ( 1 ) and ( 3 ) (with Lipshitz constants as Colombeau generalized numbers). We will need a notion of a generalized function whose growth is more restrictive than the $1/\varepsilon$  -growth (as in the definition of $\mathcal{ℰ}\text{M}$  ).
Definition 2 ([11) Let $\Omega \subset {\mathbb{R}}^{n}$  be a domain in ${\mathbb{R}}^{n}$  . Suppose we have a function $\gamma :\left(0,1\right)↦\left(0,\infty \right)$  . We say that an element $U\in \mathcal{G}\left(\Omega \right)$  ( $U\in \mathcal{G}\left(\overline{\Omega }\right)$  ) is locally of $\gamma$  -growth, if it has a representative $u\in \mathcal{ℰ}\text{M}\left(\Omega \right)$  ( $u\in \mathcal{ℰ}\text{M}\left(\overline{\Omega }\right)$  ) with the following property:
For every compact set $K\subset \Omega$  ( $K\subset \overline{\Omega }$  ) there is $N\in \mathbb{N}$  such that for every $\phi \in {\mathcal{A}}_{N}\left({\mathbb{R}}^{n}\right)$  there exist $C>0$  and $\eta >0$  with ${sup}_{x\in K}|u\left({\phi }_{\varepsilon },x\right)|\le C{\gamma }^{N}\left(\varepsilon \right)$  for $0<\varepsilon <\eta$  .
We now make assumptions on the initial data of the problem ( 1 )–( 3 ). Let $\gamma \left(\varepsilon \right)$  and ${\gamma }_{1}\left(\varepsilon \right)$  be functions from $\left(0,1\right)$  to $\left(0,\infty \right)$  such that
 $\begin{array}{c}\gamma \left(\varepsilon {\right)}^{{\gamma }^{N}\left(\varepsilon \right)}=O\left(\frac{1}{\varepsilon }\right),\gamma \left(\varepsilon {\right)}^{{\gamma }_{1}^{N}\left(\varepsilon \right)}=O\left(\frac{1}{\varepsilon }\right)\text{as}\varepsilon \to 0\end{array}$ (4)
for each $N\in \mathbb{N}$  . Assume that
• 1. $\Lambda \left(x,t\right)\in \left(\mathcal{G}\left(\overline{\Pi }\right){\right)}^{n}$  , $A\left(x\right)\in \left(\mathcal{G}\left[0,l\right]{\right)}^{n}$  .
• 2. ${\Lambda }_{i}$  for $i\le n$  are locally of $\gamma$  -growth on $\overline{\Pi }$  and invertible on $\overline{\Pi }$  .
• 3. ${\partial }_{x}{\Lambda }_{i}$  for $i\le n$  are locally of ${\gamma }_{1}$  -growth on $\overline{\Pi }$  .
• 4. $F\left(x,t,y\right)\in \left(\mathcal{G}\left(\overline{\Pi }×{\mathbb{R}}^{n}\right){\right)}^{n}$  , $H\left(t,z\right)\in \left(\mathcal{G}\left(\left[0,\infty \right)×{\mathbb{R}}^{n}\right){\right)}^{n}$  .
• 5. For every compact set $K\subset \overline{\Pi }$  and $i\le n$  the mapping $y↦{F}_{i}\left(x,t,y\right)$  and all its derivatives are polynomially bounded for all $\left(x,t\right)\in K$  with coefficients in $\mathcal{G}\left(K\right)$  .
• 6. For every compact set $K\subset \left[0,\infty \right)$  and $i\le n$  the mapping $z↦{H}_{i}\left(t,z\right)$  and all its derivatives are polynomially bounded for all $t\in K$  with coefficients in $\mathcal{G}\left(K\right)$  .
• 7. For every compact set $K\subset \overline{\Pi }$  there exists a nonnegative generalized function ${L}_{F}\left(x,t\right)\in \mathcal{G}\left(K\right)$  such that for all $\left(x,t\right)\in K$  , $i\le n$  , and ${y}^{1},{y}^{2}\in {\mathbb{R}}^{n}$  we have $|{F}_{i}\left(x,t,{y}^{1}\right)-{F}_{i}\left(x,t,{y}^{2}\right)|\le {L}_{F}\left(x,t\right){\sum }_{j=1}^{n}|{y}_{j}^{1}-{y}_{j}^{2}|.$
• 8. For every compact set $K\subset \left[0,\infty \right)$  there exists a nonnegative generalized function ${L}_{H}\left(x,t\right)\in \mathcal{G}\left(K\right)$  such that for all $t\in K$  , $i\le n$  , and ${z}^{1},{z}^{2}\in {\mathbb{R}}^{n}$  we have $|{H}_{i}\left(t,{z}^{1}\right)-{H}_{i}\left(t,{z}^{2}\right)|\le {L}_{H}\left(t\right){\sum }_{j=1}^{n}|{z}_{j}^{1}-{z}_{j}^{2}|.$
• 9. ${L}_{F}\left(x,t\right)$  and ${L}_{H}\left(t\right)$  are locally of $\gamma$  -growth on $\overline{\Pi }$  and $\left[0,\infty \right)$  , respectively.
• 10. $\text{supp}{A}_{i}\left(x\right)\subset \left(0,l\right)$  and $\text{supp}{H}_{i}\left(t,0\right)\subset \left(0,\infty \right)$  for $i\le n$  .
Assumptions imposed on ${\Lambda }_{i}$  allow them to be strongly singular and, even more, to have any desired order of singularity. Assumptions 4–6 state that, given $U\in \left(\mathcal{G}\left(\overline{\Pi }\right){\right)}^{n}$  and $V\in \left(\mathcal{G}\left[0,\infty \right){\right)}^{n}$  , $F\left(x,t,U\right)$  and $H\left(t,V\right)$  are well defined in Colombeau algebra $\mathcal{G}$  . We can interpret Assumptions 7 and 8 as the Lipshitz conditions in Colombeau sense imposed on generalized functions $F$  and $H$  . Assumption 9 allows ${L}_{F}$  and ${L}_{H}$  to be strongly singular. The last assumption ensures the compatibility of ( 2 ) and ( 3 ) of any desired order.
We now state the main result of this subsection.
Theorem 3 Under Assumptions 1–10 where the functions $\gamma$  and ${\gamma }_{1}$  are specified by ( 4 ), the problem ( 1 )–( 3 ) has a unique solution $U\in \mathcal{G}\left(\overline{\Pi }\right)$  .
Proof. We will first prove the existence of a classical smooth solution to the problem ( 1 )–( 3 ) where the initial data are smooth, satisfy Assumption 10, and the functions $F$  and $H$  have bounded gradients with respect to $U$  and $V$  , respectively, uniformly over $\left(x,t\right)$  varying in compact subsets of $\overline{\Pi }$  . In parallel, we will obtain a priori estimates for classical smooth solutions and their derivatives of any desired order. Therewith we will obtain the existence of a prospective representative $u$  of the solution $U$  in $\mathcal{G}\left(\overline{\Pi }\right)$  . To finish the existence part of the proof, we will show the moderateness of $u$  . The uniqueness of the constructed generalized solution will be proved by the same scheme.
We first reduce the problem ( 1 )–( 3 ) with smooth initial data to an equivalent integral-operator form. Denote by ${\omega }_{i}\left(\tau ;x,t\right)$  the $i$  -th characteristic of ( 1 ) passing through a point $\left(x,t\right)\in \overline{\Pi }$  . From the assumptions imposed on $\Lambda$  it follows that such characteristic exists, is smooth in $\tau ,x,t$  , and can be continued up to the boundary of $\Pi$  . The smallest value of $\tau \ge 0$  at which the characteristic $\xi ={\omega }_{i}\left(\tau ;x,t\right)$  intersects $\partial \Pi$  will be denoted by ${t}_{i}\left(x,t\right)$  . Integrating each equation of ( 1 ) along the corresponding characteristic curve, we obtain the following equivalent integral-operator form of ( 1 )–( 3 ):
 $\begin{array}{c}\begin{array}{c}{U}_{i}\left(x,t\right)=\left({R}_{i}U\right)\left(x,t\right)+{\int }_{{t}_{i}\left(x,t\right)}^{t}\left[U\left({\omega }_{i}\left(\tau ;x,t\right),\tau \right){\int }_{0}^{1}{\nabla }_{U}{F}_{i}\left({\omega }_{i}\left(\tau ;x,t\right),\tau ,\sigma U\right)d\sigma \\ +{F}_{i}\left({\omega }_{i}\left(\tau ;x,t\right),\tau ,0\right)\right]d\tau ,1\le i\le n,\end{array}\end{array}$ (5)
where $\left({R}_{i}U\right)\left(x,t\right)=\begin{cases}{V}_{i}\left({t}_{i}\left(x,t\right)\right){\int }_{0}^{1}{\nabla }_{V}{H}_{i}\left({t}_{i}\left(x,t\right),\sigma V\right)d\sigma +{H}_{i}\left({t}_{i}\left(x,t\right),0\right) & \text{if \text{if}{t}_{i}\left(x,t\right)>0\text{,}}\\ {A}_{i}\left({\omega }_{i}\left(0;x,t\right)\right) & \text{if \text{if}{t}_{i}\left(x,t\right)=0\text{.}}\end{cases}$  Given $T>0$  , denote ${\Pi }^{T}=\left\{\left(x,t\right)|0  Set ${E}_{U}\left({\alpha }_{1},{\alpha }_{2};T\right)=max\left\{|{\partial }_{x}^{{\alpha }_{1}}{\partial }_{t}^{{\alpha }_{2}}{U}_{i}\left(x,t\right)||\left(x,t\right)\in {\overline{\Pi }}^{T},1\le i\le n\right\},$  ${E}_{F}\left({\alpha }_{1},{\alpha }_{2}\right)=max\left\{|{\partial }_{x}^{{\alpha }_{1}}{\partial }_{t}^{{\alpha }_{2}}{F}_{i}\left(x,t,y\right)||\left(x,t,y\right)\in {\overline{\Pi }}^{T}×\left\{y:|y|\le {E}_{U}\left(0,0;T\right)\right\},1\le i\le n\right\},$  ${E}_{H}\left(\alpha \right)=max\left\{|{\partial }_{t}^{\alpha }{H}_{i}\left(t,z\right)||\left(t,z\right)\in \left[0,T\right]×\left\{z:|z|\le {E}_{U}\left(0,0;T\right)\right\},1\le i\le n\right\},$  ${L}_{F}^{max}=max\left\{{L}_{F}\left(x,t\right)|\left(x,t\right)\in {\overline{\Pi }}^{T}\right\},{L}_{H}^{max}=max\left\{{L}_{H}\left(t\right)|t\in \left[0,T\right]\right\}.$  Simplifying the notation, we drop the dependence of ${E}_{F}\left({\alpha }_{1},{\alpha }_{2}\right)$  , ${E}_{H}\left(\alpha \right)$  , ${L}_{F}^{max}$  , and ${L}_{H}^{max}$  on $T$  . Note that we will use these parameters for a fixed $T>0$  .
Assume that the initial data $\Lambda$  , $F$  , $A$  , and $H$  of our problem are smooth with respect to all their arguments, satisfy Assumption 10, and the functions $F$  and $H$  have bounded gradients with respect to $U$  and $V$  , uniformly over $\left(x,t\right)$  varying in compact subsets of $\overline{\Pi }$  . Fix an arbitrary $T>0$  . If $\left(x,t\right)\in {\overline{\Pi }}^{T}$  , then ${L}_{F}^{max}$  and ${L}_{H}^{max}$  are Lipshitz constants of $F$  and $H$  with respect to $U$  and $V$  , respectively. We now prove that the problem ( 5 ) has a smooth solution in ${\overline{\Pi }}^{T}$  . In parallel, we obtain global a priori estimates for smooth solutions, we will make use of for construction of a Colombeau solution. We obtain the global a priori estimates by iterating the a apriori estimates for local smooth solutions in a number of steps. The proof is split in four claims.
Claim 1. The problem ( 5 ) has a unique continuous solution in ${\overline{\Pi }}^{T}$  . We start from the local continuous solution to ( 5 ), namely, we state that there exists a unique solution $U\in \left(\text{C}\left({\overline{\Pi }}^{{t}_{0}}\right){\right)}^{n}$  to the problem ( 5 ) for some ${t}_{0}>0$  . To prove this, choose ${t}_{0}$  satisfying the condition
 $\begin{array}{c}{\omega }_{n}\left(t;0,\tau \right)<{\omega }_{1}\left(t;l,\tau \right)\forall \tau \ge 0,\forall t\in \left[\tau ,\tau +{t}_{m}\right]\end{array}$ (6)
with $m=0$  . For $t\in \left[0,{t}_{0}\right]$  we can express $V\left(t\right)$  in the form
 $\begin{array}{c}\begin{array}{c}{V}_{i}\left(t\right)={A}_{i}\left({\omega }_{i}\left(0;0,t\right)\right)+{\int }_{0}^{t}\left[U\left({\omega }_{i}\left(\tau ;0,t\right),\tau \right){\int }_{0}^{1}{\nabla }_{U}{F}_{i}\left({\omega }_{i}\left(\tau ;0,t\right),\tau ,\sigma U\right)d\sigma \\ +{F}_{i}\left({\omega }_{i}\left(\tau ;0,t\right),\tau ,0\right)\right]d\tau ,1\le i\le k,\\ {V}_{i}\left(t\right)={A}_{i}\left({\omega }_{i}\left(0;l,t\right)\right)+{\int }_{0}^{t}\left[U\left({\omega }_{i}\left(\tau ;l,t\right),\tau \right){\int }_{0}^{1}{\nabla }_{U}{F}_{i}\left({\omega }_{i}\left(\tau ;l,t\right),\tau ,\sigma U\right)d\sigma \\ +{F}_{i}\left({\omega }_{i}\left(\tau ;l,t\right),\tau ,0\right)\right]d\tau ,k+1\le i\le n.\end{array}\end{array}$ (7)
Since ( 5 ) is a system of Volterra integral equations of the second kind in ${\overline{\Pi }}^{{t}_{0}}$  , we can apply the contraction mapping principle. We apply the operator defined by the right hand side of ( 5 ) to two continuous functions ${U}^{1}$  and ${U}^{2}$  . Note that these functions have the same initial and boundary values. We cosider their difference in ${\overline{\Pi }}^{{t}_{0}}$  . Notice the estimate ${E}_{{U}^{1}-{U}^{2}}\left(0,0;{t}_{0}\right)\le {t}_{0}{q}_{0}{E}_{{U}^{1}-{U}^{2}}\left(0,0;{t}_{0}\right),$  where ${q}_{0}=n{L}_{F}^{max}\left(1+n{L}_{H}^{max}\right).$  We are able to choose ${t}_{0}$  so that the additional condition ${t}_{0}<1/{q}_{0}$  is obeyed. Then the contraction property of the operator defined by the right hand side of ( 5 ) holds with respect to ${\overline{\Pi }}^{{t}_{0}}$  . We have thus proved existence and uniqueness of a continuous solution $U$  to the problem ( 5 ) in ${\overline{\Pi }}^{{t}_{0}}$  . Furthermore, we have the following local a priori estimate:
 $\begin{array}{c}\begin{array}{c}{E}_{U}\left(0,0;{t}_{0}\right)\le \frac{1}{1-{q}_{0}{t}_{0}}\left[\left({max}_{x\in \left[0,l\right],1\le i\le n}|{A}_{i}\left(x\right)|\\ +T{max}_{\left(x,t\right)\in {\overline{\Pi }}^{T},1\le i\le n}|{F}_{i}\left(x,t,0\right)|\right)\left(1+n{L}_{H}^{max}\right)+{max}_{t\in \left[0,T\right],1\le i\le n}|{H}_{i}\left(t,0\right)|\right].\end{array}\end{array}$ (8)
Note that the value of ${q}_{0}$  depends on $T$  and does not depend on ${t}_{0}$  . This allows us to complete the proof of the claim in $⌈T/{t}_{0}⌉$  steps, iterating local existence-uniqueness result in domains $\left({\Pi }^{j{t}_{0}}\cap {\Pi }^{T}\right)\{\overline{\Pi }}^{\left(j-1\right){t}_{0}},1\le j\le ⌈T/{t}_{0}⌉.$  Moreover, using the estimate ( 8 ) $⌈T/{t}_{0}⌉$  times and each time starting with the final value of $U$  from the previous step, we derive the following global a priori estimate:
 $\begin{array}{c}\begin{array}{c}{E}_{U}\left(0,0;T\right)\le {P}_{1,0}\left(\frac{1}{1-{q}_{0}{t}_{0}},n,{L}_{H}^{max}\right)\\ ×{P}_{2,0}\left({max}_{x\in \left[0,l\right],1\le i\le n}|{A}_{i}\left(x\right)|,{max}_{\left(x,t\right)\in {\overline{\Pi }}^{T},1\le i\le n}|{F}_{i}\left(x,t,0\right)|,{max}_{t\in \left[0,T\right],1\le i\le n}|{H}_{i}\left(t,0\right)|\right),\end{array}\end{array}$ (9)
where ${P}_{1,0}$  is a polynomial of degree $3⌈T/{t}_{0}⌉$  with all coefficients identically equal to 1 and ${P}_{2,0}$  is a polynomial of the first degree with positive constant coefficients depending only on $T$  .
Claim 2. The problem ( 1 )–( 3 ) has a unique ${\text{C}}^{1}$  -solution in ${\overline{\Pi }}^{T}$  . The proof is similar to the proof of Claim 1. Let us consider the initial-boundary problem for ${\partial }_{x}U$  :
 $\begin{array}{c}\begin{array}{c}{\partial }_{x}{U}_{i}\left(x,t\right)=\left({R}_{ix}^{{}^{\prime }}U\right)\left(x,t\right)+{\int }_{{t}_{i}\left(x,t\right)}^{t}\left[{\nabla }_{U}{F}_{i}\left(\xi ,\tau ,U\right)\cdot {\partial }_{x}U\\ -\left({\partial }_{x}{\Lambda }_{i}\right)\left(\xi ,\tau \right){\partial }_{x}{U}_{i}+\left({\partial }_{x}{F}_{i}\right)\left(\xi ,\tau ,U\right)\right]{|}_{\xi ={\omega }_{i}\left(\tau ;x,t\right)}d\tau ,1\le i\le n,\end{array}\end{array}$ (10)
where $\left({R}_{ix}^{{}^{\prime }}U\right)\left(x,t\right)=\begin{cases}{\Lambda }_{i}^{-1}\left(0,\tau \right)\left[{F}_{i}\left(0,\tau ,U\right)-{\nabla }_{V}{H}_{i}\left(\tau ,V\right)\cdot {V}^{{}^{\prime }}\left(\tau \right)-\left({\partial }_{t}{H}_{i}\right)\left(\tau ,V\right)\right]{|}_{\tau ={t}_{i}\left(x,t\right)} & \text{if \text{if}{t}_{i}\left(x,t\right)>0,k+1\le i\le n\text{,}}\\ {\Lambda }_{i}^{-1}\left(l,\tau \right)\left[{F}_{i}\left(l,\tau ,U\right)-{\nabla }_{V}{H}_{i}\left(\tau ,V\right)\cdot {V}^{{}^{\prime }}\left(\tau \right)-\left({\partial }_{t}{H}_{i}\right)\left(\tau ,V\right)\right]{|}_{\tau ={t}_{i}\left(x,t\right)} & \text{if \text{if}{t}_{i}\left(x,t\right)>0,1\le i\le k\text{,}}\\ {A}_{i}^{{}^{\prime }}\left({\omega }_{i}\left(0;x,t\right)\right) & \text{if \text{if}{t}_{i}\left(x,t\right)=0\text{.}}\end{cases}$  Choose ${t}_{1}$  satisfying the condition ( 6 ) with $m=1$  . Combining ( 1 ) with ( 10 ), we get
 $\begin{array}{c}\begin{array}{c}{V}_{i}^{{}^{\prime }}\left(t\right)={F}_{i}\left(0,t,U\right)-{\Lambda }_{i}\left(0,t\right)\left({\partial }_{x}{U}_{i}\right)\left(0,t\right)\\ ={F}_{i}\left(0,t,U\right)-{\Lambda }_{i}\left(0,t\right)\left[{A}_{i}^{{}^{\prime }}\left({\omega }_{i}\left(0;0,t\right)\right)+{\int }_{0}^{t}\left[{\nabla }_{U}{F}_{i}\left(\xi ,\tau ,U\right)\cdot {\partial }_{x}U\\ -\left({\partial }_{x}{\Lambda }_{i}\right)\left(\xi ,\tau \right){\partial }_{x}{U}_{i}+\left({\partial }_{x}{F}_{i}\right)\left(\xi ,\tau ,U\right)\right]{|}_{\xi ={\omega }_{i}\left(\tau ;0,t\right)}d\tau \right],1\le i\le k,\end{array}\end{array}$ (11)
where $t\in \left[0,{t}_{1}\right]$  . The functions ${V}_{i}^{{}^{\prime }}\left(t\right)$  for $k+1\le i\le n$  can be expressed in the same form. Using the fact that $U$  is a known continuous function (see Claim 1), we now apply the operator defined by the right hand side of ( 10 ) to two continuous functions ${\partial }_{x}{U}^{1}$  and ${\partial }_{x}{U}^{2}$  . Note that these functions have the same initial and boundary values. We cosider their difference in ${\overline{\Pi }}^{{t}_{1}}$  . Notice the estimate ${E}_{{U}^{1}-{U}^{2}}\left(1,0;{t}_{1}\right)\le {t}_{1}{q}_{1}{E}_{{U}^{1}-{U}^{2}}\left(1,0;{t}_{1}\right),$  where ${q}_{1}=\left(n{L}_{F}^{max}+{E}_{\Lambda }\left(1,0;T\right)\right)\left(1+n{L}_{H}^{max}\right).$  We are able to choose ${t}_{1}$  so that the additional condition ${t}_{1}<1/{q}_{1}$  is obeyed. This shows that the operator defined by the right hand side of ( 10 ) has the contraction property with respect to the domain ${\overline{\Pi }}^{{t}_{1}}$  . Thus, we have proved the existence and the uniqueness of a solution $U\in {\text{C}}_{x,t}^{1,0}\left({\overline{\Pi }}^{t\left(1\right)}\right)$  to the problem ( 5 ). Furthermore, we have the following local a priori estimate:
 $\begin{array}{c}\begin{array}{c}{E}_{U}\left(1,0;{t}_{1}\right)\le \frac{1}{1-{q}_{1}{t}_{1}}\left[\left({max}_{x\in \left[0,l\right],1\le i\le n}|{A}_{i}^{{}^{\prime }}\left(x\right)|+T{E}_{F}\left(1,0\right)\\ +{E}_{{\Lambda }^{-1}}\left(0,0;T\right){E}_{F}\left(0,0\right)\right)\left(1+n{L}_{H}^{max}\right)+{E}_{{\Lambda }^{-1}}\left(0,0;T\right){E}_{H}\left(1\right)\right].\end{array}\end{array}$ (12)
Using the fact that the value of ${q}_{1}$  depends on $T$  and does not depend on ${t}_{1}$  and iterating the local existence-uniqueness result in domains $\left({\Pi }^{j{t}_{1}}\cap {\Pi }^{T}\right)\{\overline{\Pi }}^{\left(j-1\right){t}_{1}},1\le j\le ⌈T/{t}_{1}⌉,$  we obtain the global a priori estimate:
 $\begin{array}{c}\begin{array}{c}{E}_{U}\left(1,0;T\right)\le {P}_{1,1}\left(\frac{1}{1-{q}_{1}{t}_{1}},n,{L}_{H}^{max}\right)\\ ×{P}_{2,1}\left({max}_{x\in \left[0,l\right],1\le i\le n}|{A}_{i}^{{}^{\prime }}\left(x\right)|,{max}_{0\le \alpha \le 1}{E}_{F}\left(\alpha ,0\right),{E}_{{\Lambda }^{-1}}\left(0,0;T\right),{E}_{H}\left(1\right)\right),\end{array}\end{array}$ (13)
where ${P}_{1,1}$  is a polynomial of degree $3⌈T/{t}_{1}⌉$  with all coefficients identically equal to 1 and ${P}_{2,1}$  is a polynomial of the second degree with positive constant coefficients depending only on $T$  .
The a priori estimate for ${E}_{U}\left(0,1;T\right)$  now follows from the system ( 1 ):
${E}_{U}\left(0,1;T\right)\le {E}_{F}\left(0,0\right)+{E}_{\Lambda }\left(0,0;T\right){E}_{U}\left(1,0;T\right),$  where ${E}_{U}\left(1,0;T\right)$  satisfies the estimate ( 13 ). This finishes the proof of the claim.
Claim 3. The problem ( 1 )–( 3 ) has a unique ${\text{C}}^{2}$  -solution in ${\overline{\Pi }}^{T}$  . Following the proof of Claims 1 and 2, let us consider the following problem for ${\partial }_{x}^{2}U$  :
 $\begin{array}{c}\begin{array}{c}{\partial }_{x}^{2}{U}_{i}\left(x,t\right)=\left({R}_{ixx}^{{}^{\prime \prime }}U\right)\left(x,t\right)+{\int }_{{t}_{i}\left(x,t\right)}^{t}\left[{\nabla }_{U}{F}_{i}\left(\xi ,\tau ,U\right)\cdot {\partial }_{x}^{2}U\\ -2\left({\partial }_{x}{\Lambda }_{i}\right)\left(\xi ,\tau \right){\partial }_{x}^{2}{U}_{i}-\left({\partial }_{x}^{2}{\Lambda }_{i}\right)\left(\xi ,\tau \right){\partial }_{x}{U}_{i}+\left({\partial }_{x}^{2}{F}_{i}\right)\left(\xi ,\tau ,U\right)+2{\nabla }_{U}\left({\partial }_{x}{F}_{i}\right)\left(\xi ,\tau ,U\right)\cdot {\partial }_{x}U\\ +{\nabla }_{U}\left({\nabla }_{U}{F}_{i}\left(\xi ,\tau ,U\right)\cdot {\partial }_{x}U\right)\cdot {\partial }_{x}U\right]{|}_{\xi ={\omega }_{i}\left(\tau ;x,t\right)}d\tau ,1\le i\le n,\end{array}\end{array}$ (14)
where $\left({R}_{ixx}^{{}^{\prime \prime }}U\right)\left(x,t\right)={A}_{i}^{{}^{\prime \prime }}\left({\omega }_{i}\left(0;x,t\right)\right)\text{if}{t}_{i}\left(x,t\right)=0$  and $\left({R}_{ixx}^{{}^{\prime \prime }}U\right)\left(x,t\right)=-\left({\partial }_{t}{\Lambda }_{i}^{-1}\right)\left(0,\tau \right)\left({\partial }_{x}{U}_{i}\right)\left(0,\tau \right)-{\Lambda }_{i}^{-2}\left(0,\tau \right)\left[\left({\partial }_{t}{F}_{i}\right)\left(0,\tau ,U\right)$  $+{\nabla }_{U}{F}_{i}\left(0,\tau ,U\right)\cdot \left({\partial }_{t}U\right)\left(0,\tau \right)-{\nabla }_{V}{H}_{i}\left(\tau ,V\right)\cdot {V}^{{}^{\prime \prime }}\left(\tau \right)-{V}^{{}^{\prime }}\left(\tau \right)\cdot \left({\nabla }_{V}\left({\partial }_{t}{H}_{i}\right)\left(\tau ,V\right)$  $-\left({\partial }_{t}^{2}{H}_{i}\right)\left(\tau ,V\right)-{\nabla }_{V}\left({\partial }_{t}{H}_{i}\right)\left(\tau ,V\right)+{\nabla }_{V}\left({\nabla }_{V}{H}_{i}\left(\tau ,V\right)\cdot {V}^{{}^{\prime }}\left(\tau \right)\right)\right)\right]$  $+{\Lambda }_{i}^{-1}\left(0,\tau \right)\left[-\left({\partial }_{x}{\Lambda }_{i}\right)\left(0,\tau \right)\left({\partial }_{x}{U}_{i}\right)\left(0,\tau \right)+\left({\partial }_{x}{F}_{i}\right)\left(0,\tau ,U\right)$  $+{\nabla }_{U}{F}_{i}\left(0,\tau ,U\right)\cdot \left({\partial }_{x}U\right)\left(0,\tau \right)\right]{|}_{\tau ={t}_{i}\left(x,t\right)}\text{if}{t}_{i}\left(x,t\right)>0,1\le i\le k.$  The expressions for $\left({R}_{ixx}^{{}^{\prime \prime }}U\right)\left(x,t\right)$  if ${t}_{i}\left(x,t\right)>0$  and $k+1\le i\le n$  are similar. The expressions for ${V}_{i}^{{}^{\prime \prime }}$  are derived from the system ( 1 ), using suitable differentiations:
${V}_{i}^{{}^{\prime \prime }}={\partial }_{t}^{2}{U}_{i}=-{\partial }_{t}{\Lambda }_{i}{\partial }_{x}{U}_{i}-{\Lambda }_{i}{\partial }_{x}{\partial }_{t}{U}_{i}+{\partial }_{t}{F}_{i}+{\nabla }_{U}{F}_{i}\cdot {\partial }_{t}U$  $=-{\partial }_{t}{\Lambda }_{i}{\partial }_{x}{U}_{i}+{\Lambda }_{i}\left({\partial }_{x}{\Lambda }_{i}{\partial }_{x}{U}_{i}+{\Lambda }_{i}{\partial }_{x}^{2}{U}_{i}-{\partial }_{x}{F}_{i}-{\nabla }_{U}{F}_{i}\cdot {\partial }_{x}U\right)+{\partial }_{t}{F}_{i}+{\nabla }_{U}{F}_{i}\cdot {\partial }_{t}U,$  where ${\partial }_{x}^{2}U$  satisfies ( 14 ) and the right hand side is considered restricted to $x=0$  if $1\le i\le k$  and to $x=l$  if $k+1\le i\le n$  . Take ${t}_{2}$  satisfying the condition ( 6 ) with $m=2$  and the inequality ${t}_{2}<1/{q}_{2}$  , where ${q}_{2}=\left(n{L}_{F}^{max}+2{E}_{\Lambda }\left(1,0;T\right)\right)\left(1+n{L}_{H}^{max}\right).$  This shows that the operator defined by the right hand side of ( 14 ) has the contraction property with respect to the domain ${\overline{\Pi }}^{{t}_{2}}$  . We therefore have the following global a priori estimate:
 $\begin{array}{c}\begin{array}{c}{E}_{U}\left(2,0;T\right)\le {P}_{1,2}\left(\frac{1}{1-{q}_{2}{t}_{2}},n,{L}_{H}^{max}\right)\\ ×{P}_{2,2}\left(n,{max}_{x\in \left[0,l\right],1\le i\le n}|{A}_{i}^{{}^{\prime \prime }}\left(x\right)|,{max}_{0\le {\alpha }_{1}+{\alpha }_{2}\le 2}{E}_{\Lambda }\left({\alpha }_{1},{\alpha }_{2};T\right),{max}_{0\le {\alpha }_{1}+{\alpha }_{2}\le 1}{E}_{{\Lambda }^{-1}}\left({\alpha }_{1},{\alpha }_{2};T\right),\\ {max}_{1\le |\beta |+{\alpha }_{1}+{\alpha }_{2}\le 2}{E}_{{\partial }_{U}^{|\beta |}F}\left({\alpha }_{1},{\alpha }_{2}\right),{max}_{1\le |\beta |+{\alpha }_{1}\le 2}{E}_{{\partial }_{V}^{|\beta |}H}\left({\alpha }_{1}\right),{L}_{F}^{max},{L}_{H}^{max},{max}_{{\alpha }_{1}+{\alpha }_{2}=1}{E}_{U}\left({\alpha }_{1},{\alpha }_{2};T\right)\right),\end{array}\end{array}$ (15)
where $\beta =\left({\beta }_{1},...,{\beta }_{n}\right)$  , ${\beta }_{i}\in {\mathbb{N}}_{0}$  , $|\beta |={\beta }_{1}+...+{\beta }_{n}$  , ${P}_{1,2}$  is a polynomial of degree $3⌈T/{t}_{2}⌉$  with all coefficients identically equal to 1 and ${P}_{2,2}$  is a polynomial of the third degree with positive constant coefficients depending only on $T$  .
The estimates for ${E}_{U}\left(1,1;T\right)$  and ${E}_{U}\left(0,2;T\right)$  now follow from system ( 1 ) and its suitable differentiations. The claim is proved.
We further proceed by induction. Assume that the problem ( 5 ) has a unique ${\text{C}}^{m-1}\left({\overline{\Pi }}^{T}\right)$  -solution for an arbitrary $m\in \mathbb{N}$  .
Claim 4. The problem ( 1 )–( 3 ) has a unique ${\text{C}}^{m}\left({\overline{\Pi }}^{T}\right)$  -solution. Similarly to the proof of Claims 1–3, we consider the problem for ${\partial }_{x}^{m}U$  (using suitable differentiations and integrations of ( 1 )–( 3 )). Taking into account the induction assumption, we apply the contraction mapping principle to the operator defined by the right-hand side of the problem for ${\partial }_{x}^{m}U$  . It is not difficult to prove that the operator has the contraction property with respect to ${\overline{\Pi }}^{{t}_{m}}$  for some ${t}_{m}>0$  satisfying the condition ( 6 ) and the inequality ${t}_{m}<1/{q}_{m}$  , where
 $\begin{array}{c}{q}_{m}=\left(n{L}_{F}^{max}+m{E}_{\Lambda }\left(1,0;T\right)\right)\left(1+n{L}_{H}^{max}\right).\end{array}$ (16)
This implies the existence and the uniqueness of a ${\text{C}}_{x,t}^{m,0}\left({\overline{\Pi }}^{{t}_{m}}\right)$  -solution to the problem ( 1 )–( 3 ).
Iterating this local existence-uniqueness result in domains $\left({\Pi }^{j{t}_{m}}\cap {\Pi }^{T}\right)\{\overline{\Pi }}^{\left(j-1\right){t}_{m}},1\le j\le ⌈T/{t}_{m}⌉,$  we complete the proof of the claim in $⌈T/{t}_{m}⌉$  number of steps. In parallel, we arrive at the following global estimate:
 $\begin{array}{c}\begin{array}{c}{E}_{U}\left(m,0;T\right)\le {P}_{1,m}\left(\frac{1}{1-{q}_{m}{t}_{m}},n,{L}_{H}^{max}\right)\\ ×{P}_{2,m}\left(n,{max}_{x\in \left[0,l\right],1\le i\le n}|{A}_{i}^{\left(m\right)}\left(x\right)|,{max}_{0\le {\alpha }_{1}+{\alpha }_{2}\le m-1}{E}_{{\Lambda }^{-1}}\left({\alpha }_{1},{\alpha }_{2};T\right),\\ {max}_{0\le {\alpha }_{1}+{\alpha }_{2}\le m}{E}_{\Lambda }\left({\alpha }_{1},{\alpha }_{2};T\right),{max}_{1\le |\beta |+{\alpha }_{1}+{\alpha }_{2}\le m}{E}_{{\partial }_{U}^{|\beta |}F}\left({\alpha }_{1},{\alpha }_{2}\right),{max}_{1\le |\beta |+{\alpha }_{1}\le m}{E}_{{\partial }_{V}^{|\beta |}H}\left({\alpha }_{1}\right),\\ {L}_{F}^{max},{L}_{H}^{max},{max}_{1\le {\alpha }_{1}+{\alpha }_{2}\le m-1}{E}_{U}\left({\alpha }_{1},{\alpha }_{2};T\right)\right).\end{array}\end{array}$ (17)
Here ${P}_{1,m}$  is a polynomial of degree $3⌈T/{t}_{m}⌉$  with all coefficients identically equal to 1. Furthermore, ${P}_{2,m}$  is a polynomial whose degree depends on $m$  but neither on $T$  nor on ${t}_{m}$  and whose coefficients are positive constants depending only on $m$  and $T$  . The existence and uniqueness of a ${\text{C}}_{x,t}^{{\alpha }_{1},{\alpha }_{2}}\left({\overline{\Pi }}^{T}\right)$  -solution where ${\alpha }_{1}+{\alpha }_{2}=m$  , now follow from the system ( 1 ) and its suitable differentiations. The respective global a priori estimates for ${E}_{U}\left({\alpha }_{1},{\alpha }_{2};T\right)$  one can easily obtain from the inequality ( 17 ) and the induction assumption. The claim is proved.
The classical smooth solution to the problem ( 1 )–( 3 ) satisfying estimates ( 17 ) in ${\overline{\Pi }}^{T}$  for any $m\in {\mathbb{N}}_{0}$  can be constructed by the sequential approximation method.
We now use this solution to construct a representative of the Colombeau solution.
According to the assumptions of the theorem, we consider all the initial data as elements of the corresponding Colombeau algebras. We choose representatives $\lambda$  , $a$  , $f$  , $h$  , ${L}_{f}$  , and ${L}_{h}$  of $\Lambda$  , $A$  , $F$  , $H$  , ${L}_{F}$  , and ${L}_{H}$  , respectively, with the properties required in the theorem. Let $\phi =\phi \otimes \phi \in {\mathcal{A}}_{0}\left({\mathbb{R}}^{2}\right)$  . Consider a prospective representative $u=u\left(\phi ,x,t\right)$  of $U$  which is the classical smooth solution to the problem ( 1 )–( 3 ) with the initial data $\lambda \left(\phi ,x,t\right)$  , $a\left(\phi ,x\right)$  , $f\left(\phi ,x,t,u\left(\phi ,x,t\right)\right)$  , $h\left(\phi ,t,v\left(\phi ,t\right)\right)$  , ${L}_{f}\left(\phi ,x,t\right)$  , ${L}_{h}\left(\phi ,t\right)$  , where $v\left(\phi ,t\right)=\left({u}_{1}\left(\phi ,0,t\right),...,{u}_{k}\left(\phi ,0,t\right)$  , ${u}_{k+1}\left(\phi ,l,t\right),...,$  ${u}_{n}\left(\phi ,l,t\right)\right)$  .
For the existence part of the proof, we have to show that $u\in \mathcal{ℰ}\text{M}$  , i.e. to obtain moderate growth estimates of $u\left({\phi }_{\varepsilon },x,t\right)$  in terms of the regularization parameter $\varepsilon$  .
Fix $N\in \mathbb{N}$  to be so large that for all $\phi \in {\mathcal{A}}_{N}\left(\mathbb{R}\right)$  there exists ${\varepsilon }_{0}$  such that for all $\varepsilon <{\varepsilon }_{0}$  the following conditions are true:
a) The moderate estimate (see the definition of $\mathcal{ℰ}\text{M}$  ) holds for $a\left({\phi }_{\varepsilon },x\right)$  , $f\left({\phi }_{\varepsilon },x,t,0\right)$  , $h\left({\phi }_{\varepsilon },t,0\right)$  , ${L}_{f}\left({\phi }_{\varepsilon },x,t\right)$  , and ${L}_{h}\left({\phi }_{\varepsilon },t\right)$  .
b) The invertibility estimate (see Theorem  1 ) holds for $\lambda \left({\phi }_{\varepsilon },x,t\right)$  .
c) The local- $\gamma$  -growth estimate (see Definition  2 ) holds for $\lambda \left({\phi }_{\varepsilon },x,t\right)$  , ${L}_{f}\left({\phi }_{\varepsilon },x,t\right)$  , and ${L}_{h}\left({\phi }_{\varepsilon },t\right)$  .
d) The local ${\gamma }_{1}$  -growth estimate holds for ${\partial }_{x}\lambda \left({\phi }_{\varepsilon },x,t\right)$  .
Fix $\phi \in {\mathcal{A}}_{N}\left(\mathbb{R}\right)$  . Let ${p}_{1,m}\left(\phi \right)$  , ${p}_{2,m}\left(\phi \right)$  , and ${q}_{m}\left(\phi \right)$  denote the value of, respectively, ${P}_{1,m}$  , ${P}_{2,m}$  , and ${q}_{m}$  , where $U\left(x,t\right)$  , $\Lambda \left(x,t\right)$  , $A\left(x\right)$  , $F\left(x,t,U\left(x,t\right)\right)$  , $H\left(t,V\left(t\right)\right)$  , ${L}_{F}\left(x,t\right)$  , and ${L}_{H}\left(t\right)$  are replaced by their representatives $u\left(\phi ,x,t\right)$  , $\lambda \left(\phi ,x,t\right)$  , $a\left(\phi ,x\right)$  , $f\left(\phi ,x,t,$  $u\left(\phi ,x,t\right)\right)$  , $h\left(\phi ,t,v\left(\phi ,t\right)\right)$  , ${L}_{f}\left(\phi ,x,t\right)$  , and ${L}_{h}\left(\phi ,t\right)$  , respectively. It suffices to prove the moderate estimates for ${p}_{1,m}\left(\phi \right)$  and ${p}_{2,m}\left(\phi \right)$  for all $m\in {\mathbb{N}}_{0}$  . The expression ( 16 ) and assumptions imposed on $\Lambda$  , $F$  , and $H$  make it obvious now that ${q}_{m}\left({\phi }_{\varepsilon }\right)\le {\gamma }^{2N+1}\left(\varepsilon \right)+{\gamma }_{1}^{2N}\left(\varepsilon \right)$  for all sufficiently small $\varepsilon$  . Since ${t}_{m}\le min\left\{L/{E}_{\Lambda }\left(0,0\right),1/{q}_{m}\right\}$  and ${E}_{\Lambda }\left(0,0\right)$  $\le {\gamma }^{N+1}\left(\varepsilon \right)$  , we can choose ${t}_{m}=1/\left[2\left({\gamma }^{2N+1}\left(\varepsilon \right)+{\gamma }_{1}^{2N}\left(\varepsilon \right)\right)\right]$  . Taking into account ( 4 ), for each $m\in {\mathbb{N}}_{0}$  and for all small enough $\varepsilon$  we have
 $\begin{array}{c}{\left(\frac{n{L}_{H}^{max}}{1-{q}_{m}{t}_{m}}\right)}^{3⌈T/{t}_{m}⌉}\le \gamma \left(\varepsilon {\right)}^{6\left(N+1\right)⌈T\left({\gamma }^{2N+1}\left(\varepsilon \right)+{\gamma }_{1}^{2N}\left(\varepsilon \right)\right)⌉}\end{array}$
 $\begin{array}{c}\le \left(\gamma \left(\varepsilon {\right)}^{{\gamma }^{2N+2}\left(\varepsilon \right)}\right)\left(\gamma \left(\varepsilon {\right)}^{{\gamma }_{1}^{2N+1}\left(\varepsilon \right)}\right)=O\left(\frac{1}{{\varepsilon }^{2}}\right)\text{as}\varepsilon \to 0.\end{array}$
We conclude that for each $m\in {\mathbb{N}}_{0}$  there exists $N\in \mathbb{N}$  such that for all $\phi \in {\mathcal{A}}_{N}\left(\mathbb{R}\right)$  we have
 $\begin{array}{c}\begin{array}{c}{P}_{1,m}\left({\phi }_{\varepsilon }\right)=O\left(\frac{1}{{\varepsilon }^{2}}\right)\text{as}\varepsilon \to 0.\end{array}\end{array}$ (18)
One can easily see now that for ${\alpha }_{1}={\alpha }_{2}=0$
 $\begin{array}{c}{E}_{{u}_{\varepsilon }}\left({\alpha }_{1},{\alpha }_{2};T\right)=O\left(\frac{1}{{\varepsilon }^{N}}\right)\text{as}\varepsilon \to 0\end{array}$ (19)
for all $\phi \in {\mathcal{A}}_{N}\left(\mathbb{R}\right)$  with large enough $N\in \mathbb{N}$  , where ${u}_{\varepsilon }=u\left({\phi }_{\varepsilon },x,t\right)$  . To prove similar estimates for all derivatives of $U$  , we use induction on $\alpha ={\alpha }_{1}+{\alpha }_{2}$  . Assuming ( 19 ) to hold for $\alpha \le m-1$  , let us show that ( 19 ) is true for $\alpha =m$  as well. Indeed, let $\phi \in {\mathcal{A}}_{N}\left(\mathbb{R}\right)$  with $N$  chosen so large that for all sufficiently small $\varepsilon$  the following conditions are true:
a) The moderate estimate holds for ${\partial }_{x}^{{\alpha }_{1}}{\partial }_{t}^{{\alpha }_{2}}{u}_{i}\left({\phi }_{\varepsilon },x,t\right)$  for $0\le {\alpha }_{1}+{\alpha }_{2}\le m-1$  and $1\le i\le n$  (the induction assumption).
b) The moderate estimate holds for ${max}_{0\le {\alpha }_{1}+{\alpha }_{2}\le m}{E}_{{\lambda }_{\varepsilon }}\left({\alpha }_{1},{\alpha }_{2};T\right)$  , ${max}_{0\le {\alpha }_{1}+{\alpha }_{2}\le m-1}{E}_{{\lambda }_{\varepsilon }^{-1}}\left({\alpha }_{1},{\alpha }_{2};T\right)$  , e ${a}^{\left(m\right)}\left({\phi }_{\varepsilon },x\right)$  , where ${\lambda }_{\varepsilon }=\lambda \left({\phi }_{\varepsilon },x,t\right)$  .
c) Given $U=\left[\left(u\left({\phi }_{\varepsilon },x,t\right){\right)}_{\phi \in {\mathcal{A}}_{0}\left({\mathbb{R}}^{2}\right)}\right]\in \mathcal{G}\left({\overline{\Pi }}^{T}\right)$  satisfying the estimate ( 19 ) for ${\alpha }_{1}={\alpha }_{2}=0$  , the moderate estimate holds for ${max}_{1\le |\beta |+{\alpha }_{1}+{\alpha }_{2}\le m}{E}_{{\partial }_{U}^{|\beta |}{f}_{\varepsilon }}\left({\alpha }_{1},{\alpha }_{2}\right)$  , ${max}_{1\le |\beta |+{\alpha }_{1}\le m}{E}_{{\partial }_{V}^{|\beta |}{h}_{\varepsilon }}\left({\alpha }_{1}\right)$  , where ${f}_{\varepsilon }=f\left({\phi }_{\varepsilon },x,t,{u}_{\varepsilon }\right)$  , ${h}_{\varepsilon }=h\left({\phi }_{\varepsilon },t,v\left({\phi }_{\varepsilon },t\right)\right)$  d) the invertibility estimate holds for $\lambda \left({\phi }_{\varepsilon },x,t\right)$  .
Since ${p}_{2,m}\left({\phi }_{\varepsilon }\right)$  is a polynomial whose degree does not depend on $\varepsilon$  , the moderateness of ${p}_{2,m}\left(\phi \right)$  becomes obvious. The moderateness of ${E}_{u}\left(m,0;T\right)$  are done by ( 18 ). The moderateness property of ${E}_{u}\left({\alpha }_{1},{\alpha }_{2};T\right)$  for all other ${\alpha }_{1}$  and ${\alpha }_{2}$  such that ${\alpha }_{1}+{\alpha }_{2}=m$  is a consequence of the moderateness of ${E}_{u}\left(m,0;T\right)$  , the system ( 1 ), its suitable differentiations, and the induction assumption.
Since $T>0$  is arbitrary, the existence part of the proof is complete.
The proof of the uniqueness part follows the same scheme. The only difference is that now we consider the problem with respect to the difference $U-W$  of two Colombeau solutions $U$  and $W$  . We hence have the problem ( 1 )–( 3 ) with the right hand sides ${\int }_{0}^{1}{\nabla }_{U}F\left(x,t,\sigma U+\left(1-\sigma \right)W\right)d\sigma \cdot \left(U-W\right)+{M}_{1},$  ${\int }_{0}^{1}{\nabla }_{V}H\left(t,\sigma V+\left(1-\sigma \right){V}_{W}\right)d\sigma \cdot \left(V-{V}_{W}\right)+{M}_{3},$  and ${M}_{2}$  in ( 1 ), ( 3 ), and ( 2 ), respectively. Here ${M}_{i}\in \mathcal{N}$  and ${V}_{W}$  is equal to $V$  , where $U$  is replaced by $W$  . The analysis is even simpler since, due to [6, it suffices to check the negligibility of $U-W$  at order zero. For this purpose we rewrite the estimate ( 9 ) with respect to the function $U-W$  and use Assumptions 7–9 and the fact that $\left[\left({p}_{2,0}\left(\phi \right){\right)}_{\phi \in {\mathcal{A}}_{0}\left(\mathbb{R}\right)}\right]\in \mathcal{N}$  . This finishes the proof. $\square$

3.2 Non-Lipshitz nonlinearities

We here extend the above existence-uniqueness result to the case of non-Lipshitz nonlinearities in ( 1 ) and ( 3 ). Set ${E}_{\nabla F}\left(U\right)=max\left\{|{\nabla }_{U}{F}_{i}\left(x,t,U\left(x,t\right)\right)|:\left(x,t\right)\in {\overline{\Pi }}^{T},1\le i\le n\right\},$  ${E}_{\nabla H}\left(V\right)=max\left\{|{\nabla }_{V}{H}_{i}\left(t,V\left(t\right)\right)|:t\in \left[0,T\right],1\le i\le n\right\}.$  Simplifying the notation, we drop the dependence of ${E}_{\nabla F}\left(U\right)$  and ${E}_{\nabla H}\left(V\right)$  on $T$  .
Note that we will use these parameters for a fixed $T>0$  .
To state the main result of this section, we suppose that at least one of the following two assumptions holds.
Assumption 11.
a) $H\left(t,V\right)$  is smooth in $t,V$  and the mapping $V↦{\nabla }_{V}H\left(t,V\right)$  is globally bounded, uniformly over $t$  varying in compact subsets of $\left[0,\infty \right)$  ; b) Given $T>0$  , there exists ${C}_{F}$  such that for all $1\le i\le n$  , $\left(x,t\right)\in {\overline{\Pi }}^{T}$  , and $y\in {\mathbb{R}}^{n}$  we have $|{\nabla }_{U}{F}_{i}\left(x,t,y\right)|\le {C}_{F}loglogD\left(x,t,y\right),$  where $D\left(x,t,y\right)$  with respect to $y$  is a polynomial with coefficients in $\mathcal{G}\left({\overline{\Pi }}^{T}\right)$  .
Assumption 12.
a) Given $T>0$  , there exists ${C}_{H}$  such that for all $1\le i\le n$  , $t\in \left[0,T\right]$  , and $z\in {\mathbb{R}}^{n}$  we have $|{\nabla }_{V}{H}_{i}\left(t,z\right)|\le {C}_{H}\left(loglogB\left(t,z\right){\right)}^{1/4},$  where $B\left(t,z\right)$  with respect to $z$  is a polynomial with coefficients in $\mathcal{G}\left[0,T\right]$  .
b) Assumption 3 is true with ${\gamma }_{1}\left(\varepsilon \right)=O\left(\left(loglog1/\varepsilon {\right)}^{1/4}\right)$  ; c) Given $T>0$  , there exists ${C}_{F}$  such that for all $1\le i\le n$  , $\left(x,t\right)\in {\overline{\Pi }}^{T}$  , and $y\in {\mathbb{R}}^{n}$  we have $|{\nabla }_{U}{F}_{i}\left(x,t,y\right)|\le {C}_{F}\left(loglogD\left(x,t,y\right){\right)}^{1/4},$  where $D\left(x,t,y\right)$  with respect to $y$  is a polynomial with coefficients in $\mathcal{G}\left({\overline{\Pi }}^{T}\right)$  .
Theorem 4 Assume that Assumption 11 or 12 is true. Under Assumptions 1–6 and 10 where the functions $\gamma$  and ${\gamma }_{1}$  are specified by ( 4 ), the problem ( 1 )–( 3 ) has a unique solution $U\in \mathcal{G}\left(\overline{\Pi }\right)$  .
Proof. In the proof we will use a modified notion of $\mathcal{ℰ}\text{M}\left(\overline{\Pi }\right)$  , namely, let $u\in \mathcal{ℰ}\text{M}\left(\overline{\Pi }\right)$  iff $u\in \mathcal{ℰ}\left(\overline{\Pi }\right)$  and for every compact set $K\subset \overline{\Pi }$  there is $N\in \mathbb{N}$  such that for every $\phi \in {\mathcal{A}}_{N}\left({\mathbb{R}}^{n}\right)$  there exists $\eta >0$  with ${sup}_{x\in K}|u\left({\phi }_{\varepsilon }×{\phi }_{\varepsilon },x,t\right)|\le {\gamma }^{N}\left(\varepsilon \right)$  for all $0<\varepsilon <\eta$  .
Fix an arbitrary $T>0$  . From the proof of Theorem  3 it follows that $\left[\left({p}_{2,0}\left(\phi \right){\right)}_{\phi \in {\mathcal{A}}_{0}\left(\mathbb{R}\right)}\right]$  is a Colombeau generalized number and hence has the moderateness property. This means that there exists ${N}_{1}\in \mathbb{N}$  such that for all $\phi \in {\mathcal{A}}_{{N}_{1}}\left(\mathbb{R}\right)$  there is $\eta \left(\phi \right)>0$  with
 $\begin{array}{c}|{p}_{2,0}\left({\phi }_{\varepsilon }\right)|\le {\varepsilon }^{-{N}_{1}},0<\varepsilon <\eta \left(\phi \right).\end{array}$ (20)
Without loss of generality we can assume that ${N}_{1}$  is so large that for all $\phi \in {\mathcal{A}}_{{N}_{1}}\left(\mathbb{R}\right)$  the zero-order moderateness property holds for the coefficients of the polynomial $D\left(x,t,y\right)$  (if Assumption 11 is fulfilled) or for the coefficients of the polynomials $D\left(x,t,y\right)$  and $B\left(t,z\right)$  (if Assumption 12 is fulfilled). To simplify notation, we can suppose that, given $\phi \in {\mathcal{A}}_{{N}_{1}}\left(\mathbb{R}\right)$  , the value of $\eta \left(\phi \right)$  in ( 20 ) is so small that the zero-order moderate estimates for the coefficients of $D$  and/or $B$  are true for all $\varepsilon <\eta \left(\phi \right)$  . Note that any $U\in \mathcal{G}\left(\overline{\Pi }\right)$  has the following property: there exists ${N}_{2}\in \mathbb{N}$  such that for all $\phi \in {\mathcal{A}}_{{N}_{1}+{N}_{2}}\left(\mathbb{R}\right)$  there is ${\varepsilon }_{0}\left(\phi \right)\le \eta \left(\phi \right)$  , where the value of $\eta \left(\phi \right)$  is the same as in ( 20 ), with
 $\begin{array}{c}{sup}_{{\overline{\Pi }}^{T}}|u\left({\phi }_{\varepsilon }×{\phi }_{\varepsilon },x,t\right)|\le {\varepsilon }^{-{N}_{1}-{N}_{2}},0<\varepsilon <{\varepsilon }_{0}\left(\phi \right),\end{array}$ (21)
with the constant ${N}_{1}$  being the same as in ( 20 ). Obviously, any increase of ${N}_{2}$  and any decrease of ${\varepsilon }_{0}\left(\phi \right)$  will keep this property true. This will allow us to adjust the values of ${N}_{2}$  and ${\varepsilon }_{0}\left(\phi \right)$  according to our purposes.
Following the proof of Theorem  3 , for all $\phi \in {\mathcal{A}}_{{N}_{1}+{N}_{2}}\left(\mathbb{R}\right)$  , we arrive at the estimates ( 9 ) and ( 17 ) with ${E}_{{u}_{\varepsilon }}\left(m,0\right)$  , ${E}_{\nabla {f}_{\varepsilon }}\left({u}_{\varepsilon }\right)$  , and ${E}_{\nabla {h}_{\varepsilon }}\left({u}_{\varepsilon }\right)$  in place of ${E}_{U}\left(m,0\right)$  , ${L}_{F}$  , and ${L}_{H}$  , respectively, where $0<\varepsilon <\eta \left(\phi \right)$  and the value of $\eta \left(\phi \right)$  is the same as in ( 20 ). Recall that $u$  , $f$  , and $h$  are representatives of $U$  , $F$  , and $H$  , respectively, and ${u}_{\varepsilon }\left(x,t\right)=u\left({\phi }_{\varepsilon },x,t\right)$  . On the account of these estimates, we will obtain the existence once we prove the following assertion:
$\left(\iota \right)$  the constant ${N}_{2}\in \mathbb{N}$  can be chosen so that for all $\phi \in {\mathcal{A}}_{{N}_{1}+{N}_{2}}\left(\mathbb{R}\right)$  there exists ${\varepsilon }_{0}\left(\phi \right)$  such that
 $\begin{array}{c}{\left[2n\left(1+{E}_{\nabla {h}_{\varepsilon }}\left({u}_{\varepsilon }\right)\right)\right]}^{7T\left(1+{E}_{\nabla {h}_{\varepsilon }}\left({u}_{\varepsilon }\right)\right)\left({E}_{\nabla {f}_{\varepsilon }}\left({u}_{\varepsilon }\right)+m{E}_{\Lambda }\left(1,0;T\right)\right)}\le {\varepsilon }^{-{N}_{2}},0<\varepsilon <{\varepsilon }_{0}\left(\phi \right),\end{array}$ (22)
whatsoever $u\left(\phi ×\phi ,x,t\right)\in \mathcal{ℰ}\left({\overline{\Pi }}^{T}\right)$  satisfying the inequality ( 21 ).
Let us prove Assertion $\left(\iota \right)$  using Assumption 11. Recall that at this point ${N}_{2}$  is a constant whose exact value will be fixed below. Fix $\phi \in {\mathcal{A}}_{{N}_{1}+{N}_{2}}\left(\mathbb{R}\right)$  . By ( 21 ) and Assumption 11, there exists ${N}_{3}\in \mathbb{N}$  for which the estimate ${E}_{\nabla {f}_{\varepsilon }}\left({u}_{\varepsilon }\right)\le {C}_{F}loglogd\left({\phi }_{\varepsilon },x,t,{u}_{\varepsilon }\right)\le {C}_{F}loglog{\varepsilon }^{-{N}_{3}},0<\varepsilon <{\varepsilon }_{0}\left(\phi \right),$  is true, where $d$  is a representative of $D$  . Furthermore, there exist ${C}_{1}>1$  , ${C}_{2}>0$  , and ${k}_{1},{k}_{2}\in \mathbb{N}$  such that the left hand side of ( 22 ) is bounded from above by ${C}_{1}^{{C}_{2}\left(loglog{\varepsilon }^{-{N}_{3}}+{\gamma }_{1}\left(\varepsilon \right)\right)}\le {e}^{{C}_{3}loglog{\varepsilon }^{-{N}_{3}}}\gamma \left(\varepsilon {\right)}^{{k}_{1}{\gamma }_{1}\left(\varepsilon \right)}\le {e}^{log\left(log{\varepsilon }^{-{N}_{3}}{\right)}^{{C}_{3}}}{\varepsilon }^{-{k}_{2}}$  $\le \left(log{\varepsilon }^{-{N}_{3}}{\right)}^{{C}_{3}}{\varepsilon }^{-{k}_{2}}\le {N}_{3}^{⌈{C}_{3}⌉}{\varepsilon }^{-⌈{C}_{3}⌉-{k}_{2}},$  where ${C}_{3}={C}_{2}log{C}_{1}$  and $0<\varepsilon <{\varepsilon }_{0}\left(\phi \right)$  . It is important to note that ${C}_{3}$  and ${k}_{2}$  can be fixed so that the above estimates hold for all ${N}_{2}$  and all $\phi$  . This makes the values ${N}_{2}=2⌈{C}_{3}⌉+{k}_{2}$  and ${\varepsilon }_{0}\left(\phi \right)=min\left\{\eta \left(\phi \right),{N}_{3}^{-⌈{C}_{3}⌉}\right\}$  , which we now set up, well defined. Assertion $\left(\iota \right)$  now follows from the fact that $\phi$  is an arbitrary function from ${\mathcal{A}}_{{N}_{1}+{N}_{2}}\left(\mathbb{R}\right)$  .
Let us prove Assertion $\left(\iota \right)$  using Assumption 12. Following the same scheme as above, fix $\phi \in {\mathcal{A}}_{{N}_{1}+{N}_{2}}\left(\mathbb{R}\right)$  , where ${N}_{2}$  will be specified below. By ( 21 ) and Assumption 12, there exist ${N}_{3},{N}_{4}\in \mathbb{N}$  such that the following estimates are true:
${E}_{\nabla {f}_{\varepsilon }}\left({u}_{\varepsilon }\right)\le {C}_{F}loglogd\left({\phi }_{\varepsilon },x,t,{u}_{\varepsilon }\right)\le {C}_{F}loglog\left({\varepsilon }^{-{N}_{3}}\right),0<\varepsilon <{\varepsilon }_{0}\left(\phi \right),$  ${E}_{\nabla {h}_{\varepsilon }}\left({u}_{\varepsilon }\right)\le {C}_{H}loglogb\left({\phi }_{\varepsilon },t,{u}_{\varepsilon }\right)\le {C}_{H}loglog\left({\varepsilon }^{-{N}_{4}}\right),0<\varepsilon <{\varepsilon }_{0}\left(\phi \right),$  where $b$  is a representative of $B$  . Furthermore, there exist ${C}_{1}>1$  and ${C}_{2}>0$  such that the left hand side of ( 22 ) is bounded from above by ${\left[{C}_{1}loglog\left({\varepsilon }^{-{N}_{4}}\right)\right]}^{1/2{C}_{2}\left(loglog{\varepsilon }^{-{N}_{3}-{N}_{4}}{\right)}^{1/2}}$  $\le exp\left\{{C}_{2}log\left(log\left(log\left({\varepsilon }^{-{N}_{4}}\right){\right)}^{{C}_{1}}{\right)}^{1/2}\left(loglog{\varepsilon }^{-{N}_{3}-{N}_{4}}{\right)}^{1/2}\right\}$  $\le exp\left\{{C}_{2}log\left(log{\varepsilon }^{-{N}_{3}-{N}_{4}}{\right)}^{{C}_{1}}\right\}=\left(log{\varepsilon }^{-{N}_{3}-{N}_{4}}{\right)}^{{C}_{1}{C}_{2}}$  $={\left(\left({N}_{3}+{N}_{4}\right)log{\varepsilon }^{-1}\right)}^{{C}_{1}{C}_{2}}\le \left({N}_{3}+{N}_{4}{\right)}^{⌈{C}_{1}{C}_{2}⌉}{\varepsilon }^{-⌈{C}_{1}{C}_{2}⌉},$  where $0<\varepsilon <{\varepsilon }_{0}\left(\phi \right)$  . Note that ${C}_{1}$  and ${C}_{2}$  can be fixed so that the above estimates hold for all ${N}_{2}$  and all $\phi$  . We now set ${N}_{2}=2⌈{C}_{1}{C}_{2}⌉$  and ${\varepsilon }_{0}\left(\phi \right)=min\left\{\eta \left(\phi \right),\left({N}_{3}+{N}_{4}{\right)}^{-⌈{C}_{1}{C}_{2}⌉}\right\}$  and this value is well defined. Assertion $\left(\iota \right)$  now follows from the fact that $\phi$  is an arbitrary function in ${\mathcal{A}}_{{N}_{1}+{N}_{2}}\left(\mathbb{R}\right)$  .
Since $T>0$  is arbitrary, the existence part of the proof is complete.
The proof of the uniqueness part follows the same scheme (cf. also the proof of Theorem  3 ). We apply the estimate ( 9 ) to the difference of two generalized solutions to the problem ( 1 )–( 3 ). From the existence part of the proof we see that the first factor in the right-hand side of ( 9 ) has the moderateness propery. Since the second factor is negligible, the uniqueness follows. $\square$
Example 5 Consider $n=1$  and $F\left(x,t,U\right)=\left({G}_{1}^{2}\left(x,t\right)+{G}_{2}^{2}\left(x,t\right){U}^{2}{\right)}^{1/2}loglog\left({G}_{3}^{2}\left(x,t\right)+{G}_{4}^{2}\left(x,t\right){U}^{2}{\right)}^{1/2},$  where ${G}_{i}\left(x,t\right)\in \mathcal{G}\left(\overline{\Pi }\right)$  . Then ${\partial }_{U}F\left(x,t,U\right)=\frac{{G}_{2}^{2}U}{\left({G}_{1}^{2}+{G}_{2}^{2}{U}^{2}{\right)}^{1/2}}loglog\left({G}_{3}^{2}+{G}_{4}^{2}{U}^{2}{\right)}^{1/2}$  $+\frac{{G}_{4}^{2}U\left({G}_{1}^{2}+{G}_{2}^{2}{U}^{2}{\right)}^{1/2}}{log\left({G}_{3}^{2}+{G}_{4}^{2}{U}^{2}{\right)}^{1/2}\left({G}_{3}^{2}+{G}_{4}^{2}{U}^{2}\right)}.$  The function $F\left(x,t,U\right)$  is non-Lipshitz and satisfies Assumption 11(b).
Remark 6 The theorem states that, whatsoever singularity of the initial data of our problem and whatsoever nonlinearities of $F$  and $H$  allowed by Assumption 11 (or 12), the problem ( 1 )–( 3 ) has a unique solution in the Colombeau algebra $\mathcal{G}\left(\overline{\Pi }\right)$  .
References

1. H. Biagioni. Lecture Notes in Mathematics: A Nonlinear Theory of Generalized Functions. Springer-Verlag, Berlin-Heidelberg-New York 1990.
2. Colombeau, J. F.: Elementary Introduction to New Generalized Functions (North-Holland Mathematics Studies: Vol. 113). Elsevier Science Publishers 1985.
3. Demengel, F. and Rauch, J.: Measure valued solutions of asymptotically homogeneous semilinear hyperbolic systems in one space dimension. Proc. Edinburgh Math. Soc. 33 (1990), 443–460.
4. Gramchev, T.: Semilinear hyperbolic systems with singular initial data. Monatshefte Math. 112 (1991), 99–113.
5. Gramchev, T.: Nonlinear maps in spaces of distributions. Math. Zeitschr. 209(1992), 101 114.
6. M. Grosser, On the foundations of nonlinear generalized functions II, X-Archives: math. FA/9912215, ESI Preprint 812. To appear as a single issue of Memoirs of the AMS.
7. M. Grosser, M. Kunzinger, M. Oberguggenberger, R. Steinbauer. Geometric theory of generalized functions. Kluwer Academic Publishers, Dordrecht, 2001.
8. F. Jochmann, L. Recke. Well-posedness of an initial boundary value problem from laser dynamics, Math. Models and Methods in Applied Sciences 12 (1999), No 4, 593–606.
9. G. Hörmann and M. V. de Hoop. Microlocal analysis and global solutions of some hyperbolic equations with discontinuous coefficients. Acta Appl. Math. 67 (2001), 173–224.
10. I. Kmit. Delta waves for a strongly singular initial-boundary hyperbolic problem with integral boundary condition. Accepted for publication in J.of Analysis and its Applications (2003), 42 pages.
11. I. Kmit. Initial-Boundary Problems for Semilinear Hyperbolic Systems with Singular Coefficients. Proceedings of the International Conference on Nonlinear Differential Equations. Nonlinear Boundary Value Problems (2004) Alushta, Ukraine, to appear.
12. I. Kmit and G. Hörmann. Semilinear hyperbolic systems with nonlocal boundary conditions: reflection of singularities and delta waves, J. for Analysis and its Applications 20 (2001), No. 3, 637–659.
13. F. Lafon and M. Oberguggenberger. Generalized solutions to symmetric hyperbolic systems with discontinuous coefficients: the multidimensional case, J. Math. Anal. Appl. 160 (1991), 93–106.
14. M. Nedeljkov and S. Pilipovic̆. A note on a semilinear hyperbolic system with generalized functions as coefficients. Nonlin. Analysis, Theory, Methods and Appl. 30 (1997), 41–46.
15. M. Nedeljkov and S. Pilipovic̆. Generalized solution to a semilinear hyperbolic system with a non-Lipshitz nonlinearity, Preprint.
16. M. Oberguggenberger. Multiplication of Distributions and Applications to Partial Differential Equations, volume 259 of Pitman Research Notes in Mathematics, Longman, 1992.
17. M. Oberguggenberger. Hyperbolic systems with discontinuous coefficients: generalized solutions and a transmission problem in acoustics, J. Math. Anal. Appl. 142 (1989), 452–467.
18. Oberguggenberger, M.: Weak limits of solutions to semilinear hyperbolic systems. Math. Ann. 274 (1986), 599–607.
19. M. Radziunas, H.-J. Wünsche, B. Sartorius, O. Brox, D. Hoffmann, K. Schneider, Modelling of self-pulsating DFB lasers, WIAS preprint 516 (1999), 14 pages.
20. Rauch, J. and Reed, M.: Nonlinear superposition and absorption of delta waves in one space dimension. Funct. Anal. 73 (1987), 152–178.
21. J. Sieber, U. Bandelow, H. Wenzel, M. Eolfrum, H.-J. Wünsche, Travelling wave equations for semiconductor lasers with gain dispersion, WIAS preprint 459 (1998).
22. B. Tromborg, H. E. Lassen, H. Olesen, Travelling wave analysis of semiconductor lasers, IEEE J. of Quant. El. 30 (1994), No. 5, 939–956.
23. T. I. Zelenjak, On stationary solutions of mixed problems arising in studing of some chemical processes, Differential Equations 2 (1966), No. 2, 205–213.